A Systems Analysis of Phenotype Heterogeneity in APOE*3Leiden.CETP Mice Induced by Long-Term High-Fat High-Cholesterol Diet Feeding

Within the human population, considerable variability exists between individuals in their susceptibility to develop obesity and dyslipidemia. In humans, this is thought to be caused by both genetic and environmental variation. APOE*3-Leiden.CETP mice, as part of an inbred mouse model in which mice develop the metabolic syndrome upon being fed a high-fat high-cholesterol diet, show large inter-individual variation in the parameters of the metabolic syndrome, despite a lack of genetic and environmental variation. In the present study, we set out to resolve what mechanisms could underlie this variation. We used measurements of glucose and lipid metabolism from a six-month longitudinal study on the development of the metabolic syndrome. Mice were classified as mice with either high plasma triglyceride (responders) or low plasma triglyceride (non-responders) at the baseline. Subsequently, we fitted the data to a dynamic computational model of whole-body glucose and lipid metabolism (MINGLeD) by making use of a hybrid modelling method called Adaptations in Parameter Trajectories (ADAPT). ADAPT integrates longitudinal data, and predicts how the parameters of the model must change through time in order to comply with the data and model constraints. To explain the phenotypic variation in plasma triglycerides, the ADAPT analysis suggested a decreased cholesterol absorption, higher energy expenditure and increased fecal fatty acid excretion in non-responders. While decreased cholesterol absorption and higher energy expenditure could not be confirmed, the experimental validation demonstrated that the non-responders were indeed characterized by increased fecal fatty acid excretion. Furthermore, the amount of fatty acids excreted strongly correlated with bile acid excretion, in particular deoxycholate. Since bile acids play an important role in the solubilization of lipids in the intestine, these results suggest that variation in bile acid homeostasis may in part drive the phenotypic variation in the APOE*3-Leiden.CETP mice.

Computational modelling of energy balance in individuals with Metabolic Syndrome

Background

A positive energy balance is considered to be the primary cause of the development of obesity-related diseases. Treatment often consists of a combination of reducing energy intake and increasing energy expenditure. Here we use an existing computational modelling framework describing the long-term development of Metabolic Syndrome (MetS) in APOE3L.CETP mice fed a high-fat diet containing cholesterol with a human-like metabolic system. This model was used to analyze energy expenditure and energy balance in a large set of individual model realizations.

Results

We developed and applied a strategy to select specific individual models for a detailed analysis of heterogeneity in energy metabolism. Models were stratified based on energy expenditure. A substantial surplus of energy was found to be present during MetS development, which explains the weight gain during MetS development. In the majority of the models, energy was mainly expended in the peripheral tissues, but also distinctly different subgroups were identified.

In silico perturbation of the system to induce increased peripheral energy expenditure implied changes in lipid metabolism, but not in carbohydrate metabolism. In silico analysis provided predictions for which individual models increase of peripheral energy expenditure would be an effective treatment.

Conclusion

The computational analysis confirmed that the energy imbalance plays an important role in the development of obesity. Furthermore, the model is capable to predict whether an increase in peripheral energy expenditure – for instance by cold exposure to activate brown adipose tissue (BAT) – could resolve MetS symptoms.

Electronic supplementary material

The online version of this article (10.1186/s12918-019-0705-z) contains supplementary material, which is available to authorized users.

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.

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.

Spatiotemporal mapping of interictal epileptiform discharges in human absence epilepsy: A MEG study

PURPOSE

Although absence epilepsy is considered to be a prototypic type of generalized epilepsy, it is still under debate whether generalized 3 Hz spike-and-wave discharges (SWDs) might have a cortical focal origin. Here it is investigated whether focal interictal epileptiform discharges (IEDs), which typically occur in the electro- (EEG) and magnetoencephalogram (MEG) in case of focal epilepsy, are present in the MEG of children with absence epilepsy. Next, the location of the sources of the IEDs is established, and it is investigated whether the location is concordant to the earlier established focal cortical regions involved in the generalized SWDs of these children.

METHODS

Whole head MEG recordings of seven children with absence epilepsy were reviewed with respect to the presence of IEDs (spikes and sharp waves). These IEDs were grouped into distinct clusters, in which each contribution to a cluster yields a comparable magnetic field distribution. Source localization was then performed onto the average signal of each cluster using an equivalent current dipole model and a realistic head model of the cortical surface.

RESULTS

IEDs were detected in 6 out of 7 patients. Source reconstruction indicated most often frontal, central or parietal origins of the IED in either the left and or right hemisphere. Spatiotemporal assessment of the IEDs indicated a stable location of the averages of these discharges, indicating a single underlying cortical source.

DISCUSSION

The outcome of this pilot study shows that MEG is well suited for the detection of IEDs and suggests that their estimated sources coincide rather well with the cortical regions involved during the spikes of the SWDs. It is discussed whether the presence of IEDs, classically seen as a marker of focal epilepsies, indicate that absence epilepsy should be considered as a focal type of epilepsy, in which changes in the network are evolving rapidly.

A physiology-based model describing heterogeneity in glucose metabolism: the core of the Eindhoven Diabetes Education Simulator (E-DES)

Current diabetes education methods are costly, time-consuming, and do not actively engage the patient. Here, we describe the development and verification of the physiological model for healthy subjects that forms the basis of the Eindhoven Diabetes Education Simulator (E-DES). E-DES shall provide diabetes patients with an individualized virtual practice environment incorporating the main factors that influence glycemic control: food, exercise, and medication. The physiological model consists of 4 compartments for which the inflow and outflow of glucose and insulin are calculated using 6 nonlinear coupled differential equations and 14 parameters. These parameters are estimated on 12 sets of oral glucose tolerance test (OGTT) data (226 healthy subjects) obtained from literature. The resulting parameter set is verified on 8 separate literature OGTT data sets (229 subjects). The model is considered verified if 95% of the glucose data points lie within an acceptance range of ±20% of the corresponding model value. All glucose data points of the verification data sets lie within the predefined acceptance range. Physiological processes represented in the model include insulin resistance and β-cell function. Adjusting the corresponding parameters allows to describe heterogeneity in the data and shows the capabilities of this model for individualization. We have verified the physiological model of the E-DES for healthy subjects. Heterogeneity of the data has successfully been modeled by adjusting the 4 parameters describing insulin resistance and β-cell function. Our model will form the basis of a simulator providing individualized education on glucose control.