The Padova Type 2 Diabetes Simulator from Triple-Tracer Single-Meal Studies: In Silico Trials Also Possible in Rare but Not-So-Rare Individuals

Distribution-based sub-population selection (DSPS): A method for in-silico reproduction of clinical trials outcomes

Diabetes presents a significant challenge to healthcare due to the short- and long-term complications associated with poor blood sugar control. Computer simulation platforms have emerged as promising tools for advancing diabetes therapy by simulating patient responses to treatments in a virtual environment. The University of Virginia Virtual Lab (UVLab) is a new simulation platform engineered to mimic the metabolic behavior of individuals with type 2 diabetes (T2D) using a mathematical model of glucose homeostasis in T2D and a large population of 6062 virtual subjects. This work proposes a statistical method - the Distribution-based sub-population selection (DSPS) method - for selecting subsets of virtual subjects from this large initial pool, ensuring that the selected group possesses the desired characteristics necessary to reproduce and predict the outcomes of a clinical trial. DSPS formulates the sub-population selection as a linear programming problem, identifying the largest virtual cohort to closely resemble the statistical properties (moments) of key outcomes from real-world clinical trials. The method was applied to the insulin degludec arm of a 26-week phase 3 clinical trial, evaluating the efficacy and safety of insulin degludec and liraglutide combination therapy. DSPS selected a sub-population that mirrored clinical trial data across key metrics, including glycemic efficacy, insulin dosages, and cumulative hypoglycemia events, with a relative sum of square errors of 0.33 and a percentage error of 1.07 %. This approach bridges the gap between large population simulation platforms and clinical trials, enabling the selection of virtual sub-populations with specific properties required for targeted studies.

A new C-peptide minimal model to assess β-cell responsiveness to glucagon in individuals with and without type 2 diabetes

Glucagon regulates its own secretion indirectly by stimulating β-cells to secrete insulin. This may serve as a compensatory mechanism to enhance β-cell function in individuals with, or predisposed to, type 2 diabetes (T2D). However, tools to quantify glucagon-induced C-peptide secretion rate (SR) are lacking. To bridge this gap, we developed a novel model-based method to provide quantitative indices of β-cell function in response to a glucagon bolus in individuals without and with type 2 diabetes. Eight individuals without diabetes [3 M, age = 55 ± 9 yr, body mass index (BMI) = 32 ± 4 kg/m2] and six with T2D (1 M, age = 59 ± 5 yr, BMI = 35 ± 6 kg/m2) underwent a 210-min hyperglycemic clamp (∼9 mmol/L). After 180 min, a 1 mg bolus of glucagon was administered over 1 min, and plasma glucagon and C-peptide concentrations were frequently measured over 30 min. We tested a battery of mathematical models and selected the best one based on standard criteria. The optimal model assumes that C-peptide SR is made up of two components, one proportional to the above basal glucagon, through parameter Γs (static), and one proportional to glucagon rate of change, through parameter Γd (dynamic responsivity to the hormone). An index of total β-cell responsivity to glucagon, Γ, was also derived from Γs and Γd. The model estimated Γs and Γ were significantly higher in individuals without diabetes compared with T2D (P < 0.05), whereas Γd was not. Our findings reveal notable differences in both static and total insulin secretory response to glucagon in people with diabetes as compared with those without diabetes.NEW & NOTEWORTHY In this study, we propose a new mathematical model able to quantify C-peptide secretion and β-cell responsivity in response to a glucagon bolus in individuals without and with type 2 diabetes. We show that individuals with type 2 diabetes exhibit reduced β-cell responsivity to glucagon.

A PopPBPK-RL approach for precision dosing of benazepril in renal impaired patients

Current treatment recommendations mainly rely on rule-based protocols defined from evidence-based clinical guidelines, which are difficult to adapt for high-risk patients such as those with renal impairment. Consequently, unsuccessful therapies and the occurrence of adverse drug reactions are common. Within the context of personalized medicine, that tries to deliver the right treatment dose to maximize efficacy and minimize toxicity, the concept of model-informed precision dosing proposes the use of mechanistic models, like physiologically based pharmacokinetic (PBPK) modeling, to predict drug regimes outcomes. Nonetheless, PBPK models have limited capability when computing patients' centric optimized drug doses. Consequently, reinforcement learning (RL) has been previously used to personalize drug dosage. In this work we propose the first PBPK and RL-based precision dosing system for an orally taken drug (benazepril) considering a virtual population of patients with renal disease. Population based PBPK modeling is used in combination with RL for obtaining patient tailored dose regimes. We also perform patient stratification and feature selection to better handle dose tailoring problems. Based on patients' characteristics with best predictive capabilities, benazepril dose regimes are obtained for a population with features' diversity. Obtained regimes are evaluated based on PK parameters considered. Results show that the proof-of-concept approach herein is capable of learning good dosing regimes for most patients. The use of a PopPBPK model allowed to account for intervariability of patient characteristics and be more inclusive considering also non-frequent patients. Impact analysis of patients' features reveals that renal impairment is the main driver affecting RL capabilities.

