The Computed Sinusoid

Abstract

Hepatic sinusoids are lined with thin endothelial cells with transcellular pores, termed fenestrations. These fenestrations are open channels that connect the sinusoidal lumen to the underlying Space of Disse (SoD) and the hepatocytes of the liver parenchyma. Fenestrations range from 0.05 to 0.35 µm in diameter and cover 5–15% of the sinusoidal endothelial surface area, depending on their location along the sinusoids. The direct measurement of hemodynamic parameters, such as pressure and flow velocity, remains challenging within the narrow sinusoids. Such knowledge would increase our understanding of the physiology of the hepatic niche and possible implications in aging or diseases in which fenestrations are reduced or lost. Few simulations of liver blood flow focus on the level of the individual sinusoid, and fewer still include the transcellular pores (fenestrations) of the sinusoidal endothelium. Furthermore, none have included (i) a porosity gradient along the sinusoid wall, modeled using through-all pores rather than a porous medium, (ii) the presence of the SoD, or (iii) lymphatic drainage. Herein, computed fluid dynamics (CFD) simulations were performed using a numerical model with relevant anatomical characteristics (length, diameter, porosity, inlet/outlet pressure, and lymphatic outflow from the portal region of the SoD). The greatest contribution to luminal velocity magnitude and pressure was the overall shape of the vessel. Divergent-radius models yielded velocity magnitudes 1.5–2 times higher than constant-radius models, and pressures were 5–8% lower in the divergent-radius models compared to the constant-radius models. Porosity only modestly contributed to luminal pressure. The luminal velocity magnitude was largely unaffected by the presence or absence of lymphatic drainage. Velocity magnitudes through fenestrations were lower in higher-porosity models (20%) vs. lower-porosity models (5%) across all models (0.4–0.55-fold lower). Velocity magnitudes through the space of Disse were increased 3–4 times via the addition of lymphatic drainage to the models, while pressures were decreased by 6–12%. The flow velocity in the SoD was modified via differences in porosity, while the flow velocity in the lumens of the sinusoids was largely unaffected. The overall shape of the vessel is the single most important factor in the pressure flow behavior of the sinusoidal lumen. The flow rate over hepatocytes and the SoD is modestly affected by the distribution of porosity along the sinusoid and greatly affected by the lymphatic drainage, parameters that would be of interest for modeling the exchange of blood with the hepatic parenchyma.