Simulating Microdosimetry in a Virtual Hepatic Lobule

The liver plays a key role in removing harmful chemicals from the body and is therefore often the first tissue to suffer potentially adverse consequences. To protect public health it is necessary to quantitatively estimate the risk of long-term low dose exposure to environmental pollutants. Animal testing is the primary tool for extrapolating human risk but it is fraught with uncertainty, necessitating novel alternative approaches. Our goal is to integrate in vitro liver experiments with agent-based cellular models to simulate a spatially extended hepatic lobule. Here we describe a graphical model of the sinusoidal network that efficiently simulates portal to centrilobular mass transfer in the hepatic lobule. We analyzed the effects of vascular topology and metabolism on the cell-level distribution following oral exposure to chemicals. The spatial distribution of metabolically inactive chemicals was similar across different vascular networks and a baseline well-mixed compartment. When chemicals were rapidly metabolized, concentration heterogeneity of the parent compound increased across the vascular network. As a result, our spatially extended lobule generated greater variability in dose-dependent cellular responses, in this case apoptosis, than were observed in the classical well-mixed liver or in a parallel tubes model. The mass-balanced graphical approach to modeling the hepatic lobule is computationally efficient for simulating long-term exposure, modular for incorporating complex cellular interactions, and flexible for dealing with evolving tissues.     

An isolated rat liver model for the evaluation of thermal techniques to quantify perfusion

An isolated, thermally regulated, perfused rat liver model system is presented. The model was developed to evaluate thermal methods to quantify perfusion in small volumes of tissue. The surgically isolated rat liver is perfused with an isothermal oxygenated Krebs-Ringer bicarbonate buffer solution via the cannulated portal vein. A constant-pressure head variable-resistance scheme is utilized to control the total flow to the liver. Total flow is quantified by hepatic vein collection. The spatial distribution of perfusion within the liver is determined using two independent methods. In the first method, radio-labelled microspheres are injected into the portal vein, and the regional flow distribution is determined from the relative radioactivity of each section of tissue. In the second method, the tissue is thermally perturbed, and the time constant of the tissue temperature recovery is measured. The regional distribution is determined from the relative time constants of each section of tissue. Both methods require the measurement of total liver flow to determine the absolute perfusion at each point. Results obtained by the two methods were well correlated (0.973). The rat liver system offers a stable, controllable, and measurable perfusion model for the evaluation of new perfusion measurement techniques.     

A 3D porous media liver lobule model: the importance of vascular septa and anisotropic permeability for homogeneous perfusion

The hepatic blood circulation is complex, particularly at the microcirculatory level. Previously, 2D liver lobule models using porous media and a 3D model using real sinusoidal geometries have been developed. We extended these models to investigate the role of vascular septa (VS) and anisotropic permeability. The lobule was modelled as a hexagonal prism (with or without VS) and the tissue was treated as a porous medium (isotropic or anisotropic permeability). Models were solved using computational fluid dynamics. VS inclusion resulted in more spatially homogeneous perfusion. Anisotropic permeability resulted in a larger axial velocity component than isotropic permeability. A parameter study revealed that results are most sensitive to the lobule size and radial pressure drop. Our model provides insight into hepatic microhaemodynamics, and suggests that inclusion of VS in the model leads to perfusion patterns that are likely to reflect physiological reality. The model has potential for applications to unphysiological and pathological conditions.     

Modeling of the contrast-enhanced perfusion test in liver based on the multi-compartment flow in porous media

Journal of Mathematical Biology - The paper deals with modeling the liver perfusion intended to improve quantitative analysis of the tissue scans provided by the contrast-enhanced computed...     

Application of 128 Slice 4D CT Whole Liver Perfusion Imaging in Hepatic Tumor

