Biaxial mechanical properties of the natural and glutaraldehyde treated aortic valve cusp--Part I: Experimental results

The influence of specimen thickness and alignment on the material and failure properties of electrospun polycaprolactone nanofiber mats

Electrospinning is a versatile fabrication technique that has been recently expanded to create nanofibrous structures that mimic ECM topography. Like many materials, electrospun constructs are typically characterized on a smaller scale, and scaled up for various applications. This established practice is based on the assumption that material properties, such as toughness, failure stress and strain, are intrinsic to the material, and thus will not be influenced by specimen geometry. However, we hypothesized that the material and failure properties of electrospun nanofiber mats vary with specimen thickness. To test this, we mechanically characterized polycaprolactone (PCL) nanofiber mats of three different thicknesses in response to constant rate elongation to failure. To identify if any observed thickness-dependence could be attributed to fiber alignment, such as the effects of fiber reorientation during elongation, these tests were performed in mats with either random or aligned nanofiber orientation. Contrary to our hypothesis, the failure strain was conserved across the different thicknesses, indicating similar maximal elongation for specimens of different thickness. However, in both the aligned and randomly oriented groups, the ultimate tensile stress, short-range modulus, yield modulus, and toughness all decreased with increasing mat thickness, thereby indicating that these are not intrinsic material properties. These findings have important implications in engineered scaffolds for fibrous and soft tissue applications (e.g., tendon, ligament, muscle, and skin), where such oversights could result in unwanted laxity or reduced resistance to failure. © 2016 Wiley Periodicals, Inc. J Biomed Mater Res Part A: 104A: 2794-2800, 2016.

Characterization of three-dimensional anisotropic heart valve tissue mechanical properties using inverse finite element analysis

Computational modeling has an important role in design and assessment of medical devices. In computational simulations, considering accurate constitutive models is of the utmost importance to capture mechanical response of soft tissue and biomedical materials under physiological loading conditions. Lack of comprehensive three-dimensional constitutive models for soft tissue limits the effectiveness of computational modeling in research and development of medical devices. The aim of this study was to use inverse finite element (FE) analysis to determine three-dimensional mechanical properties of bovine pericardial leaflets of a surgical bioprosthesis under dynamic loading condition. Using inverse parameter estimation, 3D anisotropic Fung model parameters were estimated for the leaflets. The FE simulations were validated using experimental in-vitro measurements, and the impact of different constitutive material models was investigated on leaflet stress distribution. The results of this study showed that the anisotropic Fung model accurately simulated the leaflet deformation and coaptation during valve opening and closing. During systole, the peak stress reached to 3.17MPa at the leaflet boundary while during diastole high stress regions were primarily observed in the commissures with the peak stress of 1.17MPa. In addition, the Rayleigh damping coefficient that was introduced to FE simulations to simulate viscous damping effects of surrounding fluid was determined.

Impact of different aortic valve calcification patterns on the outcome of transcatheter aortic valve implantation: A finite element study

Transcatheter aortic valve implantation (TAVI) can treat symptomatic patients with calcific aortic stenosis. However, the severity and distribution of the calcification of valve leaflets can impair the TAVI efficacy. Here we tackle this issue from a biomechanical standpoint, by finite element simulation of a widely adopted balloon-expandable TAVI in three models representing the aortic root with different scenarios of calcific aortic stenosis. We developed a modeling approach realistically accounting for aortic root pressurization and complex anatomy, detailed calcification patterns, and for the actual stent deployment through balloon-expansion. Numerical results highlighted the dependency on the specific calcification pattern of the "dog-boning" of the stent. Also, local stent distortions were associated with leaflet calcifications, and led to localized gaps between the TAVI stent and the aortic tissues, with potential implications in terms of paravalvular leakage. High stresses were found on calcium deposits, which may be a risk factor for stroke; their magnitude and the extent of the affected regions substantially increased for the case of an "arc-shaped" calcification, running from commissure to commissure. Moreover, high stresses due to the interaction between the aortic wall and the leaflet calcifications were computed in the annular region, suggesting an increased risk for annular damage. Our analyses suggest a relation between the alteration of the stresses in the native anatomical components and prosthetic implant with the presence and distribution of relevant calcifications. This alteration is dependent on the patient-specific features of the calcific aortic stenosis and may be a relevant indicator of suboptimal TAVI results.

POLARIZED SPATIAL FREQUENCY DOMAIN IMAGING OF HEART VALVE FIBER STRUCTURE

