Reconstructing microvascular network skeletons from 3D images: What is the ground truth?

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.

Local Administration of Minocycline Improves Nerve Regeneration in Two Rat Nerve Injury Models.

Peripheral nerve injuries are quite common and often require a surgical intervention. However, even after surgery, patients do not often regain satisfactory sensory and motor functions. This, in turn, results in a heavy socioeconomic burden. To some extent, neurons can regenerate from the proximal nerve stump and try to reconnect to the distal stump. However, this regenerating capacity is limited, and depending on the type and size of peripheral nerve injury, this process may not lead to a positive outcome. To date, no pharmacological approach has been used to improve nerve regeneration following repair surgery. We elected to investigate the effects of local delivery of minocycline on nerve regeneration. This molecule has been studied in the central nervous system and was shown to improve the outcome in many disease models. In this study, we first tested the effects of minocycline on SCL 4.1/F7 Schwann cells in vitro and on sciatic nerve explants. We specifically focused on the Schwann cell repair phenotype, as these cells play a central role in orchestrating nerve regeneration. Finally, we delivered minocycline locally in two different rat models of nerve injury, a sciatic nerve transection and a sciatic nerve autograft, demonstrating the capacity of local minocycline treatment to improve nerve regeneration.

Micro to macro scale analysis of the intact human renal arterial tree with Synchrotron Tomography.

Background: The kidney vasculature is exquisitely structured to orchestrate renal function. Structural profiling of the vasculature in intact rodent kidneys, has provided insights into renal haemodynamics and oxygenation, but has never been extended to the human kidney beyond a few vascular generations. We hypothesised that synchrotron-based imaging of a human kidney would enable assessment of vasculature across the whole organ. Methods: An intact kidney from a 63-year-old male was scanned using hierarchical phase-contrast tomography (HiP-CT), followed by semi-automated vessel segmentation and quantitative analysis. These data were compared to published micro-CT data of whole rat kidney. Results: The intact human kidney vascular network was imaged with HiP-CT at 25 µm voxels, representing a 20-fold increase in resolution compared to clinical CT scanners. Our comparative quantitative analysis revealed the number of vessel generations, vascular asymmetry and a structural organisation optimised for minimal resistance to flow, are conserved between species, whereas the normalised radii are not. We further demonstrate regional heterogeneity in vessel geometry between renal cortex, medulla, and hilum, showing how the distance between vessels provides a structural basis for renal oxygenation and hypoxia. Conclusions: Through the application of HiP-CT, we have provided the first quantification of the human renal arterial network, with a resolution comparable to that of light microscopy yet at a scale several orders of magnitude larger than that of a renal punch biopsy. Our findings bridge anatomical scales, profiling blood vessels across the intact human kidney, with implications for renal physiology, biophysical modelling, and tissue engineering.

A combined experimental and computational framework to evaluate the behavior of therapeutic cells for peripheral nerve regeneration.

Recent studies have explored the potential of tissue-mimetic scaffolds in encouraging nerve regeneration. One of the major determinants of the regenerative success of cellular nerve repair constructs (NRCs) is the local microenvironment, particularly native low oxygen conditions which can affect implanted cell survival and functional performance. In vivo, cells reside in a range of environmental conditions due to the spatial gradients of nutrient concentrations that are established. Here we evaluate in vitro the differences in cellular behavior that such conditions induce, including key biological features such as oxygen metabolism, glucose consumption, cell death, and vascular endothelial growth factor secretion. Experimental measurements are used to devise and parameterize a mathematical model that describes the behavior of the cells. The proposed model effectively describes the interactions between cells and their microenvironment and could in the future be extended, allowing researchers to compare the behavior of different therapeutic cells. Such a combinatorial approach could be used to accelerate the clinical translation of NRCs by identifying which critical design features should be optimized when fabricating engineered nerve repair conduits.

A mathematical investigation into the uptake kinetics of nanoparticles in vitro.

