Patient-specific modeling of the heart: estimation of ventricular fiber orientations

Computational modeling of cardiac electrophysiology and arrhythmogenesis: toward clinical translation

The complexity of cardiac electrophysiology, involving dynamic changes in numerous components across multiple spatial (from ion channel to organ) and temporal (from milliseconds to days) scales, makes an intuitive or empirical analysis of cardiac arrhythmogenesis challenging. Multiscale mechanistic computational models of cardiac electrophysiology provide precise control over individual parameters, and their reproducibility enables a thorough assessment of arrhythmia mechanisms. This review provides a comprehensive analysis of models of cardiac electrophysiology and arrhythmias, from the single cell to the organ level, and how they can be leveraged to better understand rhythm disorders in cardiac disease and to improve heart patient care. Key issues related to model development based on experimental data are discussed, and major families of human cardiomyocyte models and their applications are highlighted. An overview of organ-level computational modeling of cardiac electrophysiology and its clinical applications in personalized arrhythmia risk assessment and patient-specific therapy of atrial and ventricular arrhythmias is provided. The advancements presented here highlight how patient-specific computational models of the heart reconstructed from patient data have achieved success in predicting risk of sudden cardiac death and guiding optimal treatments of heart rhythm disorders. Finally, an outlook toward potential future advances, including the combination of mechanistic modeling and machine learning/artificial intelligence, is provided. As the field of cardiology is embarking on a journey toward precision medicine, personalized modeling of the heart is expected to become a key technology to guide pharmaceutical therapy, deployment of devices, and surgical interventions.

Spatially Conserved Spiral Wave Activity During Human Atrial Fibrillation

BACKGROUND

Atrial fibrillation is the most common cardiac arrhythmia in the world and increases the risk for stroke and morbidity. During atrial fibrillation, the electric activation fronts are no longer coherently propagating through the tissue and, instead, show rotational activity, consistent with spiral wave activation, focal activity, collision, or partial versions of these spatial patterns. An unexplained phenomenon is that although simulations of cardiac models abundantly demonstrate spiral waves, clinical recordings often show only intermittent spiral wave activity.

METHODS

In silico data were generated using simulations in which spiral waves were continuously created and annihilated and in simulations in which a spiral wave was intermittently trapped at a heterogeneity. Clinically, spatio-temporal activation maps were constructed using 60 s recordings from a 64 electrode catheter within the atrium of N=34 patients (n=24 persistent atrial fibrillation). The location of clockwise and counterclockwise rotating spiral waves was quantified and all intervals during which these spiral waves were present were determined. For each interval, the angle of rotation as a function of time was computed and used to determine whether the spiral wave returned in step or changed phase at the start of each interval.

RESULTS

In both simulations, spiral waves did not come back in phase and were out of step." In contrast, spiral waves returned in step in the majority (68%; P=0.05) of patients. Thus, the intermittently observed rotational activity in these patients is due to a temporally and spatially conserved spiral wave and not due to ones that are newly created at the onset of each interval.

CONCLUSIONS

Intermittency of spiral wave activity represents conserved spiral wave activity of long, but interrupted duration or transient spiral activity, in the majority of patients. This finding could have important ramifications for identifying clinically important forms of atrial fibrillation and in guiding treatment.

Cardiac structure discontinuities revealed by ex-vivo microstructural characterization. A focus on the basal inferoseptal left ventricle region

BACKGROUND

While the microstructure of the left ventricle (LV) has been largely described, only a few studies investigated the right ventricular insertion point (RVIP). It was accepted that the aggregate cardiomyocytes organization was much more complex due to the intersection of the ventricular cavities but a precise structural characterization in the human heart was lacking even if clinical phenotypes related to right ventricular wall stress or arrhythmia were observed in this region.

METHODS

MRI-derived anatomical imaging (150 µm3) and diffusion tensor imaging (600 µm3) were performed in large mammalian whole hearts (human: N = 5, sheep: N = 5). Fractional anisotropy, aggregate cardiomyocytes orientations and tractography were compared within both species. Aggregate cardiomyocytes orientation on one ex-vivo sheep whole heart was then computed using structure tensor imaging (STI) from 21 µm isotropic acquisition acquired with micro computed tomography (MicroCT) imaging. Macroscopic and histological examination were performed. Lastly, experimental cardiomyocytes orientation distribution was then compared to the usual rule-based model using electrophysiological (EP) modeling. Electrical activity was modeled with the monodomain formulation.

RESULTS