Our group previously introduced Polarized Spatial Frequency Domain Imaging (PSFDI), a wide-field, reflectance imaging technique which we used to empirically map fiber direction in porcine pulmonary heart valve leaflets (PHVL) without optical clearing or physical sectioning of the sample. Presented is an extended analysis of our PSFDI results using an inverse Mueller matrix model of polarized light scattering that allows additional maps of fiber orientation distribution, along with instrumentation permitting increased imaging speed for dynamic PHVL fiber measurements. We imaged electrospun fiber phantoms with PSFDI, and then compared these measurements to SEM data collected for the same phantoms. PHVL was then imaged and compared to results of the same leaflets optically cleared and imaged with small angle light scattering (SALS). The static PHVL images showed distinct regional variance of fiber orientation distribution, matching our SALS results. We used our improved imaging speed to observe bovine tendon subjected to dynamic loading using a biaxial stretching device. Our dynamic imaging experiment showed trackable changes in the fiber microstructure of biological tissue under loading. Our new PSFDI analysis model and instrumentation allows characterization of fiber structure within heart valve tissues (as validated with SALS measurements), along with imaging of dynamic fiber remodeling. The experimental data will be used as inputs to our constitutive models of PHVL tissue to fully characterize these tissues' elastic behavior, and has immediate application in determining the mechanisms of structural and functional failure in PHVLs used as bio-prosthetic implants.

A meso-scale layer-specific structural constitutive model of the mitral heart valve leaflets

Fundamental to developing a deeper understanding of pathophysiological remodeling in mitral valve (MV) disease is the development of an accurate tissue-level constitutive model. In the present work, we developed a novel meso-scale (i.e. at the level of the fiber, 10-100 μm in length scale) structural constitutive model (MSSCM) for MV leaflet tissues. Due to its four-layer structure, we focused on the contributions from the distinct collagen and elastin fiber networks within each tissue layer. Requisite collagen and elastin fibrous structural information for each layer were quantified using second harmonic generation microscopy and conventional histology. A comprehensive mechanical dataset was also used to guide model formulation and parameter estimation. Furthermore, novel to tissue-level structural constitutive modeling approaches, we allowed the collagen fiber recruitment function to vary with orientation. Results indicated that the MSSCM predicted a surprisingly consistent mean effective collagen fiber modulus of 162.72 MPa, and demonstrated excellent predictive capability for extra-physiological loading regimes. There were also anterior-posterior leaflet-specific differences, such as tighter collagen and elastin fiber orientation distributions (ODF) in the anterior leaflet, and a thicker and stiffer atrialis in the posterior leaflet. While a degree of angular variance was observed, the tight valvular tissue ODF also left little room for any physically meaningful angular variance in fiber mechanical responses. Finally, a novel fibril-level (0.1-1 μm) validation approach was used to compare the predicted collagen fiber/fibril mechanical behavior with extant MV small angle X-ray scattering data. Results demonstrated excellent agreement, indicating that the MSSCM fully captures the tissue-level function. Future utilization of the MSSCM in computational models of the MV will aid in producing highly accurate simulations in non-physiological loading states that can occur in repair situations, as well as guide the form of simplified models for real-time simulation tools.

Evaluation of bioprosthetic heart valve failure using a matrix-fibril shear stress transfer approach

A matrix-fibril shear stress transfer approach is devised and developed in this paper to analyse the primary biomechanical factors which initiate the structural degeneration of the bioprosthetic heart valves (BHVs). Using this approach, the critical length of the collagen fibrils l c and the interface shear acting on the fibrils in both BHV and natural aortic valve (AV) tissues under physiological loading conditions are calculated and presented. It is shown that the required critical fibril length to provide effective reinforcement to the natural AV and the BHV tissue is l c  = 25.36 µm and l c  = 66.81 µm, respectively. Furthermore, the magnitude of the required shear force acting on fibril interface to break a cross-linked fibril in the BHV tissue is shown to be 38 µN, while the required interfacial force to break the bonds between the fibril and the surrounding extracellular matrix is 31 µN. Direct correlations are underpinned between these values and the ultimate failure strength and the failure mode of the BHV tissue compared with the natural AV, and are verified against the existing experimental data. The analyses presented in this paper explain the role of fibril interface shear and critical length in regulating the biomechanics of the structural failure of the BHVs, for the first time. This insight facilitates further understanding into the underlying causes of the structural degeneration of the BHVs in vivo.

Efficient computational simulation of actin stress fiber remodeling

Understanding collagen and stress fiber remodeling is essential for the development of engineered tissues with good functionality. These processes are complex, highly interrelated, and occur over different time scales. As a result, excessive computational costs are required to computationally predict the final organization of these fibers in response to dynamic mechanical conditions. In this study, an analytical approximation of a stress fiber remodeling evolution law was derived. A comparison of the developed technique with the direct numerical integration of the evolution law showed relatively small differences in results, and the proposed method is one to two orders of magnitude faster.

Age-dependent changes of stress and strain in the human heart valve and their relation with collagen remodeling

In order to create tissue-engineered heart valves with long-term functionality, it is essential to fully understand collagen remodeling during neo-tissue formation. Collagen remodeling is thought to maintain mechanical tissue homeostasis. Yet, the driving factor of collagen remodeling remains unidentified. In this study, we determined the collagen architecture and the geometric and mechanical properties of human native semilunar heart valves of fetal to adult age using confocal microscopy, micro-indentation and inverse finite element analysis. The outcomes were used to predict age-dependent changes in stress and stretch in the heart valves via finite element modeling. The results indicated that the circumferential stresses are different between the aortic and pulmonary valve, and, moreover, that the stress increases considerably over time in the aortic valve. Strikingly, relatively small differences were found in stretch with time and between the aortic and pulmonary valve, particularly in the circumferential direction, which is the main determinant of the collagen fiber stretch. Therefore, we suggest that collagen remodeling in the human heart valve maintains a stretch-driven homeostasis. Next to these novel insights, the unique human data set created in this study provides valuable input for the development of numerical models of collagen remodeling and optimization of tissue engineering.