To investigate the clinical significance of 128 slice whole liver four dimensional computed tomography (4D CT) in diagnosis and differential diagnosis of hepatic disease, by characterizing and comparing perfusion maps in two common hepatic tumors: hepatocellular carcinoma (HCC) and liver hemangioma. 45 patients with HCC and 40 patients with liver hemangioma were subjected to 128 slice 4D CT of the whole liver perfusion scan, perfusion images were obtained, and data were processed by the perfusion software. Four perfusion parameters generated automatically were used to characterize and compare the perfusion of tumor tissue and surrounding hepatic parenchyma: blood flow perfusion (BF), arterial liver perfusion (ALP), portal venous perfusion (PVP), and hepatic perfusion index (HPI). Volumetric CT perfusion data then reconstructed to yield 4D CT angiography. Morphological observation was made regarding to the blood supply of tumor, intrahepatic vasculature. (1) In both HCC and hepatic hemangioma, BF, ALP, HPI were higher (P < 0.01), whereas PVP were lower (P < 0.01) in tumor tissue than the surrounding hepatic parenchyma (within 1 cm of lesion). Compared with liver hemangioma tumor tissue, BF, ALP, PVP were lower in HCC tumor tissue (P < 0.05; 0.01; 0.01), but HPI is higher (P < 0.05). For the perfusion of the surrounding parenchyma, BF and ALP were higher (P < 0.001), PVP was lower (P < 0.001) in HCC, while HPI was unchanged. (2) Among 45 cases with HCC, cancer feeding artery was found in 28 cases. In 20 cases feeding artery was shown as thickening, rigid, or distorted. Tumor thrombus in portal vein was found in 14 cases. For total of 40 cases with liver hemangioma, in 23 cases blood vessels are shifted due to compression from tumor mass, the rest 17 cases show normal vasculature. With application of 128 slice 4D CT, whole liver perfusion scan can reliably reflect the hemodynamic characteristics of HCC and hepatic hemangioma, proving to be a valuable adjunct to conventional imaging techniques of liver for early detection, differential diagnosis, and determining surgical resection range as well as estimating prognosis for hepatic tumors.     

Effects of temperature on the degradation and biliary secretion of asialoorosomucoid by the perfused rat liver: evidence for two intracellular pathways

We have utilized the in situ perfused rat liver under nonrecirculating conditions to examine the effect of temperature on the metabolism and biliary secretion of [125I]-asialoorosomucid (ASOR). In this manner we were able to follow the fate of a single round of internalized ligand. In control livers perfused at 37 degrees C, approximately 50% of [125I]-ASOR injected into the portal vein was extracted on first pass. Five minutes after the injection, radioactivity, which had been extracted initially, began to appear in the hepatic venous effluent. Within 25 min, 50% of the initially extracted radioactivity was released into the perfusion medium; the bulk of this radioactivity (greater than 95%) was soluble in trichloroacetic acid. In livers perfused at temperatures slightly less than 37 degrees C (30-35 degrees C), first-pass extraction of [125I]-ASOR was similar to that observed at 37 degrees C. However, a severalfold decrease in the rate of release of radioactivity from the liver into the perfusion medium was noted at the lower perfusion temperatures; whereas greater than 50% of the initially extracted radioactivity was released within 30 min from livers perfused at 37 degrees C, only 5% was released at 30 degrees C. At the lower perfusion temperature, a larger proportion of the released radioactivity was acid precipitable (24% vs. 5%). Some radioactivity also was recovered in the bile; of the total amount of radioactivity released from the liver in 30 min at 37 degrees C, approximately 5% was directed into the bile. At lower temperatures of perfusion, a greater fraction of the radioactivity that was released from the liver was directed into the bile (20% at 30 degrees C vs. 5% at 37 degrees C). The data imply that the endosomal pathway to the lysosome is highly sensitive to slight reductions in temperature while the transcytotic route into bile is less sensitive. Lower temperatures might prolong the residence time of ASOR in the prelysosomal endosomal compartments, and thereby increase the likelihood that undegraded ligand will be returned to the blood or be missorted into bile.     

Assessment of Bioavailability after In Vitro Digestion and First Pass Metabolism of Bioactive Peptides from Collagen Hydrolysates

Collagen hydrolysates (CHs) are composed of bioactive peptides (BAPs), which possess health enhancing properties. There is a knowledge gap regarding the bioavailability of these BAPs that involves intestinal transport and hepatic first pass effects. A simulated gastrointestinal model was used to generate digesta from two CHs (CH-GL and CH-OPT), which were applied to a novel transwell co-culture of human intestinal epithelium cell line-6 (HIEC-6) and hepatic (HepG2) cells to simulate in vivo conditions of absorption and first pass metabolism. Peptide transport, hepatic first pass effects, and bioavailability were determined by measuring BAPs (Gly-Pro, Hyp-Gly, Ala-Hyp, Pro-Hyp, Gly-Pro-Hyp) using an innovative capillary electrophoresis method. All peptides were transported across the intestinal cell layer to varying degrees with both CHs; however, Gly-Pro-Hyp was transported only with CH-GL, but not CH-OPT. Notable hepatic production was observed for Ala-Hyp with both CH treatments, and for Pro-Hyp and Gly-Pro with CH-GL only. All peptides were bioavailable (>10%), except for Gly-Pro-Hyp after CH-OPT. Overall, a high degree of transport and hepatic first pass effects on CH-derived BAPs were observed. Further research is needed to explore the hepatic mechanisms related to the production of BAPs and the bifunctional effects of the bioavailable BAPs noted in this study.     

