In Silico Head-to-Head Comparison of Insulin Glargine 300 U/mL and Insulin Degludec 100 U/mL in Type 1 Diabetes

Analysis of Glucose Responsive Glucagon Therapeutics using Computational Models of the Glucoregulatory System

Glucose-responsive glucagon (GRG) therapeutics are a promising technology for reducing the risk of severe hypoglycemia as a complication of diabetes mellitus. Herein, the performance of candidate GRGs in the literature by modeling the kinetics of activation and connecting them as input into physiological glucoregulatory models is evaluated and projected the two distinct GRG designs, experimental results reported in Wu et al. (GRG-I) and Webber et al. (GRG-II) is considered. Both are evaluated using a multi-compartmental glucoregulatory model (IMPACT) and used to compare in-vivo experimental data of therapeutic performance in rats and mice. For GRG-I and GRG-II, the total integrated glucose material balances are overestimated by 41.5% ± 14% and underestimated by 24.8% ± 16% compared to in-vivo time-course data, respectively. These large differences to the relatively simple computational descriptions of glucagon dynamics in the model, which underscores the urgent need for improved glucagon models is attributed. Additionally, therapeutic insulin and glucagon infusion pumps are modeled for type 1 diabetes mellitus (T1DM) human subjects to extend the results to additional datasets. These observations suggest that both the representative physiological and non-physiological models considered in this work require additional refinement to successfully describe clinical data that involve simultaneous, coupled insulin, glucose, and glucagon dynamics.

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.

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.

Developing the UVA/Padova Type 1 Diabetes Simulator: Modeling, Validation, Refinements, and Utility

Arguably, diabetes mellitus is one of the best quantified human conditions. In the past 50 years, the metabolic monitoring technologies progressed from occasional assessment of average glycemia via HbA1c, through episodic blood glucose readings, to continuous glucose monitoring (CGM) producing data points every few minutes. The high-temporal resolution of CGM data enabled increasingly intensive treatments, from decision support assisting insulin injection or oral medication, to automated closed-loop control, known as the "artificial pancreas." Throughout this progress, mathematical models and computer simulation of the human metabolic system became indispensable for the technological progress of diabetes treatment, enabling every step, from assessment of insulin sensitivity via the now classic Minimal Model of Glucose Kinetics, to in silico trials replacing animal experiments, to automated insulin delivery algorithms. In this review, we follow these developments, beginning with the Minimal Model, which evolved through the years to become large and comprehensive and trigger a paradigm change in the design of diabetes optimization strategies: in 2007, we introduced a sophisticated model of glucose-insulin dynamics and a computer simulator equipped with a "population" of N = 300 in silico "subjects" with type 1 diabetes. In January 2008, in an unprecedented decision, the Food and Drug Administration (FDA) accepted this simulator as a substitute to animal trials for the pre-clinical testing of insulin treatment strategies. This opened the field for rapid and cost-effective development and pre-clinical testing of new treatment approaches, which continues today. Meanwhile, animal experiments for the purpose of designing new insulin treatment algorithms have been abandoned.

A New Index of Insulin Sensitivity from Glucose Sensor and Insulin Pump Data: In Silico and In Vivo Validation in Youths with Type 1 Diabetes

Background: Estimation of insulin sensitivity (SI) and its daily variation are key for optimizing insulin therapy in patients with type 1 diabetes (T1D). We recently developed a method for SI estimation from continuous glucose monitoring (CGM) and continuous subcutaneous insulin infusion (CSII) data in adults with T1D (SISP) and validated it under restrained experimental conditions. Herein, we validate in vivo a new version of SISP performing well in daily life unrestrained conditions. Methods: The new SISP was tested in both simulated and real data. The simulated dataset consists of 100 virtual adults of the UVa/Padova T1D Simulator monitored during an open-loop experiment, whereas the real dataset consists of 10 youths with T1D monitored during a hybrid closed-loop meal study. In both datasets, participants underwent two consecutive meals (breakfast and lunch, at 7 and 11 am) with the same carbohydrate content (70 g). Plasma glucose and insulin were measured during each meal to estimate the oral glucose minimal model SI (SIMM). CGM and CSII data were used for SISP calculation, which was then validated against the gold standard SIMM. Results: SISP was estimated with good precision (median coefficient of variation <20%) in 100% of the real and 91% of the simulated meals. SISP and SIMM were highly correlated, both in the simulated and real datasets (R = 0.82 and R = 0.83, P < 0.001), and exhibited a similar intraday pattern. Conclusions: SISP is suitable for estimating SI in both closed- and open-loop settings, provided that the subject wears a CGM sensor and a subcutaneous insulin pump.

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.