STATEMENT OF SIGNIFICANCE

Annually, over 280,000 heart valve replacements are performed worldwide. Tissue engineering has the potential to provide valvular disease patients with living valve substitutes that can last a lifetime. Valve functionality is mainly determined by the collagen architecture. Hence, understanding collagen remodeling is crucial for creating tissue-engineered valves with long-term functionality. In this study, we determined the structural and material properties of human native heart valves of fetal to adult age to gain insight into the mechanical stimuli responsible for collagen remodeling. The age-dependent evolutionary changes in mechanical state of the native valve suggest that collagen remodeling in heart valves is a stretch-driven process.

Multiscale modelling of solid tumour growth: the effect of collagen micromechanics

Here we introduce a model of solid tumour growth coupled with a multiscale biomechanical description of the tumour microenvironment, which facilitates the explicit simulation of fibre-fibre and tumour-fibre interactions. We hypothesise that such a model, which provides a purely mechanical description of tumour-host interactions, can be used to explain experimental observations of the effect of collagen micromechanics on solid tumour growth. The model was specified to mouse tumour data, and numerical simulations were performed. The multiscale model produced lower stresses than an equivalent continuum-like approach, due to a more realistic remodelling of the collagen microstructure. Furthermore, solid tumour growth was found to cause a passive mechanical realignment of fibres at the tumour boundary from a random to a circumferential orientation. This is in accordance with experimental observations, thus demonstrating that such a response can be explained as purely mechanical. Finally, peritumoural fibre network anisotropy was found to produce anisotropic tumour morphology. The dependency of tumour morphology on the peritumoural microstructure was reduced by adding a load-bearing non-collagenous component to the fibre network constitutive equation.

Material modeling of cardiac valve tissue: Experiments, constitutive analysis and numerical investigation

A key element of the cardiac cycle of the human heart is the opening and closing of the four valves. However, the material properties of the leaflet tissues, which fundamentally contribute to determine the mechanical response of the valves, are still an open field of research. The main contribution of the present study is to provide a complete experimental data set for porcine heart valve samples spanning all valve and leaflet types under tensile loading. The tests show a fair degree of reproducibility and are clearly indicative of a number of fundamental tissue properties, including a progressively stiffening response with increasing elongation. We then propose a simple anisotropic constitutive model, which is fitted to the experimental data set, showing a reasonable interspecimen variability. Furthermore, we present a dynamic finite element analysis of the aortic valve to show the direct usability of the obtained material parameters in computational simulations.

Mechanical stress is associated with right ventricular response to pulmonary valve replacement in patients with repaired tetralogy of Fallot

OBJECTIVE

Patients with repaired tetralogy of Fallot account for a substantial proportion of cases with late-onset right ventricular failure. The current surgical approach, which includes pulmonary valve replacement/insertion, has yielded mixed results. Therefore, it may be clinically useful to identify parameters that can be used to predict right ventricular function response to pulmonary valve replacement.

METHODS

Cardiac magnetic resonance data before and 6 months after pulmonary valve replacement were obtained from 16 patients with repaired tetralogy of Fallot (8 male, 8 female; median age, 42.75 years). Right ventricular ejection fraction change from pre- to postpulmonary valve replacement was used as the outcome. The patients were divided into group 1 (n = 8, better outcome) and group 2 (n = 8, worst outcome). Cardiac magnetic resonance-based patient-specific computational right ventricular/left ventricular models were constructed, and right ventricular mechanical stress and strain, wall thickness, curvature, and volumes were obtained for analysis.

RESULTS

Our results indicated that right ventricular wall stress was the best single predictor for postpulmonary valve replacement outcome with an area under the receiver operating characteristic curve of 0.819. Mean values of stress, strain, wall thickness, and longitudinal curvature differed significantly between the 2 groups with right ventricular wall stress showing the largest difference. Mean right ventricular stress in group 2 was 103% higher than in group 1.

CONCLUSIONS

Computational modeling and right ventricular stress may be used as tools to identify right ventricular function response to pulmonary valve replacement. Large-scale clinical studies are needed to validate these preliminary findings.

Imaging analysis of collagen fiber networks in cusps of porcine aortic valves: effect of their local distribution and alignment on valve functionality

The cusps of native aortic valve (AV) are composed of collagen bundles embedded in soft tissue, creating a heterogenic tissue with asymmetric alignment in each cusp. This study compares native collagen fiber networks (CFNs) with a goal to better understand their influence on stress distribution and valve kinematics. Images of CFNs from five porcine tricuspid AVs are analyzed and fluid-structure interaction models are generated based on them. Although the valves had similar overall kinematics, the CFNs had distinctive influence on local mechanics. The regions with dilute CFN are more prone to damage since they are subjected to higher stress magnitudes.

