FM-Track: A fiducial marker tracking software for studying cell mechanics in a three-dimensional environment

Nonlinear Elastic Bottlebrush Polymer Hydrogels Modulate Actomyosin Mediated Protrusion Formation in Mesenchymal Stromal Cells

The nonlinear elasticity of many tissue-specific extracellular matrices is difficult to recapitulate without the use of fibrous architectures, which couple strain-stiffening with stress relaxation. Herein, bottlebrush polymers are synthesized and crosslinked to form poly(ethylene glycol)-based hydrogels and used to study how strain-stiffening behavior affects human mesenchymal stromal cells (hMSCs). By tailoring the bottlebrush polymer length, the critical stress associated with the onset of network stiffening is systematically varied, and a unique protrusion-rich hMSC morphology emerges only at critical stresses within a biologically accessible stress regime. Local cell-matrix interactions are quantified using 3D traction force microscopy and small molecule inhibitors are used to identify cellular machinery that plays a critical role in hMSC mechanosensing of the engineered, strain-stiffening microenvironment. Collectively, this study demonstrates how covalently crosslinked bottlebrush polymer hydrogels can recapitulate strain-stiffening biomechanical cues at biologically relevant stresses and be used to probe how nonlinear elastic matrix properties regulate cellular processes.

An Inverse Modeling Approach to Estimate Three-Dimensional Aortic Valve Interstitial Cell Stress Fiber Force Levels

Within the aortic valve (AV) leaflet exists a population of interstitial cells (AVICs) that maintain the constituent tissues by extracellular matrix (ECM) secretion, degradation, and remodeling. AVICs can transition from a quiescent, fibroblast-like phenotype to an activated, myofibroblast phenotype in response to growth or disease. AVIC dysfunction has been implicated in AV disease processes, yet our understanding of AVIC function remains quite limited. A major characteristic of the AVIC phenotype is its contractile state, driven by contractile forces generated by the underlying stress fibers (SF). However, direct assessment of the AVIC SF contractile state and structure within physiologically mimicking three-dimensional environments remains technically challenging, as the size of single SFs are below the resolution of light microscopy. Therefore, in the present study, we developed a three-dimensional (3D) computational approach of AVICs embedded in 3D hydrogels to estimate their SF local orientations and contractile forces. One challenge with this approach is that AVICs will remodel the hydrogel, so that the gel moduli will vary spatially. We thus utilized our previous approach (Khang et al. 2023, "Estimation of Aortic Valve Interstitial Cell-Induced 3D Remodeling of Poly (Ethylene Glycol) Hydrogel Environments Using an Inverse Finite Element Approach," Acta Biomater., 160, pp. 123-133) to define local hydrogel mechanical properties. The AVIC SF model incorporated known cytosol and nucleus mechanical behaviors, with the cell membrane assumed to be perfectly bonded to the surrounding hydrogel. The AVIC SFs were first modeled as locally unidirectional hyperelastic fibers with a contractile force component. An adjoint-based inverse modeling approach was developed to estimate local SF orientation and contractile force. Substantial heterogeneity in SF force and orientations were observed, with the greatest levels of SF alignment and contractile forces occurring in AVIC protrusions. The addition of a dispersed SF orientation to the modeling approach did not substantially alter these findings. To the best of our knowledge, we report the first fully 3D computational contractile cell models which can predict locally varying stress fiber orientation and contractile force levels.

Functional differences in human aortic valve interstitial cells from patients with varying calcific aortic valve disease

Calcific aortic valve disease (CAVD) is characterized by progressive stiffening of aortic valve (AV) tissues, inducing stenosis and insufficiency. Bicuspid aortic valve (BAV) is a common congenital defect in which the AV has two leaflets rather than three, with BAV patients developing CAVD decades years earlier than in the general population. Current treatment for CAVD remains surgical replacement with its continued durability problems, as there are no pharmaceutical therapies or other alternative treatments available. Before such therapeutic approaches can be developed, a deeper understanding of CAVD disease mechanisms is clearly required. It is known that AV interstitial cells (AVICs) maintain the AV extracellular matrix and are typically quiescent in the normal state, transitioning into an activated, myofibroblast-like state during periods of growth or disease. One proposed mechanism of CAVD is the subsequent transition of AVICs into an osteoblast-like phenotype. A sensitive indicator of AVIC phenotypic state is enhanced basal contractility (tonus), so that AVICs from diseased AV will exhibit a higher basal tonus level. The goals of the present study were thus to assess the hypothesis that different human CAVD states lead to different biophysical AVIC states. To accomplish this, we characterized AVIC basal tonus behaviors from diseased human AV tissues embedded in 3D hydrogels. Established methods were utilized to track AVIC-induced gel displacements and shape changes after the application of Cytochalasin D (an actin polymerization inhibitor) to depolymerize the AVIC stress fibers. Results indicated that human diseased AVICs from the non-calcified region of TAVs were significantly more activated than AVICs from the corresponding calcified region. In addition, AVICs from the raphe region of BAVs were more activated than from the non-raphe region. Interestingly, we observed significantly greater basal tonus levels in females compared to males. Furthermore, the overall AVIC shape changes after Cytochalasin suggested that AVICs from TAVs and BAVs develop different stress fiber architectures. These findings are the first evidence of sex-specific differences in basal tonus state in human AVICs in varying disease states. Future studies are underway to quantify stress fiber mechanical behaviors to further elucidate CAVD disease mechanisms.

