Computational Multi-Scale Modeling of Drug Delivery into an Anti-Angiogenic Therapy-Treated Tumor

The present study develops a numerical model, which is the most complex one, in comparison to previous research to investigate drug delivery accompanied by the anti-angiogenesis effect. This paper simulates intravascular blood flow and interstitial fluid flow using a dynamic model. The model accounts for the non-Newtonian behavior of blood and incorporates the adaptation of the diameter of a heterogeneous microvascular network derived from modeling the evolution of endothelial cells toward a circular tumor sprouting from two-parent vessels, with and without imposing the inhibitory effect of angiostatin on a modified discrete angiogenesis model. The average solute exposure and its uniformity in solid tumors of different sizes are studied by numerically solving the convection-diffusion equation. Three different methodologies are considered for simulating anti-angiogenesis: modifying the capillary network, updating the transport properties, and considering both microvasculature and transport properties modifications. It is shown that anti-angiogenic therapy decreases drug wash-out in the periphery of the tumor. Results show the decisive role of microvascular structure, particularly its distribution, and interstitial transport properties modifications induced via vascular normalization on the quality of drug delivery, such that it is improved by 39% in uniformity by the second approach in R = 0.2 cm.

Myeloid Derived Suppressor Cells Migrate in Response to Flow and Lymphatic Endothelial Cell Interaction in the Breast Tumor Microenvironment

At the site of the tumor, myeloid derived suppressor cells (MDSCs) infiltrate and interact with elements of the tumor microenvironment in complex ways. Within the invading tumor, MDSCs are exposed to interstitial fluid flow (IFF) that exists within the chronic inflammatory tumor microenvironment at the tumor-lymphatic interface. As drivers of cell migration and invasion, the link between interstitial fluid flow, lymphatics, and MDSCs have not been clearly established. Here, we hypothesized that interstitial fluid flow and cells within the breast tumor microenvironment modulate migration of MDSCs. We developed a novel 3D model to mimic the breast tumor microenvironment and incorporated MDSCs harvested from 4T1-tumor bearing mice. Using live imaging, we found that sorted GR1+ splenocytes had reduced chemotactic index compared to the unsorted population, but their speed and displacement were similar. Using our adapted tissue culture insert assay, we show that interstitial fluid flow promotes MDSC invasion, regardless of absence or presence of tumor cells. Coordinating with lymphatic endothelial cells, interstitial fluid flow further enhanced invasion of MDSCs in the presence of 4T1 cells. We also show that VEGFR3 inhibition reduced both MDSC and 4T1 flow response. Together, these findings indicate a key role of interstitial fluid flow in MDSC migration as well as describe a tool to explore the immune microenvironment in breast cancer.

Study on the poroelastic behaviors of the defected articular cartilage

This article presented the possible mechanism of arthritis damaged changes in cartilage's interstitial fluid flowing behavior. Firstly, the analytical solutions for the pore fluid pressure and velocity in the idealized cartilage defect model were obtained, which are employed to validate the finite element (FE) method. Then according to the MRI data, an articular cartilage FE model was developed to study the effects of defect characteristics on its poroelastic behaviors. The results showed the interstitial fluid pressure and velocity in defected articular cartilage is diminished, moreover, this trend is even more severe as the defect radius or thickness increased. As the development of osteoarthritis goes, the fluid velocity is decreased and cause the even serious nutrients loss.

Computational Simulation of Exosome Transport in Tumor Microenvironment

Cellular exosome-mediated crosstalk in tumor microenvironment (TME) is a critical component of anti-tumor immune responses. In addition to particle size, exosome transport and uptake by target cells is influenced by physical and physiological factors, including interstitial fluid pressure, and exosome concentration. These variables differ under both normal and pathological conditions, including cancer. The transport of exosomes in TME is governed by interstitial flow and diffusion. Based on these determinants, mathematical models were adapted to simulate the transport of exosomes in the TME with specified exosome release rates from the tumor cells. In this study, the significance of spatial relationship in exosome-mediated intercellular communication was established by treating their movement in the TME as a continuum using a transport equation, with advection due to interstitial flow and diffusion due to concentration gradients. To quantify the rate of release of exosomes by biomechanical forces acting on the tumor cells, we used a transwell platform with confluent triple negative breast cancer cells 4T1.2 seeded in BioFlex plates exposed to an oscillatory force. Exosome release rates were quantified from 4T1.2 cells seeded at the bottom of the well following the application of either no force or an oscillatory force, and these rates were used to model exosome transport in the transwell. The simulations predicted that a larger number of exosomes reached the membrane of the transwell for 4T1.2 cells exposed to the oscillatory force when compared to controls. Additionally, we simulated the interstitial fluid flow and exosome transport in a 2-dimensional TME with macrophages, T cells, and mixtures of these two populations at two different stages of a tumor growth. Computational simulations were carried out using the commercial computational simulation package, ANSYS/Fluent. The results of this study indicated higher exosome concentrations and larger interstitial fluid pressure at the later stages of the tumor growth. Quantifying the release of exosomes by cancer cells, their transport through the TME, and their concentration in TME will afford a deeper understanding of the mechanisms of these interactions and aid in deriving predictive models for therapeutic intervention.

