A Pilot Study on Linking Tissue Mechanics with Load-Dependent Collagen Microstructures in Porcine Tricuspid Valve Leaflets

An investigation of how specimen dimensions affect biaxial mechanical characterizations with CellScale BioTester and constitutive modeling of porcine tricuspid valve leaflets

Biaxial mechanical characterizations are the accepted approach to determine the mechanical response of many biological soft tissues. Although several computational and experimental studies have examined how experimental factors (e.g., clamped vs. suture mounting) affect the acquired tissue mechanical behavior, little is known about the role of specimen dimensions in data acquisition and the subsequent modeling. In this study, we combined our established mechanical characterization framework with an iterative size-reduction protocol to test the hypothesis that specimen dimensions affect the observed mechanical behavior of biaxial characterizations. Our findings indicated that there were non-significant differences in the peak equibiaxial stretches of tricuspid valve leaflets across four specimen dimensions ranging from 4.5×4.5mm to 9 × 9mm. Further analyses revealed that there were significant differences in the low-tensile modulus of the circumferential tissue direction. These differences resulted in significantly different constitutive model parameters for the Tong-Fung model between different specimen dimensions of the posterior and septal leaflets. Overall, our findings demonstrate that specimen dimensions play an important role in experimental characterizations, but not necessarily in constitutive modeling of soft tissue mechanical behavior during biaxial testing with the commercial CellScale BioTester.

Biomechanics Assist Measurement, Modeling, Engineering Applications, and Clinical Decision Making in Medicine

Biomechanical studies of surgeries and medical devices are usually performed with human or animal models [...].

A Computational Framework for Atrioventricular Valve Modeling Using Open-Source Software

Atrioventricular valve regurgitation is a significant cause of morbidity and mortality in patients with acquired and congenital cardiac valve disease. Image-derived computational modeling of atrioventricular valves has advanced substantially over the last decade and holds particular promise to inform valve repair in small and heterogeneous populations, which are less likely to be optimized through empiric clinical application. While an abundance of computational biomechanics studies has investigated mitral and tricuspid valve disease in adults, few studies have investigated its application to vulnerable pediatric and congenital heart populations. Further, to date, investigators have primarily relied upon a series of commercial applications that are neither designed for image-derived modeling of cardiac valves nor freely available to facilitate transparent and reproducible valve science. To address this deficiency, we aimed to build an open-source computational framework for the image-derived biomechanical analysis of atrioventricular valves. In the present work, we integrated an open-source valve modeling platform, SlicerHeart, and an open-source biomechanics finite element modeling software, FEBio, to facilitate image-derived atrioventricular valve model creation and finite element analysis. We present a detailed verification and sensitivity analysis to demonstrate the fidelity of this modeling in application to three-dimensional echocardiography-derived pediatric mitral and tricuspid valve models. Our analyses achieved an excellent agreement with those reported in the literature. As such, this evolving computational framework offers a promising initial foundation for future development and investigation of valve mechanics, in particular collaborative efforts targeting the development of improved repairs for children with congenital heart disease.

Benchtop characterization of the tricuspid valve leaflet pre-strains

The pre-strains of biological soft tissues are important when relating their in vitro and in vivo mechanical behaviors. In this study, we present the first-of-its-kind experimental characterization of the tricuspid valve leaflet pre-strains. We use 3D photogrammetry and the reproducing kernel method to calculate the pre-strains within the central 10×10 mm region of the tricuspid valve leaflets from n=8 porcine hearts. In agreement with previous pre-strain studies for heart valve leaflets, our results show that all the three tricuspid valve leaflets shrink after explant from the ex vivo heart. These calculated strains are leaflet-specific and the septal leaflet experiences the most compressive changes. Furthermore, the strains observed after dissection of the central 10×10 mm region of the leaflet are smaller than when the valve is explanted, suggesting that our computed pre-strains are mainly due to the release of in situ annulus and chordae connections. The leaflets are then mounted on a biaxial testing device and preconditioned using force-controlled equibiaxial loading. We show that the employed preconditioning protocol does not 100% restore the leaflet pre-strains as removed during tissue dissection, and future studies are warranted to explore alternative preconditioning methods. Finally, we compare the calculated biomechanically oriented metrics considering five stress-free reference configurations. Interestingly, the radial tissue stretches and material anisotropies are significantly smaller compared to the post-preconditioning configuration. Extensions of this work can further explore the role of this unique leaflet-specific leaflet pre-strains on in vivo valve behavior via high-fidelity in-silico models.

Simultaneous Wide-Field Planar Strain-Fiber Orientation Distribution Measurement Using Polarized Spatial Domain Imaging