Three-dimensional analysis of hydrogel-imbedded aortic valve interstitial cell shape and its relation to contractile behavior

Cell-shape is a conglomerate of mechanical, chemical, and biological mechanisms that reflects the cell biophysical state. In a specific application, we consider aortic valve interstitial cells (AVICs), which maintain the structure and function of aortic heart valve leaflets. Actomyosin stress fibers help determine AVIC shape and facilitate processes such as adhesion, contraction, and mechanosensing. However, detailed 3D assessment of stress fiber architecture and function is currently impractical. Herein, we assessed AVIC shape and contractile behaviors using hydrogel-based 3D traction force microscopy to intuit the orientation and behavior of AVIC stress fibers. We utilized spherical harmonics (SPHARM) to quantify AVIC geometries through three days of incubation, which demonstrated a shift from a spherical shape to forming substantial protrusions. Furthermore, we assessed changes in post-three day AVIC shape and contractile function within two testing regimes: (1) normal contractile level to relaxation (cytochalasin D), and (2) normal contractile level to hyper-contraction (endothelin-1). In both scenarios, AVICs underwent isovolumic shape changes and produced complex displacement fields within the hydrogel. AVICs were more elongated when relaxed and more spherical in hyper-contraction. Locally, AVIC protrusions contracted along their long axis and expanded in their circumferential direction, indicating predominately axially aligned stress fibers. Furthermore, the magnitude of protrusion displacements was correlated with protrusion length and approached a consistent displacement plateau at a similar critical length across all AVICs. This implied that stress fiber behavior is conserved, despite great variations in AVIC shapes. We anticipate our findings will bolster future investigations into AVIC stress fiber architecture and function. STATEMENT OF SIGNIFICANCE: Within the aortic valve there exists a population of aortic valve interstitial cells, which orchestrate the turnover, secretion, and remodeling of its extracellular matrix, maintaining tissue integrity and ultimately sustaining the proper mechanical function. Alterations in these processes are thought to underlie diseases of the aortic valve, which affect hundreds of thousands domestically and world-wide. Yet, to date, there are no non-surgical treatments for aortic heart valve disease, in part due to our limited understanding of the underlying disease processes. In the present study, we built upon our previous study to include a full 3D analysis of aortic valve interstitial cell shapes at differing contractile levels. The resulting detailed shape and deformation analysis provided insight into the underlying stress-fiber structures and mechanical behaviors.

Estimation of aortic valve interstitial cell-induced 3D remodeling of poly(ethylene glycol) hydrogel environments using an inverse finite element approach

Aortic valve interstitial cells (AVICs) reside within the leaflet tissues of the aortic valve and maintain and remodel its extracellular matrix components. Part of this process is a result of AVIC contractility brought about by underlying stress fibers whose behaviors can change in various disease states. Currently, it is challenging to directly investigate AVIC contractile behaviors within dense leaflet tissues. As a result, optically clear poly (ethylene glycol) hydrogel matrices have been used to study AVIC contractility via 3D traction force microscopy (3DTFM). However, the local stiffness of the hydrogel is difficult to measure directly and is further confounded by the remodeling activity of the AVIC. Ambiguity in hydrogel mechanics can lead to large errors in computed cellular tractions. Herein, we developed an inverse computational approach to estimate AVIC-induced remodeling of the hydrogel material. The model was validated with test problems comprised of an experimentally measured AVIC geometry and prescribed modulus fields containing unmodified, stiffened, and degraded regions. The inverse model estimated the ground truth data sets with high accuracy. When applied to AVICs assessed via 3DTFM, the model estimated regions of significant stiffening and degradation in the vicinity of the AVIC. We observed that stiffening was largely localized at AVIC protrusions, likely a result of collagen deposition as confirmed by immunostaining. Degradation was more spatially uniform and present in regions further away from the AVIC, likely a result of enzymatic activity. Looking forward, this approach will allow for more accurate computation of AVIC contractile force levels. STATEMENT OF SIGNIFICANCE: The aortic valve (AV), positioned between the left ventricle and the aorta, prevents retrograde flow into the left ventricle. Within the AV tissues reside a resident population of aortic valve interstitial cells (AVICs) that replenish, restore, and remodel extracellular matrix components. Currently, it is technically challenging to directly investigate AVIC contractile behaviors within the dense leaflet tissues. As a result, optically clear hydrogels have been used to study AVIC contractility through means of 3D traction force microscopy. Herein, we developed a method to estimate AVIC-induced remodeling of PEG hydrogels. This method was able to accurately estimate regions of significant stiffening and degradation induced by the AVIC and allows a deeper understanding of AVIC remodeling activity, which can differ in normal and disease conditions.