In Vitro Model of a Fibrosa Layer of a Heart Valve

The fibrosa layer of a cardiac aortic valve is composed mostly of a dense network of type I collagen fibers oriented in circumferential direction. This main layer bears the tensile load and responds to the high stress on a leaflet. The inner fibrosa layer is also the site of pathophysiologic changes that result in valvular dysfunction, including stenosis and regurgitation. In vitro studies of these changes are limited by the absence of a substrate that mimics the circumferentially oriented structure of the fibrosa layer. In heart valve tissue engineering, generation of this layer is challenging. This study aimed to develop an artificial fibrosa layer of a native aortic leaflet. A unique morphologically biomimicked, pliable, but standalone substrate with circumferentially oriented nanofibers was fabricated by electrospinning on a novel collector designed for this study. The substrate had low-bulk tensile stiffness and ultimate strength; thus, cultured valvular interstitial cells (VICs) showed a fibroblast phenotype that is generally observed in a healthy aortic leaflet. Furthermore, gene and protein expression and morphology of VICs in substrates were close to those in the fibrosa layer of a native aortic leaflet. This artificial fibrosa layer can be useful for in vitro studies of valvular dysfunctions.

Emerging Trends in Heart Valve Engineering: Part IV. Computational Modeling and Experimental Studies

In this final portion of an extensive review of heart valve engineering, we focus on the computational methods and experimental studies related to heart valves. The discussion begins with a thorough review of computational modeling and the governing equations of fluid and structural interaction. We then move onto multiscale and disease specific modeling. Finally, advanced methods related to in vitro testing of the heart valves are reviewed. This section of the review series is intended to illustrate application of computational methods and experimental studies and their interrelation for studying heart valves.

Load-dependent extracellular matrix organization in atrioventricular heart valves: differences and similarities

The extracellular matrix of the atrioventricular (AV) valves' leaflets has a key role in the ability of these valves to properly remodel in response to constantly varying physiological loads. While the loading on mitral and tricuspid valves is significantly different, no information is available on how collagen fibers change their orientation in response to these loads. This study delineates the effect of physiological loading on AV valves' leaflets microstructures using Second Harmonic Generation (SHG) microscopy. Fresh natural porcine tricuspid and mitral valves' leaflets ( n = 12/valve type) were cut and prepared for the experiments. Histology and immunohistochemistry were performed to compare the microstructural differences between the valves. The specimens were imaged live during the relaxed, loading, and unloading phases using SHG microscopy. The images were analyzed with Fourier decomposition to mathematically seek changes in collagen fiber orientation. Despite the similarities in both AV valves as seen in the histology and immunohistochemistry data, the microstructural arrangement, especially the collagen fiber distribution and orientation in the stress-free condition, were found to be different. Uniaxial loading was dependent on the arrangement of the fibers in their relaxed mode, which led the fibers to reorient in-line with the load throughout the depth of the mitral leaflet but only to reorient in-line with the load in deeper layers of the tricuspid leaflet. Biaxial loading arranged the fibers in between the two principal axes of the stresses independently from their relaxed states. Unlike previous findings, this study concludes that the AV valves' three-dimensional extracellular fiber arrangement is significantly different in their stress-free and uniaxially loaded states; however, fiber rearrangement in response to the biaxial loading remains similar.

Three-dimensional macro-scale assessment of regional and temporal wall shear stress characteristics on aortic valve leaflets

The aortic valve (AV) achieves unidirectional blood flow between the left ventricle and the aorta. Although hemodynamic stresses have been shown to regulate valvular biology, the native wall shear stress (WSS) experienced by AV leaflets remains largely unknown. The objective of this study was to quantify computationally the macro-scale leaflet WSS environment using fluid-structure interaction modeling. An arbitrary Lagrangian-Eulerian approach was implemented to predict valvular flow and leaflet dynamics in a three-dimensional AV geometry subjected to physiologic transvalvular pressure. Local WSS characteristics were quantified in terms of temporal shear magnitude (TSM), oscillatory shear index (OSI) and temporal shear gradient (TSG). The dominant radial WSS predicted on the leaflets exhibited high amplitude and unidirectionality on the ventricularis (TSM>7.50 dyn/cm(2), OSI < 0.17, TSG>325.54 dyn/cm(2) s) but low amplitude and bidirectionality on the fibrosa (TSM < 2.73 dyn/cm(2), OSI>0.38, TSG < 191.17 dyn/cm(2) s). The radial WSS component computed in the leaflet base, belly and tip demonstrated strong regional variability (ventricularis TSM: 7.50-22.32 dyn/cm(2), fibrosa TSM: 1.26-2.73 dyn/cm(2)). While the circumferential WSS exhibited similar spatially dependent magnitude (ventricularis TSM: 1.41-3.40 dyn/cm(2), fibrosa TSM: 0.42-0.76 dyn/cm(2)) and side-specific amplitude (ventricularis TSG: 101.73-184.43 dyn/cm(2) s, fibrosa TSG: 41.92-54.10 dyn/cm(2) s), its temporal variations were consistently bidirectional (OSI>0.25). This study provides new insights into the role played by leaflet-blood flow interactions in valvular function and critical hemodynamic stress data for the assessment of the hemodynamic theory of AV disease.