In the present study, we demonstrate that soft tissue fiber architectural maps captured using polarized spatial frequency domain imaging (pSFDI) can be utilized as an effective texture source for DIC-based planar surface strain analyses. Experimental planar biaxial mechanical studies were conducted using pericardium as the exemplar tissue, with simultaneous pSFDI measurements taken. From these measurements, the collagen fiber preferred direction [Formula: see text] was determined at the pixel level over the entire strain range using established methods ( https://doi.org/10.1007/s10439-019-02233-0 ). We then utilized these pixel-level [Formula: see text] maps as a texture source to quantify the deformation gradient tensor [Formula: see text] as a function of spatial position [Formula: see text] within the specimen at time t. Results indicted that that the pSFDI approach produced accurate deformation maps, as validated using both physical markers and conventional particle based method derived from the DIC analysis of the same specimens. We then extended the pSFDI technique to extract the fiber orientation distribution [Formula: see text] as a function of [Formula: see text] from the pSFDI intensity signal. This was accomplished by developing a calibration procedure to account for the optical behavior of the constituent fibers for the soft tissue being studied. We then demonstrated that the extracted [Formula: see text] was accurately computed in both the referential (i.e. unloaded) and deformed states. Moreover, we noted that the measured [Formula: see text] agreed well with affine kinematic deformation predictions. We also demonstrated this calibration approach could also be effectively used on electrospun biomaterials, underscoring the general utility of the approach. In a final step, using the ability to simultaneously quantify [Formula: see text] and [Formula: see text], we examined the effect of deformation and collagen structural measurements on the measurement region size. For pericardial tissues, we determined a critical length of [Formula: see text] 8 mm wherein the regional variations sufficiently dissipated. This result has immediate potential in the identification of optimal length scales for meso-scale strain measurement in soft tissues and fibrous biomaterials.

Determination of a Strain Energy Density Function for the Tricuspid Valve Leaflets Using Constant Invariant-Based Mechanical Characterizations

The tricuspid valve (TV) regulates the blood flow within the right side of the heart. Despite recent improvements in understanding TV mechanical and microstructural properties, limited attention has been devoted to the development of TV-specific constitutive models. The objective of this work is to use the first-of-its-kind experimental data from constant invariant-based mechanical characterizations to determine a suitable invariant-based strain energy density function (SEDF). Six specimens for each TV leaflet are characterized using constant invariant mechanical testing. The data is then fit with three candidate SEDF forms: (i) a polynomial model-the transversely isotropic version of the Mooney-Rivlin model, (ii) an exponential model, and (iii) a combined polynomial-exponential model. Similar fitting capabilities were found for the exponential and the polynomial forms (R2=0.92-0.99 versus 0.91-0.97) compared to the combined polynomial-exponential SEDF (R2=0.65-0.95). Furthermore, the polynomial form had larger Pearson's correlation coefficients than the exponential form (0.51 versus 0.30), indicating a more well-defined search space. Finally, the exponential and the combined polynomial-exponential forms had notably smaller but more eccentric model parameter's confidence regions than the polynomial form. Further evaluations of invariant decoupling revealed that the decoupling of the invariant terms within the exponential form leads to a less satisfactory performance. From these results, we conclude that the exponential form is better suited for the TV leaflets owing to its superb fitting capabilities and smaller parameter's confidence regions.

Linking collagen fiber architecture to tissue-level biaxial mechanical behaviors of porcine semilunar heart valve cusps

The semilunar heart valves regulate the blood flow from the ventricles to the major arteries through the opening and closing of the scallop shaped cusps. These cusps are composed of collagen fibers that act as the primary loading-bearing component. The load-dependent collagen fiber architecture has been previously examined in the existing literature; however, these studies relied on chemical clearing and tissue modifications to observe the underlying changes in response to mechanical loads. In the present study, we address this gap in knowledge by quantifying the collagen fiber orientations and alignments of the aortic and pulmonary cusps through a multi-scale, non-destructive experimental approach. This opto-mechanical approach, which combines polarized spatial frequency domain imaging and biaxial mechanical testing, provides a greater field of view (10-25mm) and faster imaging time (45-50s) than other traditional collagen imaging techniques. The birefringent response of the collagen fibers was fit with a von Mises distribution, while the biaxial mechanical testing data was implemented into a modified full structural model for further analysis. Our results showed that the semilunar heart valve cusps are more extensible in the tissue's radial direction than the circumferential direction under all the varied biaxial testing protocols, together with greater material anisotropy among the pulmonary valve cusps compared to the aortic valve cusps. The collagen fibers were shown to reorient towards the direction of the greatest applied loading and incrementally realign with the increased applied stress. The collagen fiber architecture within the aortic valve cusps were found to be more homogeneous than the pulmonary valve counterparts, reflecting the differences in the physiological environments experienced by these two semilunar heart valves. Further, the von Mises distribution fitting highlighted the presence and contribution of two distinct fiber families for each of the two semilunar heart valves. The results from this work would provide valuable insight into connecting tissue-level mechanics to the underlying collagen fiber architecture-an essential information for the future development of high-fidelity aortic/pulmonary valve computational models.

Effects of enzyme-based removal of collagen and elastin constituents on the biaxial mechanical responses of porcine atrioventricular heart valve anterior leaflets

The leaflets of the atrioventricular heart valves (AHVs) regulate the one-directional flow of blood through a coordination of the extracellular matrix components, including the collagen fibers, elastin, and glycosaminoglycans. Dysfunction of the AHVs, such as those caused by unfavorable microstructural remodeling, lead to valvular heart diseases and improper blood flow, which can ultimately cause heart failure. In order to better understand the mechanics and remodeling of the AHV leaflets and how therapeutics can inadvertently cause adverse microstructural changes, a systematic characterization of the role of each constituent in the biomechanical properties is appropriate. Previous studies have quantified the contributions of the individual microstructural components to tissue-level behavior for the semilunar valve cusps, but not for the AHV leaflets. In this study, for the first time, we quantify the relationships between microstructure and mechanics of the AHV leaflet using a three-step experimental procedure: (i) biaxial tension and stress relaxation testing of control (untreated) porcine AHV anterior leaflet specimens; (ii) enzyme treatment to remove a portion of either the collagen or elastin constituent; and (iii) biaxial tensile and stress relaxation testing of the constituent-removed (treated) specimens. We have observed that the removal of ∼100% elastin resulted in a ∼10% decrease in the tissue extensibility with biaxial tension and a ∼10% increase in the overall stress reduction with stress relaxation. In contrast, removal of 46% of the collagen content insignificantly affected tissue extensibility with biaxial tension and significantly increased stress decay (10%) with stress relaxation. These findings provide an insight into the microstructure-mechanics relationship of the AHVs and will be beneficial for future developments and refinements of microstructurally informed constitutive models for the simulation of diseased and surgically intervened AHV function. STATEMENT OF SIGNIFICANCE: This study presents, for the first time, a thorough mechanical characterization of the atrioventricular heart valve leaflets before and after enzymatic removal of elastin and collagen. We found that the biaxial tensile properties of elastin-deficient tissues and collagen-deficient are stiffer. The fact of elastin supporting low-stress valve function and collagen as the main load-bearing component was evident in a decrease in the low-tension modulus for elastin-deficient tissues and in the high-tension modulus for collagen-deficient tissues. Our quantification and experimental technique could be useful in predicting the disease-related changes in heart valve mechanics. The information obtained from this work is valuable for refining the constitutive models that describe the essential microstructure-mechanics relationship.

Parameterization, geometric modeling, and isogeometric analysis of tricuspid valves

Approximately 1.6 million patients in the United States are affected by tricuspid valve regurgitation, which occurs when the tricuspid valve does not close properly to prevent backward blood flow into the right atrium. Despite its critical role in proper cardiac function, the tricuspid valve has received limited research attention compared to the mitral and aortic valves on the left side of the heart. As a result, proper valvular function and the pathologies that may cause dysfunction remain poorly understood. To promote further investigations of the biomechanical behavior and response of the tricuspid valve, this work establishes a parameter-based approach that provides a template for tricuspid valve modeling and simulation. The proposed tricuspid valve parameterization presents a comprehensive description of the leaflets and the complex chordae tendineae for capturing the typical three-cusp structural deformation observed from medical data. This simulation framework develops a practical procedure for modeling tricuspid valves and offers a robust, flexible approach to analyze the performance and effectiveness of various valve configurations using isogeometric analysis. The proposed methods also establish a baseline to examine the tricuspid valve's structural deformation, perform future investigations of native valve configurations under healthy and disease conditions, and optimize prosthetic valve designs.

Evaluation of affine fiber kinematics in porcine tricuspid valve leaflets using polarized spatial frequency domain imaging and planar biaxial testing

Collagen fibers are the primary load-bearing microstructural constituent of bodily soft tissues, and, when subjected to external loading, the collagen fibers reorient, uncrimp, and elongate. Specific to the atrioventricular heart valve leaflets, the collagen fiber kinematics form the basis of many constitutive models; however, some researchers claim that modeling the affine fiber kinematics (AFK) are sufficient for accurately predicting the macroscopic tissue deformations, while others state that modeling the non-affine kinematics (i.e., fiber uncrimping together with elastic elongation) is required. Experimental verification of the AFK theory has been previously performed for the mitral valve leaflets in the left-side heart; however, this same evaluation has yet to be performed for the morphologically distinct tricuspid valve (TV) leaflets in the right-side heart. In this work, we, for the first time, evaluated the AFK theory for the TV leaflets using an integrated biaxial testing-polarized spatial frequency domain imaging device to experimentally quantify the load-dependent collagen fiber reorientations for comparison to the AFK theory predictions. We found that the AFK theory generally underpredicted the fiber reorientations by 3.1°, on average, under the applied equibiaxial loading with greater disparity when the tissue was subjected to the applied non-equibiaxial loading. Furthermore, increased AFK errors were observed with increasing collagen fiber reorientations (Pearson coefficient r = -0.36, equibiaxial loading), suggesting the AFK theory is better suited for relatively smaller reorientations. Our findings suggest the AFK theory may require modification for more accurate predictions of the collagen fiber kinematics in the TV leaflets, which will be useful in refining modeling efforts for more accurate TV simulations.