Diameter-dependent axial prestretch of porcine coronary arteries and veins

Biomechanical effects of hemin and sildenafil treatments on the aortic wall of chronic-hypoxic lambs

Introduction: Gestation under chronic hypoxia causes pulmonary hypertension, cardiovascular remodeling, and increased aortic stiffness in the offspring. To mitigate the neonatal cardiovascular risk, pharmacological treatments (such as hemin and sildenafil) have been proposed to improve pulmonary vasodilation. However, little is known about the effects of these treatments on the aorta. Therefore, we studied the effect of hemin and sildenafil treatments in the aorta of lambs gestated and raised at highlands, thereby subjected to chronic hypoxia. Methods: Several biomechanical tests were conducted in the descending thoracic aorta (DTA) and the distal abdominal aorta (DAA), assessing 3 groups of study of hypoxic animals: non-treated (Control) and treated either with hemin or sildenafil. Based on them, the stiffness level has been quantified in both zones, along with the physiological strain in the unloaded aortic duct. Furthermore, a morphological study by histology was conducted in the DTA. Results: Biomechanical results indicate that treatments trigger an increment of axial pre-stress and circumferential residual stress levels in DTA and DAA of lambs exposed to high-altitude chronic hypoxia, which reveals a vasodilatation improvement along with an anti-hypertensive response under this characteristic environmental condition. In addition, histological findings do not reveal significant differences in either structure or microstructural content. Discussion: The biomechanics approach emerges as a valuable study perspective, providing insights to explain the physiological mechanisms of vascular function. According to established results, alterations in the function of the aortic wall may not necessarily be explained by morphostructural changes, but rather by the characteristic mechanical state of the microstructural components that are part of the studied tissue. In this sense, the reported biomechanical changes are beneficial in mitigating the adverse effects of hypobaric hypoxia exposure during gestation and early postnatal life.

Guided wave elastography of jugular veins: Theory, method and in vivo experiment

Testing the mechanical properties of veins is important for diagnosing some cardiovascular diseases such as deep venous thrombosis. Additionally, it plays a crucial role in designing body protective products such as head protective gear, where simulations are necessary to predict the mechanical responses of bridging veins during head impacts. The data on venous mechanical properties reported in the literature have mainly been obtained from ex vivo experiments, and inferring the material parameters of veins in vivo is challenging. Here, we address this issue by proposing a guided wave elastography method in which guided waves are generated in the jugular veins with focused acoustic radiation force and tracked by an ultrafast ultrasound imaging system. Then, a mechanical model considering the effects of the perivascular soft tissues and prestresses in the veins was applied to analyze the wave motions in the jugular veins. Our model enables the development of an inverse method to infer the elastic properties of the veins from measured guided waves. Phantom experiments were performed to validate the theory, and in vivo experiments were carried out to demonstrate the usefulness of the inverse method in practice.

Linking the region-specific tissue microstructure to the biaxial mechanical properties of the porcine left anterior descending artery

Coronary atherosclerosis is the main cause of death worldwide. Advancing the understanding of coronary microstructure-based mechanics is fundamental for the development of therapeutic tools and surgical procedures. Although the passive biaxial properties of the coronary arteries have been extensively explored, their regional differences and the relationship between tissue microstructure and mechanics have not been fully characterized. In this study, we characterized the passive biaxial mechanical properties and microstructural properties of the proximal, medial, and distal regions of the porcine left anterior descending artery (LADA). We also attempted to relate the biaxial stress-stretch response of the LADA and its respective birefringent responses to the polarized light for obtaining information about the load-dependent microstructural variations. We found that the LADA extensibility is reduced in the proximal-to-distal direction and that the medial region exhibits more heterogeneous mechanical behavior than the other two regions. We have also observed highly dynamic microstructural behavior where fiber families realign themselves depending on loading. In addition, we found that the microstructure of the distal region exhibited highly aligned fibers along the longitudinal axis of the artery. To verify this microstructural feature, we imaged the LADA specimens with multi-photon microscopy and observed that the adventitia microstructure transitioned from a random fiber network in the proximal region to highly aligned fibers in the distal region. Our findings could offer new perspectives for understanding coronary mechanics and aid in the development of tissue-engineered vascular grafts, which are currently limited due to their mismatch with native tissue in terms of mechanical properties and microstructural features. STATEMENT OF SIGNIFICANCE: The tissue biomechanics of coronary arteries is fundamental for the development of revascularization techniques such as coronary artery bypass. These therapeutics require a deep understanding of arterial mechanics, microstructure, and mechanobiology to prevent graft failure and reoperation. The present study characterizes the unique regional mechanical and microstructural properties of the porcine left anterior descending artery using biaxial testing, polarized-light imaging, and confocal microscopy. This comprehensive characterization provides an improved understanding of the collagen/elastin architecture in response to mechanical loads using a region-specific approach. The unique tissue properties obtained from this study will provide guidance for the selection of anastomotic sites in coronary artery bypass grafting and for the design of tissue-engineered vascular grafts.