Sarc-Graph: Automated segmentation, tracking, and analysis of sarcomeres in hiPSC-derived cardiomyocytes

A better fundamental understanding of human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) has the potential to advance applications ranging from drug discovery to cardiac repair. Automated quantitative analysis of beating hiPSC-CMs is an important and fast developing component of the hiPSC-CM research pipeline. Here we introduce “Sarc-Graph,” a computational framework to segment, track, and analyze sarcomeres in fluorescently tagged hiPSC-CMs. Our framework includes functions to segment z-discs and sarcomeres, track z-discs and sarcomeres in beating cells, and perform automated spatiotemporal analysis and data visualization. In addition to reporting good performance for sarcomere segmentation and tracking with little to no parameter tuning and a short runtime, we introduce two novel analysis approaches. First, we construct spatial graphs where z-discs correspond to nodes and sarcomeres correspond to edges. This makes measuring the network distance between each sarcomere (i.e., the number of connecting sarcomeres separating each sarcomere pair) straightforward. Second, we treat tracked and segmented components as fiducial markers and use them to compute the approximate deformation gradient of the entire tracked population. This represents a new quantitative descriptor of hiPSC-CM function. We showcase and validate our approach with both synthetic and experimental movies of beating hiPSC-CMs. By publishing Sarc-Graph, we aim to make automated quantitative analysis of hiPSC-CM behavior more accessible to the broader research community.

Heart disease is the leading cause of death worldwide. Because of this, many researchers are studying heart cells in the lab and trying to create artificial heart tissue. Recently, there has been a growing focus on human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs). These are cells that are safely sampled from living humans, for example from the blood or skin, that are then transformed into human heart muscle cells. One active research goal is to use these cells to repair the damaged heart. Another active research goal is to test new drugs on these cells before testing them in animals and humans. However, one major challenge is that hiPSC-CMs often have an irregular internal structure that is difficult to analyze. At present, their behavior is far from fully understood. To address this, we have created software to automatically analyze movies of beating hiPSC-CMs. With our software, it is possible to quantify properties such as the amount and direction of beating cell contraction, and the variation in behavior across different parts of each cell. These tools will enable further quantitative analysis of hiPSC-CMs. With these tools, it will be easier to understand, control, and optimize artificial heart tissue created with hiPSC-CMs, and quantify the effects of drugs on hiPSC-CM behavior.

On the Three-Dimensional Correlation Between Myofibroblast Shape and Contraction

Myofibroblasts are responsible for wound healing and tissue repair across all organ systems. In periods of growth and disease, myofibroblasts can undergo a phenotypic transition characterized by an increase in extracellular matrix (ECM) deposition rate, changes in various protein expression (e.g., alpha-smooth muscle actin (αSMA)), and elevated contractility. Cell shape is known to correlate closely with stress-fiber geometry and function and is thus a critical feature of cell biophysical state. However, the relationship between myofibroblast shape and contraction is complex, even as well in regards to steady-state contractile level (basal tonus). At present, the relationship between myofibroblast shape and basal tonus in three-dimensional (3D) environments is poorly understood. Herein, we utilize the aortic valve interstitial cell (AVIC) as a representative myofibroblast to investigate the relationship between basal tonus and overall cell shape. AVICs were embedded within 3D poly(ethylene glycol) (PEG) hydrogels containing degradable peptide crosslinkers, adhesive peptide sequences, and submicron fluorescent microspheres to track the local displacement field. We then developed a methodology to evaluate the correlation between overall AVIC shape and basal tonus induced contraction. We computed a volume averaged stretch tensor ⟨U⟩ for the volume occupied by the AVIC, which had three distinct eigenvalues (λ1,2,3=1.08,0.99, and 0.89), suggesting that AVIC shape is a result of anisotropic contraction. Furthermore, the direction of maximum contraction correlated closely with the longest axis of a bounding ellipsoid enclosing the AVIC. As gel-imbedded AVICs are known to be in a stable state by 3 days of incubation used herein, this finding suggests that the overall quiescent AVIC shape is driven by the underlying stress-fiber directional structure and potentially contraction level.