Nanoparticles have the potential to increase the efficacy of anticancer drugs whilst reducing off-target side effects. However, there remain uncertainties regarding the cellular uptake kinetics of nanoparticles which could have implications for nanoparticle design and delivery. Polymersomes are nanoparticle candidates for cancer therapy which encapsulate chemotherapy drugs. Here we develop a mathematical model to simulate the uptake of polymersomes via endocytosis, a process by which polymersomes bind to the cell surface before becoming internalised by the cell where they then break down, releasing their contents which could include chemotherapy drugs. We focus on two in vitro configurations relevant to the testing and development of cancer therapies: a well-mixed culture model and a tumour spheroid setup. Our mathematical model of the well-mixed culture model comprises a set of coupled ordinary differential equations for the unbound and bound polymersomes and associated binding dynamics. Using a singular perturbation analysis we identify an optimal number of ligands on the polymersome surface which maximises internalised polymersomes and thus intracellular chemotherapy drug concentration. In our mathematical model of the spheroid, a multiphase system of partial differential equations is developed to describe the spatial and temporal distribution of bound and unbound polymersomes via advection and diffusion, alongside oxygen, tumour growth, cell proliferation and viability. Consistent with experimental observations, the model predicts the evolution of oxygen gradients leading to a necrotic core. We investigate the impact of two different internalisation functions on spheroid growth, a constant and a bond dependent function. It was found that the constant function yields faster uptake and therefore chemotherapy delivery. We also show how various parameters, such as spheroid permeability, lead to travelling wave or steady-state solutions.

A four-compartment multiscale model of fluid and drug distribution in vascular tumours.

The subtle relationship between vascular network structure and mass transport is vital to predict and improve the efficacy of anticancer treatments. Here, mathematical homogenisation is used to derive a new multiscale continuum model of blood and chemotherapy transport in the vasculature and interstitium of a vascular tumour. This framework enables information at a range of vascular hierarchies to be fed into an effective description on the length scale of the tumour. The model behaviour is explored through a demonstrative case study of a simplified representation of a dorsal skinfold chamber, to examine the role of vascular network architecture in influencing fluid and drug perfusion on the length scale of the chamber. A single parameter, P, is identified that relates tumour-scale fluid perfusion to the permeability and density of the capillary bed. By fixing the topological and physiological properties of the arteriole and venule networks, an optimal value for P is identified, which maximises tumour fluid transport and is thus hypothesised to benefit chemotherapy delivery. We calculate the values for P for eight explicit network structures; in each case, vascular intervention by either decreasing the permeability or increasing the density of the capillary network would increase fluid perfusion through the cancerous tissue. Chemotherapeutic strategies are compared and indicate that single injection is consistently more successful compared with constant perfusion, and the model predicts optimal timing of a second dose. These results highlight the potential of computational modelling to elucidate the link between vascular architecture and fluid, drug distribution in tumours.

Modelling the transport of fluid through heterogeneous, whole tumours in silico.

Cancers exhibit spatially heterogeneous, unique vascular architectures across individual samples, cell-lines and patients. This inherently disorganised collection of leaky blood vessels contribute significantly to suboptimal treatment efficacy. Preclinical tools are urgently required which incorporate the inherent variability and heterogeneity of tumours to optimise and engineer anti-cancer therapies. In this study, we present a novel computational framework which incorporates whole, realistic tumours extracted ex vivo to efficiently simulate vascular blood flow and interstitial fluid transport in silico for validation against in vivo biomedical imaging. Our model couples Poiseuille and Darcy descriptions of vascular and interstitial flow, respectively, and incorporates spatially heterogeneous blood vessel lumen and interstitial permeabilities to generate accurate predictions of tumour fluid dynamics. Our platform enables highly-controlled experiments to be performed which provide insight into how tumour vascular heterogeneity contributes to tumour fluid transport. We detail the application of our framework to an orthotopic murine glioma (GL261) and a human colorectal carcinoma (LS147T), and perform sensitivity analysis to gain an understanding of the key biological mechanisms which determine tumour fluid transport. Finally we mimic vascular normalization by modifying parameters, such as vascular and interstitial permeabilities, and show that incorporating realistic vasculatures is key to modelling the contrasting fluid dynamic response between tumour samples. Contrary to literature, we show that reducing tumour interstitial fluid pressure is not essential to increase interstitial perfusion and that therapies should seek to develop an interstitial fluid pressure gradient. We also hypothesise that stabilising vessel diameters and permeabilities are not key responses following vascular normalization and that therapy may alter interstitial hydraulic conductivity. Consequently, we suggest that normalizing the interstitial microenvironment may provide a more effective means to increase interstitial perfusion within tumours.

