Liver fibrosis emulation: Impact of the vascular fibrotic alterations on hemodynamics

The liver circulatory system comprises two blood supply vascular trees (the hepatic artery and portal venous networks), microcirculation through the hepatic capillaries (the sinusoids), and a blood drainage vascular tree (the hepatic vein network). Vasculature changes due to fibrosis -located predominantly at the microcirculation level- lead to a marked increase in resistance to flow causing an increase in portal pressure (portal hypertension). Here, we present a liver fibrosis/cirrhosis model. We build on our 1D model of the healthy hepatic circulation, which considers the elasticity of the vessels walls and the pulsatile character of blood flow and pressure, and recreate the deteriorated liver vasculature due to fibrosis. We emulate altered sinusoids by fibrous tissue (stiffened, compressed and splitting) and propose boundary conditions to investigate the impact of fibrosis on hemodynamic variables within the organ. We obtain that the sinusoids stiffness leads to changes in the amplitude and shape of the blood flow and pressure waveforms but not in their mean value. For the compressed and splitting sinusoids, we observe significant increases in the mean value and amplitude of the pressure waveform in the altered sinusoids and in the portal venous network. In other words, we obtain the portal hypertension clinically observed in fibrotic/cirrhotic patients. We also study the extent of the spreading fibrosis by performing the structural fibrotic changes in an increasingly number of sinusoids. Finally, we calculate the portal pressure gradient (PPG) in the model and obtain values in agreement with those reported in the literature for fibrotic/cirrhotic patients.

1D-model of the human liver circulatory system

BACKGROUND AND OBJECTIVE

Blood flow rate and pressure can be measured in vivo by invasive and non-invasive techniques in the large vessels of the hepatic vasculature, but it is not possible to do so along the entire liver circulatory system. Here, we develop a novel 1D model of the liver circulatory system to obtain the hemodynamic signals from macrocirculation to microcirculation with a very low computational cost.

METHODS

The model considers structurally well-defined elements that constitute the entire hepatic circulatory system, the hemodynamics (the temporal-dependence of the blood flow rate and pressure), and the elasticity of the vessel walls.

RESULTS

Using flow rate signals from in vivo measurements as inputs in the model, we obtain pressure signals within their physiological range of values. Furthermore, the model allows to get and analyze the blood flow rate and pressure signals along any vessel of the hepatic vasculature. The impact of the elasticity of the different model components on the inlet pressures is also tested.

CONCLUSIONS

A 1D model of the entire blood vascular system of the human liver is presented for the first time. The model allows to obtain the hemodynamic signals along the hepatic vasculature at a low computational cost. The amplitude and shape of the flow and pressure signals has hardly been studied in the small liver vessels. In this sense, the proposed model is a useful non-invasive exploration tool of the characteristics of the hemodynamic signals. In contrast to models that partially address the hepatic vasculature or those using an electrical analogy, the model presented here is made entirely of structurally well-defined elements. Future works will allow to directly emulate structural vascular alterations due to hepatic diseases and studying their impact on pressure and blood flow signals at key locations of the vasculature.

Hierarchical Modeling of the Liver Vascular System

The liver plays a key role in the metabolic homeostasis of the whole organism. To carry out its functions, it is endowed with a peculiar circulatory system, made of three main dendritic flow structures and lobules. Understanding the vascular anatomy of the liver is clinically relevant since various liver pathologies are related to vascular disorders. Here, we develop a novel liver circulation model with a deterministic architecture based on the constructal law of design over the entire scale range (from macrocirculation to microcirculation). In this framework, the liver vascular structure is a combination of superimposed tree-shaped networks and porous system, where the main geometrical features of the dendritic fluid networks and the permeability of the porous medium, are defined from the constructal viewpoint. With this model, we are able to emulate physiological scenarios and to predict changes in blood pressure and flow rates throughout the hepatic vasculature due to resection or thrombosis in certain portions of the organ, simulated as deliberate blockages in the blood supply to these sections. This work sheds light on the critical impact of the vascular network on mechanics-related processes occurring in hepatic diseases, healing and regeneration that involve blood flow redistribution and are at the core of liver resilience.

Obstructions in Vascular Networks: Relation Between Network Morphology and Blood Supply

We relate vascular network structure to hemodynamics after vessel obstructions. We consider tree-like networks with a viscoelastic fluid with the rheological characteristics of blood. We analyze the network hemodynamic response, which is a function of the frequencies involved in the driving, and a measurement of the resistance to flow. This response function allows the study of the hemodynamics of the system, without the knowledge of a particular pressure gradient. We find analytical expressions for the network response, which explicitly show the roles played by the network structure, the degree of obstruction, and the geometrical place in which obstructions occur. Notably, we find that the sequence of resistances of the network without occlusions strongly determines the tendencies that the response function has with the anatomical place where obstructions are located. We identify anatomical sites in a network that are critical for its overall capacity to supply blood to a tissue after obstructions. We demonstrate that relatively small obstructions in such critical sites are able to cause a much larger decrease on flow than larger obstructions placed in non-critical sites. Our results indicate that, to a large extent, the response of the network is determined locally. That is, it depends on the structure that the vasculature has around the place where occlusions are found. This result is manifest in a network that follows Murray’s law, which is in reasonable agreement with several mammalian vasculatures. For this one, occlusions in early generation vessels have a radically different effect than occlusions in late generation vessels occluding the same percentage of area available to flow. This locality implies that whenever there is a tissue irrigated by a tree-like in vivo vasculature, our model is able to interpret how important obstructions are for the irrigation of such tissue.