Modeling Glucose, Insulin, C-Peptide, and Lactate Interplay in Adolescents During an Oral Glucose Tolerance Test

BACKGROUND

Lactate is not considered just a "waste product" of anaerobic glycolysis anymore. It has been proved to play a key role in several metabolic diseases, such as in the metabolic dysfunction-associated steatotic liver disease, obesity, and diabetes. The capability of simulating glucose-insulin-lactate interaction would be useful to design and test drugs targeting lactate metabolism in such pathological conditions. Minimal models are available, which describe and quantify glucose-lactate interaction but models to simulate postprandial glucose-insulin-C-peptide-lactate time courses are missing. The aim of this study is to fill this gap.

METHODS

Starting from the Padova Type 2 Diabetes Simulator (T2DS), we first added a description of glucose-lactate kinetics and then created a population of 100 in silico subjects to match glucose-insulin-C-peptide-lactate data of 44 adolescents with/without obesity who underwent a standard oral glucose tolerance test (OGTT) of 75 g.

RESULTS

The developed model accurately predicts all molecules time courses, guaranteeing precise model parameter estimates (percent coefficient of variation [CV%] median [25th-75th percentile] = 19 [9-29]%). The generated in silico population shows good agreement with the clinical data in terms of area under the curve (AUC) (P = .6, .6, .9, .6 for glucose, insulin, C-peptide, and lactate, respectively) and parameter distributions (P > .1).

CONCLUSIONS

We have developed a simulator to describe glucose, insulin, C-peptide, and lactate kinetics during an OGTT, which captures the behavior of a real population of adolescents with/without obesity both in terms of average and intersubject variability. Such simulator can be used to investigate the pharmacodynamics of drugs targeting lactate metabolic pathway in various pathological conditions.

Tailoring the Padova Type 2 Diabetes Simulator for Treatment Guidance in Target Populations

OBJECTIVE

The Padova type 2 diabetes (T2D) simulator (T2DS) has been recently proposed to optimize T2D treatments including novel long-acting insulins. It consists of a physiological model and an in silico population describing glucose dynamics, derived from early-stage T2D subjects studied with sophisticated tracer-based experimental techniques. This limits T2DS domain of validity to this specific sub-population. Conversely, running simulations in insulin-naïve or advanced T2D subjects, would be more valuable. However, it is rarely possible or cost-effective to run complex experiments in such populations. Therefore, we propose a method for tuning the T2DS to any desired T2D sub-population using published clinical data. As case study, we extended the T2DS to insulin-naïve T2D subjects, who need to start insulin therapy to compensate the reduced insulin function.

METHODS

T2DS model was identified based on literature data of the target population. The estimated parameters were used to generate a virtual cohort of insulin-naïve T2D subjects (inC1). A model of basal insulin degludec (IDeg) was also incorporated into the T2DS to enable basal insulin therapy. The resulting tailored T2DS was assessed by simulating IDeg therapy initiation and comparing simulated vs. clinical trial outcomes. For further validation, this procedure was reiterated to generate a new cohort of insulin-naïve T2D (inC2) assuming inC1 as target population.

RESULTS

No statistically significant differences were found when comparing fasting plasma glucose and IDeg dose, neither in clinical data vs. inC1, nor inC1 vs. inC2.

CONCLUSIONS

The tuned T2DS allowed reproducing the main findings of clinical studies in insulin-naïve T2D subjects.

SIGNIFICANCE

The proposed methodology makes the Padova T2DS usable for supporting treatment guidance in target T2D populations.

Digital twins and artificial intelligence in metabolic disease research