Coupled Simulation of Heart Valves: Applications to Clinical Practice

The last few decades have seen great advances in the understanding of heart valves, and consequently, in the development of novel treatment modalities and surgical procedures for valves afflicted by disease. This is due in part to the profound advancements in computing technology and noninvasive medical imaging techniques that have made it possible to numerically model the complex heart valve systems characterized by distinct features at different length scales and various interacting processes. In this article, we highlight the importance of explicitly coupling these multiple scales and diverse processes to accurately simulate the true behavior of the heart valves, in health and disease. We examine some of the computational modeling studies that have a direct consequence on clinical practice.

Simple and accurate methods for quantifying deformation, disruption, and development in biological tissues

When mechanical factors underlie growth, development, disease or healing, they often function through local regions of tissue where deformation is highly concentrated. Current optical techniques to estimate deformation can lack precision and accuracy in such regions due to challenges in distinguishing a region of concentrated deformation from an error in displacement tracking. Here, we present a simple and general technique for improving the accuracy and precision of strain estimation and an associated technique for distinguishing a concentrated deformation from a tracking error. The strain estimation technique improves accuracy relative to other state-of-the-art algorithms by directly estimating strain fields without first estimating displacements, resulting in a very simple method and low computational cost. The technique for identifying local elevation of strain enables for the first time the successful identification of the onset and consequences of local strain concentrating features such as cracks and tears in a highly strained tissue. We apply these new techniques to demonstrate a novel hypothesis in prenatal wound healing. More generally, the analytical methods we have developed provide a simple tool for quantifying the appearance and magnitude of localized deformation from a series of digital images across a broad range of disciplines.

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.

A generalized method for the analysis of planar biaxial mechanical data using tethered testing configurations

Simulation of the mechanical behavior of soft tissues is critical for many physiological and medical device applications. Accurate mechanical test data is crucial for both obtaining the form and robust parameter determination of the constitutive model. For incompressible soft tissues that are either membranes or thin sections, planar biaxial mechanical testing configurations can provide much information about the anisotropic stress-strain behavior. However, the analysis of soft biological tissue planar biaxial mechanical test data can be complicated by in-plane shear, tissue heterogeneities, and inelastic changes in specimen geometry that commonly occur during testing. These inelastic effects, without appropriate corrections, alter the stress-traction mapping and violates equilibrium so that the stress tensor is incorrectly determined. To overcome these problems, we presented an analytical method to determine the Cauchy stress tensor from the experimentally derived tractions for tethered testing configurations. We accounted for the measured testing geometry and compensate for run-time inelastic effects by enforcing equilibrium using small rigid body rotations. To evaluate the effectiveness of our method, we simulated complete planar biaxial test configurations that incorporated actual device mechanisms, specimen geometry, and heterogeneous tissue fibrous structure using a finite element (FE) model. We determined that our method corrected the errors in the equilibrium of momentum and correctly estimated the Cauchy stress tensor. We also noted that since stress is applied primarily over a subregion bounded by the tethers, an adjustment to the effective specimen dimensions is required to correct the magnitude of the stresses. Simulations of various tether placements demonstrated that typical tether placements used in the current experimental setups will produce accurate stress tensor estimates. Overall, our method provides an improved and relatively straightforward method of calculating the resulting stresses for planar biaxial experiments for tethered configurations, which is especially useful for specimens that undergo large shear and exhibit substantial inelastic effects.

On the bending properties of porcine mitral, tricuspid, aortic, and pulmonary valve leaflets

The atrioventricular valve leaflets (mitral and tricuspid) are different from the semilunar valve leaflets (aortic and pulmonary) in layered structure, ultrastructural constitution and organization, and leaflet thickness. These differences warrant a comparative look at the bending properties of the four types of leaflets. We found that the moment-curvature relationships in atrioventricular valves were stiffer than in semilunar valves, and the moment-curvature relationships of the left-side valve leaflets were stiffer than their morphological analog of the right side. These trends were supported by the moment-curvature curves and the flexural rigidity analysis (EI value decreased from mitral, tricuspid, aortic, to pulmonary leaflets). However, after taking away the geometric effect (moment of inertia I), the instantaneous effective bending modulus E showed a reversed trend. The overall trend of flexural rigidity (EI: mitral > tricuspid > aortic > pulmonary) might be correlated with the thickness variations among the four types of leaflets (thickness: mitral > tricuspid > aortic > pulmonary). The overall trend of the instantaneous effective bending modulus (E: mitral < tricuspid < aortic < pulmonary) might be correlated to the layered fibrous ultrastructures of the four types of leaflets, of which the fibers in mitral and tricuspid leaflets were less aligned, and the fibers in aortic and pulmonary leaflets were highly aligned. We also found that, for all types of leaflets, moment-curvature relationships are stiffer in against-curvature (AC) bending than in with-curvature bending (WC), which implies that leaflets tend to flex toward their natural curvature and comply with blood flow. Lastly, we observed that the leaflets were stiffer in circumferential bending compared with radial bending, likely reflecting the physiological motion of the leaflets, i.e., more bending moment and movement were experienced in radial direction than circumferential direction.