A software interface for in silico testing of type 2 diabetes treatments

BACKGROUND AND OBJECTIVE

The increasing incidence of diabetes continuously stimulates the research on new antidiabetic drugs. Computer simulation can save time and costs, alleviating the need of animal trials and providing useful information for optimal experiment design and drug dosing. We recently presented a type 2 diabetes (T2D) simulator as tool for in silico testing of new molecules and guiding treatment optimization. Here we present a user-friendly interface aimed to increase the usability of the simulator.

METHOD

The simulator, based on a large-scale glucose, insulin, and C-peptide model and equipped with 100 virtual subjects well describing system dynamics in a real T2D population, is extended to incorporate pharmacokinetics/pharmacodynamics (PK/PD) of a drug of interest. A graphical interface is developed on top of the simulator, allowing an easy design of in silico experiments: specifically, it is possible to select the population size to test, design the experiment (crossover or parallel), its duration and the sampling grid, choose glucose and insulin doses, and define treatment PK/PD and dose administered. The simulator also provides the outcome metrics requested by the user, and performs statistical comparisons among treatments and/or placebo.

RESULTS

To illustrate the potential of the simulator, we provided a case study using metformin and liraglutide. Literature-based PK/PD models of metformin and liraglutide have been incorporated in the simulator, by modulating key drug-sensitive model parameters. An in silico placebo-controlled trial has been done by simulating a three-arm meal tolerance test with subjects receiving placebo, metformin 850 mg, liraglutide 1.80 mg, respectively. The obtained results are in agreement with the clinical evidences, in terms of main glucose, insulin, and C-peptide outcome metrics.

CONCLUSIONS

We developed a user-friendly software interface for the T2D simulator to support the design and test of new antidiabetic drugs and treatments. This increases the simulator usability, making it suitable also for users who have low experience with computer programming.

A simulator with realistic and challenging scenarios for virtual T1D patients undergoing CSII and MDI therapy

In silico simulations have become essential for the development of diabetes treatments. However, currently available simulators are not challenging enough and often suffer from limitations in insulin and meal absorption variability, which is unable to realistically reflect the dynamics of people with type 1 diabetes (T1D). Additionally, T1D simulators are mainly designed for the testing of continuous subcutaneous insulin infusion (CSII) therapies. In this work, a simulator is presented that includes a generated virtual patient (VP) cohort and both fast- and long-acting Glargine-100 U/ml (Gla-100), Glargine-300 U/ml (Gla-300), and Degludec-100 U/ml (Deg-100) insulin models. Therefore, in addition to CSII therapies, multiple daily injections (MDI) therapies can also be tested. The Hovorka model and its published parameter probability distributions were used to generate cohorts of VPs that represent a T1D population. Valid patients are filtered through restrictions that guarantee that they are physiologically acceptable. To obtain more realistic scenarios, basal insulin profile patterns from the literature have been used to identify variability in insulin sensitivity. A library of mixed meals identified from real data has also been included. This work presents and validates a methodology for the creation of realistic VP cohorts that include physiological variability and a simulator that includes challenging and realistic scenarios for in silico testing. A cohort of 47 VPs has been generated and in silico simulations of both CSII and MDI therapies were performed in open-loop. The simulation outcome metrics were contrasted with literature results.

