Regional structural and biomechanical alterations of the ovine main pulmonary artery during postnatal growth

Biaxial testing system for characterization of mechanical and rupture properties of small samples

The study of damage and rupture of soft tissues using a tensile testing system is essential to understand the limits of mechanical behavior and loss of function in diseased tissues. However, commercial material testing systems are often expensive and may not be fully suitable for rupture tests of small samples. While several research laboratories have developed custom, less expensive, uniaxial or biaxial devices, there is a need for an open source, inexpensive, accurate and easy to customize biaxial material testing system to perform rupture tests in small soft samples. We designed a testing system (BiMaTS) that (a) was shown able to perform uniaxial and biaxial tests, (b) offers a large travel range for rupture tests of small samples, (c) maintains a centered field of view for effective strain mapping using digital image correlation, (d) provides a controlled temperature environment, (e) utilize many off-the-shelve components for easy manufacture and customization, and it is cost effective (∼$15 K). The instrument performance was characterized using 80%-scaled down, ASTM D412-C shaped PDMS samples. Our results demonstrate the ability of this open source, customizable, low-cost, biaxial materials testing system to successfully characterize the mechanical and rupture properties of small samples with high repeatability and accuracy.

Differences in Collagen Fiber Diameter and Waviness between Healthy and Aneurysmal Abdominal Aortas

Collagen plays a key role in the strength of aortic walls, so studying micro-structural changes during disease development is critical to better understand collagen reorganization. Second-harmonic generation microscopy is used to obtain images of human aortic collagen in both healthy and diseased states. Methods are being developed in order to efficiently determine the waviness, that is, tortuosity and amplitude, as well as the diameter, orientation, and dispersion of collagen fibers, and bundles in healthy and aneurysmal tissues. The results show layer-specific differences in the collagen of healthy tissues, which decrease in samples of aneurysmal aortic walls. In healthy tissues, the thick collagen bundles of the adventitia are characterized by greater waviness, both in the tortuosity and in the amplitude, compared to the relatively thin and straighter collagen fibers of the media. In contrast, most aneurysmal tissues tend to have a more uniform structure of the aortic wall with no significant difference in collagen diameter between the luminal and abluminal layers. An increase in collagen tortuosity compared to the healthy media is also observed in the aneurysmal luminal layer. The data set provided can help improve related material and multiscale models of aortic walls and aneurysm formation.

The Mechanobiology of Vascular Remodeling in the Aging Lung

Aging is accompanied by declining lung function and increasing susceptibility to lung diseases. The role of endothelial dysfunction and vascular remodeling in these changes is supported by growing evidence, but underlying mechanisms remain elusive. In this review we summarize functional, structural, and molecular changes in the aging pulmonary vasculature and explore how interacting aging and mechanobiological cues may drive progressive vascular remodeling in the lungs.

Ex Vivo Regional Mechanical Characterization of Porcine Pulmonary Arteries

BACKGROUND

Regional mechanical characterization of pulmonary arteries can be useful in the development of computational models of pulmonary arterial mechanics.

OBJECTIVE

We performed a biomechanical and microstructural characterization study of porcine pulmonary arteries, inclusive of the main, left, and right pulmonary arteries (MPA, LPA, and RPA, respectively).

METHODS

The specimens were initially stored at -20°C and allowed to thaw for 12-24 hours prior to testing. Each artery was further subdivided into proximal, middle, and distal regions, leading to ten location-based experimental groups. Planar equibiaxial tensile testing was performed to evaluate the mechanical behavior of the specimens, from which we calculated the stress at the maximum strain (S ₅₅), tensile modulus (TM), anisotropy index (AI), and strain energy in terms of area under the stress-strain curve (AUC). Histological quantification was performed to evaluate the area fraction of elastin and collagen content, intima-media thickness (IMT), and adventitial thickness (AT). The constitutive material behavior of each group was represented by a five-constant Holzapfel-Gasser-Ogden model.

RESULTS

The specimens exhibited non-linear stress-strain characteristics across all groups. The MPA exhibited the highest mean wall stress and TM in the longitudinal and circumferential directions, while the bifurcation region yielded the highest values of AI and AUC. All regions revealed a higher stiffness in the longitudinal direction compared to the circumferential direction, suggesting a degree of anisotropy that is believed to be within the margin of experimental uncertainty. Collagen content was found to be the highest in the MPA and decreased significantly at the bifurcation, LPA and RPA. Elastin content did not yield such significant differences amongst the ten groups. The MPA had the highest IMT, which decreased concomitantly to the distal LPA and RPA. No significant differences were found in the AT amongst the ten groups.

CONCLUSION

The mechanical properties of porcine pulmonary arteries exhibit strong regional dissimilarities, which can be used to inform future studies of high fidelity finite element models.

3D printed biaxial stretcher compatible with live fluorescence microscopy

Mechanical characterization and tensile testing of biological samples is important when determining the material properties of a tissue; however, performing tensile testing and tissue stretching while monitoring cellular changes via fluorescence microscopy is often challenging. Additionally, commercially available cell/tissue stretchers are often expensive, hard to customize, and limited in their fluorescence imaging compatibility. We have developed a 3D printed Open source Biaxial Stretcher (OBS) to be a low-cost stage top mountable biaxial stretching system for use with live cell fluorescence microscopy in both upright and inverted microscope configurations. Our OBS takes advantage of readily available open source desktop 3D printer hardware and software to deliver a fully motorized high precision (10 ± 0.5 μm movement accuracy) low cost biaxial stretching device capable of 4.5 cm of XY travel with a touch screen control panel, and an integrated heated platform with sample bath to maintain cell and tissue viability. Further, we designed a series of tissue mounts and clamps to accommodate varying samples from synthetic materials to biological tissue. By creating a low-profile design, we can directly mount the stretcher onto a microscope stage, and through coordinated biaxial stretching we maintain a constant field of view facilitating real-time sample tracking and time-lapse fluorescence imaging.