Limits to flow detection in phase contrast MRI

Pulsed gradient spin echo (PGSE) complex signal behavior becomes dominated by attenuation rather than oscillation when displacements due to flow are similar or less than diffusive displacements. In this "slow-flow" regime, the optimal displacement encoding parameter q for phase contrast velocimetry depends on the diffusive length scale q s l o w = 1 / l D = 1 / 2 D Δ rather than the velocity encoding parameter v _enc = π/(qΔ). The minimum detectable mean velocity using the difference between the phase at +q _slow and -q _slow is 〈 v m i n 〉 = 1 / SNR D / Δ . These theories are then validated and applied to MRI by performing PGSE echo planar imaging experiments on water flowing through a column with a bulk region and a beadpack region at controlled flow rates. Velocities as slow as 6 μm/s are detected with velocimetry. Theories, MRI experimental protocols, and validation on a controlled phantom help to bridge the gap between porous media NMR and pre-clinical phase contrast and diffusion MRI.

Investigating Low-Velocity Fluid Flow in Tumors with Convection-MRI

Several distinct fluid flow phenomena occur in solid tumors, including intravascular blood flow and interstitial convection. Interstitial fluid pressure is often raised in solid tumors, which can limit drug delivery. To probe low-velocity flow in tumors resulting from raised interstitial fluid pressure, we developed a novel MRI technique named convection-MRI, which uses a phase-contrast acquisition with a dual-inversion vascular nulling preparation to separate intra- and extravascular flow. Here, we report the results of experiments in flow phantoms, numerical simulations, and tumor xenograft models to investigate the technical feasibility of convection-MRI. We observed a significant correlation between estimates of effective fluid pressure from convection-MRI with gold-standard, invasive measurements of interstitial fluid pressure in mouse models of human colorectal carcinoma. Our results show how convection-MRI can provide insights into the growth and responsiveness to vascular-targeting therapy in colorectal cancers.Significance: A noninvasive method for measuring low-velocity fluid flow caused by raised fluid pressure can be used to assess changes caused by therapy. Cancer Res; 78(7); 1859-72. ©2018 AACR.

[Effect of artery pulse on the osteonal interstitial fluid flow behavior]

There are two main types of fluid in bone tissue: blood and interstitial fluid. The metabolism of cells mainly relies on the microenvironment of the interstitial fluid. Researches of osteonal fluid seepage behavior based on the microstructure of bone tissue have become a hot point. The aim of the present research work is to assess the effect of blood pressure oscillation on the osteonal interstitial fluid seepage behavior. We established finite element osteon models for a hollow and that considering blood pressure oscillation, respectively, with COMSOL Multiphysics software in order to compare their fluid flow behavior under the axial loading. The results predicted that the interstitial fluid pressure field was enlarged considering the blood pressure oscillation, while the velocity filed changed little. Specifically, the increase of blood pressure oscillatory amplitude could result in the increase of osteonal interstitial fluid pressure, while the blood pressure oscillatory frequency had limited effects on the osteonal pore fluid pressure. Moreover, the blood pressure oscillatory amplitude and frequency had no effect on the osteonal interstitial fluid velocity. The finite element model can be used for the study of the poroelastic behaviors of the osteon under non-axisymmetric loads and microcracks, and can also be a new way to study the mechanism of bone mechanotransduction and electromechanotransduction.

Micromotion-induced peri-prosthetic fluid flow around a cementless femoral stem

