# Desynchronization Strain Patterns and Contractility in Left Bundle Branch Block through Computer Model Simulation

Kimi P Owashi, Marion Taconné, Nicolas Courtial, Antoine Simon, Mireille

Garreau, Alfredo Hernandez, Erwan Donal, Virginie Le Rolle, Elena Galli

<figure><img src="image_1.png" /></figure>

► To cite this version:

Kimi P Owashi, Marion Taconné, Nicolas Courtial, Antoine Simon, Mireille Garreau, et al. Desynchronization Strain Patterns and Contractility in Left Bundle Branch Block through Computer Model Simulation. Journal of Cardiovascular Development and Disease, 2022, 9 (2), pp.53. 10.3390/jcdd9020053. [https://hal.science/hal-03593832](https://hal.science/hal-03593832)

HAL Id: hal-03593832

https://hal.science/hal-03593832

Submitted on 2 Mar 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

Article

# Desynchronization Strain Patterns and Contractility in Left Bundle Branch Block through Computer Model Simulation

Kimi Owashi †, Marion Taconné † ID, Nicolas Courtial, Antoine Simon, Mireille Garreau, Alfredo Hernandez, Erwan Donal ID, Virginie Le Rolle ‡ ID and Elena Galli *,‡

INSERM, CHU Rennes, University of Rennes, LTSI-UMR 1099, F-35000 Rennes, France;
kimiowashi@gmail.com (K.O.); marion.taconne@univ-rennes1.fr (M.T.); nicolas.courtial@univ-rennes1.fr (N.C.);
antoine.simon@univ-rennes1.fr (A.S.); mireille.garreau@univ-rennes1.fr (M.G.);
alfredo.hernardes@univ-rennes1.fr (A.H.); erwan.donal@chu-rennes.fr (E.D.);
virginie.lerolle@univ-rennes1.fr (V.L.R.)

* Correspondence: gallelena@gmail.com

† These authors contributed equally to this work.

‡ These authors contributed equally to this work.

**Abstract:** Left bundle branch block (LBBB) is associated with specific septal-to-lateral wall activation patterns which are strongly influenced by the intrinsic left ventricular (LV) contractility and myocardial scar localization. The objective of this study was to propose a computational-model-based interpretation of the different patterns of LV contraction observed in the case of LBBB and preserved contractility or myocardial scarring. Two-dimensional transthoracic echocardiography was used to obtain LV volumes and deformation patterns in three patients with LBBB: (1) a patient with non-ischemic dilated cardiomyopathy, (2) a patient with antero-septal myocardial scar, and (3) a patient with lateral myocardial scar. Scar was confirmed by the distribution of late gadolinium enhancement with cardiac magnetic resonance imaging (cMRI). Model parameters were evaluated manually to reproduce patient-derived data such as strain curves obtained from echocardiographic apical views. The model was able to reproduce the specific strain patterns observed in patients. A typical septal flash with pre-ejection shortening, rebound stretch, and delayed lateral wall activation was observed in the case of non-ischemic cardiomyopathy. In the case of lateral scar, the contractility of the lateral wall was significantly impaired and septal flash was absent. In the case of septal scar, septal flash and rebound stretch were also present as previously described in the literature. Interestingly, the model was also able to simulate the specific contractile properties of the myocardium, providing an excellent localization of LV scar in ischemic patients. The model was able to simulate the electromechanical delay and specific contractility patterns observed in patients with LBBB of ischemic and non-ischemic etiology. With further improvement and validation, this technique might be a useful tool for the diagnosis and treatment planning of heart failure patients needing CRT.

**Keywords:** cardiac resynchronization therapy; left bundle branch block; left ventricular deformation; personalized computational model; myocardial scar

<figure><img src="image_4.png" /></figure>

### check for updates

**Citation:** Owashi, K.; Taconné, M.; Courtial, N.; Simon, A.; Garreau, M.; Hernandez, A.; Donal, E.; Le Rolle, V.; Galli, E. Desynchronization Strain Patterns and Contractility in Left Bundle Branch Block through Computer Model Simulation. *J. Cardiovasc. Dev. Dis.* **2022**, *9*, 53.
https://doi.org/10.3390/jcdd9020053

Academic Editor: Enrico Ferrari

Received: 29 December 2021

Accepted: 25 January 2022

Published: 6 February 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_5.png" /></figure>

**Copyright:** © 2022 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

## 1. Introduction

Cardiac resynchronization therapy (CRT) is an established therapy for patients with systolic heart failure and widened QRS (>130 ms) who remain symptomatic despite optimized medical therapy [1]. Despite this, 30–40% of patients receiving CRT according to current recommendations do not experience significant left ventricular (LV) reverse remodeling after CRT [2]. In recent years, several measurements of LV dyssynchrony have been proposed, but none of them have been shown to significantly improve responders’ rate [1]. This might be because only specific mechanical substrates associated with LV electrical delay are amenable to be responsive to CRT [3]. The identification of these substrates is particularly challenging because it should consider the specific electro-mechanical properties of
the myocardium, which are highly patient-specific. In this context, computational modeling might be an efficient tool to integrate cardiac electrical activation, mechanical properties, and hemodynamic conditions in data processing. Indeed, various computational models have been proposed to underline the different physiological mechanisms. Most of them are based on the finite element method [4–6], including a complete representation of the geometry, which requires high computational resources and is difficult to personalize. Alternative approaches have been proposed to overcome this computational cost by simplifying the patient anatomy [7,8]. Although these models successfully represented a patient diversity in LV dyssynchrony and proposed keys to understand CRT response, patient specification with real and non-invasive clinical data must be made to propose a personalized medicine.

Various models of the cardiovascular system have been previously developed by our team [9,10]. We proposed the first patient-specific model of the left ventricle that successfully simulates myocardial strains in dilated cardiomyopathy [11]. In this paper, a model of the cardiovascular system was adapted from [12], integrating four main sub-models: (1) the cardiac electrical system, (2) right and left atria, (3) multi-segment representation of the right and the left ventricles, and (4) systemic and pulmonary circulations. The model considered not only the electrical activation delay of each segment, but also the differences in regional contractility, which are known to largely contribute to strain morphology and global cardiac mechanics.

The aim of this study was to propose a model-based interpretation of the different patterns of LV contraction observed in patients with LBBB in case of myocardial scarring and non-ischemic cardiomyopathy.

## 2. Materials and Methods

### 2.1. *Experimental Data*

#### 2.1.1. Study Population

We prospectively included three CRT candidates with typical LBBB (QRS width > 130 ms) participating in ongoing CRT studies at the Department of Cardiology of the University Hospital of Rennes, France. One patient had non-ischemic dilated cardiomyopathy, one patient had ischemic cardiomyopathy with lateral scar, and one patient had ischemic cardiomyopathy with antero-septal scar. All patients underwent transthoracic echocardiography before CRT implantation. We performed cardiac magnetic resonance imaging (cMRI) with late gadolinium enhancement (LGE) to exclude previous myocardial infarction and to localize myocardial scar when needed.

Table 1 summarizes patients’ clinical characteristics. The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the local Institutional Ethics Committee (French Brittany Regional Ethic Committee validation number: 35RC14-9767).

<table><caption>**Table 1.** Main patient characteristics.</caption><thead><tr><td></td><td>Age (years)</td><td>Etiology</td><td>QRS Width (ms)</td><td>LVEF (%)</td></tr></thead><tbody><tr><td>Patient 1</td><td>77</td><td>Non-ischemic cardiomyopathy</td><td>136</td><td>33</td></tr><tr><td>Patient 2</td><td>80</td><td>Ischemic cardiomyopathy—lateral scar</td><td>157</td><td>24</td></tr><tr><td>Patient 3</td><td>82</td><td>Ischemic cardiomyopathy—antero-septal scar</td><td>152</td><td>30</td></tr></tbody></table>

LVEF, left ventricular ejection fraction.

#### 2.1.2. Echocardiography

All patients underwent standard transthoracic echocardiography using a Vivid E9 or E95 ultrasound system (GE Healthcare, Horten, Norway) equipped with an M5S 3.5-MHz transducer before CRT implantation.

LV volumes were measured according to current recommendations, and LV ejection fraction (LVEF) was calculated using the Simpson biplane method [13].
Two-dimensional grayscale images were acquired in the apical 4-, 3-, and 2-chamber views at a frame rate ≥60 frames/s. The recordings were processed offline using an acoustic-tracking dedicated software (EchoPAC version 202, GE Healthcare, Horten, Norway) to estimate LV global longitudinal strain (GLS). Numerized strain traces (.txt files) were exported for computational analysis.

### 2.1.3. Cardiac Magnetic Resonance Imaging

The cMRI was performed before CRT implantation using a 3-T clinical magnetic resonance system (Ingenia, Philips Medical Systems, Best, The Netherlands) with a 32-channel cardiovascular array coil. LGE images were acquired 10–15 min after intravenous administration of 0.2 mmol/kg of gadolinium (Gadoterate meglumine, Dotarem, Guerbet, Aulnay-sous-Bois, France), using 2D breath-hold inversion-recovery and phase-sensitive inversion-recovery sequences in the short-axis plane (spoiled gradient-echo, slice thickness 8 mm, repetition time 6.1 ms, echo time 2.9 ms, flip angle 25°, inversion time adjusted to null normal myocardium, typical breath-hold 11 s). The presence and localization of myocardial scar were assessed by a trained radiologist on a 16-segment LV model and the regional LGE extent was determined qualitatively and assessed on a per-segment basis.

Concerning two ischemic LBBB patients, both inner and outer contours were manually delineated on LGE-SAX images to have a precise estimation of the fibrotic burden. Then, the fuzzy c-means method proposed by [14] was used to delineate myocardial fibrosis. This is an unsupervised method classifying the voxels in the myocardium as belonging to one of two possible classes: late-gadolinium-enhanced voxel, or non-LGE voxel. This method proved appropriate correlations with biochemical myocardial infarct (scar) quantification as well as LV function parameters [15].

## 2.2. Computational Model

The computational model included four components: (1) cardiac electrical system, (2) right and left atria, (3) multi-segment representation of the right and the left ventricles, and (4) systemic and pulmonary circulations. The model was adapted from previous works of our team [11,12,16–18]. Model parameters were adjusted manually to fit clinical data obtained from patients (Figure 1). To experimentally reproduce the lateral and septal strain traces of representative ischemic and non-ischemic patients with typical LBBB, we varied the value of one model parameter at a time (electro-mechanical driving, contractility, and depolarization delay) until we were able to simulate the pattern of electro-mechanical activation observed in patients.

### 2.2.1. Electrical Activity

The cardiac electrical activity was represented by a set of interconnected cellular automata, adapted from [11,19]. Each automaton was associated with the main electrophysiological activation periods: slow diastolic depolarization (SDD), upstroke depolarization period (UDP), absolute refractory period (ARP), and relative refractory period (RRP). Each automaton cycled between these four stages/phases, sending out a stimulation signal to its neighbors when the UDP ended. The model consisted of 26 automata: the sinoatrial node (NSA), right and left atria (RA and LA), the atrioventricular node (NAV), upper bundle of His (UH), right and left bundle branches (RBB and LBB), 3 segments of right ventricle (RV) and 16 segments of left ventricle (LV). Each automaton created an electrical signal based on a template and tuned by the activation period timings. These electrical signals were gathered to create an electrocardiogram signal and access to each activation delay of the LV segments.
<figure><img src="image_6.png" /><figcaption>**Figure 1.** Exemplification of the exploitation of ECG and echocardiographic data obtained from patients for the construction of an electro-mechanical heart model and characteristics of the model. ((Left) panel): ECG cycle-time and strain-derived curves obtained from each patient were compared to the outputs of the model and used to adjust the model's outputs to increase the patient's specification of the model. ((Middle) panel): Closed-loop model of the cardiovascular system including a multi-segment representation of the left and right ventricles. Numbers (1 to 19) represent the different cardiac segments: (1) basal anterior; (2) basal anteroseptal; (3) basal inferoseptal; (4) basal inferior; (5) basal inferolateral; (6) basal anterolateral; (7) mid anterior; (8) mid anteroseptal; (9) mid inferoseptal; (10) mid inferior; (11) mid inferolateral; (12) mid anterolateral; (13) apical anterior; (14) apical septal; (15) apical inferior; (16) apical lateral; (17) right ventricular basal; (18) right ventricular median; (19) right ventricular apical. The systemic (sys) and pulmonary (pul) circulations include: aorta (ao), vena cava (vc), pulmonary artery (pa), and pulmonary veins (pu). ((Right) panel): State diagram of the cellular automaton that shows the correspondence of the transition parameters with the myocardial action potential dynamic.</figcaption></figure>

### 2.2.2. Right and Left Atria

The right and left atrial pressures were defined as linear functions of instantaneous volumes and elastances representing the elastic properties of the atrial wall. This linear relation was driven by a Gaussian function that cycles between atrial diastole and systole influenced by each atria activation time, established from the stimulation of the atrial automaton [16,20].

### 2.2.3. Multi-Segment Representation of the Right and the Left Ventricles

The multi-segment model was adapted from previous work of our team [11,12]. The left ventricle wall was divided into 16 segments according to the standardized segmentation of the AHA [21]. The right ventricle wall was divided into three layers: base, medium, and apex. The left ventricular basal and medium walls were separated into six components: anterior, anteroseptal, inferoseptal, inferior, inferolateral, and anterolateral walls. The apical layer was divided into four components: anterior, septal, inferior, and lateral.

The cardiac mechanical activity was represented by active and passive components. The active tension is related with the myofiber contractility and muscle kinetic properties, and the passive tension depends on the power conservation law between LV pressure and wall tension. Longitudinal strain was calculated from current and reference fibers' length. Each segment was activated by an electro-mechanical driving function triggered by the corresponding LV and RV automaton. The hydraulic behavior of the blood volume in contact with the wall segment was represented by its inertial and resistive effects. Ventricular flow accounted for the contribution of each segment and intra-ventricular
cavity flow. Septal segment pressures depended directly on the pressure gradient across the septum.

### 2.2.4. Systemic and Pulmonary Circulations

The cardiovascular system model integrated the right and left atria and ventricles, the cardiac valves, the pulmonary arteries, capillaries and veins, and the systemic arteries and veins [16,22]. Concerning each cavity, the pressure was related to volume using a linear relationship defined by an elastance. Then, blood flow between two chambers was calculated from the pressure gradient between the chambers and the corresponding resistance. The cardiac valves, located at the inlet and outlet of the ventricles, were represented as perfect diodes.

Figure 1 shows how ECG cycle time and strain curves obtained at transthoracic echocardiography were used to improve the patient specification of the model.

## 3. Results

### 3.1. *Baseline Model Simulations*

Figure 2 illustrates the simulation results from the proposed computational model at baseline. Systolic LV pressure was equal to 130 mmHg, the aortic pressure varied between 60 and 130 mmHg, and LV volume varied between 40 and 90 mL; simulated aortic and mitral flows were like those observed in healthy subjects.

<figure><img src="image_7.png" /><figcaption>**Figure 2.** Hemodynamic simulations in baseline conditions.</figcaption></figure>

Figure 3 illustrates the strain signals obtained at baseline. Strain traces presented similar morphologies and timing across the 16 LV segments, due to the absence of LV mechanical dyssynchrony and the presence of a preserved, homogeneous LV contractility.

### 3.2. *LBBB Strain Patterns*

Figure 4 illustrates the experimental and model-based strain traces obtained in the septal and lateral wall for patients with non-ischemic cardiomyopathy, as well as lateral and anterior ischemia. In the case of non-ischemic LBBB, the model was able to simulate a typical septo-to-lateral activation pattern, which is characterized by a pre-ejection contraction of the septal wall followed by an immediate re-lengthening (septal rebound stretch) of the wall. In the case of lateral scar, the LV activation pattern simulated by the model was quite different and characterized by a significant reduction in lateral wall strain, loss of the septal
rebound stretch, and a marked increase in septal deformation. Finally, the model simulated LV activation in the case of LBBB and anteroseptal scar. In this case, the septal deformation pattern was still characterized by an early septal shortening followed by rebound stretch. The root-mean-square errors (RMSEs) calculated between estimated and observed strain curves of the septal and lateral walls were equal to (i) 3.99% and 1.86% for the non-ischemic patient, (ii) 3.39% and 1.86% for the lateral ischemia, and (iii) 1.72% and 4.42% for the anteroseptal ischemia.

<figure><img src="image_8.png" /><figcaption>**Figure 3.** Strain simulations in baseline conditions. The presented strain curves correspond to the longitudinal deformation of each of the 16 left ventricular segments and 3 right ventricular segments during one cardiac cycle obtained after simulation of a healthy patient.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 4.** LBBB strain patterns obtained from simulated and clinical data for septal and lateral walls.</figcaption></figure>

The septal wall curve was obtained as the average of the five strain curves of septal segments (apical septum, mid inferoseptum, basal inferoseptum, mid anteroseptum, and basal anteroseptum).

The lateral wall curve was obtained as the average of the five strain curves of the lateral segments (apical lateral, mid inferolateral, basal inferolateral, mid anterolateral, and basal anterolateral).
Figure 5 describes how the model was able to reproduce the LV electrical delay and localized alteration in contractility. In all cases, we observed a septal-to-lateral wall delay typical of LBBB (lower panel). Furthermore, the model was able to reproduce regional modifications in LV contractility which are due to the LBBB, but also to local scarring. In the case of isolated LBBB, we observed increased contractility of the lateral wall compared to the septum. In the case of lateral scar, we observed a significant impairment in lateral contractility. In the case of anteroseptal scar, a higher reduction in contractility was observed in the septal and apical segments. As depicted in Figure 5, reduced contractility in ischemic patients corresponded to the areas of transmural distribution of late gadolinium enhancement observed with cMRI (Figure 5). A higher percentage of transmurality translates into larger fibrotic areas, which are associated with low contractility. Therefore, regional contractility levels allow distinction between ischemic and non-ischemic cases, where reduced contractility could be associated with damaged tissues.

<figure><img src="image_10.png" /><figcaption>**Figure 5.** Contractility levels (%), electrical delays (ms) estimated with the model-based approach, and transmurality degree (%) in ischemic patients obtained by cMRI.</figcaption></figure>

**4. Discussion**

In this paper we propose a model-based approach which replicates the intrinsic electromechanical properties of the myocardium observed in three real-life patients with typical LBBB undergoing CRT. Our model is also able to reproduce local alteration in LV contractility, which corresponds exactly to the areas of myocardial scar observed with cMRI.

These results were obtained by adjusting parameters such as the electrical periods of some automata or the active and passive parameters of the electro-mechanical driving function in order to simulate the LV deformation curves observed in patients with LBBB. In the case of ischemic patients, we varied the active and passive parameters of the electro-mechanical driving function to fit the experimental strain curves. At the end, we noticed
that the ischemic segments also presented reduced contractility, as expressed by the smaller active parameters.

## 4.1. Strain Patterns and Scar Localization in LBBB

The presence of abnormal LV contraction patterns in patients with LBBB has been largely described and referred to as septal flash and apical rocking [23]. The visual identification of septal flash and apical rocking has been associated with a better response to CRT and better survival in CRT recipients [24]. Nevertheless, the visual identification of these contraction patterns is largely influenced by the experience of the echo-reader and might be associated with poor reproducibility. Moreover, in recent years, some pathophysiological studies have shown that these specific contraction patterns are related to both the electrical delay due to the LBBB and specific LV mechanical properties [25,26]. Namely, the presence of the septal rebound stretch, which is a known predictor of CRT response [27,28], is influenced by the presence of lateral ischemia, which abolishes the septal rebound stretch and increases septal contractility. On the other hand, anteroseptal ischemia is associated with the persistence of septal rebound stretch and a decrease in septal contractility [25,29]. Our computational model was able to disclose the characteristics of the LV deformation patterns observed in patients with typical LBBB of both ischemic and non-ischemic etiology.

Another important property of our model is its capability to reproduce LV contractility. The area of reduced contractility in ischemic patients corresponded to the localization of LV scar in patients with previous myocardial infarction. This has pivotal importance because it is well known that the localization of scar in the lateral wall can be associated with a lower rate of CRT response because of the inefficient stimulation of an ischemic territory [30]. On the other hand, Aalen et al. have recently shown that the assessment of septal viability is also important because septal recruitment through CRT delivery significantly contributes to LV reverse remodeling after CRT [31]. This seems particularly important because some patients with septal scar might have a mechanical activation pattern similar to those observed in the case of non-ischemic LBBB, which might complicate the identification of CRT responders [32].

## 4.2. Modelling and Personalized Medicine

In our study, all patients had typical LBBB of similar duration and had similar LVEF, so they were all good CRT candidates according to current recommendations. Nevertheless, the proposed model-based approach was able to reproduce the electro-mechanical delay typical of each patient and simulate regional alterations in LV contractility by the simple analysis of ECG and strain traces. The corresponding local alteration in LV performance largely corresponded with the area of scar identified by cMRI in patients with ischemic cardiomyopathy. This process might be of use in the pre-implantation phase, to aid interpretation of the pathophysiological substrate associated with the development of LBBB and potentially improve the identification of CRT responders. Therefore, the application of a model-based approach could provide additional information on the regional electrical and mechanical function of the LV and can help to disclose the intrinsic complexity of LV mechanics in CRT candidates, representing a step forward in the development of personalized LV modelling in the field of CRT.

Although previous models [7,8] successfully represented a patient diversity in LV dyssynchrony and proposed keys to understand CRT response [33], our approach provides interesting advantages and originalities. This model proposes a sub-model of the right and left ventricle, which allows not only the analysis of septal and lateral strain curves but also a complete analysis of all the ventricular regions. Contrary to the models based on the finite element method, our approach requires low computational resources, facilitating personalized therapy. Therefore, the level of detail of the model facilitates translation to clinical practice. More interestingly, our approach seems able to simulate the electrical and mechanical properties observed in patients with both ischemic and non-ischemic LBBB, in a unique and reliable way. This might improve the interpretation of diagnostic
measures, provide interesting insight to the comprehension of the pathophysiological
mechanisms associated with LBBB, and provide an alternative to multimodality assessment
in CRT candidates.

4.3. *Limitations and Perspectives*

The main aims of computer modelling in the field of CRT are to (1) ease the com-
prehension of the pathophysiological mechanisms related to LBBB, (2) identify the patho-
physiological substrates that are associated with CRT response, and (3) eventually help
in the planification of therapy. The ideal model will be able to fully reproduce the LV
electromechanical properties observed in a CRT candidate and correctly simulate the effect
of CRT delivery before its effective release.

According to this goal, the main limitation of our model is the lack of parameter
identification, which prevents us from proposing a patient-specific model, as developed in
previous work of our team [11].

Moreover, the development of computational models for biological structures should
always balance the computational costs that are necessary to obtain a detailed model
with the necessary simplification process. This is also evident in our work because the
simplification and assumptions made to propose a computational model for LBBB imply
several limitations. At the tissue level, the segmental modelling of the ventricles induces
a less-accurate representation of the electrical and mechanical behaviors such as fiber
torsion or a complete mechanical continuity. Despite these simplifications, we were able
to reproduce distinct strain patterns of very different patients. To consider a clinical
application for our model, we should develop an automatic and robust identification
process. As such, the next step of our research should be parameter identification and
the validation of our results in a larger clinical dataset. After this phase, the information
provided by the model will be investigated to propose keys to better understand patients'
profiles and support the planification of therapy by a patient-specific characterization.

**5. Conclusions**

The present computational model of typical LBBB is tailored to patients’ characteristics and can simulate LV electrical, mechanical, and contractility patterns encountered in real-life patients. This might provide significant help in the interpretation of septal contraction patterns and facilitate the identification of patient’s characteristics that are associated with CRT efficacy.

**Author Contributions:** Conceptualization, V.L.R. and E.G.; data curation, K.O., V.L.R. and E.G.; formal analysis, M.T., K.O. and V.L.R.; funding acquisition, A.H. and V.L.R.; investigation, M.T., K.O., V.L.R. and E.G.; methodology, M.T., K.O., N.C., A.S., M.G., A.H., V.L.R. and E.G.; project administration, V.L.R. and E.G.; resources, E.D. and E.G.; software, M.T., K.O. and N.C.; supervision, A.S., M.G., A.H., E.D., V.L.R. and E.G.; validation, A.H., E.D., V.L.R. and E.G.; visualization, V.L.R. and E.G.; writing—original draft, M.T., K.O., V.L.R. and E.G.; writing—review and editing, V.L.R. and E.G. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by the French National Research Agency (ANR) (ANR-16-CE19-0008-01) (project MAESTRO) and the French Brittany Council (ADvICE project).

**Institutional Review Board Statement:** The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Ethics Committee (French Brittany Regional Ethic Committee validation number: 35RC14-9767.

**Informed Consent Statement:** Informed consent was obtained from all subjects involved in the study.

**Data Availability Statement:** Anonymized strain traces are available if requested.

**Conflicts of Interest:** The authors declare no conflict of interest.
## References

1. Glikson, M.; Nielsen, J.C.; Kronborg, M.B.; Michowitz, Y.; Auricchio, A.; Barbash, I.M.; Barrabés, J.A.; Boriani, G.; Braunschweig, F.; Brignole, M.; et al. 2021 ESC Guidelines on cardiac pacing and cardiac resynchronization therapy. *Eur. Heart J.* **2021**, *42*, 3427–3520. [CrossRef] [PubMed]

2. Chung, E.S.; Leon, A.R.; Tavazzi, L.; Sun, J.P.; Nihoyannopoulos, P.; Merlino, J.; Abraham, W.T.; Ghio, S.; Leclercq, C.; Bax, J.J.; et al. Results of the predictors of response to crt (prospect) trial. *Circulation* **2008**, *117*, 2608–2616. [CrossRef] [PubMed]

3. Lumens, J.; Tayal, B.; Walmsley, J.; Delgado-Montero, A.; Huntjens, P.R.; Schwartzman, D.; Althouse, A.D.; Delhaas, T.; Prinzen, F.W.; Gorcsan, J. Differentiating Electromechanical from Non-Electrical Substrates of Mechanical Discoordination to Identify Responders to Cardiac Resynchronization Therapy. *Circ. Cardiovasc. Imaging* **2015**, *8*, e003744. [CrossRef]

4. Le Rolle, V.; Galli, E.; Danan, D.; El Houari, K.; Hubert, A.; Donal, E.; Hernández, A.I. Sensitivity Analysis of a Left Ventricle Model in the Context of Intraventricular Dyssynchrony. *Acta Biotheor.* **2020**, *68*, 45–59. [CrossRef]

5. Constantino, J.; Hu, Y.; Trayanova, N.A. A computational approach to understanding the cardiac electromechanical activation sequence in the normal and failing heart, with translation to the clinical practice of CRT. *Prog. Biophys. Mol. Biol.* **2012**, *110*, 372–379. [CrossRef]

6. Crozier, A.; Blazevic, B.; Lamata, P.; Plank, G.; Ginks, M.; Duckett, S.; Sohal, M.; Shetty, A.; Rinaldi, C.A.; Razavi, R.; et al. The relative role of patient physiology and device optimisation in cardiac resynchronisation therapy: A computational modelling study. *J. Mol. Cell. Cardiol.* **2016**, *96*, 93–100. [CrossRef]

7. Van Everdingen, W.M.; Walmsley, J.; Cramer, M.J.; van Hagen, I.; De Boeck, B.W.L.; Meine, M.; Delhaas, T.; Doevendans, P.A.; Prinzen, F.W.; Lumens, J.; et al. Echocardiographic Prediction of Cardiac Resynchronization Therapy Response Requires Analysis of Both Mechanical Dyssynchrony and Right Ventricular Function: A Combined Analysis of Patient Data and Computer Simulations. *J. Am. Soc. Echocardiogr.* **2017**, *30*, 1012–1020. [CrossRef]

8. Lumens, J.; Leenders, G.E.; Cramer, M.J.; De Boeck, B.W.L.; Doevendans, P.A.; Prinzen, F.W.; Delhaas, T. Mechanistic evaluation of echocardiographic dyssynchrony indices patient data combined with multiscale computer simulations. *Circ. Cardiovasc. Imaging* **2012**, *5*, 491–499. [CrossRef]

9. Guerrero, G.; Lerolle, V.; Loiodice, C.; Amblard, A.; Pepin, J.-L.; Hernandez, A.I. Modeling patient-specific desaturation patterns in sleep apnea. *IEEE Trans. Biomed. Eng.* **2021**. [CrossRef]

10. Gallard, A.; Hubert, A.; Smiseth, O.; Voigt, J.U.; Le Rolle, V.; Leclercq, C.; Bidaut, A.; Galli, E.; Donal, E.; Hernandez, A.I. Prediction of response to cardiac resynchronization therapy using a multi-feature learning method. *Int. J. Cardiovasc. Imaging* **2020**, *37*, 989–998. [CrossRef]

11. Le Rolle, V.; Hernández, A.I.; Richard, P.Y.; Donal, E.; Carrault, G. Model-based analysis of myocardial strain data acquired by tissue Doppler imaging. *Artif. Intell. Med.* **2008**, *44*, 201–219. [CrossRef]

12. Owashi, K.; Galli, E.; Hubert, A.; Donal, E.; Hernandez, A.; Le Rolle, V. Patients-Specific Analysis of Myocardial Strains in Left Bundle Branch Block Based on Computational Models. In Proceedings of the 2020 Computing in Cardiology, Rimini, Italy, 13–16 September 2020. [CrossRef]

13. Lang, R.M.; Badano, L.P.; Victor, M.A.; Afilalo, J.; Armstrong, A.; Ernande, L.; Flachskampf, F.A.; Foster, E.; Goldstein, S.A.; Kuznetsova, T.; et al. Recommendations for cardiac chamber quantification by echocardiography in adults: An update from the American Society of Echocardiography and the European Association of Cardiovascular Imaging. *J. Am. Soc. Echocardiogr.* **2015**, *16*, 233–271. [CrossRef]

14. Kachenoura, N.; Redheuil, A.; Herment, A.; Mousseaux, E.; Frouin, F. Robust assessment of the transmural extent of myocardial infarction in late gadolinium-enhanced MRI studies using appropriate angular and circumferential subdivision of the myocardium. *Eur. Radiol.* **2008**, *18*, 2140–2147. [CrossRef] [PubMed]

15. Baron, N.; Kachenoura, N.; Cluzel, P.; Frouin, F.; Herment, A.; Grenier, P.; Montalescot, G.; Beygui, F. Comparison of various methods for quantitative evaluation of myocardial infarct volume from magnetic resonance delayed enhancement data. *Int. J. Cardiol.* **2013**, *167*, 739–744. [CrossRef] [PubMed]

16. Calvo, M.; Le Rolle, V.; Romero, D.; Béhar, N.; Gomis, P.; Mabo, P.; Hernández, A.I. Model-based analysis of the autonomic response to head-up tilt testing in Brugada syndrome. *Comput. Biol. Med.* **2018**, *103*, 82–92. [CrossRef]

17. Owashi, K.P.; Hubert, A.; Galli, E.; Donal, E.; Hernández, A.I.; Rolle, V.L. Model-based estimation of left ventricular pressure and myocardial work in aortic stenosis. *PLoS ONE* **2020**, *15*, e0229609. [CrossRef]

18. Rolle, V.L.; Hernández, A.I.; Richard, P.-Y.; Donal, E.; Carrault, G. A Tissue-Level Model of the Left Ventricle for the Analysis of Regional Myocardial Function. In Proceedings of the 2007 29th Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Lyon, France, 22–26 August 2007. [CrossRef]

19. Hernández, A.I.; Carrault, G.; Mora, F. Model-based interpretation of cardiac beats by evolutionary algorithms. *Comput. Cardiol.* **2002**, *29*, 93–96. [CrossRef]

20. Ve Koon, K.T.; Thebault, C.; Le Rolle, V.; Donal, E.; Hernández, A.I. Atrioventricular delay optimization in cardiac resynchronization therapy assessed by a computer model. *Comput. Cardiol.* **2010**, *37*, 333–336.

21. Cerqueira, M.D.; Weissman, N.J.; Dilsizian, V.; Jacobs, A.K.; Kaul, S.; Laskey, W.K.; Pennell, D.J.; Rumberger, J.A.; Ryan, T.J.; Verani, M.S. Standardized myocardial segmentation and nomenclature for tomographic imaging of the heart. *J. Cardiovasc. Magn. Reson.* **2002**, *4*, 203–210. [CrossRef]
22. Rolle, V.L.; Ojeda, D.; Hernández, A.I. Embedding a cardiac pulsatile model into an integrated model of the cardiovascular regulation for heart failure followup. *IEEE Trans. Biomed. Eng.* **2011**, *58*, 2982–2986. [CrossRef]

23. Parsai, C.; Bijnens, B.; Sutherland, G.R.; Baltabaeva, A.; Claus, P.; Marciniak, M.; Paul, V.; Scheffer, M.; Donal, E.; Derumeaux, G.; et al. Toward understanding response to cardiac resynchronization therapy: Left ventricular dyssynchrony is only one of multiple mechanisms. *Eur. Heart J.* **2009**, *30*, 940–949. [CrossRef] [PubMed]

24. Stankovic, I.; Prinz, C.; Ciarka, A.; Daraban, A.M.; Kotrc, M.; Aarones, M.; Szulik, M.; Winter, S.; Belmans, A.; Neskovic, A.N.; et al. Relationship of visually assessed apical rocking and septal flash to response and long-term survival following cardiac resynchronization therapy (PREDICT-CRT). *Eur. Heart J. Cardiovasc. Imaging* **2016**, *17*, 262–269. [CrossRef]

25. Leenders, G.E.; Lumens, J.; Cramer, M.J.; Delhaas, T.; Prinzen, F.W. Septal Deformation Patterns Delineate Mechanical Dyssynchrony and Regional Differences in Contractility Analysis of Patient Data Using a Computer Model. *Circ. Heart Fail.* **2011**, *5*, 87–96. [CrossRef]

26. Suga, H.; Sagawa, K.; Shoukas, A.A. Load independence of the instantaneous pressure-volume ratio of the canine left ventricle and effects of epinephrine and heart rate on the ratio. *Circ. Res.* **1973**, *32*, 314–322. [CrossRef]

27. Gorcsan, J.; Anderson, C.P.; Tayal, B.; Sugahara, M.; Walmsley, J.; Starling, R.C.; Lumens, J. Systolic Stretch Characterizes the Electromechanical Substrate Responsive to Cardiac Resynchronization Therapy. *JACC Cardiovasc. Imaging* **2019**, *12*, 1741–1752. [CrossRef]

28. Bohbot, Y.; De Meester de Ravenstein, C.; Chadha, G.; Rusinaru, D.; Belkhir, K.; Trouillet, C.; Pasquet, A.; Marechaux, S.; Vanoverschelde, J.L.; Tribouilloy, C. Relationship Between Left Ventricular Ejection Fraction and Mortality in Asymptomatic and Minimally Symptomatic Patients With Severe Aortic Stenosis. *JACC Cardiovasc. Imaging* **2019**, *12*, 38–48. [CrossRef]

29. Aalen, J.M.; Remme, E.W.; Larsen, C.K.; Andersen, O.S.; Krogh, M.; Duchenne, J.; Hopp, E.; Ross, S.; Beela, A.S.; Kongsgaard, E.; et al. Mechanism of Abnormal Septal Motion in Left Bundle Branch Block: Role of Left Ventricular Wall Interactions and Myocardial Scar. *JACC Cardiovasc. Imaging* **2019**, *12*, 2402–2413. [CrossRef]

30. Khan, F.Z.; Virdee, M.S.; Palmer, C.R.; Pugh, P.J.; O'Halloran, D.; Elsik, M.; Read, P.A.; Begley, D.; Fynn, S.P.; Dutka, D.P. Targeted left ventricular lead placement to guide cardiac resynchronization therapy: The TARGET study: A randomized, controlled trial. *J. Am. Coll. Cardiol.* **2012**, *59*, 1509–1518. [CrossRef]

31. Aalen, J.M.; Donal, E.; Larsen, C.K.; Duchenne, J.; Lederlin, M.; Cvijic, M.; Hubert, A.; Voros, G.; Leclercq, C.; Bogaert, J.; et al. Imaging predictors of response to cardiac resynchronization therapy: Left ventricular work asymmetry by echocardiography and septal viability by cardiac magnetic resonance. *Eur. Heart J.* **2020**, *41*, 3813–3823. [CrossRef]

32. Gjesdal, O.; Remme, E.W.; Opdahl, A.; Skulstad, H.; Russell, K.; Kongsgaard, E.; Edvardsen, T.; Smiseth, O.A. Mechanisms of abnormal systolic motion of the interventricular septum during left bundle-branch block. *Circ. Cardiovasc. Imaging* **2011**, *4*, 264–273. [CrossRef]

33. Willemen, E.; Schreurs, R.; Huntjens, P.R.; Strik, M.; Plank, G.; Vigmond, E.; Walmsley, J.; Vernooy, K.; Delhaas, T.; Prinzen, F.W.; et al. The left and right ventricles respond differently to variation of pacing delays in cardiac resynchronization therapy: A combined experimental-computational approach. *Front. Physiol.* **2019**, *10*, 17. [CrossRef] [PubMed]