Prediction of matrix-to-cell stress transfer in heart valve tissues

Non-linear and anisotropic heart valve leaflet tissue mechanics manifest principally from the stratification, orientation, and inhomogeneity of their collagenous microstructures. Disturbance of the native collagen fiber network has clear consequences for valve and leaflet tissue mechanics and presumably, by virtue of their intimate embedment, on the valvular interstitial cell stress-strain state and concomitant phenotype. In the current study, a set of virtual biaxial stretch experiments were conducted on porcine pulmonary valve leaflet tissue photomicrographs via an image-based finite element approach. Stress distribution evolution during diastolic valve closure was predicted at both the tissue and cellular levels. Orthotropic material properties consistent with distinct stages of diastolic loading were applied. Virtual experiments predicted tissue- and cellular-level stress fields, providing insight into how matrix-to-cell stress transfer may be influenced by the inhomogeneous collagen fiber architecture, tissue anisotropic material properties, and the cellular distribution within the leaflet tissue. To the best of the authors' knowledge, this is the first study reporting on the evolution of stress fields at both the tissue and cellular levels in valvular tissue and thus contributes toward refining our collective understanding of valvular tissue micromechanics while providing a computational tool enabling the further study of valvular cell-matrix interactions.

DIRECTIONAL BIOMECHANICAL PROPERTIES OF PORCINE SKIN TISSUE

Skin is a multilayered composite material and composed principally of the proteins collagen, elastic fibers and fibroblasts. The direction-dependent material properties of skin tissue is important for physiological functions like skin expansion. The current study has developed methods to characterize the directional biomechanical properties of porcine skin tissues as studies have shown that pigs represent a useful animal model due to similarities between porcine and human skin. It is observed that skin tissue has a nonlinear anisotropy biomechanical behavior, where the parameters of material modulus is 378 ± 160 kPa in the preferred-fiber direction and 65.96 ± 40.49 kPa in the cross-fiber direction when stretching above 30% strain equibiaxially. The result from the study provides methods of characterizing biaxial mechanical properties of skin tissue, as the collagen fiber direction appears to be one of the primary determinants of tissue anisotropy.

ECM hydrogel coating mitigates the chronic inflammatory response to polypropylene mesh

Polypropylene has been used as a surgical mesh material for several decades. This non-degradable synthetic polymer provides mechanical strength, a predictable host response, and its use has resulted in reduced recurrence rates for ventral hernia and pelvic organ prolapse. However, polypropylene and similar synthetic materials are associated with a chronic local tissue inflammatory response and dense fibrous tissue deposition. These outcomes have prompted variations in mesh design to minimize the surface area interface and increase integration with host tissue. In contrast, biologic scaffold materials composed of extracellular matrix (ECM) are rapidly degraded in-vivo and are associated with constructive tissue remodeling and minimal fibrosis. The objective of the present study was to assess the effects of an ECM hydrogel coating on the long-term host tissue response to polypropylene mesh in a rodent model of abdominal muscle injury. At 14 days post implantation, the ECM coated polypropylene mesh devices showed a decreased inflammatory response as characterized by the number and distribution of M1 macrophages (CD86+/CD68+) around mesh fibers when compared to the uncoated mesh devices. At 180 days the ECM coated polypropylene showed decreased density of collagen and amount of mature type I collagen deposited between mesh fibers when compared to the uncoated mesh devices. This study confirms and extends previous findings that an ECM coating mitigates the chronic inflammatory response and associated scar tissue deposition characteristic of polypropylene.

Computational modeling of cardiac valve function and intervention

In the past two decades, major advances have been made in the clinical evaluation and treatment of valvular heart disease owing to the advent of noninvasive cardiac imaging modalities. In clinical practice, valvular disease evaluation is typically performed on two-dimensional (2D) images, even though most imaging modalities offer three-dimensional (3D) volumetric, time-resolved data. Such 3D data offer researchers the possibility to reconstruct the 3D geometry of heart valves at a patient-specific level. When these data are integrated with computational models, native heart valve biomechanical function can be investigated, and preoperative planning tools can be developed. In this review, we outline the advances in valve geometry reconstruction, tissue property modeling, and loading and boundary definitions for the purpose of realistic computational structural analysis of cardiac valve function and intervention.

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.

The living aortic valve: From molecules to function

