FiberSim: A flexible open-source model of myofilament-level contraction.

FiberSim is a flexible open-source model of myofilament-level contraction. The code uses a spatially explicit technique, meaning that it tracks the position and status of each contractile molecule within the lattice framework. This allows the model to simulate some of the mechanical effects modulated by myosin-binding protein C, as well as the dose dependence of myotropes and the effects of varying isoform expression levels. This paper provides a short introduction to FiberSim and presents simulations of tension-pCa curves with and without regulation of thick and thin filament activation by myosin-binding protein C. A myotrope dose-dependent response as well as slack/re-stretch maneuvers to assess rates of tension recovery are also presented. The software was designed to be flexible (the user can define their own model and/or protocol) and computationally efficient (simulations can be performed on a regular laptop). We hope that other investigators will use FiberSim to explore myofilament level mechanisms and to accelerate research focusing on the contractile properties of sarcomeres.

A Matched-Filter-Based Algorithm for Subcellular Classification of T-System in Cardiac Tissues.

In mammalian ventricular cardiomyocytes, invaginations of the surface membrane form the transverse tubular system (T-system), which consists of transverse tubules (TTs) that align with sarcomeres and Z-lines as well as longitudinal tubules (LTs) that are present between Z-lines in some species. In many cardiac disease etiologies, the T-system is perturbed, which is believed to promote spatially heterogeneous, dyssynchronous Ca2+ release and inefficient contraction. In general, T-system characterization approaches have been directed primarily at isolated cells and do not detect subcellular T-system heterogeneity. Here, we present MatchedMyo, a matched-filter-based algorithm for subcellular T-system characterization in isolated cardiomyocytes and millimeter-scale myocardial sections. The algorithm utilizes "filters" representative of TTs, LTs, and T-system absence. Application of the algorithm to cardiomyocytes isolated from rat disease models of myocardial infarction (MI), dilated cardiomyopathy induced via aortic banding, and sham surgery confirmed and quantified heterogeneous T-system structure and remodeling. Cardiomyocytes from post-MI hearts exhibited increasing T-system disarray as proximity to the infarct increased. We found significant (p < 0.05, Welch's t-test) increases in LT density within cardiomyocytes proximal to the infarct (12 ± 3%, data reported as mean ± SD, n = 3) versus sham (4 ± 2%, n = 5), but not distal to the infarct (7 ± 1%, n = 3). The algorithm also detected decreases in TTs within 5° of the myocyte minor axis for isolated aortic banding (36 ± 9%, n = 3) and MI cardiomyocytes located intermediate (37 ± 4%, n = 3) and proximal (34 ± 4%, n = 3) to the infarct versus sham (57 ± 12%, n = 5). Application of bootstrapping to rabbit MI tissue revealed distal sections comprised 18.9 ± 1.0% TTs, whereas proximal sections comprised 10.1 ± 0.8% TTs (p < 0.05), a 46.6% decrease. The matched-filter approach therefore provides a robust and scalable technique for T-system characterization from isolated cells through millimeter-scale myocardial sections.

3D dSTORM imaging reveals novel detail of ryanodine receptor localization in rat cardiac myocytes.

KEY POINTS: Using 3D direct stochastic optical reconstruction microscopy (dSTORM), we developed novel approaches to quantitatively describe the nanoscale, 3D organization of ryanodine receptors (RyRs) in cardiomyocytes. Complex arrangements of RyR clusters were observed in 3D space, both at the cell surface and within the cell interior, with allocation to dyadic and non-dyadic pools. 3D imaging importantly allowed discernment of clusters overlapping in the z-axis, for which detection was obscured by conventional 2D imaging techniques. Thus, RyR clusters were found to be significantly smaller than previous 2D estimates. Ca2+ release units (CRUs), i.e. functional groupings of neighbouring RyR clusters, were similarly observed to be smaller than earlier reports. Internal CRUs contained more RyRs in more clusters than CRUs on the cell surface, and yielded longer duration Ca2+ sparks. ABSTRACT: Cardiomyocyte contraction is dependent on Ca2+ release from ryanodine receptors (RyRs). However, the precise localization of RyRs remains unknown, due to shortcomings of imaging techniques which are diffraction limited or restricted to 2D. We aimed to determine the 3D nanoscale organization of RyRs in rat cardiomyocytes by employing direct stochastic optical reconstruction microscopy (dSTORM) with phase ramp technology. Initial observations at the cell surface showed an undulating organization of RyR clusters, resulting in their frequent overlap in the z-axis and obscured detection by 2D techniques. Non-overlapping clusters were imaged to create a calibration curve for estimating RyR number based on recorded fluorescence blinks. Employing this method at the cell surface and interior revealed smaller RyR clusters than 2D estimates, as erroneous merging of axially aligned RyRs was circumvented. Functional groupings of RyR clusters (Ca2+ release units, CRUs), contained an average of 18 and 23 RyRs at the surface and interior, respectively, although half of all CRUs contained only a single 'rogue' RyR. Internal CRUs were more tightly packed along z-lines than surface CRUs, contained larger and more numerous RyR clusters, and constituted ∼75% of the roughly 1 million RyRs present in an average cardiomyocyte. This complex internal 3D geometry was underscored by correlative imaging of RyRs and t-tubules, which enabled quantification of dyadic and non-dyadic RyR populations. Mirroring differences in CRU size and complexity, Ca2+ sparks originating from internal CRUs were of longer duration than those at the surface. These data provide novel, nanoscale insight into RyR organization and function across cardiomyocytes.