1
|
Walsh CL, Berg M, West H, Holroyd NA, Walker-Samuel S, Shipley RJ. Reconstructing microvascular network skeletons from 3D images: What is the ground truth? Comput Biol Med 2024; 171:108140. [PMID: 38422956 DOI: 10.1016/j.compbiomed.2024.108140] [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: 10/18/2023] [Revised: 01/29/2024] [Accepted: 02/12/2024] [Indexed: 03/02/2024]
Abstract
Structural changes to microvascular networks are increasingly highlighted as markers of pathogenesis in a wide range of disease, e.g. Alzheimer's disease, vascular dementia and tumour growth. This has motivated the development of dedicated 3D imaging techniques, alongside the creation of computational modelling frameworks capable of using 3D reconstructed networks to simulate functional behaviours such as blood flow or transport processes. Extraction of 3D networks from imaging data broadly consists of two image processing steps: segmentation followed by skeletonisation. Much research effort has been devoted to segmentation field, and there are standard and widely-applied methodologies for creating and assessing gold standards or ground truths produced by manual annotation or automated algorithms. The Skeletonisation field, however, lacks widely applied, simple to compute metrics for the validation or optimisation of the numerous algorithms that exist to extract skeletons from binary images. This is particularly problematic as 3D imaging datasets increase in size and visual inspection becomes an insufficient validation approach. In this work, we first demonstrate the extent of the problem by applying 4 widely-used skeletonisation algorithms to 3 different imaging datasets. In doing so we show significant variability between reconstructed skeletons of the same segmented imaging dataset. Moreover, we show that such a structural variability propagates to simulated metrics such as blood flow. To mitigate this variability we introduce a new, fast and easy to compute super metric that compares the volume, connectivity, medialness, bifurcation point identification and homology of the reconstructed skeletons to the original segmented data. We then show that such a metric can be used to select the best performing skeletonisation algorithm for a given dataset, as well as to optimise its parameters. Finally, we demonstrate that the super metric can also be used to quickly identify how a particular skeletonisation algorithm could be improved, becoming a powerful tool in understanding the complex implication of small structural changes in a network.
Collapse
Affiliation(s)
- Claire L Walsh
- Department of Mechanical Engineering, University College London, United Kingdom
| | - Maxime Berg
- Department of Mechanical Engineering, University College London, United Kingdom.
| | - Hannah West
- Department of Mechanical Engineering, University College London, United Kingdom
| | - Natalie A Holroyd
- Centre for Computational Medicine, Division of Medicine, University College London, United Kingdom
| | - Simon Walker-Samuel
- Centre for Computational Medicine, Division of Medicine, University College London, United Kingdom
| | - Rebecca J Shipley
- Department of Mechanical Engineering, University College London, United Kingdom; Centre for Computational Medicine, Division of Medicine, University College London, United Kingdom
| |
Collapse
|
2
|
Fletcher AG, Osborne JM. Seven challenges in the multiscale modeling of multicellular tissues. WIREs Mech Dis 2022; 14:e1527. [PMID: 35023326 DOI: 10.1002/wsbm.1527] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2020] [Revised: 11/23/2020] [Accepted: 03/25/2021] [Indexed: 11/11/2022]
Abstract
The growth and dynamics of multicellular tissues involve tightly regulated and coordinated morphogenetic cell behaviors, such as shape changes, movement, and division, which are governed by subcellular machinery and involve coupling through short- and long-range signals. A key challenge in the fields of developmental biology, tissue engineering and regenerative medicine is to understand how relationships between scales produce emergent tissue-scale behaviors. Recent advances in molecular biology, live-imaging and ex vivo techniques have revolutionized our ability to study these processes experimentally. To fully leverage these techniques and obtain a more comprehensive understanding of the causal relationships underlying tissue dynamics, computational modeling approaches are increasingly spanning multiple spatial and temporal scales, and are coupling cell shape, growth, mechanics, and signaling. Yet such models remain challenging: modeling at each scale requires different areas of technical skills, while integration across scales necessitates the solution to novel mathematical and computational problems. This review aims to summarize recent progress in multiscale modeling of multicellular tissues and to highlight ongoing challenges associated with the construction, implementation, interrogation, and validation of such models. This article is categorized under: Reproductive System Diseases > Computational Models Metabolic Diseases > Computational Models Cancer > Computational Models.
Collapse
Affiliation(s)
- Alexander G Fletcher
- School of Mathematics and Statistics, University of Sheffield, Sheffield, UK.,Bateson Centre, University of Sheffield, Sheffield, UK
| | - James M Osborne
- School of Mathematics and Statistics, University of Melbourne, Parkville, Victoria, Australia
| |
Collapse
|
3
|
Kurdyaeva T, Milias-Argeitis A. Uncertainty propagation for deterministic models of biochemical networks using moment equations and the extended Kalman filter. J R Soc Interface 2021; 18:20210331. [PMID: 34343452 DOI: 10.1098/rsif.2021.0331] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
Differential equation models of biochemical networks are frequently associated with a large degree of uncertainty in parameters and/or initial conditions. However, estimating the impact of this uncertainty on model predictions via Monte Carlo simulation is computationally demanding. A more efficient approach could be to track a system of low-order statistical moments of the state. Unfortunately, when the underlying model is nonlinear, the system of moment equations is infinite-dimensional and cannot be solved without a moment closure approximation which may introduce bias in the moment dynamics. Here, we present a new method to study the time evolution of the desired moments for nonlinear systems with polynomial rate laws. Our approach is based on solving a system of low-order moment equations by substituting the higher-order moments with Monte Carlo-based estimates from a small number of simulations, and using an extended Kalman filter to counteract Monte Carlo noise. Our algorithm provides more accurate and robust results compared to traditional Monte Carlo and moment closure techniques, and we expect that it will be widely useful for the quantification of uncertainty in biochemical model predictions.
Collapse
Affiliation(s)
- Tamara Kurdyaeva
- Molecular Systems Biology, Groningen Biomolecular Sciences & Biotechnology Institute, University of Groningen, Groningen, Netherlands
| | - Andreas Milias-Argeitis
- Molecular Systems Biology, Groningen Biomolecular Sciences & Biotechnology Institute, University of Groningen, Groningen, Netherlands
| |
Collapse
|
4
|
McCarthy GD, Drewell RA, Dresch JM. Analyzing the stability of gene expression using a simple reaction-diffusion model in an early Drosophila embryo. Math Biosci 2019; 316:108239. [PMID: 31454629 DOI: 10.1016/j.mbs.2019.108239] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2018] [Revised: 08/20/2019] [Accepted: 08/22/2019] [Indexed: 11/28/2022]
Abstract
In all complex organisms, the precise levels and timing of gene expression controls vital biological processes. In higher eukaryotes, including the fruit fly Drosophila melanogaster, the complex molecular control of transcription (the synthesis of RNA from DNA) and translation (the synthesis of proteins from RNA) events driving this gene expression are not fully understood. In particular, for Drosophila melanogaster, there is a plethora of experimental data, including quantitative measurements of both RNA and protein concentrations, but the precise mechanisms that control the dynamics of gene expression during early development and the processes which lead to steady-state levels of certain proteins remain elusive. This study analyzes a current mathematical modeling approach in an attempt to better understand the long-term behavior of gene regulation. The model is a modified reaction-diffusion equation which has been previously employed in predicting gene expression levels and studying the relative contributions of transcription and translation events to protein abundance [10,11,24]. Here, we use Matrix Algebra and Analysis techniques to study the stability of the gene expression system and analyze equilibria, using very general assumptions regarding the parameter values incorporated into the model. We prove that, given realistic biological parameter values, the system will result in a unique, stable equilibrium solution. Additionally, we give an example of this long-term behavior using the model alongside actual experimental data obtained from Drosophila embryos.
Collapse
Affiliation(s)
- Gregory D McCarthy
- School of Natural Science, Hampshire College, Amherst, MA 01002, United States.
| | - Robert A Drewell
- Biology Department, Clark University, Worcester, MA 01610, United States.
| | - Jacqueline M Dresch
- Department of Mathematics and Computer Science, Clark University, Worcester, MA 01610, United States.
| |
Collapse
|
5
|
Wolff HB, Davidson LA, Merks RMH. Adapting a Plant Tissue Model to Animal Development: Introducing Cell Sliding into VirtualLeaf. Bull Math Biol 2019; 81:3322-3341. [PMID: 30927191 PMCID: PMC6677868 DOI: 10.1007/s11538-019-00599-9] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2018] [Accepted: 03/11/2019] [Indexed: 11/16/2022]
Abstract
Cell-based, mathematical modeling of collective cell behavior has become a prominent tool in developmental biology. Cell-based models represent individual cells as single particles or as sets of interconnected particles and predict the collective cell behavior that follows from a set of interaction rules. In particular, vertex-based models are a popular tool for studying the mechanics of confluent, epithelial cell layers. They represent the junctions between three (or sometimes more) cells in confluent tissues as point particles, connected using structural elements that represent the cell boundaries. A disadvantage of these models is that cell-cell interfaces are represented as straight lines. This is a suitable simplification for epithelial tissues, where the interfaces are typically under tension, but this simplification may not be appropriate for mesenchymal tissues or tissues that are under compression, such that the cell-cell boundaries can buckle. In this paper, we introduce a variant of VMs in which this and two other limitations of VMs have been resolved. The new model can also be seen as on off-the-lattice generalization of the Cellular Potts Model. It is an extension of the open-source package VirtualLeaf, which was initially developed to simulate plant tissue morphogenesis where cells do not move relative to one another. The present extension of VirtualLeaf introduces a new rule for cell-cell shear or sliding, from which cell rearrangement (T1) and cell extrusion (T2) transitions emerge naturally, allowing the application of VirtualLeaf to problems of animal development. We show that the updated VirtualLeaf yields different results than the traditional vertex-based models for differential adhesion-driven cell sorting and for the neighborhood topology of soft cellular networks.
Collapse
Affiliation(s)
- Henri B Wolff
- Centrum Wiskunde and Informatica, Science Park 123, 1098 XG, Amsterdam, The Netherlands
- Departments of Bioengineering, Developmental Biology, and Computational and Systems Biology, University of Pittsburgh, Bioscience Tower 3-5059 3501 Fifth Avenue, Pittsburgh, PA, USA
- Department of Epidemiology and Biostatistics, Decision Modeling Center VUmc, Amsterdam UMC location VUmc, PO Box 7057, 1007 MB, Amsterdam, The Netherlands
| | - Lance A Davidson
- Departments of Bioengineering, Developmental Biology, and Computational and Systems Biology, University of Pittsburgh, Bioscience Tower 3-5059 3501 Fifth Avenue, Pittsburgh, PA, USA.
| | - Roeland M H Merks
- Centrum Wiskunde and Informatica, Science Park 123, 1098 XG, Amsterdam, The Netherlands.
- Mathematical Institute, University Leiden, P.O. Box 9512, 2300 RA, Leiden, The Netherlands.
- Mathematical Institute and Institute of Biology, Leiden University, P.O. Box 9505, 2300 RA, Leiden, The Netherlands.
| |
Collapse
|
6
|
Hsieh NH, Reisfeld B, Bois FY, Chiu WA. Applying a Global Sensitivity Analysis Workflow to Improve the Computational Efficiencies in Physiologically-Based Pharmacokinetic Modeling. Front Pharmacol 2018; 9:588. [PMID: 29937730 PMCID: PMC6002508 DOI: 10.3389/fphar.2018.00588] [Citation(s) in RCA: 40] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2018] [Accepted: 05/16/2018] [Indexed: 11/13/2022] Open
Abstract
Traditionally, the solution to reduce parameter dimensionality in a physiologically-based pharmacokinetic (PBPK) model is through expert judgment. However, this approach may lead to bias in parameter estimates and model predictions if important parameters are fixed at uncertain or inappropriate values. The purpose of this study was to explore the application of global sensitivity analysis (GSA) to ascertain which parameters in the PBPK model are non-influential, and therefore can be assigned fixed values in Bayesian parameter estimation with minimal bias. We compared the elementary effect-based Morris method and three variance-based Sobol indices in their ability to distinguish “influential” parameters to be estimated and “non-influential” parameters to be fixed. We illustrated this approach using a published human PBPK model for acetaminophen (APAP) and its two primary metabolites APAP-glucuronide and APAP-sulfate. We first applied GSA to the original published model, comparing Bayesian model calibration results using all the 21 originally calibrated model parameters (OMP, determined by “expert judgment”-based approach) vs. the subset of original influential parameters (OIP, determined by GSA from the OMP). We then applied GSA to all the PBPK parameters, including those fixed in the published model, comparing the model calibration results using this full set of 58 model parameters (FMP) vs. the full set influential parameters (FIP, determined by GSA from FMP). We also examined the impact of different cut-off points to distinguish the influential and non-influential parameters. We found that Sobol indices calculated by eFAST provided the best combination of reliability (consistency with other variance-based methods) and efficiency (lowest computational cost to achieve convergence) in identifying influential parameters. We identified several originally calibrated parameters that were not influential, and could be fixed to improve computational efficiency without discernable changes in prediction accuracy or precision. We further found six previously fixed parameters that were actually influential to the model predictions. Adding these additional influential parameters improved the model performance beyond that of the original publication while maintaining similar computational efficiency. We conclude that GSA provides an objective, transparent, and reproducible approach to improve the performance and computational efficiency of PBPK models.
Collapse
Affiliation(s)
- Nan-Hung Hsieh
- Department of Veterinary Integrative Biosciences, College of Veterinary Medicine and Biomedical Sciences, Texas A&M University, College Station, TX, United States
| | - Brad Reisfeld
- Chemical and Biological Engineering and School of Biomedical Engineering, Colorado State University, Fort Collins, CO, United States
| | | | - Weihsueh A Chiu
- Department of Veterinary Integrative Biosciences, College of Veterinary Medicine and Biomedical Sciences, Texas A&M University, College Station, TX, United States
| |
Collapse
|
7
|
Palm MM, Dallinga MG, van Dijk E, Klaassen I, Schlingemann RO, Merks RMH. Computational Screening of Tip and Stalk Cell Behavior Proposes a Role for Apelin Signaling in Sprout Progression. PLoS One 2016; 11:e0159478. [PMID: 27828952 PMCID: PMC5102492 DOI: 10.1371/journal.pone.0159478] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2016] [Accepted: 05/24/2016] [Indexed: 12/30/2022] Open
Abstract
Angiogenesis involves the formation of new blood vessels by sprouting or splitting of existing blood vessels. During sprouting, a highly motile type of endothelial cell, called the tip cell, migrates from the blood vessels followed by stalk cells, an endothelial cell type that forms the body of the sprout. To get more insight into how tip cells contribute to angiogenesis, we extended an existing computational model of vascular network formation based on the cellular Potts model with tip and stalk differentiation, without making a priori assumptions about the differences between tip cells and stalk cells. To predict potential differences, we looked for parameter values that make tip cells (a) move to the sprout tip, and (b) change the morphology of the angiogenic networks. The screening predicted that if tip cells respond less effectively to an endothelial chemoattractant than stalk cells, they move to the tips of the sprouts, which impacts the morphology of the networks. A comparison of this model prediction with genes expressed differentially in tip and stalk cells revealed that the endothelial chemoattractant Apelin and its receptor APJ may match the model prediction. To test the model prediction we inhibited Apelin signaling in our model and in an in vitro model of angiogenic sprouting, and found that in both cases inhibition of Apelin or of its receptor APJ reduces sprouting. Based on the prediction of the computational model, we propose that the differential expression of Apelin and APJ yields a "self-generated" gradient mechanisms that accelerates the extension of the sprout.
Collapse
Affiliation(s)
- Margriet M. Palm
- Life Sciences Group, Centrum Wiskunde & Informatica, Amsterdam, the Netherlands
| | | | - Erik van Dijk
- Life Sciences Group, Centrum Wiskunde & Informatica, Amsterdam, the Netherlands
| | - Ingeborg Klaassen
- Ocular Angiogenesis Group, Academic Medical Center, Amsterdam, the Netherlands
| | | | - Roeland M. H. Merks
- Life Sciences Group, Centrum Wiskunde & Informatica, Amsterdam, the Netherlands
- Mathematical Institute, Leiden University, Leiden, the Netherlands
| |
Collapse
|
8
|
Jeanquartier F, Jean-Quartier C, Cemernek D, Holzinger A. In silico modeling for tumor growth visualization. BMC SYSTEMS BIOLOGY 2016; 10:59. [PMID: 27503052 PMCID: PMC4977902 DOI: 10.1186/s12918-016-0318-8] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/10/2016] [Accepted: 07/12/2016] [Indexed: 12/18/2022]
Abstract
BACKGROUND Cancer is a complex disease. Fundamental cellular based studies as well as modeling provides insight into cancer biology and strategies to treatment of the disease. In silico models complement in vivo models. Research on tumor growth involves a plethora of models each emphasizing isolated aspects of benign and malignant neoplasms. Biologists and clinical scientists are often overwhelmed by the mathematical background knowledge necessary to grasp and to apply a model to their own research. RESULTS We aim to provide a comprehensive and expandable simulation tool to visualizing tumor growth. This novel Web-based application offers the advantage of a user-friendly graphical interface with several manipulable input variables to correlate different aspects of tumor growth. By refining model parameters we highlight the significance of heterogeneous intercellular interactions on tumor progression. Within this paper we present the implementation of the Cellular Potts Model graphically presented through Cytoscape.js within a Web application. The tool is available under the MIT license at https://github.com/davcem/cpm-cytoscape and http://styx.cgv.tugraz.at:8080/cpm-cytoscape/ . CONCLUSION In-silico methods overcome the lack of wet experimental possibilities and as dry method succeed in terms of reduction, refinement and replacement of animal experimentation, also known as the 3R principles. Our visualization approach to simulation allows for more flexible usage and easy extension to facilitate understanding and gain novel insight. We believe that biomedical research in general and research on tumor growth in particular will benefit from the systems biology perspective.
Collapse
Affiliation(s)
- Fleur Jeanquartier
- Holzinger Group, Research Unit HCI-KDD, Institute for Medical Informatics, Statistics and Documentation, Medical University Graz, Auenbruggerplatz 2/V, 8036, AT, Graz, Austria
| | - Claire Jean-Quartier
- Holzinger Group, Research Unit HCI-KDD, Institute for Medical Informatics, Statistics and Documentation, Medical University Graz, Auenbruggerplatz 2/V, 8036, AT, Graz, Austria
| | - David Cemernek
- Holzinger Group, Research Unit HCI-KDD, Institute for Medical Informatics, Statistics and Documentation, Medical University Graz, Auenbruggerplatz 2/V, 8036, AT, Graz, Austria
| | - Andreas Holzinger
- Holzinger Group, Research Unit HCI-KDD, Institute for Medical Informatics, Statistics and Documentation, Medical University Graz, Auenbruggerplatz 2/V, 8036, AT, Graz, Austria. .,Institute of Information Systems and Computer Media, Graz University of Technology, Inffeldgasse 16c, Graz, 8010, AT, Austria.
| |
Collapse
|