Discrete-to-continuum models of pre-stressed cytoskeletal filament networks

We introduce a mathematical model for the mechanical behaviour of the eukaryotic cell cytoskeleton. This discrete model involves a regular array of pre-stressed protein filaments that exhibit resistance to enthalpic stretching, joined at crosslinks to form a network. Assuming that the inter-crosslink distance is much shorter than the lengthscale of the cell, we upscale the discrete force balance to form a continuum system of governing equations and deduce the corresponding macroscopic stress tensor. We use these discrete and continuum models to analyse the imposed displacement of a bead placed in the domain, characterising the cell rheology through the force-displacement curve. We further derive an analytical approximation to the stress and strain fields in the limit of small bead radius, predicting the net force required to generate a given deformation and elucidating the dependency on the microscale properties of the filaments. We apply these models to networks of the intermediate filament vimentin and demonstrate good agreement between predictions of the discrete, continuum and analytical approaches. In particular, our model predicts that the network stiffness increases sublinearly with the filament pre-stress and scales logarithmically with the bead size.Comment: 38 pages, 12 figure

Optimising dead-end cake filtration using poroelasticity theory

Understanding the operation of filters used to remove particulates from fluids is important in many practical industries. Typically the particles are larger than the pores in the filter so a cake layer of particles forms on the filter surface. Here we extend existing models for filter blocking to account for deformation of the filter material and the cake layer due to the applied pressure that drives the fluid. These deformations change the permeability of the filter and the cake and hence the flow. We develop a new theory of compressible-cake filtration based on a simple poroelastic model in which we assume that the permeability depends linearly on local deformation. This assumption allows us to derive an explicit filtration law. The model predicts the possible shutdown of the filter when the imposed pressure difference is sufficiently large to reduce the permeability at some point to zero. The theory is applied to industrially relevant operating conditions, namely constant flux, maximising flux and constant pressure drop. Under these conditions, further analytical results are obtained, which yield predictions for optimal filter design with respect to given properties of the filter materials and the particles

Enhanced perfusion following exposure to radiotherapy: a theoretical investigation

Tumour angiogenesis leads to the formation of blood vessels that are structurally and spatially heterogeneous. Poor blood perfusion, in conjunction with increased hypoxia and oxygen heterogeneity, impairs a tumour’s response to radiotherapy. The optimal strategy for enhancing tumour perfusion remains unclear, preventing its regular deployment in combination therapies. In this work, we first identify vascular architectural features that correlate with enhanced perfusion following radiotherapy, using in vivo imaging data from vascular tumours. Then, we present a novel computational model to determine the relationship between these architectural features and blood perfusion in silico. If perfusion is defined to be the proportion of vessels that support blood flow, we find that vascular networks with small mean diameters and large numbers of angiogenic sprouts show the largest increases in perfusion post-irradiation for both biological and synthetic tumours. We also identify cases where perfusion increases due to the pruning of hypoperfused vessels, rather than blood being rerouted. These results indicate the importance of considering network composition when determining the optimal irradiation strategy. In the future, we aim to use our findings to identify tumours that are good candidates for perfusion enhancement and to improve the efficacy of combination therapies

Enhanced perfusion following exposure to radiotherapy: a theoretical investigation

Tumour angiogenesis leads to the formation of blood vessels that are structurally and spatially heterogeneous. Poor blood perfusion, in conjunction with increased hypoxia and oxygen heterogeneity, impairs a tumour’s response to radiotherapy. The optimal strategy for enhancing tumour perfusion remains unclear, preventing its regular deployment in combination therapies. In this work, we first identify vascular architectural features that correlate with enhanced perfusion following radiotherapy, using in vivo imaging data from vascular tumours. Then, we present a novel computational model to determine the relationship between these architectural features and blood perfusion in silico. If perfusion is defined to be the proportion of vessels that support blood flow, we find that vascular networks with small mean diameters and large numbers of angiogenic sprouts show the largest increases in perfusion post-irradiation for both biological and synthetic tumours. We also identify cases where perfusion increases due to the pruning of hypoperfused vessels, rather than blood being rerouted. These results indicate the importance of considering network composition when determining the optimal irradiation strategy. In the future, we aim to use our findings to identify tumours that are good candidates for perfusion enhancement and to improve the efficacy of combination therapies

Abnormal morphology biases haematocrit distribution in tumour vasculature and contributes to heterogeneity in tissue oxygenation

Oxygen heterogeneity in solid tumors is recognized as a limiting factor for therapeutic efficacy. This heterogeneity arises from the abnormal vascular structure of the tumor, but the precise mechanisms linking abnormal structure and compromised oxygen transport are only partially understood. In this paper, we investigate the role that red blood cell (RBC) transport plays in establishing oxygen heterogeneity in tumor tissue. We focus on heterogeneity driven by network effects, which are challenging to observe experimentally due to the reduced fields of view typically considered. Motivated by our findings of abnormal vascular patterns linked to deviations from current RBC transport theory, we calculated average vessel lengths L⎯⎯ and diameters d⎯⎯ from tumor allografts of three cancer cell lines and observed a substantial reduction in the ratio λ=L⎯⎯/d⎯⎯ compared to physiological conditions. Mathematical modeling reveals that small values of the ratio λ (i.e., λ<6 ) can bias hematocrit distribution in tumor vascular networks and drive heterogeneous oxygenation of tumor tissue. Finally, we show an increase in the value of λ in tumor vascular networks following treatment with the antiangiogenic cancer agent DC101. Based on our findings, we propose λ as an effective way of monitoring the efficacy of antiangiogenic agents and as a proxy measure of perfusion and oxygenation in tumor tissue undergoing antiangiogenic treatment

Dependence of cell-free-layer width on rheological parameters: Combining empirical data on flow separation at microvascular bifurcations with geometrical considerations

When blood flows through vessel networks, red blood cells (RBCs) are typically concentrated close to the vessel center line, creating a lubrication layer near the vessel wall. As RBCs bind oxygen, the width of this cell-free layer (CFL) impacts not only the blood rheology inside the vasculature, but also oxygen delivery to the tissues they perfuse and, hence, their function. Existing attempts to relate the width of the CFL to the rheological properties of the blood and the geometrical properties of the vessel are based on an analysis of the forces acting on RBCs suspended in the blood. However, the complexity of interactions in the blood makes this a challenging task. Here, we propose an alternative, two-step approach to derive such a functional relationship. First, we extend widely accepted empirical fits describing the minimum flow fraction needed for RBCs to enter a daughter vessel downstream of a microvascular bifurcation so that it depends not only on the diameter and discharge haematocrit of the parent vessel, but also on its average shear rate. Second, we propose a simple geometrical model for the minimum flow fraction based on the cross-sectional blood flow profile within the parent vessel upstream of the bifurcation—considering uniform, parabolic, and blunt velocity profiles—and derive the leading-order approximation to this model for small CFL widths. By equating the functional relationships obtained using these two approaches, we derive expressions relating the CFL width to the vessel diameter, discharge haematocrit, and mean shear rate. The resulting expressions are in good agreement with available in vivo data and represent a promising basis for future research

Enhanced perfusion following exposure to radiotherapy: a theoretical investigation (revised manuscript data)

<p>Dataset and software supporting the revised manuscript: "Enhanced perfusion following exposure to radiotherapy: a theoretical investigation." See the enclosed README and manuscript for further information.</p&gt