Developing an In Vitro Model to Screen Drugs for Nerve Regeneration.

Peripheral nerve injuries (PNI) have a high prevalence and can be debilitating, resulting in life-long loss or disturbance in end-organ function, which compromises quality of life for patients. Current therapies use microsurgical approaches but there is the potential for enhancing recovery through other therapeutic modalities such as; cell-based conduits, gene therapy and small molecules. A number of molecular targets and drugs which have the potential to improve nerve regeneration have been identified, however, there are challenges associated with moving therapies toward clinical translation. Due to the lack of detailed knowledge about the pro-regenerative effect of potential drug treatments, there is a need for effective in vitro models to screen compounds to inform future pre-clinical and clinical studies. The interaction between regenerating neurites and supporting Schwann cells is a key feature of the nerve environment, therefore, in vitro models that mimic this cellular association are useful tools. In this study, we have investigated various cell culture models, including simple monolayer systems and more complex 3D-engineered co-cultures, as models for use in PNI drug development. Anat Rec, 301:1628-1637, 2018. © 2018 The Authors. The Anatomical Record published by Wiley Periodicals, Inc. on behalf of American Association of Anatomists.

The Mechanics of Brow-Suspension Ptosis Repair: A Comparative Study of Fox Pentagon and Crawford Triangle Techniques.

PURPOSE: To perform quantitative analysis of the most commonly used brow-suspension configurations. METHODS: The inflection positions for Fox pentagon and Crawford triangle configurations were marked on 49 healthy volunteers (male and female) and photographs taken in 3 states: "normal," "closed," and "raised." The skin marks were measured vectorially with respect to the medial canthus, and displacement changes were evaluated for "normal-to-closed" ("blinking") and from "closed-to-raised" ("eye-opening") states. The distance between a pair of inflection marks, representing the approximate path of sling configurations, were also measured and analyzed in relation to the mechanical properties of a variety of synthetic brow-suspension materials. RESULTS: "Blinking" resulted in the greatest displacement in the medial eyelid incision, resulting in the greatest strain on the line connecting the medial eyelid and medial brow inflections. No significant differences in the strains for individual lines were found between the Fox and Crawford techniques, although the former shows a significantly lower overall strain in the whole loop than the latter. The displacements of some inflections and of the strains of a few lines differed significantly in men and women. CONCLUSIONS: Within the scope of this study, the blinking action was shown to result in the maximum strain of ~40%, which lies within the elastic region of stress-strain curves for some commonly used synthetic brow-suspension materials. No one method was statistically superior, although the Fox pentagon gave a significantly lower overall strain when the sling material was assumed to move somewhat around the inflections within a closed loop.

Regulation of O2 consumption by the PI3K and mTOR pathways contributes to tumor hypoxia.

BACKGROUND: Inhibitors of the phosphatidylinositol 3-kinase (PI3K) and the mammalian target of rapamycin (mTOR) pathway are currently in clinical trials. In addition to antiproliferative and proapoptotic effects, these agents also diminish tumor hypoxia. Since hypoxia is a major cause of resistance to radiotherapy, we sought to understand how it is regulated by PI3K/mTOR inhibition. METHODS: Whole cell, mitochondrial, coupled and uncoupled oxygen consumption were measured in cancer cells after inhibition of PI3K (Class I) and mTOR by pharmacological means or by RNAi. Mitochondrial composition was assessed by immunoblotting. Hypoxia was measured in spheroids, in tumor xenografts and predicted with mathematical modeling. RESULTS: Inhibition of PI3K and mTOR reduced oxygen consumption by cancer cell lines is predominantly due to reduction of mitochondrial respiration coupled to ATP production. Hypoxia in tumor spheroids was reduced, but returned after removal of the drug. Murine tumors had increased oxygenation even in the absence of average perfusion changes or tumor necrosis. CONCLUSIONS: Targeting the PI3K/mTOR pathway substantially reduces mitochondrial oxygen consumption thereby reducing tumor hypoxia. These alterations in tumor hypoxia should be considered in the design of clinical trials using PI3K/mTOR inhibitors, particularly in conjunction with radiotherapy.

