Mathematical modeling improves EC50 estimations from classical dose-response curves

Synthetic aporphine alkaloids are potential therapeutics for Leigh syndrome

Mitochondrial diseases are mainly caused by dysfunction of mitochondrial respiratory chain complexes and have a variety of genetic variants or phenotypes. There are only a few approved treatments, and fundamental therapies are yet to be developed. Leigh syndrome (LS) is the most severe type of progressive encephalopathy. We previously reported that apomorphine, an anti- "off" agent for Parkinson's disease, has cell-protective activity in patient-derived skin fibroblasts in addition to strong dopamine agonist effect. We obtained 26 apomorphine analogs, synthesized 20 apomorphine derivatives, and determined their anti-cell death effect, dopamine agonist activity, and effects on the mitochondrial function. We found three novel apomorphine derivatives with an active hydroxy group at position 11 of the aporphine framework, with a high anti-cell death effect without emetic dopamine agonist activity. These synthetic aporphine alkaloids are potent therapeutics for mitochondrial diseases without emetic side effects and have the potential to overcome the low bioavailability of apomorphine. Moreover, they have high anti-ferroptotic activity and therefore have potential as a therapeutic agent for diseases related to ferroptosis.

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.

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.

Pharmacokinetic-pharmacodynamic modelling of atazanavir in hair among adolescents on antiretroviral treatment in Zimbabwe

BACKGROUND

Drug potency is a pharmacological parameter defining dose or concentration of drug required to obtain 50% of the drug's maximal effect. Pharmacokinetic-pharmacodynamic modelling and simulation allows estimation of potency and evaluate strategies improving treatment outcome. The objective of our study is to determine potency of atazanavir in hair, defined as atazanavir level in hair associated with 50% probability of failing to achieve viral load below 1000 copies/ml among adolescents, and explore the effect of participant specific variables on potency.

METHODS

A secondary analysis was performed on data from a previous study conducted in HIV-infected adolescents failing 2nd line ART from Harare central hospital, Zimbabwe, between 2015 and 2016. We simulated atazanavir concentrations in hair using NONMEM (version 7.3) ADVAN 13, based on a previously established pharmacokinetic model. Logistic regression methods were used for PKPD analysis. Simulations utilising PKPD model focused on estimation of potency and exploring the effect of covariates.

RESULTS

The potency of atazanavir in hair was found to be 4.5 ng/mg hair before adjusting for covariate effects. Participants at three months follow-up, reporting adequate adherence, having normal BMI-for-age, and cared for by mature guardians had increased potency of atazanavir in hair of 2.6 ng/mg, however the follow-up event was the only statistically significant factor at 5% level.

CONCLUSION

Atazanavir in hair in the range 2.6 to 4.5 ng/mg is associated with above 50% probability of early viral load suppression. Adherence monitoring to adolescents with lower potency of atazanavir is recommended. The effect self-reported adherence level, BMI-for-age, and caregiver status require further evaluation.

Properties of cell signaling pathways and gene expression systems operating far from steady-state

Ligand-receptor systems, covalent modification cycles, and transcriptional networks are basic units of signaling systems and their steady-state properties are well understood. However, the behavior of such systems before steady-state is poorly characterized. Here, we analyzed the properties of input-output curves for each of these systems as they approach steady-state. In ligand-receptor systems, the EC₅₀ (concentration of the ligand that occupies 50% of the receptors) is higher before the system reaches steady-state. Based on this behavior, we have previously defined PRESS (for pre-equilibrium sensing and signaling), a general “systems level” mechanism cells may use to overcome input saturation. Originally, we showed that, given a step stimulation, PRESS operates when the kinetics of ligand-receptor binding are slower than the downstream signaling steps. Now, we show that, provided the input increases slowly, it is not essential for the ligand binding reaction itself to be slow. In addition, we demonstrate that covalent modification cycles and gene expression systems may also operate in PRESS mode. Thus, nearly all biochemical processes may operate in PRESS mode, suggesting that this mechanism may be ubiquitous in cell signaling systems.

Systems biology reveals uncoupling beyond UCP1 in human white fat-derived beige adipocytes