Digital twin technology is emerging as a transformative paradigm for personalized medicine in the management of chronic conditions. In this article, we explore the concept and key characteristics of a digital twin and its applications in chronic non-communicable metabolic disease management, with a focus on diabetes case studies. We cover various types of digital twin models, including mechanistic models based on ODEs, data-driven ML algorithms, and hybrid modeling strategies that combine the strengths of both approaches. We present successful case studies demonstrating the potential of digital twins in improving glucose outcomes for individuals with T1D and T2D, and discuss the benefits and challenges of translating digital twin research applications to clinical practice.

Smart Algorithms for Efficient Insulin Therapy Initiation in Individuals With Type 2 Diabetes: An in Silico Study

BACKGROUND

Insulin-naive subjects with type 2 diabetes (T2D) start basal insulin titration from a low initial insulin dose (IID), which is adjusted weekly or twice per week based on fasting plasma glucose (FPG) measurement as recommended by the American Diabetes Association (ADA). The procedure to reach the optimal insulin dose (OID) is time-consuming, especially in subjects with high insulin needs (HIN). The aim of this study is to provide a fast and effective, but still safe, insulin titration algorithm in insulin-naive T2D subjects with HIN.

METHOD

To do that, we in silico cloned 300 subjects, matching a real population of insulin-naive T2D and used a logistic regression model to classify them as subjects with HIN or subjects with low insulin needs (LIN). Then, we applied to the subjects with HIN both a more aggressive insulin dose initiation (SMART-IID) and two newly developed titration algorithms (continuous glucose monitoring [CGM]-BASED and SMART-CGM-BASED) in which CGM was used to guide the decision-making process.

RESULTS

The new titration algorithm applied to HIN-classified individuals guaranteed a faster reaching of OID, with significant improvements in time in range (TIR) and reduction in time above range (TAR) in the first months of the trial, without any clinically significant increase in the risk of hypoglycemia.

CONCLUSIONS

Smart basal insulin titration algorithms enable insulin-naive T2D individuals to achieve OID and improve their glycemic control faster than standard guidelines, without jeopardizing patient safety.

Short- and Long-Term Effects on Glucose Control of Nonadherence to Insulin Therapy in People With Type 2 Diabetes An In Silico Study

BACKGROUND

Strict adherence to multiple daily insulin (MDI) therapy is a cornerstone for the achievement of good glucose control in people with advanced type 2 diabetes (T2D). Here, we aim to in silico assess glucose control in T2D subjects with poor adherence to MDI therapy.

METHODS

We tuned the Padova T2D Simulator, originally describing early-stage T2D physiology, around advanced T2D people. One hundred in silico advanced T2D subjects were generated and equipped with optimal MDI therapy: specifically, basal and bolus insulin amounts and injection times were individualized for each subject by applying titration algorithms that iteratively update insulin dose based on glucose deviation from its target. Then, the effect of nonadhering to MDI therapy was assessed using standard glucose control metrics calculated in two 6-month 3-meal/day in silico scenarios: in Scenario 1, subjects received the optimal basal and prandial insulin bolus at each meal; in Scenario 2, subjects received optimal basal insulin and randomly delayed or skipped the prandial insulin bolus in 3 lunches during working days and 1 dinner during weekends.

RESULTS

A statistically significant degradation was found in all glucose control outcome metrics in Scenario 2 versus Scenario 1: e.g., percent time above 180 mg/dL increased by 22.2% and glucose management index by 0.2%.

CONCLUSIONS

Impaired adherence to MDI therapy in T2D leads to glucose control deteriorations in both short and long terms. Interestingly, short-term hyperglycemia seems being contrasted by residual endogenous insulin secretion, which statistically increased by 3-fold after delayed/skipped insulin boluses compared with optimal ones.

Socioeconomic Level and the Relationship in Glycemic Behavior in the Mexican Population. A Nutritional Alternative Focused on Vulnerable Populations

In this study, numerical approximations were generated to analyze the behavior of glycemic variations in the non-diabetic population of the Mexican republic. The main objective of this work is to obtain an overview of the glycemic variations in the non-diabetic population from different socioeconomic statuses in Mexico (Medium-high, medium, and low). Thus, evaluating the effect on the glucose level under a healthy diet considering the socioeconomic capabilities of the population. Through the national health and nutrition survey of Mexico 2020 and the Mexican food base, 1420 virtual patients were proposed (522 low status, 485 medium status and 413 Medium-High status) focused on simulating the glycemic behavior in each of the survey participants. Considering that the average food expenditure of the Mexican population is $107.00 MXN, and the cost of a healthy diet is $66.50 MXN, the economic sustainability of the Mexican population to adopt a healthy diet is revealed. The particularity of this work is focused on obtaining diverse data that are difficult to access in the development of population analyses. Such is the case of the approach proposed for different socioeconomic statuses. In this way, the proposed methodology provides a framework for complementary research contributions to the subject.