Microstructure and mechanical properties of synthetic brow-suspension materials.

Levator palpebrae superioris (LPS) is a muscle responsible for lifting the upper eyelid and its malfunction leads to a condition called "ptosis", resulting in disfigurement and visual impairment. Severe ptosis is generally treated with "brow-suspension" surgery, whereby the eyelid is cross-connected to the mobile tissues above the eyebrow using a cord-like material, either natural (e.g. fascia lata harvested from the patient) or a synthetic cord. Synthetic brow-suspension materials are widely used, due to not requiring the harvesting of fascia lata that can be associated with pain and donor-site complications. The mechanical properties of some commonly-used synthetic brow-suspension materials were investigated--namely, monofilament polypropylene (Prolene®), sheathed braided polyamide (Supramid Extra® II), silicone frontalis suspension rod (Visitec® Seiff frontalis suspension set), woven polyester (Mersilene® mesh), and expanded polytetrafluoroethylene (Ptose-Up). Each material underwent a single tensile loading to the failure of the material, at three different displacement rates (1, 750 and 1500 mm/min). All the materials exhibited elastic-plastic tensile stress-strain behaviour with considerable differences in elastic modulus, ultimate tensile strength, elastic limit and work of fracture. The results suggest that, as compared to other materials, the silicone brow-suspension rod (Visitec® SFSS) might be the most suitable, providing relatively long-lasting stability and desirable performance. These findings, together with other factors such as commercial availability, cost and clinical outcomes, will provide clinicians with a more rational basis for selection of brow-suspension materials.

Combined mathematical modelling and experimentation to predict polymersome uptake by oral cancer cells.

This study is motivated by understanding and controlling the key physical properties underlying internalisation of nano drug delivery. We consider the internalisation of specific nanometre size delivery vehicles, comprised of self-assembling amphiphilic block copolymers, called polymersomes that have the potential to specifically deliver anticancer therapeutics to tumour cells. The possible benefits of targeted polymersome drug delivery include reduced off-target toxic effects in healthy tissue and increased drug uptake by diseased tissue. Through a combination of in vitro experimentation and mathematical modelling, we develop a validated model of nanoparticle uptake by cells via the clathrin-mediated endocytotic pathway, incorporating receptor binding, clustering and recycling. The model predicts how the characteristics of receptor targeting, and the size and concentration of polymersomes alter uptake by tumour cells. The number of receptors per cell was identified as being the dominant mechanism accounting for the difference between cell types in polymersome uptake rate. FROM THE CLINICAL EDITOR: This article reports on a validated model developed through a combination of in vitro experimentation and mathematical modeling of nanoparticle uptake by cells via the clathrin-mediated endocytotic pathway. The model incorporates receptor binding, clustering, and recycling and predicts how the characteristics of receptor targeting, the size and concentration alter polymersome uptake by cancer cells.

Multiscale modelling of fluid and drug transport in vascular tumours.

A model for fluid and drug transport through the leaky neovasculature and porous interstitium of a solid tumour is developed. The transport problems are posed on a micro-scale characterized by the inter-capillary distance, and the method of multiple scales is used to derive the continuum equations describing fluid and drug transport on the length scale of the tumour (under the assumption of a spatially periodic microstructure). The fluid equations comprise a double porous medium, with coupled Darcy flow through the interstitium and vasculature, whereas the drug equations comprise advection-reaction equations; in each case the dependence of the transport coefficients on the vascular geometry is determined by solving micro-scale cell problems.

Multiscale modeling of fluid transport in tumors.

A model for fluid flow through the leaky neovasculature and porous interstitium of a solid tumor is developed. A network of isolated capillaries is analyzed in the limit of small capillary radius, and analytical expressions for the hydraulic conductivities and fractional leakage coefficients derived. This model is then homogenized to give a continuum description in terms of the vascular density. The resulting equations comprise a double porous medium with coupled Darcy flow through the interstitium and vasculature.