A computational multi-scale approach to investigate mechanically-induced changes in tricuspid valve anterior leaflet microstructure

The tricuspid valve is an atrioventricular valve that prevents blood backflow from the right ventricle into the right atrium during ventricular contractions. It is important to study mechanically induced microstructural alterations in the tricuspid valve leaflets, as this aids both in understanding valvular diseases and in the development of new engineered tissue replacements. The structure and composition of the extracellular matrix (ECM) fiber networks are closely tied to an overall biomechanical function of the tricuspid valve. In this study, we conducted experiments and implemented a multiscale modeling approach to predict ECM microstructural changes to tissue-level mechanical responses in a controlled loading environment. In particular, we characterized a sample of a porcine anterior leaflet at a macroscale using a biaxial mechanical testing method. We then generated a three-dimensional finite element model, to which computational representations of corresponding fiber networks were incorporated based on properties of the microstructural architecture obtained from small angle light scattering. Using five different biaxial boundary conditions, we performed iterative simulations to obtain model parameters with an overall R² value of 0.93. We observed that mechanical loading could markedly alter the underlying ECM architecture. For example, a relatively isotropic fiber network (with an anisotropy index value α of 28%) became noticeably more anisotropic (with an α of 40%) when it underwent mechanical loading. We also observed that the mechanical strain was distributed in a different manner at the ECM/fiber level as compared to the tissue level. The approach presented in this study has the potential to be implemented in pathophysiologically altered biomechanical and structural conditions and to bring insights into the mechanobiology of the tricuspid valve. STATEMENT OF SIGNIFICANCE: Quantifying abnormal cellar/ECM-level deformation of tricuspid valve leaflets subjected to a modified loading environment is of great importance, as it is believed to be linked to valvular remodeling responses. For example, developing surgical procedures or engineered tissue replacements that maintain/mimic ECM-level mechanical homeostasis could lead to more durable outcomes. To quantify leaflet deformation, we built a multiscale framework encompassing the contributions of disorganized ECM components and organized fibers, which can predict the behavior of the tricuspid valve leaflets under physiological loading conditions both at the tissue level and at the ECM level. In addition to future in-depth studies of tricuspid valve pathologies, our model can be used to characterize tissues in other valves of the heart.

Evaluation of an elastic decellularized tendon-derived scaffold for the vascular tissue engineering application

Due to the limited success rate of currently available vascular replacements, tissue engineering has received tremendous attention in recent years. A main challenge in the field of regenerative medicine is creating a mechanically functional tissue with a well-organized extracellular matrix, particularly of collagen and elastin. In this study, the native collagen scaffold derived from decellularized tendon sections, as a scaffold having the potential to be used for vascular tissue engineering applications, was studied. We showed that the elasticity of the scaffolds was improved when crosslinked with the bovine elastin. The effect of different concentrations of elastin on mechanical properties of the collagen scaffolds was evaluated of which 15% elastin concentration was selected for further analysis based on the results. Addition of 15% elastin to collagen scaffolds significantly decreased Young's modulus and the tensile stress at the maximum load and increased the tensile strain at the maximum load of the constructs as compared to those of the collagen scaffolds or control samples. Moreover, tubular elastin modified collagen scaffolds showed significantly higher burst pressure compared to the control samples. Smooth muscle cells and endothelial cells cultured on the elastin modified collagen scaffolds showed high viability (>80%) after 1, 3, and 7 days. Furthermore, the cells showed a high tendency to align with the collagen fibers within the scaffold and produced their own extracellular matrix over time. In conclusion, the results show that the decellularized tendon sections have a great potential to be used as scaffolds for vascular tissue engineering applications. © 2019 Wiley Periodicals, Inc. J Biomed Mater Res Part A: 107A: 1225-1234, 2019.

Collagen fiber recruitment: A microstructural basis for the nonlinear response of the posterior pole of the eye to increases in intraocular pressure

Our goal was to quantify and characterize how the collagen fiber crimp waviness of the lamina cribrosa (LC) and peripapillary sclera (PPS) changes with intraocular pressure (IOP). Thirteen sheep (ovine) eyes were immersion and perfusion fixed while maintaining IOP at 0, 10, 15, 20, or 50 mmHg. Coronal optic nerve head (ONH) sections (30 µm) were imaged with polarized light microscopy (PLM) and analyzed for collagen fiber orientation and waviness (SD of fiber orientation). In the LC, the waviness of every LC beam was measured. In the PPS, at least 900 collagen bundles were measured per eye. Using the waviness at 50 mmHg IOP, we defined tissue-specific thresholds to determine the fraction of loaded or recruited fibers. We found that fiber waviness decreased with IOP (P < 0.001). At every IOP, the waviness of the collagen fibers, and the fraction of fibers recruited in the LC were smaller or equal than those of the PPS (P < 0.001). At 15 mmHg IOP, both LC and PPS had ¾ recruited fibers and ¼ left in reserve. The decreased waviness with IOP and associated fiber recruitment is experimental evidence of fiber-based nonlinear biomechanical behavior of the ONH. At all IOPs the PPS had an equal or larger fraction of fibers recruited than the LC. That both LC and PPS had the same fraction of recruited and reserve fibers at normal IOP suggests that this may be an optimal fraction of recruitment for healthy eyes. Whether this extends to human eyes remains unknown.

STATEMENT OF SIGNIFICANCE