Evaluation of Reference Metrics for Continuous Glucose Monitoring in Persons Without Diabetes and Prediabetes

BACKGROUND

Recent guidelines have been developed for continuous glucose monitoring (CGM) metrics in persons with diabetes. To understand what glucose profiles should be judged as normal in clinical practice and glucose-lowering trials, we examined the glucose profile of healthy individuals using CGM.

METHODS

Persons without diabetes or prediabetes were included after passing a normal oral glucose tolerance test, two-hour value <8.9 mmol/L (160 mg/dL), fasting glucose <6.1 mmol/L (110 mg/dL), and HbA1c <6.0% (<42 mmol/mol). CGM metrics were evaluated using the Dexcom G4 Platinum.

RESULTS

In total, 60 persons were included, mean age was 43.0 years, 70.0% were women, mean HbA1c was 5.3% (34 mmol/mol), and mean body mass index was 25.7 kg/m². Median and mean percent times in hypoglycemia <3.9 mmol/L (70 mg/dL) were 1.6% (IQR 0.6-3.2), and 3.2% (95% CI 2.0; 4.3), respectively. For glucose levels <3.0 mmol/L (54 mg/dL), the corresponding estimates were 0.0% (IQR 0.0-0.4) and 0.5% (95% CI 0.2; 0.8). Median and mean time-in-range (3.9-10.0 mmol/L [70-180 mg/dL]) was 97.3% (IQR 95.4-98.7) and 95.4% (95% CI 94.0; 96.8), respectively. Median and mean standard deviations were 1.04 mmol/L (IQR 0.92-1.29) and 1.15 mmol/L (95% CI 1.05; 1.24), respectively. Measures of glycemic variability (standard deviation, coefficient of variation, mean amplitude of glycemic excursions) were significantly greater during daytime compared with nighttime, whereas others did not differ.

CONCLUSIONS

People without prediabetes or diabetes show a non-negligible % time in hypoglycemia, median 1.6% and mean 3.2%, which needs to be accounted for in clinical practice and glucose-lowering trials. Glycemic variability measures differ day and night in this population.

Methods to Develop an in silico Clinical Trial: Computational Head-to-Head Comparison of Lisdexamfetamine and Methylphenidate

Regulatory agencies encourage computer modeling and simulation to reduce the time and cost of clinical trials. Although still not classified in formal guidelines, system biology-based models represent a powerful tool for generating hypotheses with great molecular detail. Herein, we have applied a mechanistic head-to-head in silico clinical trial (ISCT) between two treatments for attention-deficit/hyperactivity disorder, to wit lisdexamfetamine (LDX) and methylphenidate (MPH). The ISCT was generated through three phases comprising (i) the molecular characterization of drugs and pathologies, (ii) the generation of adult and children virtual populations (vPOPs) totaling 2,600 individuals and the creation of physiologically based pharmacokinetic (PBPK) and quantitative systems pharmacology (QSP) models, and (iii) data analysis with artificial intelligence methods. The characteristics of our vPOPs were in close agreement with real reference populations extracted from clinical trials, as did our PBPK models with in vivo parameters. The mechanisms of action of LDX and MPH were obtained from QSP models combining PBPK modeling of dosing schemes and systems biology-based modeling technology, i.e., therapeutic performance mapping system. The step-by-step process described here to undertake a head-to-head ISCT would allow obtaining mechanistic conclusions that could be extrapolated or used for predictions to a certain extent at the clinical level. Altogether, these computational techniques are proven an excellent tool for hypothesis-generation and would help reach a personalized medicine.

Modeling Intraperitoneal Insulin Absorption in Patients with Type 1 Diabetes