Influence of tauroursodeoxycholic and taurodeoxycholic acids on hepatic metabolism and biliary secretion of phosphatidylcholine in the isolated rat liver

Studies were carried out using an isolated rat liver system to define: the contribution of exogenous phosphatidylcholine (PC) to biliary phospholipid secretion; and its hepatic metabolism during perfusion of the livers with conjugated bile salts with different hydrophilic/hydrophobic properties. A tracer dose of sn-1-palmitoyl-sn-2-[14C]linoleoylPC was injected as a bolus into the recirculating liver perfusate, under constant infusion of 0.75 mumol/min of tauroursodeoxycholate or taurodeoxycholate. The effects on bile flow, biliary lipid secretion, 14C disappearance from the perfusate and its appearance in bile, as well as hepatic and biliary biotransformation were determined. With both the bile salts, about 40% of the [14C]PC was taken up by the liver from the perfusate over 100 min. During the same period less than 2% of the given radioactivity was secreted into bile. More than 95% of the 14C recovered in bile was located within the identical injected PC molecular species. The biliary secretion of labeled as well as unlabeled PC, however, was significantly higher in livers perfused with taurodeoxycholate than tauroursodeoxycholate, while the reverse was observed with respect to bile flow and total bile salt secretion. The exogenous PC underwent extensive hepatic metabolization which appeared to be influenced by the type of bile salt perfusing the liver. After 2 h perfusion, the liver radioactivity was found, in decreasing order, in PC, triacylglycerol, phosphatidylethanolamine and diacylglycerol. In addition, the specific activity of triacylglycerol was significantly higher in tauroursodeoxycholate than in taurodeoxycholate-perfused livers (P less than 0.025), while the reverse was true for the specific activity of hepatic PC (P less than 0.01). Because taurodeoxycholate and tauroursodeoxycholate showed opposite effects on both biliary lipid secretion and hepatic PC biotransformations, we conclude that the hepatic metabolism of glycerolipids is influenced by the physiochemical properties of bile salts.     

Mathematical analysis of cancer chemotherapy

This paper presents a mathematical analysis of a tumor model first proposed by Skipper and Zubrod. The tumor model is comprised of three compartments, a proliferative compartment, a nonproliferative but viable compartment, and a dead compartment. By the suitable selection of functions describing loss of cells from the proliferative and nonproliferative compartments, the model is capable of describing tumor behavior during periods of growth and drug treatment. The loss functions during treatment are related to pharmacokinetic functions and may be chosen according to known drug properties. Tumor properties may be simulated by the appropriate choice of cell cycle parameters. It therefore seems feasible to simulate tumor behavior for scheduled treatment with chemotherapeutic agents. Another important result of this analysis is the derivation of a fraction labelled mitosis function which incorporates the nonproliferative compartments.     

Liver bioengineering: From the stage of liver decellularized matrix to the multiple cellular actors and bioreactor special effects

Liver bioengineering has been a field of intense research and popular excitement in the past decades. It experiences great interest since the introduction of whole liver acellular scaffolds generated by perfusion decellularization^1–3. Nevertheless, the different strategies developed so far have failed to generate hepatic tissue in vitro bioequivalent to native liver tissue. Even notable novel strategies that rely on iPSC-derived liver progenitor cells potential to self-organize in association with endothelial cells in hepatic organoids are lacking critical components of the native tissue (e.g., bile ducts, functional vascular network, hepatic microarchitecture, etc)⁴. Hence, it is vital to understand the strengths and short comes of our current strategies in this quest to re-create liver organogenesis in vitro. To shed some light into these issues, this review describes the different actors that play crucial roles in liver organogenesis and highlights the steps still missing to successfully generate whole livers and hepatic organoids in vitro for multiple applications.     

Liver bioengineering