Maximizing the Flap Inflow in a Foot Reconstruction: Ultrasonographic Evaluation of Artery Flow in Accordance with the Angle of the Ankle

SUMMARY

Adequate arterial inflow is mandatory for a successful free tissue transfer. Maximizing this inflow appears to be especially important during distal lower extremity reconstruction, where perfusion pressures can be significantly lower than in the trunk or head and neck. In the current study, the authors used color and pulsed-wave Doppler ultrasonography to analyze blood flow volumes in the foot in accordance with the ankle position. Twenty feet of 10 healthy subjects were included in the analyses. All subjects were examined by Doppler ultrasonography at three different ankle angles (dorsiflexion at 90 degrees, relaxed position, and complete plantar flexion). The average flow volumes of the dorsalis pedis artery varied dramatically at the different ankle positions, as follows: dorsiflexion, 19.81 ± 3.88 ml/minute; relaxed position, 11.71 ± 4.72 ml/minute; and plantar flexion, 3.47 ± 1.41 ml/minute. The average flow volume of the posterior tibial artery also varied dramatically at these different ankle positions: dorsiflexion, 9.08 ± 8.25 ml/minute; relaxed position, 14.78 ± 2.45 ml/minute; and plantar flexion, 19.03 ± 5.03 ml/minute. The ankle position dramatically affects the flow velocity in the anterior and posterior tibial vessels, which is an important consideration when performing free tissue transfer reconstruction of the foot. Ankle dorsiflexion would be recommended when using the anterior tibial artery as a recipient vessel, whereas plantar flexion of the ankle would be advised when the posterior tibial artery is used as a recipient vessel.

Branching Pattern of the Cerebral Arterial Tree

Quantitative data on branching patterns of the human cerebral arterial tree are lacking in the 1.0–0.1 mm radius range. We aimed to collect quantitative data in this range, and to study if the cerebral artery tree complies with the principle of minimal work (Law of Murray).

To enable easy quantification of branching patterns a semi‐automatic method was employed to measure 1,294 bifurcations and 2,031 segments on 7 T‐MRI scans of two corrosion casts embedded in a gel. Additionally, to measure segments with a radius smaller than 0.1 mm, 9.4 T‐MRI was used on a small cast section to characterize 1,147 bifurcations and 1,150 segments. Besides MRI, traditional methods were employed. Seven hundred thirty‐three bifurcations were manually measured on a corrosion cast and 1,808 bifurcations and 1,799 segment lengths were manually measured on a fresh dissected cerebral arterial tree. Data showed a large variation in branching pattern parameters (asymmetry‐ratio, area‐ratio, length‐radius‐ratio, tapering). Part of the variation may be explained by the variation in measurement techniques, number of measurements and location of measurement in the vascular tree.

This study confirms that the cerebral arterial tree complies with the principle of minimum work. These data are essential in the future development of more accurate mathematical blood flow models. Anat Rec, 302:1434–1446, 2019. © 2018 The Authors. The Anatomical Record published by Wiley Periodicals, Inc. on behalf of American Association of Anatomists.

A General Approach to Derive Stress and Elasticity Tensors for Hyperelastic Isotropic and Anisotropic Biomaterials

