Ventricular anatomical complexity and sex differences impact predictions from electrophysiological computational models

The aim of this work was to analyze the influence of sex hormones and anatomical details (trabeculations and false tendons) on the electrophysiology of healthy human hearts. Additionally, sex- and anatomy-dependent effects of ventricular tachycardia (VT) inducibility are presented. To this end, four anatomically normal, human, biventricular geometries (two male, two female), with identifiable trabeculations, were obtained from high-resolution, ex-vivo MRI and represented by detailed and smoothed geometrical models (with and without the trabeculations). Additionally one model was augmented by a scar. The electrophysiology finite element model (FEM) simulations were carried out, using O’Hara-Rudy human myocyte model with sex phenotypes of Yang and Clancy. A systematic comparison between detailed vs smooth anatomies, male vs female normal hearts was carried out. The heart with a myocardial infarction was subjected to a programmed stimulus protocol to identify the effects of sex and anatomical detail on ventricular tachycardia inducibility. All female hearts presented QT-interval prolongation however the prolongation interval in comparison to the male phenotypes was anatomy-dependent and was not correlated to the size of the heart. Detailed geometries showed QRS fractionation and increased T-wave magnitude in comparison to the corresponding smoothed geometries. A variety of sustained VTs were obtained in the detailed and smoothed male geometries at different pacing locations, which provide evidence of the geometry-dependent differences regarding the prediction of the locations of reentry channels. In the female phenotype, sustained VTs were induced in both detailed and smooth geometries with RV apex pacing, however no consistent reentry channels were identified. Anatomical and physiological cardiac features play an important role defining risk in cardiac disease. These are often excluded from cardiac electrophysiology simulations. The assumption that the cardiac endocardium is smooth may produce inaccurate predictions towards the location of reentry channels in in-silico tachycardia inducibility studiesJA-S, FS, GH and MV are supported by the European Union’s Horizon 2020 research and innovation programme under grant agreements No 675451 (Compbiomed project phase 1) and No 823712 (CompBioMed project, phase 2) and project No 777204 (SilicoFCM project). Part of the simulation computing hours were provided by the CompBioMed project phase 1. JA-S was awarded computation time from Red Espanola de Supercomputacion (RES). (Activity IDs: FI-2018-2-0049 and BCV-2019-2-0014) JA-S is funded by a Ramon y Cajal fellowship (RYC-2017-22532), Ministerio de Ciencia e Innovacion, Spain; and by Plan Estatal de Investigacion Cientifica y Tecnica y de Innovacion 2017-2020 from the Ministerio de Ciencia e Innovacion y Universidades (PID2019-104356RBC41/AEI/10.13039/501100011033): meHeart ME PID2019-104356RB-C44. CB is funded by the Torres Quevedo Program (PTQ2018-010290), Ministerio de Ciencia e Innovacion, Spain. MV, GH and CB are funded by the Spanish Neotec project EXP - 00123159/SNEO-20191113 Generador de corazones virtuales. LKGM was funded by Fundacion Carolina-BBVA. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. There was no additional external funding received for this study.Peer Reviewed"Article signat per 11 autors/es: Pablo Gonzalez-Martin,Federica Sacco,Constantine Butakoff,Ruben Doste,Carlos Bederian,Lilian K. Gutierrez Espinosa de los Monteros,Guillaume Houzeaux,Paul A. Iaizzo,Tinen L. Iles,Mariano Vazquez,Jazmin Aguado-Sierra"Postprint (published version

Mechanism based therapies enable personalised treatment of hypertrophic cardiomyopathy

Cardiomyopathies have unresolved genotype–phenotype relationships and lack disease-specific treatments. Here we provide a framework to identify genotype-specific pathomechanisms and therapeutic targets to accelerate the development of precision medicine. We use human cardiac electromechanical in-silico modelling and simulation which we validate with experimental hiPSC-CM data and modelling in combination with clinical biomarkers. We select hypertrophic cardiomyopathy as a challenge for this approach and study genetic variations that mutate proteins of the thick (MYH7R403Q/+) and thin filaments (TNNT2R92Q/+, TNNI3R21C/+) of the cardiac sarcomere. Using in-silico techniques we show that the destabilisation of myosin super relaxation observed in hiPSC-CMs drives disease in virtual cells and ventricles carrying the MYH7R403Q/+ variant, and that secondary effects on thin filament activation are necessary to precipitate slowed relaxation of the cell and diastolic insufficiency in the chamber. In-silico modelling shows that Mavacamten corrects the MYH7R403Q/+ phenotype in agreement with hiPSC-CM experiments. Our in-silico model predicts that the thin filament variants TNNT2R92Q/+ and TNNI3R21C/+ display altered calcium regulation as central pathomechanism, for which Mavacamten provides incomplete salvage, which we have corroborated in TNNT2R92Q/+ and TNNI3R21C/+ hiPSC-CMs. We define the ideal characteristics of a novel thin filament-targeting compound and show its efficacy in-silico. We demonstrate that hybrid human-based hiPSC-CM and in-silico studies accelerate pathomechanism discovery and classification testing, improving clinical interpretation of genetic variants, and directing rational therapeutic targeting and design

Digital twinning of the human ventricular activation sequence to clinical 12-lead ECGs and magnetic resonance imaging using realistic Purkinje networks for in silico clinical trials

Cardiac in silico clinical trials can virtually assess the safety and efficacy of therapies using human-based modelling and simulation. These technologies can provide mechanistic explanations for clinically observed pathological behaviour. Designing virtual cohorts for in silico trials requires exploiting clinical data to capture the physiological variability in the human population. The clinical characterisation of ventricular activation and the Purkinje network is challenging, especially non-invasively. Our study aims to present a novel digital twinning pipeline that can efficiently generate and integrate Purkinje networks into human multiscale biventricular models based on subject-specific clinical 12-lead electrocardiogram and magnetic resonance recordings. Essential novel features of the pipeline are the human-based Purkinje network generation method, personalisation considering ECG R wave progression as well as QRS morphology, and translation from reduced-order Eikonal models to equivalent biophysically-detailed monodomain ones. We demonstrate ECG simulations in line with clinical data with clinical image-based multiscale models with Purkinje in four control subjects and two hypertrophic cardiomyopathy patients (simulated and clinical QRS complexes with Pearson's correlation coefficients > 0.7). Our methods also considered possible differences in the density of Purkinje myocardial junctions in the Eikonal-based inference as regional conduction velocities. These differences translated into regional coupling effects between Purkinje and myocardial models in the monodomain formulation. In summary, we demonstrate a digital twin pipeline enabling simulations yielding clinically consistent ECGs with clinical CMR image-based biventricular multiscale models, including personalised Purkinje in healthy and cardiac disease conditions

Multiscale Modelling of β-Adrenergic Stimulation in Cardiac Electromechanical Function

β-adrenergic receptor stimulation (β-ARS) is a physiological mechanism that regulates cardiovascular function under stress conditions or physical exercise. Triggered during the so-called “fight-or-flight” response, the activation of the β-adrenergic receptors located on the cardiomyocyte membrane initiates a phosphorylation cascade of multiple ion channel targets that regulate both cellular excitability and recovery and of different proteins involved in intracellular calcium handling. As a result, β-ARS impacts both the electrophysiological and the mechanical response of the cardiomyocyte. β-ARS also plays a crucial role in several cardiac pathologies, greatly modifying cardiac output and potentially causing arrhythmogenic events. Mathematical patient-specific models are nowadays envisioned as an important tool for the personalised study of cardiac disease, the design of tailored treatments, or to inform risk assessment. Despite that, only a reduced number of computational studies of heart disease have incorporated β-ARS modelling. In this review, we describe the main existing multiscale frameworks to equip cellular models of cardiac electrophysiology with a β-ARS response. We also outline various applications of these multiscale frameworks in the study of cardiac pathology. We end with a discussion of the main current limitations and the future steps that need to be taken to adapt these models to a clinical environment and to incorporate them in organ-level simulations

Remodelling of potassium currents underlies arrhythmic action potential prolongation under beta-adrenergic stimulation in hypertrophic cardiomyopathy

Hypertrophic cardiomyopathy (HCM) patients often present an enhanced arrhythmogenicity that can lead to lethal arrhythmias, especially during exercise. Recent studies have indicated an abnormal response of HCM cardiomyocytes to β-adrenergic receptor stimulation (β-ARS), with prolongation of their action potential rather than shortening. The mechanisms underlying this aberrant response to sympathetic stimulation and its possible proarrhythmic role remain unknown. The aims of this study are to investigate the key ionic mechanisms underlying the HCM abnormal response to β-ARS and the resultant repolarisation abnormalities using human-based experimental and computational methodologies. We integrated and calibrated the latest models of human ventricular electrophysiology and β-ARS using experimental measurements of human adult cardiomyocytes from control and HCM patients. Our major findings include: (1) the developed in silico models of β-ARS capture the behaviour observed in the experimental data, including the aberrant response of HCM cardiomyocytes to β-ARS; (2) the reduced increase of potassium currents under β-ARS was identified as the main mechanism of action potential prolongation in HCM, rather than a more sustained inward calcium current; (3) action potential duration differences between healthy and HCM cardiomyocytes were increased upon β-ARS, while endocardial to epicardial differences in HCM cardiomyocytes were reduced; (4) models presenting repolarisation abnormalities were characterised by downregulation of the rapid delayed rectifier potassium current and the sodium‑potassium pump, while inward currents were upregulated. In conclusion, our results identify causal relationships between the HCM phenotype and its arrhythmogenic response to β-ARS through the downregulation of potassium currents

Pharmacological management of hypertrophic cardiomyopathy: From bench to bedside

Hypertrophic cardiomyopathy (HCM), the most common inherited heart disease, is still orphan of a specific drug treatment. The erroneous consideration of HCM as a rare disease has hampered the design and conduct of large, randomized trials in the last 50 years, and most of the indications in the current guidelines are derived from small non-randomized studies, case series, or simply from the consensus of experts. Guideline-directed therapy of HCM includes non-selective drugs such as disopyramide, non-dihydropyridine calcium channel blockers, or β-adrenergic receptor blockers, mainly used in patients with symptomatic obstruction of the outflow tract. Following promising preclinical studies, several drugs acting on potential HCM-specific targets were tested in patients. Despite the huge efforts, none of these studies was able to change clinical practice for HCM patients, because tested drugs were proven to be scarcely effective or hardly tolerated in patients. However, novel compounds have been developed in recent years specifically for HCM, addressing myocardial hypercontractility and altered energetics in a direct manner, through allosteric inhibition of myosin. In this paper, we will critically review the use of different classes of drugs in HCM patients, starting from “old” established agents up to novel selective drugs that have been recently trialed in patients

Training machine learning models with synthetic data improves the prediction of ventricular origin in outflow tract ventricular arrhythmias

In order to determine the site of origin (SOO) in outflow tract ventricular arrhythmias (OTVAs) before an ablation procedure, several algorithms based on manual identification of electrocardiogram (ECG) features, have been developed. However, the reported accuracy decreases when tested with different datasets. Machine learning algorithms can automatize the process and improve generalization, but their performance is hampered by the lack of large enough OTVA databases. We propose the use of detailed electrophysiological simulations of OTVAs to train a machine learning classification model to predict the ventricular origin of the SOO of ectopic beats. We generated a synthetic database of 12-lead ECGs (2,496 signals) by running multiple simulations from the most typical OTVA SOO in 16 patient-specific geometries. Two types of input data were considered in the classification, raw and feature ECG signals. From the simulated raw 12-lead ECG, we analyzed the contribution of each lead in the predictions, keeping the best ones for the training process. For feature-based analysis, we used entropy-based methods to rank the obtained features. A cross-validation process was included to evaluate the machine learning model. Following, two clinical OTVA databases from different hospitals, including ECGs from 365 patients, were used as test-sets to assess the generalization of the proposed approach. The results show that V2 was the best lead for classification. Prediction of the SOO in OTVA, using both raw signals or features for classification, presented high accuracy values (>0.96). Generalization of the network trained on simulated data was good for both patient datasets (accuracy of 0.86 and 0.84, respectively) and presented better values than using exclusively real ECGs for classification (accuracy of 0.84 and 0.76 for each dataset). The use of simulated ECG data for training machine learning-based classification algorithms is critical to obtain good SOO predictions in OTVA compared to real data alone. The fast implementation and generalization of the proposed methodology may contribute towards its application to a clinical routine.This work has been funded by Generalitat Valenciana Grant AICO/2021/318 (Consolidables 2021) and Grant PID2020-114291RB-I00 funded by MCIN/10.13039/501100011033 and by “ERDF A way of making Europe”

Ventricular anatomical complexity and sex differences impact predictions from electrophysiological computational models

The aim of this work was to analyze the influence of sex hormones and anatomical details (trabeculations and false tendons) on the electrophysiology of healthy human hearts. Additionally, sex- and anatomy-dependent effects of ventricular tachycardia (VT) inducibility are presented. To this end, four anatomically normal, human, biventricular geometries (two male, two female), with identifiable trabeculations, were obtained from high-resolution, ex-vivo MRI and represented by detailed and smoothed geometrical models (with and without the trabeculations). Additionally one model was augmented by a scar. The electrophysiology finite element model (FEM) simulations were carried out, using O’Hara-Rudy human myocyte model with sex phenotypes of Yang and Clancy. A systematic comparison between detailed vs smooth anatomies, male vs female normal hearts was carried out. The heart with a myocardial infarction was subjected to a programmed stimulus protocol to identify the effects of sex and anatomical detail on ventricular tachycardia inducibility. All female hearts presented QT-interval prolongation however the prolongation interval in comparison to the male phenotypes was anatomy-dependent and was not correlated to the size of the heart. Detailed geometries showed QRS fractionation and increased T-wave magnitude in comparison to the corresponding smoothed geometries. A variety of sustained VTs were obtained in the detailed and smoothed male geometries at different pacing locations, which provide evidence of the geometry-dependent differences regarding the prediction of the locations of reentry channels. In the female phenotype, sustained VTs were induced in both detailed and smooth geometries with RV apex pacing, however no consistent reentry channels were identified. Anatomical and physiological cardiac features play an important role defining risk in cardiac disease. These are often excluded from cardiac electrophysiology simulations. The assumption that the cardiac endocardium is smooth may produce inaccurate predictions towards the location of reentry channels in in-silico tachycardia inducibility studies.Fil: Gonzalez Martin, Pablo. Barcelona Supercomputing Center - Centro Nacional de Supercomputacion; EspañaFil: Sacco, Federica. Universitat Pompeu Fabra; España. Barcelona Supercomputing Center - Centro Nacional de Supercomputacion; EspañaFil: Butakoff, Constantine. Barcelona Supercomputing Center - Centro Nacional de Supercomputacion; EspañaFil: Doste, Ruben. University of Oxford; Reino UnidoFil: Bederian, Carlos Sergio. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Física Enrique Gaviola. Universidad Nacional de Córdoba. Instituto de Física Enrique Gaviola; ArgentinaFil: Gutierrez Espinosa de los Monteros, Lilian K.. Instituto de Salud Carlos III; EspañaFil: Houzeaux, Guillaume. Barcelona Supercomputing Center - Centro Nacional de Supercomputacion; EspañaFil: Iaizzo, Paul A.. University of Minnesota; Estados UnidosFil: Iles, Tinen L.. University of Minnesota; Estados UnidosFil: Vazquez, Mariano. Barcelona Supercomputing Center - Centro Nacional de Supercomputacion; EspañaFil: Aguado Sierra, Jazmin. Barcelona Supercomputing Center - Centro Nacional de Supercomputacion; Españ

Ventricular anatomical complexity and sex differences impact predictions from electrophysiological computational models.

The aim of this work was to analyze the influence of sex hormones and anatomical details (trabeculations and false tendons) on the electrophysiology of healthy human hearts. Additionally, sex- and anatomy-dependent effects of ventricular tachycardia (VT) inducibility are presented. To this end, four anatomically normal, human, biventricular geometries (two male, two female), with identifiable trabeculations, were obtained from high-resolution, ex-vivo MRI and represented by detailed and smoothed geometrical models (with and without the trabeculations). Additionally one model was augmented by a scar. The electrophysiology finite element model (FEM) simulations were carried out, using O'Hara-Rudy human myocyte model with sex phenotypes of Yang and Clancy. A systematic comparison between detailed vs smooth anatomies, male vs female normal hearts was carried out. The heart with a myocardial infarction was subjected to a programmed stimulus protocol to identify the effects of sex and anatomical detail on ventricular tachycardia inducibility. All female hearts presented QT-interval prolongation however the prolongation interval in comparison to the male phenotypes was anatomy-dependent and was not correlated to the size of the heart. Detailed geometries showed QRS fractionation and increased T-wave magnitude in comparison to the corresponding smoothed geometries. A variety of sustained VTs were obtained in the detailed and smoothed male geometries at different pacing locations, which provide evidence of the geometry-dependent differences regarding the prediction of the locations of reentry channels. In the female phenotype, sustained VTs were induced in both detailed and smooth geometries with RV apex pacing, however no consistent reentry channels were identified. Anatomical and physiological cardiac features play an important role defining risk in cardiac disease. These are often excluded from cardiac electrophysiology simulations. The assumption that the cardiac endocardium is smooth may produce inaccurate predictions towards the location of reentry channels in in-silico tachycardia inducibility studies