Evolutionary algorithm for the optimization of meal intake and insulin administration in patients with type 2 diabetes

The optimal management of type 2 diabetes (T2DM) is complex and involves an appropriate combination of diet, exercise, and different pharmacological treatments. Artificial intelligence-based tools have been shown to be very useful for the diagnosis and treatment of diverse pathologies, including diabetes. In the present study, we present a proof of concept of the potential of an evolutionary algorithm to optimize the meal size, timing and insulin dose for the control of glycemia. We found that an appropriate distribution of food intake throughout the day permits a reduction in the insulin dose required to maintain glycemia within the range recommended by the American Diabetes Association for patients with T2DM of a range of severities. Furthermore, the effects of restrictions to both the timing and amount of food ingested were assessed, and we found that an increase in the amount of insulin was required to control glycemia as dietary intake became more restricted. In the near future, the use of these computational tools should permit patients with T2DM to optimize their personal meal schedule and insulin dose, according to the severity of their diabetes.

Gluclas: A software for computer-aided modulation of glucose infusion in glucose clamp experiments

BACKGROUND AND OBJECTIVE

The glucose clamp (GC) is an experimental technique for assessing several aspects of glucose metabolism. In these experiments, investigators face the non-trivial challenge of accurately adjusting the rate of intravenous glucose infusion to drive subjects' blood glucose (BG) concentration towards a desired plateau level. In this work we present Gluclas, an open-source software to support researchers in the modulation of glucose infusion rate (GIR) during GC experiments.

METHODS

Gluclas uses a proportional-integrative-derivative controller to provide GIR suggestions based on BG measurements. The controller embeds an anti-wind-up scheme to account for actuator physical limits and suitable corrections of control action to accommodate for possible sampling jitter due to manual measurement and actuation. The software also provides a graphic user interface to increase its usability. A preliminary validation of the controller is performed for different clamp scenarios (hyperglycemic, euglycemic, hypoglycemic) on a simulator of glucose metabolism in healthy subjects, which also includes models of measurement error and sampling delay for increased realism. In silico trials are performed on 50 virtual subjects. We also report the results of the first in-vivo application of the software in three subjects undergoing a hypoglycemic clamp.

RESULTS

In silico, during the plateau period, the coefficient of variation (CV) is in median below 5% for every protocol, with 5% being considered the threshold for sufficient quality. In terms of median [5th percentile, 95th percentile], average BG level during the plateau period is 12.18 [11.58 - 12.53] mmol/l in the hyperglycemic clamp (target: 12.4 mmol/), 4.92 [4.51 - 5.14] mmol/l in the euglycemic clamp (target: 5.5 mmol/) and 2.38 [2.33 - 2.64] in the hypoglycemic clamp (target: 2.5 mmol/). Results in vivo are consistent with those obtained in silico during the plateau period: average BG levels are between 2.56 and 2.68 mmol/l (target: 2.5 mmol/l); CV is below 5% for all three experiments.

CONCLUSIONS

Gluclas offered satisfactory control for tested GC protocols. Although its safety and efficacy need to be further validated in vivo, this preliminary validation suggest that Gluclas offers a reliable and non-expensive solution for reducing investigator bias and improving the quality of GC experiments.

Minimal and Maximal Models to Quantitate Glucose Metabolism: Tools to Measure, to Simulate and to Run in Silico Clinical Trials