Liver bioengineering has been a field of intense research and popular excitement in the past decades. It experiences great interest since the introduction of whole liver acellular scaffolds generated by perfusion decellularization^1–3. Nevertheless, the different strategies developed so far have failed to generate hepatic tissue in vitro bioequivalent to native liver tissue. Even notable novel strategies that rely on iPSC-derived liver progenitor cells potential to self-organize in association with endothelial cells in hepatic organoids are lacking critical components of the native tissue (e.g., bile ducts, functional vascular network, hepatic microarchitecture, etc)⁴. Hence, it is vital to understand the strengths and short comes of our current strategies in this quest to re-create liver organogenesis in vitro. To shed some light into these issues, this review describes the different actors that play crucial roles in liver organogenesis and highlights the steps still missing to successfully generate whole livers and hepatic organoids in vitro for multiple applications.     

A real time hyperelastic tissue model

Real-time soft tissue modeling has a potential application in medical training, procedure planning and image-guided therapy. This paper characterizes the mechanical properties of organ tissue using a hyperelastic material model, an approach which is then incorporated into a real-time finite element framework. While generalizable, in this paper we use the published mechanical properties of pig liver to characterize an example application. Specifically, we calibrate the parameters of an exponential model, with a least-squares method (LSM) using the assumption that the material is isotropic and incompressible in a uniaxial compression test. From the parameters obtained, the stress–strain curves generated from the LSM are compared to those from the corresponding computational model solved by ABAQUS and also to experimental data, resulting in mean errors of 1.9 and 4.8%, respectively, which are considerably better than those obtained when employing the Neo-Hookean model. We demonstrate our approach through the simulation of a biopsy procedure, employing a tetrahedral mesh representation of human liver generated from a CT image. Using the material properties along with the geometric model, we develop a nonlinear finite element framework to simulate the behaviour of liver during an interventional procedure with a real-time performance achieved through the use of an interpolation approach.     

Portal branch ligation induces a hepatic arterial buffer response, microvascular remodeling, normoxygenation, and cell proliferation in portal blood-deprived liver tissue

Portal branch ligation (PBL) may prevent liver failure after extended hepatic resection. However, clinical studies indicate that tumors within the ligated lobe develop accelerated growth. Although it is well known that tumor growth depends on the host's microvascularization, there is no information about how PBL affects the hepatic microcirculation. Our aims were to determine hepatic artery response, liver microcirculation, tissue oxygenation, and cell proliferation after PBL. Therefore, we used intravital multifluorescence microscopy, laser-Doppler flowmetry, immunohistochemistry, and biochemical techniques to examine microcirculatory responses, microvascular remodeling, and cellular consequences after left lateral PBL in BALB/c mice. During the first 7 days, PBL induced a reduction of left hilar blood flow by approximately 50%. This resulted in 80% sinusoidal perfusion failure, significant parenchymal hypoxia, and liver atrophy. After 14 days, however, left hilar blood flow was found to be restored. However, remodeling of the microvasculature included a rarefaction of the sinusoidal network, however, without substantial perfusion failure, compensated by a hepatic arterial buffer response and significant sinusoidal dilatation. This resulted in normalization of tissue oxygenation, indicating arterialization of the ligated lobe. Interestingly, late microvascular remodeling was associated with increased endothelial nitric oxide synthase expression, significant hepatocellular proliferation, and weight gain of the ligated lobe. Thus PBL induces only an initial microcirculatory failure with liver atrophy, followed by a hepatic arterial buffer response, microvascular remodeling, normoxygenation, and hepatocellular proliferation. This may explain the accelerated tumor progression occasionally observed in patients after PBL.     

A real time hyperelastic tissue model

Real-time soft tissue modeling has a potential application in medical training, procedure planning and image-guided therapy. This paper characterizes the mechanical properties of organ tissue using a hyperelastic material model, an approach which is then incorporated into a real-time finite element framework. While generalizable, in this paper we use the published mechanical properties of pig liver to characterize an example application. Specifically, we calibrate the parameters of an exponential model, with a least-squares method (LSM) using the assumption that the material is isotropic and incompressible in a uniaxial compression test. From the parameters obtained, the stress-strain curves generated from the LSM are compared to those from the corresponding computational model solved by ABAQUS and also to experimental data, resulting in mean errors of 1.9 and 4.8%, respectively, which are considerably better than those obtained when employing the Neo-Hookean model. We demonstrate our approach through the simulation of a biopsy procedure, employing a tetrahedral mesh representation of human liver generated from a CT image. Using the material properties along with the geometric model, we develop a nonlinear finite element framework to simulate the behaviour of liver during an interventional procedure with a real-time performance achieved through the use of an interpolation approach.     