The aortic valve lies in a unique hemodynamic environment, one characterized by a range of stresses (shear stress, bending forces, loading forces and strain) that vary in intensity and direction throughout the cardiac cycle. Yet, despite its changing environment, the aortic valve opens and closes over 100,000 times a day and, in the majority of human beings, will function normally over a lifespan of 70–90 years. Until relatively recently heart valves were considered passive structures that play no active role in the functioning of a valve, or in the maintenance of its integrity and durability. However, through clinical experience and basic research the aortic valve can now be characterized as a living, dynamic organ with the capacity to adapt to its complex mechanical and biomechanical environment through active and passive communication between its constituent parts. The clinical relevance of a living valve substitute in patients requiring aortic valve replacement has been confirmed. This highlights the importance of using tissue engineering to develop heart valve substitutes containing living cells which have the ability to assume the complex functioning of the native valve.

Effect of glutaraldehyde fixation on the frictional response of immature bovine articular cartilage explants

This study examined functional properties and biocompatibility of glutaraldehyde-fixed bovine articular cartilage over several weeks of incubation at body temperature to investigate its potential use as a resurfacing material in joint arthroplasty. In the first experiment, treated cartilage disks were fixed using 0.02, 0.20 and 0.60% glutaraldehyde for 24h then incubated, along with an untreated control group, in saline for up to 28d at 37°C. Both the equilibrium compressive and tensile moduli increased nearly twofold in treated samples compared to day 0 control, and remained at that level from day 1 to 28; the equilibrium friction coefficient against glass rose nearly twofold immediately after fixation (day 1) but returned to control values after day 7. Live explants co-cultured with fixed explants showed no quantitative difference in cell viability over 28d. In general, no significant differences were observed between 0.20 and 0.60% groups, so 0.20% was deemed sufficient for complete fixation. In the second experiment, cartilage-on-cartilage frictional measurements were performed under a migrating contact configuration. In the treated group, one explant was fixed using 0.20% glutaraldehyde while the apposing explant was left untreated; in the control group both explants were left untreated. From day 1 to 28, the treated group exhibited either no significant difference or slightly lower friction coefficient than the untreated group. These results suggest that a properly titrated glutaraldehyde treatment can reproduce the desired functional properties of native articular cartilage and maintain these properties for at least 28d at body temperature.

Interlayer micromechanics of the aortic heart valve leaflet

While the mechanical behaviors of the fibrosa and ventricularis layers of the aortic valve (AV) leaflet are understood, little information exists on their mechanical interactions mediated by the GAG-rich central spongiosa layer. Parametric simulations of the interlayer interactions of the AV leaflets in flexure utilized a tri-layered finite element (FE) model of circumferentially oriented tissue sections to investigate inter-layer sliding hypothesized to occur. Simulation results indicated that the leaflet tissue functions as a tightly bonded structure when the spongiosa effective modulus was at least 25 % that of the fibrosa and ventricularis layers. Novel studies that directly measured transmural strain in flexure of AV leaflet tissue specimens validated these findings. Interestingly, a smooth transmural strain distribution indicated that the layers of the leaflet indeed act as a bonded unit, consistent with our previous observations (Stella and Sacks in J Biomech Eng 129:757-766, 2007) of a large number of transverse collagen fibers interconnecting the fibrosa and ventricularis layers. Additionally, when the tri-layered FE model was refined to match the transmural deformations, a layer-specific bimodular material model (resulting in four total moduli) accurately matched the transmural strain and moment-curvature relations simultaneously. Collectively, these results provide evidence, contrary to previous assumptions, that the valve layers function as a bonded structure in the low-strain flexure deformation mode. Most likely, this results directly from the transverse collagen fibers that bind the layers together to disable physical sliding and maintain layer residual stresses. Further, the spongiosa may function as a general dampening layer while the AV leaflets deforms as a homogenous structure despite its heterogeneous architecture.

Material properties of aged human mitral valve leaflets

This study aimed to characterize the mechanical properties of aged human anterior mitral leaflets (AML) and posterior mitral leaflets (PML). The AML and PML samples from explanted human hearts (n = 21, mean age of 82.62 ± 8.77-years-old) were subjected to planar biaxial mechanical tests. The material stiffness, extensibility, and degree of anisotropy of the leaflet samples were quantified. The microstructure of the samples was assessed through histology. Both the AML and PML samples exhibited a nonlinear and anisotropic behavior with the circumferential direction being stiffer than the radial direction. The AML samples were significantly stiffer than the PML samples in both directions, suggesting that they should be modeled with separate sets of material properties in computational studies. Histological analysis indicated the changes in the tissue elastic constituents, including the fragmented and disorganized elastin network, the presence of fibrosis and proteoglycan/glycosaminoglycan infiltration and calcification, suggesting possible valvular degenerative characteristics in the aged human leaflet samples. Overall, stiffness increased and areal strain decreased with calcification severity. In addition, leaflet tissues from hypertensive individuals also exhibited a higher stiffness and low areal strain than normotensive individuals. There are significant differences in the mechanical properties of the two human mitral valve leaflets from this advanced age group. The morphologic changes in the tissue composition and structure also infer the structural and functional difference between aged human valves and those of animals.

Real-time strain mapping via biaxial stretching in heart valve tissues