Several models have been proposed to describe the glucose system at whole-body, organ/tissue and cellular level, designed to measure non-accessible parameters (minimal models), to simulate system behavior and run in silico clinical trials (maximal models). Here, we will review the authors' work, by putting it into a concise historical background. We will discuss first the parametric portrait provided by the oral minimal models-building on the classical intravenous glucose tolerance test minimal models-to measure otherwise non-accessible key parameters like insulin sensitivity and beta-cell responsivity from a physiological oral test, the mixed meal or the oral glucose tolerance tests, and what can be gained by adding a tracer to the oral glucose dose. These models were used in various pathophysiological studies, which we will briefly review. A deeper understanding of insulin sensitivity can be gained by measuring insulin action in the skeletal muscle. This requires the use of isotopic tracers: both the classical multiple-tracer dilution and the positron emission tomography techniques are discussed, which quantitate the effect of insulin on the individual steps of glucose metabolism, that is, bidirectional transport plasma-interstitium, and phosphorylation. Finally, we will present a cellular model of insulin secretion that, using a multiscale modeling approach, highlights the relations between minimal model indices and subcellular secretory events. In terms of maximal models, we will move from a parametric to a flux portrait of the system by discussing the triple tracer meal protocol implemented with the tracer-to-tracee clamp technique. This allows to arrive at quasi-model independent measurement of glucose rate of appearance (Ra), endogenous glucose production (EGP), and glucose rate of disappearance (Rd). Both the fast absorbing simple carbs and the slow absorbing complex carbs are discussed. This rich data base has allowed us to build the UVA/Padova Type 1 diabetes and the Padova Type 2 diabetes large scale simulators. In particular, the UVA/Padova Type 1 simulator proved to be a very useful tool to safely and effectively test in silico closed-loop control algorithms for an artificial pancreas (AP). This was the first and unique simulator of the glucose system accepted by the U.S. Food and Drug Administration as a substitute to animal trials for in silico testing AP algorithms. Recent uses of the simulator have looked at glucose sensors for non-adjunctive use and new insulin molecules.

Dynamic of Glucose Homeostasis in Virtual Patients: A Comparison between Different Behaviors

This work presents a mathematical model of homeostasis dynamics in healthy individuals, focusing on the generation of conductive data on glucose homeostasis throughout the day under dietary and physical activity factors. Two case studies on glucose dynamics for populations under conditions of physical activity and sedentary lifestyle were developed. For this purpose, two types of virtual populations were generated, the first population was developed according to the data of a total of 89 physical persons between 20 and 75 years old and the second was developed using the Monte Carlo approach, obtaining a total of 200 virtual patients. In both populations, each participant was classified as an active or sedentary person depending on the physical activity performed. The results obtained demonstrate the capacity of virtual populations in the generation of in-silico approximations similar to those obtained from in-vivo studies. Obtaining information that is only achievable through specific in-vivo experiments. Being a tool that generates information for the approach of alternatives in the prevention of the development of type 2 Diabetes.

Modeling Between-Subject Variability in Subcutaneous Absorption of a Long-Acting Insulin Glargine 100 U/mL by a Nonlinear Mixed Effects Approach

Subcutaneous insulin absorption is well-known to vary significantly both between and within subjects (BSV and WSV, respectively). This variability considerably obstacles the establishing of a reproducible and effective insulin therapy. Some models exist to describe the subcutaneous kinetics of both fast and long-acting insulin analogues; however, none of them account for the BSV. The aim of this study is to develop a nonlinear mixed effects model able to describe the BSV observed in the subcutaneous absorption of a long-acting insulin glargine 100 U/mL. Four stochastic models of the BSV were added to a previously validated model of subcutaneous absorption of insulin glargine 100 U/mL. These were assessed on a database of 47 subjects with type 1 diabetes. The best model was selected based on residual analysis, precision of the estimates and parsimony criteria. The selected model provided good fit of individual data, precise population parameter estimates and allowed quantifying the BSV of the insulin glargine 100 U/mL pharmacokinetics. Future model development will include the description of the WSV of long- acting insulin absorption.

A Novel Method for Generation of In Silico Subjects with Type 2 Diabetes

A type 2 diabetes (T2D) simulator has been recently proposed for supporting drug development and treatment optimization. This tool consists of a physiological model of glucose/insulin/C-peptide dynamics and a virtual cohort of T2D subjects (i.e., random extractions of model parameterizations from a joint parameter distribution) well describing both average and variability realistic T2D dynamics . However, the state-of-art procedure to get a reliable virtual population requires some post-processing after subject extraction, in order to discard implausible behaviors. We propose an improved method for virtual subjects' generation to overcome this burdensome task. To do so, we first assessed a refined joint parameter distribution, from which extracting a number of subjects, greater than the target population size. Then, target-size subsets are undersampled from the large cohort. The final virtual population is selected among the subsets as the one maximizing the similarity with T2D data and model parameter distribution, by means of measurement' outcome metrics and Euclidian distance (Δ), respectively. In the final population, almost all the outcome metrics are statistically identical to the clinical counterparts (p-value>0.05) and model parameters' distribution differs by ~5-10% from that derived from data. The methodology described here is flexible, thus resulting suitable for different T2D stages and type 1 diabetes, as well.Clinical Relevance- A straightforward subjects' generation would ease the availability of tailored in silico trials for testing diabetes treatment in a specific population.