空间直观景观模型的验证方法

空间直观景观模型已是当前景观生态学研究的一大热点。空间景观模型模拟空间格局变化。其模拟结果包含非空间数据和空间数据。空间直观景观模型的验证除进行非空间数据的验证外,还需要进行空间数据的验证。本文回顾了空间直观模型发展历程,总结现有的空间直观模型验证方法。包括主观评价、图形比较、偏差分析、回归分析、假设检验、多尺度拟合度分析和景观指数分析,同时提出今后空间直观景观模型验证方法研究的重点方向。     

The Kinetics of Distribution of the Fat-Soluble Inert Gas Cyclopropane in the Body

A kinetic analysis is made of the experimentally measured time course of respiratory uptake of the highly fat-soluble, inert gas cyclopropane by normal human subjects. The analysis is based on the well-known perfusion-limited model in which a number of body compartments are arranged in parallel with the lungs via the circulating blood. Three distinct body compartments are derived from the data. These are tentatively identified as: (a) adipose tissue (b) fat-poor tissue of low perfusion such as resting muscle, skin, and connective tissue (c) fat-poor tissue of high perfusion such as brain, heart, gut, liver, and kidney. Blood flow rates to the several compartments are also derived from the data. The rates to compartments (a) and (b) are each approximately 10 per cent of the estimated total cardiac output. The derived perfusion (blood flow rate/compartment weight) of the three compartments are in the range, respectively, (a) 2 to 4, (b) 1 to 2.5, (c) 25 to 75 ml/min/100 gm. Uncertainties arising from the experimental data and from simplifications of the model (neglect of lung fill-up phase of uptake and gross diffusion of cyclopropane from one tissue into another) are discussed. The present type of uptake experiment is significant for the problems of total body fat determination, of gross body composition in relation to weight change, of gross shunting of blood flow from one compartment to another, of anesthesia by fat-soluble substances, and of decompression sickness.     

A nonlinear kinetic model of γ-aminobutyric acid metabolism in a rabbit model of acute liver failure

To investigate the mechanisms by which serum levels of γ-aminobutyric acid (GABA) become elevated in experimental acute liver failure, a multicompartmental model of GABA metabolism has been constructed and used to simulate previously generated data on the kinetics of ³H-GABA uptake by isolated hepatocytes from normal rats and the kinetics of ³H-GABA in the plasma of normal rabbits, rabbits with galactosamine-induced acute liver failure, and rabbits with divascularized livers. Modeling analysis revealed that acute liver failure was associated with values for the mean fractional catabolic rate of GABA, plasma volume, and hepatic extraction of GABA that were 29%, 12%, and 49% less, respectively, than the corresponding control values. The defect in hepatic tissue extraction of GABA was sufficient to account for only 60% of the 10-fold increase in serum GABA levels that occurs in acute liver failure. Furthermore the 10-fold increase in serum GABA levels occured in acute liver failure before the onset of overt hepatic encephalopathy when hepatic extraction of GABA was not appreciably different from that found in normal rabbits. Thus the increase in serum GABA levels that occurs in acute liver failure cannot be attributed to a defect in hepatic extraction of GABA alone. Indeed, the modeling analysis indicated that in acute liver failure there is a 3—8-fold increase in the rate of delivery of GABA to the systemic circulation, but did not indicate its source.     

Comparaison de deux modèles hybrides simulant la propagation de la lumière dans les tissus biologiques

In this paper, two numerical hybrid methods to model photon transport phenomena in biological tissues are compared. The coupled radiative transfer–diffusion model is based on the finite element solution of the radiative transfer equation and its approximation. The hybrid Monte-Carlo–diffusion consists in modeling the propagation of laser light in turbid media with the pure statistical Monte-Carlo method in the vicinity of the source and the boundaries and the diffusion approximation elsewhere in the domain. We apply these codes to calculate the spatially resolved reflectance amplitude and phase resulting from an intensity modulated laser beam. The results show that the hybrid methods can be used to simulate light propagation with good accuracy and speed.