Modeling tumor transport and growth with poroelastic biopolymer networks

The mechanical properties of the extracellular matrix (ECM) regulate tumor growth and invasion in the tumor microenvironment. Models of biopolymer networks have been used to investigate the impact of elasticity and viscoelasticity of ECM on tumor behavior. Under tumor compression, these networks also show poroelastic behavior that is governed by the resistance to water flow through their pores. This work investigates the hypothesis that poroelastic properties regulate tumor growth. Here, alginate hydrogels with tunable ionic and hybrid ionic/covalent crosslinking are used as a model biopolymer system. Hydrogel stiffness, viscoelasticity, and stress relaxation behavior were characterized using stepwise axial compression. Among these properties, we find poroelastic fluid outflow dominates ECM stress relaxation, as the measured water flux was significantly affected under compression. Continuum mechanics-based modeling was developed to formulate and calculate the chemical potential gradients of water (solvent) in the hydrogels under compression. This framework was extended into an advection-diffusion framework to quantify growth factor (solute) distribution under varying strengths of stress and diffusion indexed by the relative strength of convective to diffusive transport, characterized by the Péclet number. An agent-based computational simulation showed that tumor growth was affected by Péclet number. Together, these results highlight the role of the poroelastic properties of ECM on water flux and transport in the tumor microenvironment.