Collagen fibers exhibit a natural waviness or crimp that largely determine the nonlinear biomechanics of soft tissue. Experimental measurements of crimp morphology in the sheep eye, and how it changes with intraocular pressure (IOP), however, are exceedingly sparse. We quantified how posterior eye crimp changes with increasing IOP. We found that the lamina cribrosa and peripapillary sclera have fundamentally different crimp, and with increasing IOP, have different proportions of fibers that straighten, or get recruited, versus remaining wavy, or in reserve. Interestingly, at physiologic IOP of 15 mmHg, both tissues had about ¾ fibers recruited and ¼ fibers in reserve, indicating there may be an optimal fraction of fibers.

Structural and Functional Differences Between Porcine Aorta and Vena Cava

Elastin and collagen fibers are the major load-bearing extracellular matrix (ECM) constituents of the vascular wall. Arteries function differently than veins in the circulatory system; however as a result from several treatment options, veins are subjected to sudden elevated arterial pressure. It is thus important to recognize the fundamental structure and function differences between a vein and an artery. Our research compared the relationship between biaxial mechanical function and ECM structure of porcine thoracic aorta and inferior vena cava. Our study suggests that aorta contains slightly more elastin than collagen due to the cyclical extensibility, but vena cava contains almost four times more collagen than elastin to maintain integrity. Furthermore, multiphoton imaging of vena cava showed longitudinally oriented elastin and circumferentially oriented collagen that is recruited at supraphysiologic stress, but low levels of strain. However in aorta, elastin is distributed uniformly, and the primarily circumferentially oriented collagen is recruited at higher levels of strain than vena cava. These structural observations support the functional finding that vena cava is highly anisotropic with the longitude being more compliant and the circumference stiffening substantially at low levels of strain. Overall, our research demonstrates that fiber distributions and recruitment should be considered in addition to relative collagen and elastin contents. Also, the importance of accounting for the structural and functional differences between arteries and veins should be taken into account when considering disease treatment options.

Experimental Study of Anisotropic Stress/Strain Relationships of Aortic and Pulmonary Artery Homografts and Synthetic Vascular Grafts

Homografts and synthetic grafts are used in surgery for congenital heart disease (CHD). Determining these materials' mechanical properties will aid in understanding tissue behavior when subjected to abnormal CHD hemodynamics. Homograft tissue samples from anterior/posterior aspects, of ascending/descending aorta (AA, DA), innominate artery (IA), left subclavian artery (LScA), left common carotid artery (LCCA), main/left/right pulmonary artery (MPA, LPA, RPA), and synthetic vascular grafts, were obtained in three orientations: circumferential, diagonal (45 deg relative to circumferential direction), and longitudinal. Samples were subjected to uniaxial tensile testing (UTT). True strain-Cauchy stress curves were individually fitted for each orientation to calibrate Fung model. Then, they were used to calibrate anisotropic Holzapfel–Gasser model (R2 > 0.95). Most samples demonstrated a nonlinear hyperelastic strain–stress response to UTT. Stiffness (measured by tangent modulus at different strains) in all orientations were compared and shown as contour plots. For each vessel segment at all strain levels, stiffness was not significantly different among aspects and orientations. For synthetic grafts, stiffness was significantly different among orientations (p < 0.042). Aorta is significantly stiffer than pulmonary artery at 10% strain, comparing all orientations, aspects, and regions (p = 0.0001). Synthetic grafts are significantly stiffer than aortic and pulmonary homografts at all strain levels (p < 0.046). Aortic, pulmonary artery, and synthetic grafts exhibit hyperelastic biomechanical behavior with anisotropic effect. Differences in mechanical properties among vascular grafts may affect native tissue behavior and ventricular/arterial mechanical coupling, and increase the risk of deformation due to abnormal CHD hemodynamics.

Transmural variation in elastin fiber orientation distribution in the arterial wall

The complex three-dimensional elastin network is a major load-bearing extracellular matrix (ECM) component of an artery. Despite the reported anisotropic behavior of arterial elastin network, it is usually treated as an isotropic material in constitutive models. Our recent multiphoton microscopy study reported a relatively uniform elastin fiber orientation distribution in porcine thoracic aorta when imaging from the intima side (Chow et al., 2014). However it is questionable whether the fiber orientation distribution obtained from a small depth is representative of the elastin network structure in the arterial wall, especially when developing structure-based constitutive models. To date, the structural basis for the anisotropic mechanical behavior of elastin is still not fully understood. In this study, we examined the transmural variation in elastin fiber orientation distribution in porcine thoracic aorta and its association with elastin anisotropy. Using multi-photon microscopy, we observed that the elastin fibers orientation changes from a relatively uniform distribution in regions close to the luminal surface to a more circumferential distribution in regions that dominate the media, then to a longitudinal distribution in regions close to the outer media. Planar biaxial tensile test was performed to characterize the anisotropic behavior of elastin network. A new structure-based constitutive model of elastin network was developed to incorporate the transmural variation in fiber orientation distribution. The new model well captures the anisotropic mechanical behavior of elastin network under both equi- and nonequi-biaxial loading and showed improvements in both fitting and predicting capabilities when compared to a model that only considers the fiber orientation distribution from the intima side. We submit that the transmural variation in fiber orientation distribution is important in characterizing the anisotropic mechanical behavior of elastin network and should be considered in constitutive modeling of an artery.

The Coupled Bio-Chemo-Electro-Mechanical Behavior of Glucose Exposed Arterial Elastin

Elastin, the principle protein component of the elastic fiber, is a critical extracellular matrix (ECM) component of the arterial wall providing structural resilience and biological signaling essential in vascular morphogenesis and maintenance of mechanical homeostasis. Pathogenesis of many cardiovascular diseases have been associated with alterations of elastin. As a long-lived ECM protein that is deposited and organized before adulthood, elastic fibers can suffer from cumulative effects of biochemical exposure encountered during aging and/or disease, which greatly compromise their mechanical function. This review article covers findings from recent studies of the mechanical and structural contribution of elastin to vascular function, and the effects of biochemical degradation. Results from diverse experimental methods including tissue-level mechanical characterization, fiber-level nonlinear optical imaging, piezoelectric force microscopy, and nuclear magnetic resonance are reviewed. The intriguing coupled bio-chemo-electro-mechanical behavior of elastin calls for a multi-scale and multi-physical understanding of ECM mechanics and mechanobiology in vascular remodeling.

