Multiscale analysis of the coupled effects governing the movement of interstitial fluid in cortical bone

A review on computer modeling of bone piezoelectricity and its application to bone adaptation and regeneration

Bone is a hierarchical, multiphasic and anisotropic structure which in addition possess piezoelectric properties. The generation of piezoelectricity in bone is a complex process which has been shown to play a key role both in bone adaptation and regeneration. In order to understand the complex biological, mechanical and electrical interactions that take place during these processes, several computer models have been developed and used to test hypothesis on potential mechanisms behind experimental observations. This paper aims to review the available literature on computer modeling of bone piezoelectricity and its application to bone adaptation and healing. We first provide a brief overview of the fundamentals of piezoelectricity and bone piezoelectric effects. We then review how these properties have been used in computational models of bone adaptation and electromechanical behaviour of bone. In addition, in the last section, we summarize current limitations and potential directions for future work.

The Virtual Physiological Human: Ten Years After

Biomedical research and clinical practice are struggling to cope with the growing complexity that the progress of health care involves. The most challenging diseases, those with the largest socioeconomic impact (cardiovascular conditions; musculoskeletal conditions; cancer; metabolic, immunity, and neurodegenerative conditions), are all characterized by a complex genotype-phenotype interaction and by a "systemic" nature that poses a challenge to the traditional reductionist approach. In 2005 a small group of researchers discussed how the vision of computational physiology promoted by the Physiome Project could be translated into clinical practice and formally proposed the term Virtual Physiological Human. Our knowledge about these diseases is fragmentary, as it is associated with molecular and cellular processes on the one hand and with tissue and organ phenotype changes (related to clinical symptoms of disease conditions) on the other. The problem could be solved if we could capture all these fragments of knowledge into predictive models and then compose them into hypermodels that help us tame the complexity that such systemic behavior involves. In 2005 this was simply not possible-the necessary methods and technologies were not available. Now, 10 years later, it seems the right time to reflect on the original vision, the results achieved so far, and what remains to be done.

In silico bone mechanobiology: modeling a multifaceted biological system

Mechanobiology, the study of the influence of mechanical loads on biological processes through signaling to cells, is fundamental to the inherent ability of bone tissue to adapt its structure in response to mechanical stimulation. The immense contribution of computational modeling to the nascent field of bone mechanobiology is indisputable, having aided in the interpretation of experimental findings and identified new avenues of inquiry. Indeed, advances in computational modeling have spurred the development of this field, shedding new light on problems ranging from the mechanical response to loading by individual cells to tissue differentiation during events such as fracture healing. To date, in silico bone mechanobiology has generally taken a reductive approach in attempting to answer discrete biological research questions, with research in the field broadly separated into two streams: (1) mechanoregulation algorithms for predicting mechanobiological changes to bone tissue and (2) models investigating cell mechanobiology. Future models will likely take advantage of advances in computational power and techniques, allowing multiscale and multiphysics modeling to tie the many separate but related biological responses to loading together as part of a larger systems biology approach to shed further light on bone mechanobiology. Finally, although the ever‐increasing complexity of computational mechanobiology models will inevitably move the field toward patient‐specific models in the clinic, the determination of the context in which they can be used safely for clinical purpose will still require an extensive combination of computational and experimental techniques applied to in vitro and in vivo applications. WIREs Syst Biol Med 2016, 8:485–505. doi: 10.1002/wsbm.1356

For further resources related to this article, please visit the WIREs website.

Entropic electrokinetics: recirculation, particle separation, and negative mobility

We show that when particles are suspended in an electrolyte confined between corrugated charged surfaces, electrokinetic flows lead to a new set of phenomena such as particle separation, mixing for low-Reynolds micro- and nanometric devices, and negative mobility. Our analysis shows that such phenomena arise, for incompressible fluids, due to the interplay between the electrostatic double layer and the corrugated geometrical confinement and that they are magnified when the width of the channel is comparable to the Debye length. Our characterization allows us to understand the physical origin of such phenomena, therefore, shedding light on their possible relevance in a wide variety of situations ranging from nano- and microfluidic devices to biological systems.