Previous studies show that the collagen fiber architecture is key to the heart valves tissue mechanical property. We report a real-time strain mapping approach that provides displacement vectors and principal strain directions during the mechanical characterization of heart valve tissues. The strain maps reported in the current study allows an individual to quickly identify the approximate strain imposed on a location of the sample. The result shows that when samples are biaxially stretched under 18% strain, less anisotropy is observed in both aortic and pulmonary valve leaflet samples. Moreover, when samples are stretched from 28% to 35%, pulmonary valves leaflet samples exhibits a stronger anisotropic effect than aortic valve. Therefore, a higher degree of straightening is required for collagen fibers to be fully aligned. This work provides an easy approach to quantify mechanical properties with the corresponding strain maps of heart valve tissues and potentially facilitates the developments of tissue engineering heart valves.

Biaxial tension of fibrous tissue: using finite element methods to address experimental challenges arising from boundary conditions and anisotropy

Planar biaxial tension remains a critical loading modality for fibrous soft tissue and is widely used to characterize tissue mechanical response, evaluate treatments, develop constitutive formulas, and obtain material properties for use in finite element studies. Although the application of tension on all edges of the test specimen represents the in situ environment, there remains a need to address the interpretation of experimental results. Unlike uniaxial tension, in biaxial tension the applied forces at the loading clamps do not transmit fully to the region of interest (ROI), which may lead to improper material characterization if not accounted for. In this study, we reviewed the tensile biaxial literature over the last ten years, noting experimental and analysis challenges. In response to these challenges, we used finite element simulations to quantify load transmission from the clamps to the ROI in biaxial tension and to formulate a correction factor that can be used to determine ROI stresses. Additionally, the impact of sample geometry, material anisotropy, and tissue orientation on the correction factor were determined. Large stress concentrations were evident in both square and cruciform geometries and for all levels of anisotropy. In general, stress concentrations were greater for the square geometry than the cruciform geometry. For both square and cruciform geometries, materials with fibers aligned parallel to the loading axes reduced stress concentrations compared to the isotropic tissue, resulting in more of the applied load being transferred to the ROI. In contrast, fiber-reinforced specimens oriented such that the fibers aligned at an angle to the loading axes produced very large stress concentrations across the clamps and shielding in the ROI. A correction factor technique was introduced that can be used to calculate the stresses in the ROI from the measured experimental loads at the clamps. Application of a correction factor to experimental biaxial results may lead to more accurate representation of the mechanical response of fibrous soft tissue.

Age-related analysis of structural, biochemical and mechanical properties of the porcine mitral heart valve leaflets

The aim of this study was to assess structural and biochemical differences in the extracellular matrix of the fetal and adult porcine mitral heart valves in relation to their mechanical characteristics. Using tensile tests it was demonstrated that the material properties of porcine mitral heart valves progressively change with age. The collagen content of the adult heart valve, as estimated by hydroxyproline assay, increases three times as compared with fetal heart valves. Transmission electron microscopy demonstrated that the diameter of collagen fibrils increased in adult heart valves compared with fetal heart valves. The level of collagen cross-linking is lower in the fetal heart valve than the adult heart valve. The reported age differences in the material properties of fetal and adult porcine heart valves were associated with increases in collagen content, the diameter of collagen fibrils and the level of collagen cross-linking. These data lay a foundation for systematic elucidation of the structural determinants of material properties of heart valves during embryonic and postnatal valvulogenesis. They are also essential to define the desirable level of tissue maturation in heart valve tissue engineering.

Polypropylene surgical mesh coated with extracellular matrix mitigates the host foreign body response

Surgical mesh devices composed of synthetic materials are commonly used for ventral hernia repair. These materials provide robust mechanical strength and are quickly incorporated into host tissue; factors that contribute to reduced hernia recurrence rates. However, such mesh devices cause a foreign body response with the associated complications of fibrosis and patient discomfort. In contrast, surgical mesh devices composed of naturally occurring extracellular matrix (ECM) are associated with constructive tissue remodeling, but lack the mechanical strength of synthetic materials. A method for applying a porcine dermal ECM hydrogel coating to a polypropylene mesh is described herein with the associated effects upon the host tissue response and biaxial mechanical behavior. Uncoated and ECM coated heavy-weight BARD™ Mesh were compared to the light-weight ULTRAPRO™ and BARD™ Soft Mesh devices in a rat partial thickness abdominal defect overlay model. The ECM coated mesh attenuated the pro-inflammatory response compared to all other devices, with a reduced cell accumulation and fewer foreign body giant cells. The ECM coating degraded by 35 days, and was replaced with loose connective tissue compared to the dense collagenous tissue associated with the uncoated polypropylene mesh device. Biaxial mechanical characterization showed that all of the mesh devices were of similar isotropic stiffness. Upon explanation, the light-weight mesh devices were more compliant than the coated or uncoated heavy-weight devices. This study shows that an ECM coating alters the default host response to a polypropylene mesh, but not the mechanical properties in an acute in vivo abdominal repair model.