Microstructure and mechanics of human resistance arteries

Vascular diseases such as diabetes and hypertension cause changes to the vasculature that can lead to vessel stiffening and the loss of vasoactivity. The microstructural bases of these changes are not presently fully understood. We present a new methodology for stain-free visualization, at a microscopic scale, of the morphology of the main passive components of the walls of unfixed resistance arteries and their response to changes in transmural pressure. Human resistance arteries were dissected from subcutaneous fat biopsies, mounted on a perfusion myograph, and imaged at varying transmural pressures using a multimodal nonlinear microscope. High-resolution three-dimensional images of elastic fibers, collagen, and cell nuclei were constructed. The honeycomb structure of the elastic fibers comprising the internal elastic layer became visible at a transmural pressure of 30 mmHg. The adventitia, comprising wavy collagen fibers punctuated by straight elastic fibers, thinned under pressure as the collagen network straightened and pulled taut. Quantitative measurements of fiber orientation were made as a function of pressure. A multilayer analytical model was used to calculate the stiffness and stress in each layer. The adventitia was calculated to be up to 10 times as stiff as the media and experienced up to 8 times the stress, depending on lumen diameter. This work reveals that pressure-induced reorganization of fibrous proteins gives rise to very high local strain fields and highlights the unique mechanical roles of both fibrous networks. It thereby provides a basis for understanding the micromechanical significance of structural changes that occur with age and disease.

Biomechanical Behavior of Bioprosthetic Heart Valve Heterograft Tissues: Characterization, Simulation, and Performance

The use of replacement heart valves continues to grow due to the increased prevalence of valvular heart disease resulting from an ageing population. Since bioprosthetic heart valves (BHVs) continue to be the preferred replacement valve, there continues to be a strong need to develop better and more reliable BHVs through and improved the general understanding of BHV failure mechanisms. The major technological hurdle for the lifespan of the BHV implant continues to be the durability of the constituent leaflet biomaterials, which if improved can lead to substantial clinical impact. In order to develop improved solutions for BHV biomaterials, it is critical to have a better understanding of the inherent biomechanical behaviors of the leaflet biomaterials, including chemical treatment technologies, the impact of repetitive mechanical loading, and the inherent failure modes. This review seeks to provide a comprehensive overview of these issues, with a focus on developing insight on the mechanisms of BHV function and failure. Additionally, this review provides a detailed summary of the computational biomechanical simulations that have been used to inform and develop a higher level of understanding of BHV tissues and their failure modes. Collectively, this information should serve as a tool not only to infer reliable and dependable prosthesis function, but also to instigate and facilitate the design of future bioprosthetic valves and clinically impact cardiology.

Multi-scale Modeling of the Cardiovascular System: Disease Development, Progression, and Clinical Intervention

Cardiovascular diseases (CVDs) are the leading cause of death in the western world. With the current development of clinical diagnostics to more accurately measure the extent and specifics of CVDs, a laudable goal is a better understanding of the structure-function relation in the cardiovascular system. Much of this fundamental understanding comes from the development and study of models that integrate biology, medicine, imaging, and biomechanics. Information from these models provides guidance for developing diagnostics, and implementation of these diagnostics to the clinical setting, in turn, provides data for refining the models. In this review, we introduce multi-scale and multi-physical models for understanding disease development, progression, and designing clinical interventions. We begin with multi-scale models of cardiac electrophysiology and mechanics for diagnosis, clinical decision support, personalized and precision medicine in cardiology with examples in arrhythmia and heart failure. We then introduce computational models of vasculature mechanics and associated mechanical forces for understanding vascular disease progression, designing clinical interventions, and elucidating mechanisms that underlie diverse vascular conditions. We conclude with a discussion of barriers that must be overcome to provide enhanced insights, predictions, and decisions in pre-clinical and clinical applications.

Microstructure and Mechanical Property of Glutaraldehyde-Treated Porcine Pulmonary Ligament

There is a significant need for fixed biological tissues with desired structural and material constituents for tissue engineering applications. Here, we introduce the lung ligament as a fixed biological material that may have clinical utility for tissue engineering. To characterize the lung tissue for potential clinical applications, we studied glutaraldehyde-treated porcine pulmonary ligament (n = 11) with multiphoton microscopy (MPM) and conducted biaxial planar experiments to characterize the mechanical property of the tissue. The MPM imaging revealed that there are generally two families of collagen fibers distributed in two distinct layers: The first family largely aligns along the longitudinal direction with a mean angle of θ = 10.7 ± 9.3 deg, while the second one exhibits a random distribution with a mean θ = 36.6 ± 27.4. Elastin fibers appear in some intermediate sublayers with a random orientation distribution with a mean θ = 39.6 ± 23 deg. Based on the microstructural observation, a microstructure-based constitutive law was proposed to model the elastic property of the tissue. The material parameters were identified by fitting the model to the biaxial stress-strain data of specimens, and good fitting quality was achieved. The parameter e0 (which denotes the strain beyond which the collagen can withstand tension) of glutaraldehyde-treated tissues demonstrated low variability implying a relatively consistent collagen undulation in different samples, while the stiffness parameters for elastin and collagen fibers showed relatively greater variability. The fixed tissues presented a smaller e0 than that of fresh specimen, confirming that glutaraldehyde crosslinking increases the mechanical strength of collagen-based biomaterials. The present study sheds light on the biomechanics of glutaraldehyde-treated porcine pulmonary ligament that may be a candidate for tissue engineering.