Modeling Between-Subject Variability in Subcutaneous Absorption of a Fast-Acting Insulin Analogue by a Nonlinear Mixed Effects Approach

Despite the great progress made in insulin preparation and titration, many patients with diabetes are still experiencing dangerous fluctuations in their blood glucose levels. This is mainly due to the large between- and within-subject variability, which considerably hampers insulin therapy, leading to defective dosing and timing of the administration process. In this work, we present a nonlinear mixed effects model describing the between-subject variability observed in the subcutaneous absorption of fast-acting insulin. A set of 14 different models was identified on a large and frequently-sampled database of lispro pharmacokinetic data, collected from 116 subjects with type 1 diabetes. The tested models were compared, and the best one was selected on the basis of the ability to fit the data, the precision of the estimated parameters, and parsimony criteria. The selected model was able to accurately describe the typical trend of plasma insulin kinetics, as well as the between-subject variability present in the absorption process, which was found to be related to the subject’s body mass index. The model provided a deeper understanding of the insulin absorption process and can be incorporated into simulation platforms to test and develop new open- and closed-loop treatment strategies, allowing a step forward toward personalized insulin therapy.

Mathematical Modeling for the Physiological and Clinical Investigation of Glucose Homeostasis and Diabetes

Mathematical modeling in the field of glucose metabolism has a longstanding tradition. The use of models is motivated by several reasons. Models have been used for calculating parameters of physiological interest from experimental data indirectly, to provide an unambiguous quantitative representation of pathophysiological mechanisms, to determine indices of clinical usefulness from simple experimental tests. With the growing societal impact of type 2 diabetes, which involves the disturbance of the glucose homeostasis system, development and use of models in this area have increased. Following the approaches of physiological and clinical investigation, the focus of the models has spanned from representations of whole body processes to those of cells, i.e., from in vivo to in vitro research. Model-based approaches for linking in vivo to in vitro research have been proposed, as well as multiscale models merging the two areas. The success and impact of models has been variable. Two kinds of models have received remarkable interest: those widely used in clinical applications, e.g., for the assessment of insulin sensitivity and β-cell function and some models representing specific aspects of the glucose homeostasis system, which have become iconic for their efficacy in describing clearly and compactly key physiological processes, such as insulin secretion from the pancreatic β cells. Models are inevitably simplified and approximate representations of a physiological system. Key to their success is an appropriate balance between adherence to reality, comprehensibility, interpretative value and practical usefulness. This has been achieved with a variety of approaches. Although many models concerning the glucose homeostasis system have been proposed, research in this area still needs to address numerous issues and tackle new opportunities. The mathematical representation of the glucose homeostasis processes is only partial, also because some mechanisms are still only partially understood. For in vitro research, mathematical models still need to develop their potential. This review illustrates the problems, approaches and contribution of mathematical modeling to the physiological and clinical investigation of glucose homeostasis and diabetes, focusing on the most relevant and stimulating models.

In Silico Cloning of Target Type 2 Diabetes Population for Treatments Development and Decision Support. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2020;2020:5111-5114. [PMID: 33019136 DOI: 10.1109/embc44109.2020.9175271]

Therapies for treatment of type 2 diabetes (T2D) involve a variety of medications, depending on the stage of T2D progression. It is now an accepted knowledge that in silico trials can help to accelerate drug development and support treatment optimization. A T2D simulator (T2DS), consisting of a model of the glucose-insulin system and an in silico population describing glucose-insulin dynamics in T2D subjects, has been recently developed based on early-stage T2D data, studied with sophisticated experimental techniques. This limits the domain of validity of the simulator to this specific sub-population of T2D. Here we proposed a method for tuning the T2DS to any desired T2D target population, e.g. insulin-naïve (i.e., not experienced with insulin) patients, without the need to resort to complex and expensive clinical studies. This will allow to use the T2DS for testing treatments in the target population. To illustrate the methodology, we used a case study: extending the T2DS to reproduce the behavior of insulin-naïve T2D subjects. The methodology described here can be extended to other stages of T2D, allowing an extensive in silico testing phase of different treatments before human trials.