The RVIP at the level of the inferior ventricular septum presented a unique arrangement of aggregate cardiomyocytes. An abrupt, mid-myocardial change in cardiomyocytes orientation was observed, delimiting a triangle-shaped region, present in both sheep and human hearts. FA's histogram distribution (mean ± std: 0.29 ± 0.06) of the identified region as well as the main dimension (22.2 mm ± 5.6 mm) was found homogeneous across samples and species. Averaged volume is 0.34 cm3 ± 0.15 cm3. Both local activation time (LAT) and morphology of pseudo-ECGs were strongly impacted with delayed LAT and change in peak-to-peak amplitude in the simulated wedge model.

CONCLUSION

The study was the first to describe the 3D cardiomyocytes architecture of the basal inferoseptal left ventricle region in human hearts and identify the presence of a well-organized aggregate cardiomyocytes arrangement and cardiac structural discontinuities. The results might offer a better appreciation of clinical phenotypes like RVIP-late gadolinium enhancement or uncommon idiopathic ventricular arrhythmias (VA) originating from this region.

Design and Evaluation of Personalized Percutaneous Coronary Intervention Surgery Simulation System

In recent years, medical simulators have been widely applied to a broad range of surgery training tasks. However, most of the existing surgery simulators can only provide limited immersive environments with a few pre-processed organ models, while ignoring the instant modeling of various personalized clinical cases, which brings substantive differences between training experiences and real surgery situations. To this end, we present a virtual reality (VR) based surgery simulation system for personalized percutaneous coronary intervention (PCI). The simulation system can directly take patient-specific clinical data as input and generate virtual 3D intervention scenarios. Specially, we introduce a fiber-based patient-specific cardiac dynamic model to simulate the nonlinear deformation among the multiple layers of the cardiac structure, which can well respect and correlate the atriums, ventricles and vessels, and thus gives rise to more effective visualization and interaction. Meanwhile, we design a tracking and haptic feedback hardware, which can enable users to manipulate physical intervention instruments and interact with virtual scenarios. We conduct quantitative analysis on deformation precision and modeling efficiency, and evaluate the simulation system based on the user studies from 16 cardiologists and 20 intervention trainees, comparing it to traditional desktop intervention simulators. The results confirm that our simulation system can provide a better user experience, and is a suitable platform for PCI surgery training and rehearsal.

Evaluation of the Differences of Myocardial Fibers between Acute and Chronic Myocardial Infarction: Application of Diffusion Tensor Magnetic Resonance Imaging in a Rhesus Monkey Model

Objective

To understand microstructural changes after myocardial infarction (MI), we evaluated myocardial fibers of rhesus monkeys during acute or chronic MI, and identified the differences of myocardial fibers between acute and chronic MI.

Materials and Methods

Six fixed hearts of rhesus monkeys with left anterior descending coronary artery ligation for 1 hour or 84 days were scanned by diffusion tensor magnetic resonance imaging (MRI) to measure apparent diffusion coefficient (ADC), fractional anisotropy (FA) and helix angle (HA).

Results

Comparing with acute MI monkeys (FA: 0.59 ± 0.02; ADC: 5.0 ± 0.6 × 10^-4 mm²/s; HA: 94.5 ± 4.4°), chronic MI monkeys showed remarkably decreased FA value (0.26 ± 0.03), increased ADC value (7.8 ± 0.8 × 10^-4mm²/s), decreased HA transmural range (49.5 ± 4.6°) and serious defects on endocardium in infarcted regions. The HA in infarcted regions shifted to more components of negative left-handed helix in chronic MI monkeys (-38.3 ± 5.0°–11.2 ± 4.3°) than in acute MI monkeys (-41.4 ± 5.1°–53.1 ± 3.7°), but the HA in remote regions shifted to more components of positive right-handed helix in chronic MI monkeys (-43.8 ± 2.7°–66.5 ± 4.9°) than in acute MI monkeys (-59.5 ± 3.4°–64.9 ± 4.3°).

Conclusion

Diffusion tensor MRI method helps to quantify differences of mechanical microstructure and water diffusion of myocardial fibers between acute and chronic MI monkey's models.

Computer-aided design of the human aortic root

The development of computer-based 3D models of the aortic root is one of the most important problems in constructing the prostheses for transcatheter aortic valve implantation. In the current study, we analyzed data from 117 patients with and without aortic valve disease and computed tomography data from 20 patients without aortic valvular diseases in order to estimate the average values of the diameter of the aortic annulus and other aortic root parameters. Based on these data, we developed a 3D model of human aortic root with unique geometry. Furthermore, in this study we show that by applying different material properties to the aortic annulus zone in our model, we can significantly improve the quality of the results of finite element analysis. To summarize, here we present four 3D models of human aortic root with unique geometry based on computational analysis of ECHO and CT data. We suggest that our models can be utilized for the development of better prostheses for transcatheter aortic valve implantation.