Collagen Matrix Remodeling in Stented Pulmonary Arteries after Transapical Heart Valve Replacement

The use of valved stents for minimally invasive replacement of semilunar heart valves is expected to change the extracellular matrix and mechanical function of the native artery and may thus impair long-term functionality of the implant. Here we investigate the impact of the stent on matrix remodeling of the pulmonary artery in a sheep model, focusing on matrix composition and collagen (re)orientation of the host tissue. Ovine native pulmonary arteries were harvested 8 (n = 2), 16 (n = 4) and 24 (n = 2) weeks after transapical implantation of self-expandable stented heart valves. Second harmonic generation (SHG) microscopy was used to assess the collagen (re)orientation of fresh tissue samples. The collagen and elastin content was quantified using biochemical assays. SHG microscopy revealed regional differences in collagen organization in all explants. In the adventitial layer of the arterial wall far distal to the stent (considered as the control tissue), we observed wavy collagen fibers oriented in the circumferential direction. These circumferential fibers were more straightened in the adventitial layer located behind the stent. On the luminal side of the wall behind the stent, collagen fibers were aligned along the stent struts and randomly oriented between the struts. Immediately distal to the stent, however, fibers on both the luminal and the adventitial side of the wall were oriented in the axial direction, demonstrating the stent impact on the collagen structure of surrounding arterial tissues. Collagen orientation patterns did not change with implantation time, and biochemical analyses showed no changes in the trend of collagen and elastin content with implantation time or location of the vascular wall. We hypothesize that the collagen fibers on the adventitial side of the arterial wall and behind the stent straighten in response to the arterial stretch caused by oversizing of the stent. However, the collagen organization on the luminal side suggests that stent-induced remodeling is dominated by contact guidance.

Microstructural quantification of collagen fiber orientations and its integration in constitutive modeling of the porcine carotid artery

BACKGROUND

Mechanical characteristics of vascular tissue may play a role in different arterial pathologies, which, amongst others, requires robust constitutive descriptions to capture the vessel wall's anisotropic and non-linear properties.Specifically, the complex 3D network of collagen and its interaction with other structural elements has a dominating effect of arterial properties at higher stress levels.The aim of this study is to collect quantitative collagen organization as well as mechanical properties to facilitate structural constitutive models for the porcine carotid artery.This helps the understanding of the mechanics of swine carotid arteries, being a standard in clinical hypothesis testing, in endovascular preclinical trials for example.

METHOD

Porcine common carotid arteries (n=10) were harvested and used to (i) characterize the collagen fiber organization with polarized light microscopy, and (ii) the biaxial mechanical properties by inflation testing.The collagen organization was quantified by the Bingham orientation density function (ODF), which in turn was integrated in a structural constitutive model of the vessel wall.A one-layered and thick-walled model was used to estimate mechanical constitutive parameters by least-square fitting the recorded in vitro inflation test results.Finally, uniaxial data published elsewhere were used to validate the mean collagen organization described by the Bingham ODF.

RESULTS

Thick collagen fibers, i.e.the most mechanically relevant structure, in the common carotid artery are dispersed around the circumferential direction.In addition, almost all samples showed two distinct families of collagen fibers at different elevation, but not azimuthal, angles.Collagen fiber organization could be accurately represented by the Bingham ODF (κ1,2,3=[13.5,0.0,25.2] and κ1,2,3=[14.7,0.0,26.6]; average error of about 5%), and their integration into a structural constitutive model captured the inflation characteristics of individual carotid artery samples.Specifically, only four mechanical parameters were required to reasonably (average error from 14% to 38%) cover the experimental data over a wide range of axial and circumferential stretches.However, it was critical to account for fibrilar links between thick collagen fibers.Finally, the mean Bingham ODF provide also good approximation to uniaxial experimental data.

CONCLUSIONS

The applied structural constitutive model, based on individually measured collagen orientation densities, was able to capture the biaxial properties of the common carotid artery. Since the model required coupling amongst thick collagen fibers, the collagen fiber orientations measured from polarized light microscopy, alone, seem to be insufficient structural information. Alternatively, a larger dispersion of collagen fiber orientations, that is likely to arise from analyzing larger wall sections, could have had a similar effect, i.e. could have avoided coupling amongst thick collagen fibers.

STATEMENT OF SIGNIFICANCE

The applied structural constitutive model, based on individually measured collagen orientation densities, was able to capture the biaxial and uniaxial properties of the common carotid artery. Since the model required coupling amongst thick collagen fibers, an effective orientation density that accounts for cross-links between the main collagen fibers has been porposed. The model provides a good approximation to the experimental data.

Arterial mechanics considering the structural and mechanical contributions of ECM constituents

Considering the organization and engagement behavior of different extracellular matrix (ECM) constituents in the medial and adventitial layer of the arterial wall, in this study, we proposed a new constitutive model of ECM mechanics that considers the distinct structural and mechanical contributions of medial elastin, medial collagen, and adventitial collagen, to incorporate the constituent-specific fiber orientation and the sequential fiber engagement in arterial mechanics. Planar biaxial tensile testing method was used to characterize the orthotropic and hyperelastic behavior of porcine thoracic aorta. Fiber distribution functions of medial elastin, medial collagen, and adventitial collagen were incorporated into the constitutive model. Considering the sequential engagement of ECM constituents in arterial mechanics, a recruitment density function was incorporated into the model to capture the delayed engagement of adventitial collagen. A freely jointed chain model was used to capture the mechanical behavior of elastin and collagen at the fiber level. The tissue-level ECM mechanics was obtained by incorporating fiber distribution, engagement, and elastin and collagen content. The multi-scale constitutive model considering the structural and mechanical contributions of the three major ECM constituents allows us to directly incorporate information obtained from quantitative multi-photon imaging and analysis, and biochemical assay for the prediction of tissue-level mechanical response. Moreover, the model shows promises in fitting and predicting with a small set of material parameters, which has physical meanings and can be related to the structure of the ECM.