Micromotion-induced interstitial fluid flow at the bone-implant interface has been proposed to play an important role in aseptic loosening of cementless implants. High fluid velocities are thought to promote aseptic loosening through activation of osteoclasts, shear stress induced control of mesenchymal stem cells differentiation, or transport of molecules. In this study, our objectives were to characterize and quantify micromotion-induced fluid flow around a cementless femoral stem using finite element modeling. With a 2D model of the bone-implant interface and full-factorial design, we first evaluated the relative influence of material properties, and bone-implant micromotion and gap on fluid velocity. Transverse sections around a femoral stem were built from computed tomography images, while boundary conditions were obtained from experimental measurements on the same femur. In a second step, a 3D model was built from the same data-set to estimate the shear stress experienced by cells hosted in the peri-implant tissues. The full-factorial design analysis showed that local micromotion had the most influence on peak fluid velocity at the interface. Remarkable variations in fluid velocity were observed in the macrostructures at the surface of the implant in the 2D transverse sections of the stem. The 3D model predicted peak fluid velocities extending up to 2.2 mm/s in the granulation tissue and to 3.9 mm/s in the trabecular bone. Peak shear stresses on the cells hosted in these tissues ranged from 0.1 to 12.5 Pa. These results offer insight into mechanical stimuli encountered at the bone-implant interface.

Biomechanical aspects of the muscle-bone interaction

There is growing interest in the interaction between skeletal muscle and bone, particularly at the genetic and molecular levels. However, the genetic and molecular linkages between muscle and bone are achieved only within the context of the essential mechanical coupling of the tissues. This biomechanical and physiological linkage is readily evident as muscles attach to bone and induce exposure to varied mechanical stimuli via functional activity. The responsiveness of bone cells to mechanical stimuli, or their absence, is well established. However, questions remain regarding how muscle forces applied to bone serve to modulate bone homeostasis and adaptation. Similarly, the contributions of varied, but unique, stimuli generated by muscle to bone (such as low-magnitude, high-frequency stimuli) remains to be established. The current article focuses upon the mechanical relationship between muscle and bone. In doing so, we explore the stimuli that muscle imparts upon bone, models that enable investigation of this relationship, and recent data generated by these models.

Quantification of Lacunar-Canalicular Interstitial Fluid Flow Through Computational Modeling of Fluorescence Recovery After Photobleaching

Skeletal adaptation to mechanical loading has been widely hypothesized to involve the stimulation of osteocytes by interstitial fluid flow (IFF). However, direct investigation of this hypothesis has been difficult due in large part to the inability to directly measure IFF velocities within the lacunar-canalicular system. Measurements of fluorescence recovery after photobleaching (FRAP) within individual lacunae could be used to quantify lacunar-canalicular IFF when combined with mathematical modeling. In this study, we used a computational transport model to characterize the relationship between flow frequency (0.5-10 Hz), peak flow velocity (0-300 μm/s), tracer diffusion coefficient (100-300 μm(2)/s), and transport enhancement (i.e., (k/k(0)) - 1, where k and k(0) are the transport rates in the presence/absence of flow) during lacunar FRAP investigations. We show that this relationship is well described by a simple power law with frequency-dependent coefficients, and is relatively insensitive to variations in lacunar geometry. Using this power law relationship, we estimated peak IFF velocities in hindlimb mice subjected to intramedullary pressurization using values of k and k(0) previously obtained from ex vivo lacunar FRAP investigations. Together, our findings suggest that skeletal adaptation in hindlimb suspended mice subjected to dynamic intramedullary pressure occurred in the presence of IFF at levels associated with physiological loading.

Microfluidic enhancement of intramedullary pressure increases interstitial fluid flow and inhibits bone loss in hindlimb suspended mice

Interstitial fluid flow (IFF) has been widely hypothesized to mediate skeletal adaptation to mechanical loading. Although a large body of in vitro evidence has demonstrated that fluid flow stimulates osteogenic and antiresorptive responses in bone cells, there is much less in vivo evidence that IFF mediates loading-induced skeletal adaptation. This is due in large part to the challenges associated with decoupling IFF from matrix strain. In this study we describe a novel microfluidic system for generating dynamic intramedullary pressure (ImP) and IFF within the femurs of alert mice. By quantifying fluorescence recovery after photobleaching (FRAP) within individual lacunae, we show that microfluidic generation of dynamic ImP significantly increases IFF within the lacunocanalicular system. In addition, we demonstrate that dynamic pressure loading of the intramedullary compartment for 3 minutes per day significantly eliminates losses in trabecular and cortical bone mineral density in hindlimb suspended mice, enhances trabecular and cortical structural integrity, and increases endosteal bone formation rate. Unlike previously developed modalities for enhancing IFF in vivo, this is the first model that allows direct and dynamic modulation of ImP and skeletal IFF within mice. Given the large number of genetic tools for manipulating the mouse genome, this model is expected to serve as a powerful investigative tool in elucidating the role of IFF in skeletal adaptation to mechanical loading and molecular mechanisms mediating this process.