Pharmaceutical induction of metabolically active beige adipocytes in the normally energy storing white adipose tissue has potential to reduce obesity. Mitochondrial uncoupling in beige adipocytes, as in brown adipocytes, has been reported to occur via the uncoupling protein 1 (UCP1). However, several previous in vitro characterizations of human beige adipocytes have only measured UCP1 mRNA fold increase, and assumed a direct correlation with metabolic activity. Here, we provide an example of pharmaceutical induction of beige adipocytes, where increased mRNA levels of UCP1 are not translated into increased protein levels, and perform a thorough analysis of this example. We incorporate mRNA and protein levels of UCP1, time-resolved mitochondrial characterizations, and numerous perturbations, and analyze all data with a new fit-for-purpose mathematical model. The systematic analysis challenges the seemingly obvious experimental conclusion, i.e., that UCP1 is not active in the induced cells, and shows that hypothesis testing with iterative modeling and experimental work is needed to sort out the role of UCP1. The analyses demonstrate, for the first time, that the uncoupling capability of human beige adipocytes can be obtained without UCP1 activity. This finding thus opens the door to a new direction in drug discovery that targets obesity and its associated comorbidities. Furthermore, the analysis advances our understanding of how to evaluate UCP1-independent thermogenesis in human beige adipocytes.

Mathematical modeling of white adipocyte exocytosis predicts adiponectin secretion and quantifies the rates of vesicle exo- and endocytosis

Adiponectin is a hormone secreted from white adipocytes and takes part in the regulation of several metabolic processes. Although the pathophysiological importance of adiponectin has been thoroughly investigated, the mechanisms controlling its release are only partly understood. We have recently shown that adiponectin is secreted via regulated exocytosis of adiponectin-containing vesicles, that adiponectin exocytosis is stimulated by cAMP-dependent mechanisms, and that Ca²⁺ and ATP augment the cAMP-triggered secretion. However, much remains to be discovered regarding the molecular and cellular regulation of adiponectin release. Here, we have used mathematical modeling to extract detailed information contained within our previously obtained high-resolution patch-clamp time-resolved capacitance recordings to produce the first model of adiponectin exocytosis/secretion that combines all mechanistic knowledge deduced from electrophysiological experimental series. This model demonstrates that our previous understanding of the role of intracellular ATP in the control of adiponectin exocytosis needs to be revised to include an additional ATP-dependent step. Validation of the model by introduction of data of secreted adiponectin yielded a very close resemblance between the simulations and experimental results. Moreover, we could show that Ca²⁺-dependent adiponectin endocytosis contributes to the measured capacitance signal, and we were able to predict the contribution of endocytosis to the measured exocytotic rate under different experimental conditions. In conclusion, using mathematical modeling of published and newly generated data, we have obtained estimates of adiponectin exo- and endocytosis rates, and we have predicted adiponectin secretion. We believe that our model should have multiple applications in the study of metabolic processes and hormonal control thereof.

Facing the challenges of multiscale modelling of bacterial and fungal pathogen-host interactions

Recent and rapidly evolving progress on high-throughput measurement techniques and computational performance has led to the emergence of new disciplines, such as systems medicine and translational systems biology. At the core of these disciplines lies the desire to produce multiscale models: mathematical models that integrate multiple scales of biological organization, ranging from molecular, cellular and tissue models to organ, whole-organism and population scale models. Using such models, hypotheses can systematically be tested. In this review, we present state-of-the-art multiscale modelling of bacterial and fungal infections, considering both the pathogen and host as well as their interaction. Multiscale modelling of the interactions of bacteria, especially Mycobacterium tuberculosis, with the human host is quite advanced. In contrast, models for fungal infections are still in their infancy, in particular regarding infections with the most important human pathogenic fungi, Candida albicans and Aspergillus fumigatus. We reflect on the current availability of computational approaches for multiscale modelling of host-pathogen interactions and point out current challenges. Finally, we provide an outlook for future requirements of multiscale modelling.

Mechanistic Mathematical Modeling Tests Hypotheses of the Neurovascular Coupling in fMRI