Automatic Segmentation of Mechanically Inhomogeneous Tissues Based on Deformation Gradient Jump

Variations in properties, active behavior, injury, scarring, and/or disease can all cause a tissue's mechanical behavior to be heterogeneous. Advances in imaging technology allow for accurate full-field displacement tracking of both in vitro and in vivo deformation from an applied load. While detailed strain fields provide some insight into tissue behavior, material properties are usually determined by fitting stress-strain behavior with a constitutive equation. However, the determination of the mechanical behavior of heterogeneous soft tissue requires a spatially varying constitutive equation (i.e., one in which the material parameters vary with position). We present an approach that computationally dissects the sample domain into many homogeneous subdomains, wherein subdomain boundaries are formed by applying a betweenness based graphical analysis to the deformation gradient field to identify locations with large discontinuities. This novel partitioning technique successfully determined the shape, size and location of regions with locally similar material properties for: (1) a series of simulated soft tissue samples prescribed with both abrupt and gradual changes in anisotropy strength, prescribed fiber alignment, stiffness, and nonlinearity, (2) tissue analogs (PDMS and collagen gels) which were tested biaxially and speckle tracked (3) and soft tissues which exhibited a natural variation in properties (cadaveric supraspinatus tendon), a pathologic variation in properties (thoracic aorta containing transmural plaque), and active behavior (contracting cardiac sheet). The routine enables the dissection of samples computationally rather than physically, allowing for the study of small tissues specimens with unknown and irregular inhomogeneity.

Quantification and simulation of layer-specific mitral valve interstitial cells deformation under physiological loading

Within each of the four layers of mitral valve (MV) leaflet tissues there resides a heterogeneous population of interstitial cells that maintain the structural integrity of the MV tissue via protein biosynthesis and enzymatic degradation. There is increasing evidence that tissue stress-induced MV interstitial cell (MVIC) deformations can have deleterious effects on their biosynthetic states that are potentially related to the reduction of tissue-level maintenance and to subsequent organ-level failure. To better understand the interrelationships between tissue-level loading and cellular responses, we developed the following integrated experimental-computational approach. Since in vivo cellular deformations are not directly measurable, we quantified the in-situ layer-specific MVIC deformations for each of the four layers under a controlled biaxial tension loading device coupled to multi-photon microscopy. Next, we explored the interrelationship between the MVIC stiffness and deformation to layer-specific tissue mechanical and structural properties using a macro-micro finite element computational model. Experimental results indicated that the MVICs in the fibrosa and ventricularis layers deformed significantly more than those in the atrialis and spongiosa layers, reaching a nucleus aspect ratio of 3.3 under an estimated maximum physiological tension of 150N/m. The simulated MVIC moduli for the four layers were found to be all within a narrow range of 4.71-5.35kPa, suggesting that MVIC deformation is primarily controlled by each tissue layer's respective structure and mechanical behavior rather than the intrinsic MVIC stiffness. This novel result further suggests that while the MVICs may be phenotypically and biomechanically similar throughout the leaflet, they experience layer-specific mechanical stimulatory inputs due to distinct extracellular matrix architecture and mechanical behaviors of the four MV leaflet tissue layers. This also suggests that MVICs may behave in a layer-specific manner in response to mechanical stimuli in both normal and surgically modified MVs.

On the presence of affine fibril and fiber kinematics in the mitral valve anterior leaflet

In this study, we evaluated the hypothesis that the constituent fibers follow an affine deformation kinematic model for planar collagenous tissues. Results from two experimental datasets were utilized, taken at two scales (nanometer and micrometer), using mitral valve anterior leaflet (MVAL) tissues as the representative tissue. We simulated MVAL collagen fiber network as an ensemble of undulated fibers under a generalized two-dimensional deformation state, by representing the collagen fibrils based on a planar sinusoidally shaped geometric model. The proposed approach accounted for collagen fibril amplitude, crimp period, and rotation with applied macroscopic tissue-level deformation. When compared to the small angle x-ray scattering measurements, the model fit the data well, with an r(2) = 0.976. This important finding suggests that, at the homogenized tissue-level scale of ∼1 mm, the collagen fiber network in the MVAL deforms according to an affine kinematics model. Moreover, with respect to understanding its function, affine kinematics suggests that the constituent fibers are largely noninteracting and deform in accordance with the bulk tissue. It also suggests that the collagen fibrils are tightly bounded and deform as a single fiber-level unit. This greatly simplifies the modeling efforts at the tissue and organ levels, because affine kinematics allows a straightforward connection between the macroscopic and local fiber strains. It also suggests that the collagen and elastin fiber networks act independently of each other, with the collagen and elastin forming long fiber networks that allow for free rotations. Such freedom of rotation can greatly facilitate the observed high degree of mechanical anisotropy in the MVAL and other heart valves, which is essential to heart valve function. These apparently novel findings support modeling efforts directed toward improving our fundamental understanding of tissue biomechanics in healthy and diseased conditions.

Structure-based constitutive model can accurately predict planar biaxial properties of aortic wall tissue

