Solute transport within porous biofilms: diffusion or dispersion?

Many microorganisms live within surface-associated consortia, termed biofilms, that can form intricate porous structures interspersed with a network of fluid channels. In such systems, transport phenomena, including flow and advection, regulate various aspects of cell behaviour by controllling nutrient supply, evacuation of waste products and permeation of antimicrobial agents. This study presents multiscale analysis of solute transport in these porous biofilms. We start our analysis with a channel-scale description of mass transport and use the method of volume averaging to derive a set of homogenized equations at the biofilmscale. We show that solute transport may be described via two coupled partial differential equations for the averaged concentrations, or telegrapher’s equations. These models are particularly relevant for chemical species, such as some antimicrobial agents, that penetrate cell clusters very slowly. In most cases, especially for nutrients, solute penetration is faster, and transport can be described via an advection-dispersion equation. In this simpler case, the effective diffusion is characterised by a second-order tensor whose components depend on: (1) the topology of the channels’ network; (2) the solute’s diffusion coefficients in the fluid and the cell clusters; (3) hydrodynamic dispersion effects; and (4) an additional dispersion term intrinsic to the two-phase configuration. Although solute transport in biofilms is commonly thought to be diffusion-dominated, this analysis shows that dispersion effects may significantly contribute to transport

Chemical Power for Microscopic Robots in Capillaries

The power available to microscopic robots (nanorobots) that oxidize bloodstream glucose while aggregated in circumferential rings on capillary walls is evaluated with a numerical model using axial symmetry and time-averaged release of oxygen from passing red blood cells. Robots about one micron in size can produce up to several tens of picowatts, in steady-state, if they fully use oxygen reaching their surface from the blood plasma. Robots with pumps and tanks for onboard oxygen storage could collect oxygen to support burst power demands two to three orders of magnitude larger. We evaluate effects of oxygen depletion and local heating on surrounding tissue. These results give the power constraints when robots rely entirely on ambient available oxygen and identify aspects of the robot design significantly affecting available power. More generally, our numerical model provides an approach to evaluating robot design choices for nanomedicine treatments in and near capillaries.Comment: 28 pages, 7 figure

The role of the microvascular network structure on diffusion and consumption of anticancer drugs

We investigate the impact of microvascular geometry on the transport of drugs in solid tumors, focusing on the diffusion and consumption phenomena. We embrace recent advances in the asymptotic homogenization literature starting from a double Darcy—double advection-diffusion-reaction system of partial differential equations that is obtained exploiting the sharp length separation between the intercapillary distance and the average tumor size. The geometric information on the microvascular network is encoded into effective hydraulic conductivities and diffusivities, which are numerically computed by solving periodic cell problems on appropriate microscale representative cells. The coefficients are then injected into the macroscale equations, and these are solved for an isolated, vascularized spherical tumor. We consider the effect of vascular tortuosity on the transport of anticancer molecules, focusing on Vinblastine and Doxorubicin dynamics, which are considered as a tracer and as a highly interacting molecule, respectively. The computational model is able to quantify the treatment performance through the analysis of the interstitial drug concentration and the quantity of drug metabolized in the tumor. Our results show that both drug advection and diffusion are dramatically impaired by increasing geometrical complexity of the microvasculature, leading to nonoptimal absorption and delivery of therapeutic agents. However, this effect apparently has a minor role whenever the dynamics are mostly driven by metabolic reactions in the tumor interstitium, eg, for highly interacting molecules. In the latter case, anticancer therapies that aim at regularizing the microvasculature might not play a major role, and different strategies are to be developed

Distributed Control of Microscopic Robots in Biomedical Applications

Current developments in molecular electronics, motors and chemical sensors could enable constructing large numbers of devices able to sense, compute and act in micron-scale environments. Such microscopic machines, of sizes comparable to bacteria, could simultaneously monitor entire populations of cells individually in vivo. This paper reviews plausible capabilities for microscopic robots and the physical constraints due to operation in fluids at low Reynolds number, diffusion-limited sensing and thermal noise from Brownian motion. Simple distributed controls are then presented in the context of prototypical biomedical tasks, which require control decisions on millisecond time scales. The resulting behaviors illustrate trade-offs among speed, accuracy and resource use. A specific example is monitoring for patterns of chemicals in a flowing fluid released at chemically distinctive sites. Information collected from a large number of such devices allows estimating properties of cell-sized chemical sources in a macroscopic volume. The microscopic devices moving with the fluid flow in small blood vessels can detect chemicals released by tissues in response to localized injury or infection. We find the devices can readily discriminate a single cell-sized chemical source from the background chemical concentration, providing high-resolution sensing in both time and space. By contrast, such a source would be difficult to distinguish from background when diluted throughout the blood volume as obtained with a blood sample

Multiphase modelling of vascular tumour growth in two spatial dimensions

In this paper we present a continuum mathematical model of vascular tumour growth which is based on a multiphase framework in which the tissue is decomposed into four distinct phases and the principles of conservation of mass and momentum are applied to the normal/healthy cells, tumour cells, blood vessels and extracellular material. The inclusion of a diffusible nutrient, supplied by the blood vessels, allows the vasculature to have a nonlocal influence on the other phases. Two-dimensional computational simulations are carried out on unstructured, triangular meshes to allow a natural treatment of irregular geometries, and the tumour boundary is captured as a diffuse interface on this mesh, thereby obviating the need to explicitly track the (potentially highly irregular and ill-defined) tumour boundary. A hybrid finite volume/finite element algorithm is used to discretise the continuum model: the application of a conservative, upwind, finite volume scheme to the hyperbolic mass balance equations and a finite element scheme with a stable element pair to the generalised Stokes equations derived from momentum balance, leads to a robust algorithm which does not use any form of artificial stabilisation. The use of a matrix-free Newton iteration with a finite element scheme for the nutrient reaction-diffusion equations allows full nonlinearity in the source terms of the mathematical model. Numerical simulations reveal that this four-phase model reproduces the characteristic pattern of tumour growth in which a necrotic core forms behind an expanding rim of well-vascularised proliferating tumour cells. The simulations consistently predict linear tumour growth rates. The dependence of both the speed with which the tumour grows and the irregularity of the invading tumour front on the model parameters are investigated

Brownian motion near an elastic cell membrane: A theoretical study

Elastic confinements are an important component of many biological systems and dictate the transport properties of suspended particles under flow. In this chapter, we review the Brownian motion of a particle moving in the vicinity of a living cell whose membrane is endowed with a resistance towards shear and bending. The analytical calculations proceed through the computation of the frequency-dependent mobility functions and the application of the fluctuation-dissipation theorem. Elastic interfaces endow the system with memory effects that lead to a long-lived anomalous subdiffusive regime of nearby particles. In the steady limit, the diffusional behavior approaches that near a no-slip hard wall. The analytical predictions are validated and supplemented with boundary-integral simulations.Comment: 16 pages, 7 figures and 161 references. Contributed chapter to the flowing matter boo