Functional magnetic resonance imaging (fMRI) measures brain activity by detecting the blood-oxygen-level dependent (BOLD) response to neural activity. The BOLD response depends on the neurovascular coupling, which connects cerebral blood flow, cerebral blood volume, and deoxyhemoglobin level to neuronal activity. The exact mechanisms behind this neurovascular coupling are not yet fully investigated. There are at least three different ways in which these mechanisms are being discussed. Firstly, mathematical models involving the so-called Balloon model describes the relation between oxygen metabolism, cerebral blood volume, and cerebral blood flow. However, the Balloon model does not describe cellular and biochemical mechanisms. Secondly, the metabolic feedback hypothesis, which is based on experimental findings on metabolism associated with brain activation, and thirdly, the neurotransmitter feed-forward hypothesis which describes intracellular pathways leading to vasoactive substance release. Both the metabolic feedback and the neurotransmitter feed-forward hypotheses have been extensively studied, but only experimentally. These two hypotheses have never been implemented as mathematical models. Here we investigate these two hypotheses by mechanistic mathematical modeling using a systems biology approach; these methods have been used in biological research for many years but never been applied to the BOLD response in fMRI. In the current work, model structures describing the metabolic feedback and the neurotransmitter feed-forward hypotheses were applied to measured BOLD responses in the visual cortex of 12 healthy volunteers. Evaluating each hypothesis separately shows that neither hypothesis alone can describe the data in a biologically plausible way. However, by adding metabolism to the neurotransmitter feed-forward model structure, we obtained a new model structure which is able to fit the estimation data and successfully predict new, independent validation data. These results open the door to a new type of fMRI analysis that more accurately reflects the true neuronal activity.

Functional magnetic resonance imaging (fMRI) is a widely used technique for measuring brain activity. However, the signal registered by fMRI is not a direct measurement of the neuronal activity in the brain, but it is influenced by the interplay between the metabolism, blood flow and blood volume in the active area. This signal is called the blood-oxygen-level dependent (BOLD) response and occurs when the blood supply to the active area increases in response to neuronal activity. The mechanisms that the cells use to influence the blood supply are not fully known, and therefore it is difficult to know the true neuronal signalling only from inspection of the fMRI signal. In this article, we present a new mathematical model built on the physiological mechanisms thought to underlie the BOLD response. We could successfully fit the model to data and predict the activity caused by new stimuli. By using the validated model we investigated physiological mechanisms that cause different parts of the BOLD response.

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.

Nonlinear mixed-effects modelling for single cell estimation: when, why, and how to use it

Background

Studies of cell-to-cell variation have in recent years grown in interest, due to improved bioanalytical techniques which facilitates determination of small changes with high uncertainty. Like much high-quality data, single-cell data is best analysed using a systems biology approach. The most common systems biology approach to single-cell data is the standard two-stage (STS) approach. In STS, data from each cell is analysed in a separate sub-problem, meaning that only data from the same cell is used to calculate the parameter values within that cell. Because only parts of the data are considered, problems with parameter unidentifiability are exaggerated in STS. In contrast, a related approach to data analysis has been developed for the studies of patient-to-patient variations. This approach, called nonlinear mixed-effects modelling (NLME), makes use of all data, when estimating the patient-specific parameters. NLME would therefore be advantageous compared to STS also for the study of cell-to-cell variation. However, no such systematic evaluation of the two approaches exists.

Results

Herein, such a systematic comparison between STS and NLME has been performed. Different examples, both linear and nonlinear, and both simulated and real experimental data, have been examined. With informative data, there is no significant difference in the results for either parameter or noise estimation. However, when data becomes uninformative, NLME is significantly superior to STS. These results hold independently of whether the loss of information is due to a low signal-to-noise ratio, too few data points, or a bad input signal. The improvement is shown to come from both the consideration of a joint likelihood (JLH) function, describing all parameters and data, and from an a priori postulated form of the population parameters. Finally, we provide a small tutorial that shows how to use NLME for single-cell analysis, using the free and user-friendly software Monolix.

Conclusions

When considering uninformative single-cell data, NLME yields more accurate parameter and noise estimates, compared to more traditional approaches, such as STS and JLH.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-015-0203-x) contains supplementary material, which is available to authorized users.