Structure-based constitutive models might help in exploring mechanisms by which arterial wall histology is linked to wall mechanics. This study aims to validate a recently proposed structure-based constitutive model. Specifically, the model's ability to predict mechanical biaxial response of porcine aortic tissue with predefined collagen structure was tested. Histological slices from porcine thoracic aorta wall (n=9) were automatically processed to quantify the collagen fiber organization, and mechanical testing identified the non-linear properties of the wall samples (n=18) over a wide range of biaxial stretches. Histological and mechanical experimental data were used to identify the model parameters of a recently proposed multi-scale constitutive description for arterial layers. The model predictive capability was tested with respect to interpolation and extrapolation. Collagen in the media was predominantly aligned in circumferential direction (planar von Mises distribution with concentration parameter bM=1.03 ± 0.23), and its coherence decreased gradually from the luminal to the abluminal tissue layers (inner media, b=1.54 ± 0.40; outer media, b=0.72 ± 0.20). In contrast, the collagen in the adventitia was aligned almost isotropically (bA=0.27 ± 0.11), and no features, such as families of coherent fibers, were identified. The applied constitutive model captured the aorta biaxial properties accurately (coefficient of determination R(2)=0.95 ± 0.03) over the entire range of biaxial deformations and with physically meaningful model parameters. Good predictive properties, well outside the parameter identification space, were observed (R(2)=0.92 ± 0.04). Multi-scale constitutive models equipped with realistic micro-histological data can predict macroscopic non-linear aorta wall properties. Collagen largely defines already low strain properties of media, which explains the origin of wall anisotropy seen at this strain level. The structure and mechanical properties of adventitia are well designed to protect the media from axial and circumferential overloads.

Contribution of collagen fiber undulation to regional biomechanical properties along porcine thoracic aorta

As major extracellular matrix components, elastin, and collagen play crucial roles in regulating the mechanical properties of the aortic wall and, thus, the normal cardiovascular function. The mechanical properties of aorta, known to vary with age and multitude of diseases as well as the proximity to the heart, have been attributed to the variations in the content and architecture of wall constituents. This study is focused on the role of layer-specific collagen undulation in the variation of mechanical properties along the porcine descending thoracic aorta. Planar biaxial tensile tests are performed to characterize the hyperelastic anisotropic mechanical behavior of tissues dissected from four locations along the thoracic aorta. Multiphoton microscopy is used to image the associated regional microstructure. Exponential-based and recruitment-based constitutive models are used to account for the observed mechanical behavior while considering the aortic wall as a composite of two layers with independent properties. An elevated stiffness is observed in distal regions compared to proximal regions of thoracic aorta, consistent with sharper and earlier collagen recruitment estimated for medial and adventitial layers in the models. Multiphoton images further support our prediction that higher stiffness in distal regions is associated with less undulation in collagen fibers. Recruitment-based models further reveal that regardless of the location, collagen in the media is recruited from the onset of stretching, whereas adventitial collagen starts to engage with a delay. A parameter sensitivity analysis is performed to discriminate between the models in terms of the confidence in the estimated model parameters.

Arterial extracellular matrix: a mechanobiological study of the contributions and interactions of elastin and collagen

The complex network structure of elastin and collagen extracellular matrix (ECM) forms the primary load bearing components in the arterial wall. The structural and mechanobiological interactions between elastin and collagen are important for properly functioning arteries. Here, we examined the elastin and collagen organization, realignment, and recruitment by coupling mechanical loading and multiphoton imaging. Two-photon excitation fluorescence and second harmonic generation methods were performed with a multiphoton video-rate microscope to capture real time changes to the elastin and collagen structure during biaxial deformation. Enzymatic removal of elastin was performed to assess the structural changes of the remaining collagen structure. Quantitative analysis of the structural changes to elastin and collagen was made using a combination of two-dimensional fast Fourier transform and fractal analysis, which allows for a more complete understanding of structural changes. Our study provides new quantitative evidence, to our knowledge on the sequential engagement of different arterial ECM components in response to mechanical loading. The adventitial collagen exists as large wavy bundles of fibers that exhibit fiber engagement after 20% strain. The medial collagen is engaged throughout the stretching process, and prominent elastic fiber engagement is observed up to 20% strain after which the engagement plateaus. The fiber orientation distribution functions show remarkably different changes in the ECM structure in response to mechanical loading. The medial collagen shows an evident preferred circumferential distribution, however the fiber families of adventitial collagen are obscured by their waviness at no or low mechanical strains. Collagen fibers in both layers exhibit significant realignment in response to unequal biaxial loading. The elastic fibers are much more uniformly distributed and remained relatively unchanged due to loading. Removal of elastin produces similar structural changes in collagen as mechanical loading. Our study suggests that the elastic fibers are under tension and impart an intrinsic compressive stress on the collagen.

Insights into regional adaptations in the growing pulmonary artery using a meso-scale structural model: effects of ascending aorta impingement

As the next step in our investigations into the structural adaptations of the main pulmonary artery (PA) during postnatal growth, we utilized the extensive experimental measurements of the growing ovine PA from our previous study (Fata et al., 2013, "Estimated in vivo Postnatal Surface Growth Patterns of the Ovine Main Pulmonary Artery and Ascending Aorta," J. Biomech. Eng., 135(7), pp. 71010-71012). to develop a structural constitutive model for the PA wall tissue. Novel to the present approach was the treatment of the elastin network as a distributed fiber network rather than a continuum phase. We then utilized this model to delineate structure-function differences in the PA wall at the juvenile and adult stages. Overall, the predicted elastin moduli exhibited minor differences remained largely unchanged with age and region (in the range of 150 to 200 kPa). Similarly, the predicted collagen moduli ranged from ∼1,600 to 2700 kPa in the four regions studied in the juvenile state. Interestingly, we found for the medial region that the elastin and collagen fiber splay underwent opposite changes (collagen standard deviation juvenile = 17 deg to adult = 28 deg, elastin standard deviation juvenile = 35 deg to adult = 27 deg), along with a trend towards more rapid collagen fiber strain recruitment with age, along with a drop in collagen fiber moduli, which went from 2700 kPa for the juvenile stage to 746 kPa in the adult. These changes were likely due to the previously observed impingement of the relatively stiff ascending aorta on the growing PA medial region. Intuitively, the effects of the local impingement would be to lower the local wall stress, consistent with the observed parallel decrease in collagen modulus. These results suggest that during the postnatal somatic growth period local stresses can substantially modulate regional tissue microstructure and mechanical behaviors in the PA. We further underscore that our previous studies indicated an increase in effective PA wall stress with postnatal maturation. When taken together with the fact that the observed changes in mechanical behavior and structure in the growing PA wall were modest in the other three regions studied, our collective results suggest that the majority of the growing PA wall is subjected to increasing stress levels with age without undergoing major structural adaptations. This observation is contrary to the accepted theory of maintenance of homeostatic stress levels in the regulation of vascular function, and suggests alternative mechanisms might regulate postnatal somatic growth. Understanding the underlying mechanisms will help to improve our understanding of congenital defects of the PA and lay the basis for functional duplication in their repair and replacement.

