1
|
Kobayashi M, Miyauchi A, Jimbo EF, Oishi N, Aoki S, Watanabe M, Yoshikawa Y, Akiyama Y, Yamagata T, Osaka H. Synthetic aporphine alkaloids are potential therapeutics for Leigh syndrome. Sci Rep 2024; 14:11561. [PMID: 38773300 PMCID: PMC11109252 DOI: 10.1038/s41598-024-62445-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2024] [Accepted: 05/16/2024] [Indexed: 05/23/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Mizuki Kobayashi
- Department of Pediatrics, Division of Pediatrics, Jichi Medical University, 3311-1 Yakushiji, Shimotsuke, Tochigi, 329-0498, Japan
| | - Akihiko Miyauchi
- Department of Pediatrics, Division of Pediatrics, Jichi Medical University, 3311-1 Yakushiji, Shimotsuke, Tochigi, 329-0498, Japan
| | - Eriko F Jimbo
- Department of Pediatrics, Division of Pediatrics, Jichi Medical University, 3311-1 Yakushiji, Shimotsuke, Tochigi, 329-0498, Japan
| | - Natsumi Oishi
- Department of Pediatrics, Division of Pediatrics, Jichi Medical University, 3311-1 Yakushiji, Shimotsuke, Tochigi, 329-0498, Japan
| | - Shiho Aoki
- Department of Pediatrics, Division of Pediatrics, Jichi Medical University, 3311-1 Yakushiji, Shimotsuke, Tochigi, 329-0498, Japan
| | - Miyuki Watanabe
- Department of Pediatrics, Division of Pediatrics, Jichi Medical University, 3311-1 Yakushiji, Shimotsuke, Tochigi, 329-0498, Japan
| | - Yasushi Yoshikawa
- School of Life Science and Technology, Tokyo Institute of Technology, Yokohama, 226-8501, Japan
- Department of Computer Science, School of Computing, Tokyo Institute of Technology, Meguro-ku, Tokyo, 152-8550, Japan
- Middle-Molecule IT-Based Drug Discovery Laboratory (MIDL), Tokyo Institute of Technology, Kawasaki, Kanagawa, 210-0821, Japan
| | - Yutaka Akiyama
- School of Life Science and Technology, Tokyo Institute of Technology, Yokohama, 226-8501, Japan
- Department of Computer Science, School of Computing, Tokyo Institute of Technology, Meguro-ku, Tokyo, 152-8550, Japan
- Middle-Molecule IT-Based Drug Discovery Laboratory (MIDL), Tokyo Institute of Technology, Kawasaki, Kanagawa, 210-0821, Japan
| | - Takanori Yamagata
- Department of Pediatrics, Division of Pediatrics, Jichi Medical University, 3311-1 Yakushiji, Shimotsuke, Tochigi, 329-0498, Japan
| | - Hitoshi Osaka
- Department of Pediatrics, Division of Pediatrics, Jichi Medical University, 3311-1 Yakushiji, Shimotsuke, Tochigi, 329-0498, Japan.
| |
Collapse
|
2
|
Casas B, Vilén L, Bauer S, Kanebratt KP, Wennberg Huldt C, Magnusson L, Marx U, Andersson TB, Gennemark P, Cedersund G. Integrated experimental-computational analysis of a HepaRG liver-islet microphysiological system for human-centric diabetes research. PLoS Comput Biol 2022; 18:e1010587. [PMID: 36260620 PMCID: PMC9621595 DOI: 10.1371/journal.pcbi.1010587] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2022] [Revised: 10/31/2022] [Accepted: 09/19/2022] [Indexed: 11/05/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Belén Casas
- Drug Metabolism and Pharmacokinetics, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca, Gothenburg, Sweden
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden
| | - Liisa Vilén
- Drug Metabolism and Pharmacokinetics, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca, Gothenburg, Sweden
| | | | - Kajsa P. Kanebratt
- Drug Metabolism and Pharmacokinetics, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca, Gothenburg, Sweden
| | - Charlotte Wennberg Huldt
- Bioscience, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca, Gothenburg, Sweden
| | - Lisa Magnusson
- Bioscience, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca, Gothenburg, Sweden
| | | | - Tommy B. Andersson
- Drug Metabolism and Pharmacokinetics, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca, Gothenburg, Sweden
| | - Peter Gennemark
- Drug Metabolism and Pharmacokinetics, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca, Gothenburg, Sweden
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden
| | - Gunnar Cedersund
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden
- Center for Medical Image Science and Visualization (CMIV), Linköping University, Linköping, Sweden
- * E-mail:
| |
Collapse
|
3
|
Lövfors W, Ekström J, Jönsson C, Strålfors P, Cedersund G, Nyman E. A systems biology analysis of lipolysis and fatty acid release from adipocytes in vitro and from adipose tissue in vivo. PLoS One 2021; 16:e0261681. [PMID: 34972146 PMCID: PMC8719686 DOI: 10.1371/journal.pone.0261681] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2021] [Accepted: 12/07/2021] [Indexed: 12/03/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- William Lövfors
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden
- Department of Mathematics, Linköping University, Linköping, Sweden
| | - Jona Ekström
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden
| | - Cecilia Jönsson
- Department of Biomedical and Clinical Sciences, Linköping University, Linköping, Sweden
| | - Peter Strålfors
- Department of Biomedical and Clinical Sciences, Linköping University, Linköping, Sweden
| | - Gunnar Cedersund
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden
- Center for Medical Image Science and Visualization (CMIV), Linköping University, Linköping, Sweden
| | - Elin Nyman
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden
| |
Collapse
|
4
|
Ngara B, Zvada S, Chawana TD, Nhachi CFB, Rusakaniko S. Pharmacokinetic-pharmacodynamic modelling of atazanavir in hair among adolescents on antiretroviral treatment in Zimbabwe. BMC Pharmacol Toxicol 2021; 22:29. [PMID: 34030726 PMCID: PMC8147021 DOI: 10.1186/s40360-021-00497-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2020] [Accepted: 05/04/2021] [Indexed: 11/17/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Bernard Ngara
- Department of Community Medicine, University of Zimbabwe College of Health Sciences, Mazowe Street, Parirenyatwa Complex, P. O Box A178 Avondale, Harare, Zimbabwe.
| | - Simbarashe Zvada
- Department of Clinical Pharmacology, Stellenbosch University, Private Bag X1, Matieland, Stellenbosch, 7602, South Africa
| | - Tariro Dianah Chawana
- Department of Clinical Pharmacology, University of Zimbabwe College of Health Sciences, Mazowe Street, Parirenyatwa Complex, P. O Box A178 Avondale, Harare, Zimbabwe
| | - Charles Fungai Brian Nhachi
- Department of Clinical Pharmacology, University of Zimbabwe College of Health Sciences, Mazowe Street, Parirenyatwa Complex, P. O Box A178 Avondale, Harare, Zimbabwe
| | - Simbarashe Rusakaniko
- Department of Community Medicine, University of Zimbabwe College of Health Sciences, Mazowe Street, Parirenyatwa Complex, P. O Box A178 Avondale, Harare, Zimbabwe
| |
Collapse
|
5
|
Joyce W, Egginton S, Farrell AP, Axelsson M. Adrenergic and adenosinergic regulation of the cardiovascular system in an Antarctic icefish: Insight into central and peripheral determinants of cardiac output. Comp Biochem Physiol A Mol Integr Physiol 2019; 230:28-38. [DOI: 10.1016/j.cbpa.2018.12.012] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2018] [Revised: 12/19/2018] [Accepted: 12/19/2018] [Indexed: 01/27/2023]
|
6
|
Properties of cell signaling pathways and gene expression systems operating far from steady-state. Sci Rep 2018; 8:17035. [PMID: 30451879 PMCID: PMC6242903 DOI: 10.1038/s41598-018-34766-0] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2018] [Accepted: 10/23/2018] [Indexed: 11/10/2022] Open
Abstract
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 EC50 (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.
Collapse
|
7
|
Systems biology reveals uncoupling beyond UCP1 in human white fat-derived beige adipocytes. NPJ Syst Biol Appl 2017; 3:29. [PMID: 28983409 PMCID: PMC5626775 DOI: 10.1038/s41540-017-0027-y] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2017] [Revised: 08/04/2017] [Accepted: 08/23/2017] [Indexed: 01/08/2023] Open
Abstract
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.
Collapse
|
8
|
Brännmark C, Lövfors W, Komai AM, Axelsson T, El Hachmane MF, Musovic S, Paul A, Nyman E, Olofsson CS. Mathematical modeling of white adipocyte exocytosis predicts adiponectin secretion and quantifies the rates of vesicle exo- and endocytosis. J Biol Chem 2017; 292:20032-20043. [PMID: 28972187 DOI: 10.1074/jbc.m117.801225] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2017] [Revised: 09/22/2017] [Indexed: 11/06/2022] Open
Abstract
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 Ca2+ 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 Ca2+-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.
Collapse
Affiliation(s)
- Cecilia Brännmark
- Department of Physiology/Metabolic Physiology, Institute of Neuroscience and Physiology, The Sahlgrenska Academy at University of Gothenburg, Medicinaregatan 11, SE-405 30 Göteborg
| | - William Lövfors
- Departments of Biomedical Engineering, SE-581 83 Linköping; Mathematics, Linköping University, SE-581 83 Linköping
| | - Ali M Komai
- Department of Physiology/Metabolic Physiology, Institute of Neuroscience and Physiology, The Sahlgrenska Academy at University of Gothenburg, Medicinaregatan 11, SE-405 30 Göteborg
| | - Tom Axelsson
- Departments of Biomedical Engineering, SE-581 83 Linköping
| | - Mickaël F El Hachmane
- Department of Physiology/Metabolic Physiology, Institute of Neuroscience and Physiology, The Sahlgrenska Academy at University of Gothenburg, Medicinaregatan 11, SE-405 30 Göteborg
| | - Saliha Musovic
- Department of Physiology/Metabolic Physiology, Institute of Neuroscience and Physiology, The Sahlgrenska Academy at University of Gothenburg, Medicinaregatan 11, SE-405 30 Göteborg
| | - Alexandra Paul
- Department of Biology and Biological Engineering, Chalmers University of Technology, Kemivägen 10 SE-412 96 Göteborg
| | - Elin Nyman
- Departments of Biomedical Engineering, SE-581 83 Linköping; Cardiovascular and Metabolic Diseases iMed Biotech Unit, AstraZeneca R&D, 431 83 Gothenburg, Sweden.
| | - Charlotta S Olofsson
- Department of Physiology/Metabolic Physiology, Institute of Neuroscience and Physiology, The Sahlgrenska Academy at University of Gothenburg, Medicinaregatan 11, SE-405 30 Göteborg.
| |
Collapse
|
9
|
Schleicher J, Conrad T, Gustafsson M, Cedersund G, Guthke R, Linde J. Facing the challenges of multiscale modelling of bacterial and fungal pathogen-host interactions. Brief Funct Genomics 2017; 16:57-69. [PMID: 26857943 PMCID: PMC5439285 DOI: 10.1093/bfgp/elv064] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
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.
Collapse
Affiliation(s)
| | | | | | | | | | - Jörg Linde
- Corresponding author: Jörg Linde, Leibniz Institute for Natural Product Research and Infection Biology—Hans Knöll Institute, Jena, Germany. Tel.: +49-3641-532-1290; E-mail:
| |
Collapse
|
10
|
Lundengård K, Cedersund G, Sten S, Leong F, Smedberg A, Elinder F, Engström M. Mechanistic Mathematical Modeling Tests Hypotheses of the Neurovascular Coupling in fMRI. PLoS Comput Biol 2016; 12:e1004971. [PMID: 27310017 PMCID: PMC4911100 DOI: 10.1371/journal.pcbi.1004971] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2015] [Accepted: 05/04/2016] [Indexed: 11/22/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Karin Lundengård
- Department of Medical and Health Sciences, Linköping University, Linköping, Sweden
- Center for Medical Image Science and Visualization (CMIV), Linköping University, Linköping, Sweden
| | - Gunnar Cedersund
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden
- Department of Clinical and Experimental Medicine, Linköping University, Linköping, Sweden
| | - Sebastian Sten
- Department of Medical and Health Sciences, Linköping University, Linköping, Sweden
| | - Felix Leong
- Department of Medical and Health Sciences, Linköping University, Linköping, Sweden
| | - Alexander Smedberg
- Department of Medical and Health Sciences, Linköping University, Linköping, Sweden
| | - Fredrik Elinder
- Department of Clinical and Experimental Medicine, Linköping University, Linköping, Sweden
| | - Maria Engström
- Department of Medical and Health Sciences, Linköping University, Linköping, Sweden
- Center for Medical Image Science and Visualization (CMIV), Linköping University, Linköping, Sweden
- * E-mail:
| |
Collapse
|
11
|
Nyman E, Rozendaal YJW, Helmlinger G, Hamrén B, Kjellsson MC, Strålfors P, van Riel NAW, Gennemark P, Cedersund G. Requirements for multi-level systems pharmacology models to reach end-usage: the case of type 2 diabetes. Interface Focus 2016; 6:20150075. [PMID: 27051506 DOI: 10.1098/rsfs.2015.0075] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Elin Nyman
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden; CVMD iMed DMPK AstraZeneca R&D, Gothenburg, Sweden
| | - Yvonne J W Rozendaal
- Department of Biomedical Engineering , Eindhoven University of Technology , Eindhoven , The Netherlands
| | - Gabriel Helmlinger
- Quantitative Clinical Pharmacology, AstraZeneca , Pharmaceuticals LP, Waltham, MA , USA
| | - Bengt Hamrén
- Quantitative Clinical Pharmacology , AstraZeneca , Gothenburg , Sweden
| | - Maria C Kjellsson
- Department of Pharmaceutical Biosciences , Uppsala University , Uppsala , Sweden
| | - Peter Strålfors
- Department of Clinical and Experimental Medicine , Linköping University , Linköping , Sweden
| | - Natal A W van Riel
- Department of Biomedical Engineering , Eindhoven University of Technology , Eindhoven , The Netherlands
| | | | - Gunnar Cedersund
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden; Department of Clinical and Experimental Medicine, Linköping University, Linköping, Sweden
| |
Collapse
|
12
|
Karlsson M, Janzén DLI, Durrieu L, Colman-Lerner A, Kjellsson MC, Cedersund G. Nonlinear mixed-effects modelling for single cell estimation: when, why, and how to use it. BMC SYSTEMS BIOLOGY 2015; 9:52. [PMID: 26335227 PMCID: PMC4559169 DOI: 10.1186/s12918-015-0203-x] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/28/2015] [Accepted: 08/22/2015] [Indexed: 11/29/2022]
Abstract
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.
Collapse
Affiliation(s)
- Markus Karlsson
- Department of Biomedical Engineering, Linköping University, Linköping, SE-58185, Sweden.
| | - David L I Janzén
- Department of Biomedical Engineering, Linköping University, Linköping, SE-58185, Sweden. .,Department of Clinical and Experimental Medicine, Linköping University, Uppsala, SE-58185, Sweden. .,Current Address: Biomedical and Biological Systems Laboratory, School of Engineering, University of Warwick, Coventry, CV4 7AL, UK. .,Modeling and Simulation, AstraZeneca, Mölndal, Sweden. .,Department of Systems and Data Analysis, Fraunhofer-Chalmers Centre, Chalmers Science Park, Gothenburg, SE-412 88, Sweden.
| | - Lucia Durrieu
- Instituto de Fisiología, Biología Molecular y Neurociencias, Consejo Nacional de Investigaciones Científicas y Técnicas and Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, Argentina.
| | - Alejandro Colman-Lerner
- Instituto de Fisiología, Biología Molecular y Neurociencias, Consejo Nacional de Investigaciones Científicas y Técnicas and Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, Argentina.
| | - Maria C Kjellsson
- Pharmacometrics Group, Pharmaceutical Biosciences, Uppsala University, Uppsala, SE-75124, Sweden.
| | - Gunnar Cedersund
- Department of Biomedical Engineering, Linköping University, Linköping, SE-58185, Sweden. .,Department of Clinical and Experimental Medicine, Linköping University, Uppsala, SE-58185, Sweden. .,IKE, Linköping University, Linköping, 58185, Sweden.
| |
Collapse
|