Textural versus electrostatic exclusion-enrichment effects in the effective chemical transport within the cortical bone: a numerical investigation

Interstitial fluid within bone tissue is known to govern the remodelling signals' expression. Bone fluid flow is generated by skeleton deformation during the daily activities. Due to the presence of charged surfaces in the bone porous matrix, the electrochemical phenomena occurring in the vicinity of mechanosensitive bone cells, the osteocytes, are key elements in the cellular communication. In this study, a multiscale model of interstitial fluid transport within bone tissues is proposed. Based on an asymptotic homogenization method, our modelling takes into account the physicochemical properties of bone tissue. Thanks to this multiphysical approach, the transport of nutrients and waste between the blood vessels and the bone cells can be quantified to better understand the mechanotransduction of bone remodelling. In particular, it is shown that the electrochemical tortuosity may have stronger implications in the mass transport within the bone than the purely morphological one.

Multiphysical modelling of fluid transport through osteo-articular media

In this study, a multiphysical description of fluid transport through osteo-articular porous media is presented. Adapted from the model of Moyne and Murad, which is intended to describe clayey materials behaviour, this multiscale modelling allows for the derivation of the macroscopic response of the tissue from microscopical information. First the model is described. At the pore scale, electrohydrodynamics equations governing the electrolyte movement are coupled with local electrostatics (Gauss-Poisson equation), and ionic transport equations. Using a change of variables and an asymptotic expansion method, the macroscopic description is carried out. Results of this model are used to show the importance of couplings effects on the mechanotransduction of compact bone remodelling.

Modeling fluorescence recovery after photobleaching in loaded bone: potential applications in measuring fluid and solute transport in the osteocytic lacunar-canalicular system

Solute transport through the bone lacunar-canalicular system is essential for osteocyte viability and function, and it can be measured using fluorescence recovery after photobleaching (FRAP). The mathematical model developed here aims to analyze solute transport during FRAP in mechanically loaded bone. Combining both whole bone-level poroelasticity and cellular-level solute transport, we found that load-induced solute transport during FRAP is characterized by an exponential recovery rate, which is determined by the dimensionless Strouhal (St) number that characterizes the oscillation effects over the mean flows, and that significant transport occurs only for St values below a threshold, when the solute stroke displacement exceeds the distance between the source and sink (the canalicular length). This threshold mechanism explains the general flow behaviors such as increasing transport with increasing magnitude and decreasing frequency. Mechanical loading is predicted to enhance transport of all tracers relative to diffusion, with the greatest enhancement for medium-sized tracers and less enhancement for small and large tracers. This study provides guidelines for future FRAP experiments, based on which the model can be used to quantify bone permeability, solute-matrix interaction, and flow velocities. These studies should provide insights into bone adaptation and metabolism, and help to treat various bone diseases and conditions.

Study of the influence of fibrous pericellular matrix in the cortical interstitial fluid movement with hydroelectrochemical effects

Fluid flow within cortical bone tissue is modeled through an upscaling approach of a local description of the fluid movement. At the pore scale, the coupled phenomena (Poiseuille effect, osmosis, and electro-osmosis) governing the interstitial fluid movement are considered. Thus, actions of electro-osmotic and osmotic motions, in addition to the classical Poiseuille flow, are studied at the canaliculus scale by deriving a coupled Darcy law. The addition of a Brinkman-like term in this macroscopic result helps us to take into account the influence of the pericellular matrix on the coupled transport phenomena. At the canaliculus scale, the general trends that can be drawn from this study are as follows: (i) The presence of the fibrous matrix tends to reduce the fluid flow considerably; (ii) the role of osmotic and electro-osmotic effects is no longer negligible for dense fibrous media.