Cardiac resynchronisation therapy (CRT) is an established treatment for heart failure, however the effective selection of patients and optimisation of therapy remain controversial. While extensive research is ongoing, it remains unclear whether improvements in patient selection or therapy planning offers a greater opportunity for the improvement of clinical outcomes. This computational study investigates the impact of both physiological conditions that guide patient selection and the optimisation of pacing lead placement on CRT outcomes. A multi-scale biophysical model of cardiac electromechanics was developed and personalised to patient data in three patients. These models were separated into components representing cardiac anatomy, pacing lead location, myocardial conductivity and stiffness, afterload, active contraction and conduction block for each individual, and recombined to generate a cohort of 648 virtual patients. The effect of these components on the change in total activation time of the ventricles (ΔTAT) and acute haemodynamic response (AHR) was analysed. The pacing site location was found to have the largest effect on ΔTAT and AHR. Secondary effects on ΔTAT and AHR were found for functional conduction block and cardiac anatomy. The simulation results highlight a need for a greater emphasis on therapy optimisation in order to achieve the best outcomes for patients.

Computational cardiac physiology has great potential to improve the management of cardiovascular diseases. One of the main bottlenecks in this field is the customization of the computational model to the anatomical and physiological status of the patient. We present a fully automatic service for the geometrical personalization of cardiac ventricular meshes with high-order interpolation from segmented images. The method is versatile (able to work with different species and disease conditions) and robust (fully automatic results fulfilling accuracy and quality requirements in 87% of 255 cases). Results also illustrate the capability to minimize the impact of segmentation errors, to overcome the sparse resolution of dynamic studies and to remove the sometimes unnecessary anatomical detail of papillary and trabecular structures. The smooth meshes produced can be used to simulate cardiac function, and in particular mechanics, or can be used as diagnostic descriptors of anatomical shape by cardiologists. This fully automatic service is deployed in a cloud infrastructure, and has been made available and accessible to the scientific community.

The quality of a computational mesh is an important characteristic for stable and accurate simulations. Quality depends on the regularity of the initial mesh, and in mechanical simulations it evolves in time, with deformations causing changes in volume and distortion of mesh elements. Mesh quality metrics are therefore relevant for both mesh personalization and the monitoring of the simulation process. This work evaluates the significance, in meshes with high order interpolation, of four quality metrics described in the literature, applying them to analyse the stability of the simulation of the heart beat. It also investigates how image registration and mesh warping parameters affect the quality and stability of meshes. Jacobian-based metrics outperformed or matched the results of coarse geometrical metrics of aspect ratio or orthogonality, although they are more expensive computationally. The stability of simulations of a complete heart cycle was best predicted with a specificity of 61%, sensitivity of 85%, and only nominal differences were found changing the intra-element and per-element combination of quality values. A compromise between fitting accuracy and mesh stability and quality was found. Generic geometrical quality metrics have a limited success predicting stability, and an analysis of the simulation problem may be required for an optimal definition of quality.

This computational modelling study supports pacing in laterobasal locations, optimizing pacing site by minimizing paced QRS duration and pacing in regions activated late at sinus rhythm. Results demonstrate information content is redundant using multiple preclinical metrics. Of significance, the correlation of AHR with ΔTAT indicates that minimization of QRSd is a promising metric for optimization of lead placement.