Hyperelastic models are of particular interest in modeling biomaterials. In order to implement them, one must derive the stress and elasticity tensors from the given potential energy function explicitly. However, it is often cumbersome to do so because researchers in biomechanics may not be well-exposed to systematic approaches to derive the stress and elasticity tensors as it is vaguely addressed in literature. To resolve this, we present a framework of a general approach to derive the stress and elasticity tensors for hyperelastic models. Throughout the derivation we carefully elaborate the differences between formulas used in the displacement-based formulation and the displacement/pressure mixed formulation. Three hyperelastic models, Mooney-Rivlin, Yeoh and Holzapfel-Gasser-Ogden models that span from first-order to higher order and from isotropic to anisotropic materials, are served as examples. These detailed derivations are validated with numerical experiments that demonstrate excellent agreements with analytical and other computational solutions. Following this framework, one could implement with ease any hyperelastic model as user-defined functions in software packages or develop as an original source code from scratch.

Mechanobiological model of arterial growth and remodeling

A coupled agent-based model (ABM) and finite element analysis (FEA) computational framework is developed to study the interplay of bio-chemo-mechanical factors in blood vessels and their role in maintaining homeostasis. The agent-based model implements the power of REPAST Simphony libraries and adapts its environment for biological simulations. Coupling a continuum-level model (FEA) to a cellular-level model (ABM) has enabled this computational framework to capture the response of blood vessels to increased or decreased levels of growth factors, proteases and other signaling molecules (on the micro scale) as well as altered blood pressure. Performance of the model is assessed by simulating porcine left anterior descending artery under normotensive conditions and transient increases in blood pressure and by analyzing sensitivity of the model to variations in the rule parameters of the ABM. These simulations proved that the model is stable under normotensive conditions and can recover from transient increases in blood pressure. Sensitivity studies revealed that the model is most sensitive to variations in the concentration of growth factors that affect cellular proliferation and regulate extracellular matrix composition (mainly collagen).

Evaluation of Hemodynamics in a Prestressed and Compliant Tapered Femoral Artery Using an Optimization-Based Inverse Algorithm

The important factors that affect the arterial wall compliance are the tissue properties of the arterial wall, the in vivo pulsatile pressure, and the prestressed condition of the artery. It is necessary to obtain the load-free geometry for determining the physiological level of prestress in the arterial wall. The previously developed optimization-based inverse algorithm was improved to obtain the load-free geometry and the wall prestress of an idealized tapered femoral artery of a dog under varying arterial wall properties. The compliance of the artery was also evaluated over a range of systemic pressures (72.5–140.7 mmHg), associated blood flows, and artery wall properties using the prestressed arterial geometry. The results showed that the computed load-free outer diameter at the inlet of the tapered artery was 6.7%, 9.0%, and 12% smaller than the corresponding in vivo diameter for the 25% softer, baseline, and 25% stiffer arterial wall properties, respectively. In contrast, the variations in the prestressed geometry and circumferential wall prestress were less than 2% for variable arterial wall properties. The computed compliance at the inlet of the prestressed artery for the baseline arterial wall property was 0.34%, 0.19%, and 0.13% diameter change/mmHg for time-averaged pressures of 72.5, 104.1, and 140.7 mmHg, respectively. However, the variation in compliance due to the change in arterial wall property was less than 6%. The load-free and prestressed geometries of the idealized tapered femoral artery were accurately (error within 1.2% of the in vivo geometry) computed under variable arterial wall properties using the modified inverse algorithm. Based on the blood-arterial wall interaction results, the arterial wall compliance was influenced significantly by the change in average pressure. In contrast, the change in arterial wall property did not influence the arterial wall compliance.

Pulsatile arterial wall-blood flow interaction with wall pre-stress computed using an inverse algorithm

Background

The computation of arterial wall deformation and stresses under physiologic conditions requires a coupled compliant arterial wall-blood flow interaction model. The in-vivo arterial wall motion is constrained by tethering from the surrounding tissues. This tethering, together with the average in-vivo pressure, results in wall pre-stress. For an accurate simulation of the physiologic conditions, it is important to incorporate the wall pre-stress in the computational model. The computation of wall pre-stress is complex, as the un-loaded and un-tethered arterial shape with residual stress is unknown. In this study, the arterial wall deformation and stresses in a canine femoral artery under pulsatile pressure was computed after incorporating the wall pre-stresses. A nonlinear least square optimization based inverse algorithm was developed to compute the in-vivo wall pre-stress.