Standard insulin therapy to treat type 1 diabetes (T1D) consists of exogenous insulin administration through the subcutaneous (SC) tissue. Despite recent advances in insulin formulations, the SC route still suffers from delays and large inter/intra-subject variability that limiting optimal glucose control. Intraperitoneal (IP) insulin administration, despite its higher invasiveness, was shown to represent a valid alternative to the SC one. To date, no mathematical model describing the absorption and distribution of insulin after IP administration is available. Here, we aim to fill this gap by using data from eight patients with T1D, treated by implanted IP pump, studied in a hospitalized setting, with frequent measurements of plasma insulin and glucose concentration. A battery of models describing insulin kinetics after IP administration were tested. Model comparison and selection were performed based on model ability to predict the data, precision of parameters and parsimony criteria. The selected model assumed that the insulin absorption from the IP space was described by a linear, two-compartment model, coupled with a two-compartment model of whole-body insulin kinetics with hepatic insulin extraction controlled by hepatic insulin. Future developments include model incorporation into the UVa/Padova T1D Simulator for testing open- and closed-loop therapies with IP insulin administration.

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 Models of Meal Amount and Timing Variability With Implementation in the Type-1 Diabetes Patient Decision Simulator

BACKGROUND

In type 1 diabetes (T1D) research, in-silico clinical trials (ISCTs) have proven effective in accelerating the development of new therapies. However, published simulators lack a realistic description of some aspects of patient lifestyle which can remarkably affect glucose control. In this paper, we develop a mathematical description of meal carbohydrates (CHO) amount and timing, with the aim to improve the meal generation module in the T1D Patient Decision Simulator (T1D-PDS) published in Vettoretti et al.

METHODS

Data of 32 T1D subjects under free-living conditions for 4874 days were used. Univariate probability density function (PDF) parametric models with different candidate shapes were fitted, individually, against sample distributions of: CHO amounts of breakfast (CHOB), lunch (CHOL), dinner (CHOD), and snack (CHOS); breakfast timing (TB); and time between breakfast-lunch (TBL) and between lunch-dinner (TLD). Furthermore, a support vector machine (SVM) classifier was developed to predict the occurrence of a snack in future fixed-length time windows. Once embedded inside the T1D-PDS, an ISCT was performed.

RESULTS

Resulting PDF models were: gamma (CHOB, CHOS), lognormal (CHOL, TB), loglogistic (CHOD), and generalized-extreme-values (TBL, TLD). The SVM showed a classification accuracy of 0.8 over the test set. The distributions of simulated meal data were not statistically different from the distributions of the real data used to develop the models (α = 0.05).

CONCLUSIONS

The models of meal amount and timing variability developed are suitable for describing real data. Their inclusion in modules that describe patient behavior in the T1D-PDS can permit investigators to perform more realistic, reliable, and insightful ISCTs.

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

Background:In silico trials in type 2 diabetes (T2D) would be useful for testing diabetes treatments and accelerating the development of new antidiabetic drugs. In this study, we present a T2D simulator able to reproduce the variability observed in a T2D population. The simulator also allows to safely experiment on virtual subjects with severe (and possibly rare) pathological conditions. Methods: A meal simulation model of glucose, insulin, and C-peptide systems, made of 15 differential equations and 39 parameters, has been identified using a system decomposition and forcing function Bayesian strategy on data of 51 T2D subjects undergoing a single triple-tracer mixed meal. One hundred T2D in silico subjects have been generated from the joint distribution of estimated model parameters. A case study is presented to illustrate the simulator use for testing a virtual drug (improving insulin action and secretion) in a subpopulation of rare, extremely impaired, T2D subjects. Results: The model well fitted T2D data and parameters were estimated with precision. Simulated plasma glucose, insulin, and C-peptide well matched the data (e.g., median [25^th-75^th percentile] glucose area under the curves of 6.9 [6.1-8.5] 10⁴ mg/dL·min in silico vs. 7.0 [5.6-8.2] 10⁴ mg/dL·min in vivo). The potential use of the simulator was shown in a case study, in which the (virtual) antidiabetic drug dose was optimized for very insulin-resistant T2D subjects. Conclusions: We have developed a T2D simulator that captures the behavior of T2D population during a meal, both in terms of average and intersubject variability. The simulator represents a cost-effective way to test new antidiabetic drugs, before moving to human trials.

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.