Simulation of planar soft tissues using a structural constitutive model: Finite element implementation and validation

Computational implementation of physical and physiologically realistic constitutive models is critical for numerical simulation of soft biological tissues in a variety of biomedical applications. It is well established that the highly nonlinear and anisotropic mechanical behaviors of soft tissues are an emergent behavior of the underlying tissue microstructure. In the present study, we have implemented a structural constitutive model into a finite element framework specialized for membrane tissues. We noted that starting with a single element subjected to uniaxial tension, the non-fibrous tissue matrix must be present to prevent unrealistic tissue deformations. Flexural simulations were used to set the non-fibrous matrix modulus because fibers have little effects on tissue deformation under three-point bending. Multiple deformation modes were simulated, including strip biaxial, planar biaxial with two attachment methods, and membrane inflation. Detailed comparisons with experimental data were undertaken to insure faithful simulations of both the macro-level stress-strain insights into adaptations of the fiber architecture under stress, such as fiber reorientation and fiber recruitment. Results indicated a high degree of fidelity and demonstrated interesting microstructural adaptions to stress and the important role of the underlying tissue matrix. Moreover, we apparently resolve a discrepancy in our 1997 study (Billiar and Sacks, 1997. J. Biomech. 30 (7), 753-756) where we observed that under strip biaxial stretch the simulated fiber splay responses were not in good agreement with the experimental results, suggesting non-affine deformations may have occurred. However, by correctly accounting for the isotropic phase of the measured fiber splay, good agreement was obtained. While not the final word, these simulations suggest that affine fiber kinematics for planar collagenous tissues is a reasonable assumption at the macro level. Simulation tools such as these are imperative in the design and simulation of native and engineered tissues.

Estimated in vivo postnatal surface growth patterns of the ovine main pulmonary artery and ascending aorta

Delineating the normal postnatal development of the pulmonary artery (PA) and ascending aorta (AA) can inform our understanding of congenital abnormalities, as well as pulmonary and systolic hypertension. We thus conducted the following study to delineate the PA and AA postnatal growth deformation characteristics in an ovine model. MR images were obtained from endoluminal surfaces of 11 animals whose ages ranged from 1.5 months/15.3 kg mass (very young) to 12 months/56.6 kg mass (adult). A bicubic Hermite finite element surface representation was developed for the each artery from each animal. Under the assumption that the relative locations of surface points were retained during growth, the individual animal surface fits were subsequently used to develop a method to estimate the time-evolving local effective surface growth (relative to the youngest measured animal) in the end-diastolic state. Results indicated that the spatial and temporal surface growth deformation patterns of both arteries, especially in the circumferential direction, were heterogeneous, leading to an increase in taper and increase in cross-sectional ellipticity of the PA. The longitudinal PA growth stretch of a large segment on the posterior wall reached 2.57 ± 0.078 (mean ± SD) at the adult stage. In contrast, the longitudinal growth of the AA was smaller and more uniform (1.80 ± 0.047). Interestingly, a region of the medial wall of both arteries where both arteries are in contact showed smaller circumferential growth stretches-specifically 1.12 ± 0.012 in the PA and 1.43 ± 0.071 in the AA at the adult stage. Overall, our results indicated that contact between the PA and AA resulted in increasing spatial heterogeneity in postnatal growth, with the PA demonstrating the greatest changes. Parametric studies using simplified geometric models of curved arteries during growth suggest that heterogeneous effective surface growth deformations must occur to account for the changes in measured arterial shapes during the postnatal growth period. This result suggests that these first results are a reasonable first-approximation to the actual effective growth patterns. Moreover, this study clearly underscores how functional growth of the PA and AA during postnatal maturation involves complex, local adaptations in tissue formation. Moreover, the present results will help to lay the basis for functional replacement by defining critical geometric metrics.

Mechanical analysis of ovine and pediatric pulmonary artery for heart valve stent design

Transcatheter heart valve replacement is an attractive and promising technique for congenital as well as acquired heart valve disease. In this procedure, the replacement valve is mounted in a stent that is expanded at the aimed valve position and fixated by clamping. However, for this technique to be appropriate for pediatric patients, the material properties of the host tissue need to be determined to design stents that can be optimized for this particular application. In this study we performed equibiaxial tensile tests on four adult ovine pulmonary artery walls and compared the outcomes with one pediatric pulmonary artery. Results show that the pediatric pulmonary artery was significantly thinner (1.06 ± 0.36 mm (mean ± SD)) than ovine tissue (2.85 ± 0.40 mm), considerably stiffer for strain values that exceed the physiological conditions (beyond 50% strain in the circumferential and 60% in the longitudinal direction), more anisotropic (with a significant difference in stiffness between the longitudinal and circumferential directions beyond 60% strain) and presented stronger non-linear stress-strain behavior at equivalent strains (beyond 26% strain) compared to ovine tissue. These discrepancies suggest that stents validated and optimized using the ovine pre-clinical model might not perform satisfactorily in pediatric patients. The material parameters derived from this study may be used to develop stent designs for both applications using computational models.