Methods

First, the proposed inverse algorithm was used to obtain the un-loaded and un-tethered arterial geometry from the unstressed in-vivo geometry. Then, the un-loaded, and un-tethered arterial geometry was pre-stressed by applying a mean in-vivo pressure of 104.5 mmHg and an axial stretch of 48% from the un-tethered length. Finally, the physiologic pressure pulse was applied at the inlet and the outlet of the pre-stressed configuration to calculate the in-vivo deformation and stresses. The wall material properties were modeled with an incompressible, Mooney-Rivlin model derived from previously published experimental stress-strain data (Attinger et al., 1968).

Results

The un-loaded and un-tethered artery geometry computed by the inverse algorithm had a length, inner diameter and thickness of 35.14 mm, 3.10 mm and 0.435 mm, respectively. The pre-stressed arterial wall geometry was obtained by applying the in-vivo axial-stretch and average in-vivo pressure to the un-loaded and un-tethered geometry. The length of the pre-stressed artery, 51.99 mm, was within 0.01 mm (0.019%) of the in-vivo length of 52.0 mm; the inner diameter of 3.603 mm was within 0.003 mm (0.08%) of the corresponding in-vivo diameter of 3.6 mm, and the thickness of 0.269 mm was within 0.0015 mm (0.55%) of the in-vivo thickness of 0.27 mm. Under physiologic pulsatile pressure applied to the pre-stressed artery, the time averaged longitudinal stress was found to be 42.5% higher than the circumferential stresses. The results of this study are similar to the results reported by Zhang et al., (2005) for the left anterior descending coronary artery.

Conclusions

An inverse method was adopted to compute physiologic pre-stress in the arterial wall before conducting pulsatile hemodynamic calculations. The wall stresses were higher in magnitude in the longitudinal direction, under physiologic pressure after incorporating the effect of in-vivo axial stretch and pressure loading.

Finite Element simulation of buckling-induced vein tortuosity and influence of the wall constitutive properties

The mechanisms giving rise to vein tortuosity, which is often associated with varicosis, are poorly understood. Recent works suggest that significant biological changes in the wall of varicose veins may precede the mechanical aspects of the disease. To test the hypothesis of tortuosity being a consequence of these changes, a Finite Element model was developed based on previous experimental work on vein buckling. The model was then used to evaluate the effect of alterations of the mechanical behavior of the wall on tortuosity onset and severity. The results showed that increasing anisotropy toward the circumferential direction promotes tortuosity. An increase in wall stiffness tends to decrease the level of tortuosity but interestingly, if the vein segment is little or not pre-stretched such increase will not prevent, or it will even promote, the onset of tortuosity. These results provide additional arguments supporting the hypothesis of tortuosity being the consequence of biologically-induced changes in the varicose vein wall. Based on a 3D model of the leg and in vivo identification of the material properties of varicose veins, a clinical validation of these findings is being developed.

Slackness between vessel and myocardium is necessary for coronary flow reserve

Tone regulation in coronary microvessels has largely been studied in isolated vessels in the absence of myocardial tethering. Here, the potential effect of radial tethering and interstitial space connective tissue (ISCT) between coronary microvessels and the surrounding myocardium was studied. We hypothesized that rigid tethering between microvessels and the myocardium would constrain the active contraction of arterioles and is not compatible with the observed tone regulation. The ISCT between coronary microvessels and myocardium in five swine was found to increase exponentially from 0.22 ± 0.02 μm in capillaries (modified Strahler order 0) of the endocardium to 34.9 ± 7.1 μm in epicardial vessels (order 10). Microvessels with both soft tethering and ISCT gap were capable of significant changes in vessel resistance (up to an ∼1,600% increase), consistent with experimental measurements of high coronary flow reserve. Additionally, the mechanical energy required for myogenic contraction was estimated. The results indicate that rigid tethering requires up to four times more mechanical energy than soft tethering in the absence of a gap. Hence, the experimental measurements and model predictions suggest that effectiveness and efficiency in tone regulation can be achieved only if the vessel is both softly tethered to and separated from the myocardium in accordance with the experimental findings of ISCT gap. These results have fundamental implications on future simulations of coronary circulation.