Recent advances in magnetic resonance (MR) imaging technology have unveiled a wealth of information regarding cardiac histoanatomical complexity. However, methods to faithfully translate this level of fine-scale structural detail into computational whole ventricular models are still in their infancy, and, thus, the relevance of this additional complexity for simulations of cardiac function has yet to be elucidated. Here, we describe the development of a highly detailed finite-element computational model (resolution: ∼125 μm) of rabbit ventricles constructed from high-resolution MR data (raw data resolution: 43 × 43 × 36 μm), including the processes of segmentation (using a combination of level-set approaches), identification of relevant anatomical features, mesh generation, and myocyte orientation representation (using a rule-based approach). Full access is provided to the completed model and MR data. Simulation results were compared with those from a simplified model built from the same images but excluding finer anatomical features (vessels/endocardial structures). Initial simulations showed that the presence of trabeculations can provide shortcut paths for excitation, causing regional differences in activation after pacing between models. Endocardial structures gave rise to small-scale virtual electrodes upon the application of external field stimulation, which appeared to protect parts of the endocardium in the complex model from strong polarizations, whereas intramural virtual electrodes caused by blood vessels and extracellular cleft spaces appeared to reduce polarization of the epicardium. Postshock, these differences resulted in the genesis of new excitation wavefronts that were not observed in more simplified models. Furthermore, global differences in the stimulus recovery rates of apex/base regions were observed, causing differences in the ensuing arrhythmogenic episodes. In conclusion, structurally simplified models are well suited for a large range of cardiac modeling applications. However, important differences are seen when behavior at microscales is relevant, particularly when examining the effects of external electrical stimulation on tissue electrophysiology and arrhythmia induction. This highlights the utility of histoanatomically detailed models for investigations of cardiac function, in particular for future patient-specific modeling.

We investigate acute effects of axial stretch, applied by carbon fibers (CFs), on diastolic Ca2± spark rate in rat isolated cardiomyocytes. CFs were attached either to both cell ends (to maximize the stretched region), or to the center and one end of the cell (to compare responses in stretched and nonstretched half-cells). Sarcomere length was increased by 8.01 ± 0.94% in the stretched cell fraction, and time series of XY confocal images were recorded to monitor diastolic Ca2± spark frequency and dynamics. Whole-cell stretch causes an acute increase of Ca2± spark rate (to 130.7 ± 6.4%) within 5 seconds, followed by a return to near background levels (to 104.4±5.1%) within 1 minute of sustained distension. Spark rate increased only in the stretched cell region, without significant differences in spark amplitude, time to peak, and decay time constants of sparks in stretched and nonstretched areas. Block of stretch-activated ion channels (2 gmol/L GsMTx-4), perfusion with Na±/Ca2±-free solution, and block of nitric oxide synthesis (1 mmol/L L-NAME) all had no effect on the stretch-induced acute increase in Ca2± spark rate. Conversely, interference with cytoskeletal integrity (2 hours of 10 gmol/L colchicine) abolished the response. Subsequent electron microscopic tomography confirmed the close approximation of microtubules with the T-tubular–sarcoplasmic reticulum complex (to within · 10−8m). In conclusion, axial stretch of rat cardiomyocytes acutely and transiently increases sarcoplasmic reticulum Ca2± spark rate via a mechanism that is independent of sarcolemmal stretch-activated ion channels, nitric oxide synthesis, or availability of extracellular calcium but that requires cytoskeletal integrity. The potential of microtubule-mediated modulation of ryanodine receptor function warrants further investigation.

This paper presents methods to build histo-anatomically detailed individualized cardiac models. The models are based on high-resolution three-dimensional anatomical and/or diffusion tensor magnetic resonance images, combined with serial histological sectioning data, and are used to investigate individualized cardiac function. The current state of the art is reviewed, and its limitations are discussed. We assess the challenges associated with the generation of histo-anatomically representative individualized in silico models of the heart. The entire processing pipeline including image acquisition, image processing, mesh generation, model set-up and execution of computer simulations, and the underlying methods are described. The multifaceted challenges associated with these goals are highlighted, suitable solutions are proposed, and an important application of developed highresolution structure-function models in elucidating the effect of individual structural heterogeneity upon wavefront dynamics is demonstrated.