Translating myosin-binding protein C and titin abnormalities to whole-heart function using a novel calcium-contraction coupling model.

Mutations in cardiac myosin-binding protein C (cMyBP-C) or titin may respectively lead to hypertrophic (HCM) or dilated (DCM) cardiomyopathies. The mechanisms leading to these phenotypes remain unclear because of the challenge of translating cellular abnormalities to whole-heart and system function. We developed and validated a novel computer model of calcium-contraction coupling incorporating the role of cMyBP-C and titin based on the key assumptions: 1) tension in the thick filament promotes cross-bridge attachment mechanochemically, 2) with increasing titin tension, more myosin heads are unlocked for attachment, and 3) cMyBP-C suppresses cross-bridge attachment. Simulated stationary calcium-tension curves, isotonic and isometric contractions, and quick release agreed with experimental data. The model predicted that a loss of cMyBP-C function decreases the steepness of the calcium-tension curve, and that more compliant titin decreases the level of passive and active tension and its dependency on sarcomere length. Integrating this cellular model in the CircAdapt model of the human heart and circulation showed that a loss of cMyBP-C function resulted in HCM-like hemodynamics with higher left ventricular end-diastolic pressures and smaller volumes. More compliant titin led to higher diastolic pressures and ventricular dilation, suggesting DCM-like hemodynamics. The novel model of calcium-contraction coupling incorporates the role of cMyBP-C and titin. Its coupling to whole-heart mechanics translates changes in cellular calcium-contraction coupling to changes in cardiac pump and circulatory function and identifies potential mechanisms by which cMyBP-C and titin abnormalities may develop into HCM and DCM phenotypes. This modeling platform may help identify distinct mechanisms underlying clinical phenotypes in cardiac diseases.

How to assess and treat right ventricular electromechanical dyssynchrony in post-repair tetralogy of Fallot: insights from imaging, invasive studies, and computational modelling.

BACKGROUND AND AIMS: Right bundle branch block (RBBB) and resulting right ventricular (RV) electromechanical discoordination are thought to play a role in the disease process of subpulmonary RV dysfunction that frequently occur post-repair tetralogy of Fallot (ToF). We sought to describe this disease entity, the role of pulmonary re-valvulation, and the potential added value of RV cardiac resynchronization therapy (RV-CRT). METHODS: Two patients with repaired ToF, complete RBBB, pulmonary regurgitation, and significantly decreased RV function underwent echocardiography, cardiac magnetic resonance, and an invasive study to evaluate the potential for RV-CRT as part of the management strategy. The data were used to personalize the CircAdapt model of the human heart and circulation. Resulting Digital Twins were analysed to quantify the relative effects of RV pressure and volume overload and to predict the effect of RV-CRT. RESULTS: Echocardiography showed components of a classic RV dyssynchrony pattern which could be reversed by RV-CRT during invasive study and resulted in acute improvement in RV systolic function. The Digital Twins confirmed a contribution of electromechanical RV dyssynchrony to RV dysfunction and suggested improvement of RV contraction efficiency after RV-CRT. The one patient who underwent successful permanent RV-CRT as part of the pulmonary re-valvulation procedure carried improvements that were in line with the predictions based on his Digital Twin. CONCLUSION: An integrative diagnostic approach to RV dysfunction, including the construction of Digital Twins may help to identify candidates for RV-CRT as part of the lifetime management of ToF and similar congenital heart lesions.

Parameter subset reduction for imaging-based digital twin generation of patients with left ventricular mechanical discoordination.

BACKGROUND: Integration of a patient's non-invasive imaging data in a digital twin (DT) of the heart can provide valuable insight into the myocardial disease substrates underlying left ventricular (LV) mechanical discoordination. However, when generating a DT, model parameters should be identifiable to obtain robust parameter estimations. In this study, we used the CircAdapt model of the human heart and circulation to find a subset of parameters which were identifiable from LV cavity volume and regional strain measurements of patients with different substrates of left bundle branch block (LBBB) and myocardial infarction (MI). To this end, we included seven patients with heart failure with reduced ejection fraction (HFrEF) and LBBB (study ID: 2018-0863, registration date: 2019-10-07), of which four were non-ischemic (LBBB-only) and three had previous MI (LBBB-MI), and six narrow QRS patients with MI (MI-only) (study ID: NL45241.041.13, registration date: 2013-11-12). Morris screening method (MSM) was applied first to find parameters which were important for LV volume, regional strain, and strain rate indices. Second, this parameter subset was iteratively reduced based on parameter identifiability and reproducibility. Parameter identifiability was based on the diaphony calculated from quasi-Monte Carlo simulations and reproducibility was based on the intraclass correlation coefficient ( ICC ) obtained from repeated parameter estimation using dynamic multi-swarm particle swarm optimization. Goodness-of-fit was defined as the mean squared error ( χ 2 ) of LV myocardial strain, strain rate, and cavity volume. RESULTS: A subset of 270 parameters remained after MSM which produced high-quality DTs of all patients ( χ 2 < 1.6), but minimum parameter reproducibility was poor ( ICC min = 0.01). Iterative reduction yielded a reproducible ( ICC min = 0.83) subset of 75 parameters, including cardiac output, global LV activation duration, regional mechanical activation delay, and regional LV myocardial constitutive properties. This reduced subset produced patient-resembling DTs ( χ 2 < 2.2), while septal-to-lateral wall workload imbalance was higher for the LBBB-only DTs than for the MI-only DTs (p < 0.05). CONCLUSIONS: By applying sensitivity and identifiability analysis, we successfully determined a parameter subset of the CircAdapt model which can be used to generate imaging-based DTs of patients with LV mechanical discoordination. Parameters were reproducibly estimated using particle swarm optimization, and derived LV myocardial work distribution was representative for the patient's underlying disease substrate. This DT technology enables patient-specific substrate characterization and can potentially be used to support clinical decision making.

Electrophysiological effects of stretch-activated ion channels: a systematic computational characterization.

Cardiac electrophysiology and mechanics are strongly interconnected. Their interaction is, among others, mediated by mechano-electric feedback through stretch-activated ion channels (SACs). The electrophysiological changes induced by SACs may contribute to arrhythmogenesis, but the precise SAC-induced electrophysiological changes remain incompletely understood. Here, we provide a systematic characterization of stretch effects through three distinguished SACs on cardiac electrophysiology using computational modelling. We implemented potassium-selective, calcium-selective and non-selective SACs in the Tomek-Rodriguez-O'Hara-Rudy model of human ventricular electrophysiology. The model was calibrated to experimental data from isolated cardiomyocytes undergoing stretch, considering inter-species differences, and disease-related remodelling of SACs. SAC-mediated effects on the action potential (AP) were analysed by varying stretch amplitude, application timing and/or duration. Afterdepolarizations of different amplitudes were observed with transient 10-ms stretch stimuli of 15-18% applied during phase 4, while stretch ≥18% during phase 4 elicited triggered APs. Longer stimuli shifted the threshold of AP trigger during phase 4 to lower amplitudes, while shorter stimuli increased it. Continuous stretch provoked electrophysiological remodelling. Furthermore, stretch shortened duration or changed morphology of a subsequent electrically evoked AP, and, if applied during a vulnerable time window with sufficient amplitude, prevented its occurrence because of stretch-induced modulation of sodium and L-type calcium channel gating. These effects were more pronounced with disease-related SAC remodelling due to increased stretch sensitivity of diseased hearts. We showed that SACs may induce afterdepolarizations and triggered activities, and prevent subsequent AP generation or change its morphology. These effects depend on cardiomyocyte stretch characteristics and disease-related SACs remodelling and may contribute to cardiac arrhythmogenesis. KEY POINTS: The interplay between cardiac electrophysiology and mechanics is mediated by mechano-electric feedback through stretch-activated ion channels (SACs). These channels may be pro-arrhythmic, but their precise effect on electrophysiology remains unclear. Here we present a systematic in silico characterization of stretch effects through three SACs by implementing inter-species differences as well as disease-related remodelling of SACs in a novel computational model of human ventricular cardiomyocyte electrophysiology. Our simulations showed that, at the cellular level, SACs may provoke electrophysiological remodelling, afterdepolarizations, triggered activities, change the morphology or shorten subsequent electrically evoked action potentials. The model further suggests that a vulnerable window exists in which stretch prevents the following electrically triggered beat occurrence. The pro-arrhythmic effects of stretch strongly depend on disease-related SAC remodelling as well as on stretch characteristics, such as amplitude, time of application and duration. Our study helps in understanding the role of stretch in cardiac arrhythmogenesis and revealing the underlying cellular mechanisms.

Comparison of novel ventricular pacing strategies using an electro-mechanical simulation platform.

AIMS: Focus of pacemaker therapy is shifting from right ventricular (RV) apex pacing (RVAP) and biventricular pacing (BiVP) to conduction system pacing. Direct comparison between the different pacing modalities and their consequences to cardiac pump function is difficult, due to the practical implications and confounding variables. Computational modelling and simulation provide the opportunity to compare electrical, mechanical, and haemodynamic consequences in the same virtual heart. METHODS AND RESULTS: Using the same single cardiac geometry, electrical activation maps following the different pacing strategies were calculated using an Eikonal model on a three-dimensional geometry, which were then used as input for a lumped mechanical and haemodynamic model (CircAdapt). We then compared simulated strain, regional myocardial work, and haemodynamic function for each pacing strategy. Selective His-bundle pacing (HBP) best replicated physiological electrical activation and led to the most homogeneous mechanical behaviour. Selective left bundle branch (LBB) pacing led to good left ventricular (LV) function but significantly increased RV load. RV activation times were reduced in non-selective LBB pacing (nsLBBP), reducing RV load but increasing heterogeneity in LV contraction. LV septal pacing led to a slower LV and more heterogeneous LV activation than nsLBBP, while RV activation was similar. BiVP led to a synchronous LV-RV, but resulted in a heterogeneous contraction. RVAP led to the slowest and most heterogeneous contraction. Haemodynamic differences were small compared to differences in local wall behaviour. CONCLUSION: Using a computational modelling framework, we investigated the mechanical and haemodynamic outcome of the prevailing pacing strategies in hearts with normal electrical and mechanical function. For this class of patients, nsLBBP was the best compromise between LV and RV function if HBP is not possible.

Virtual pacing of a patient's digital twin to predict left ventricular reverse remodelling after cardiac resynchronization therapy.

AIMS: Identifying heart failure (HF) patients who will benefit from cardiac resynchronization therapy (CRT) remains challenging. We evaluated whether virtual pacing in a digital twin (DT) of the patient's heart could be used to predict the degree of left ventricular (LV) reverse remodelling post-CRT. METHODS AND RESULTS: Forty-five HF patients with wide QRS complex (≥130 ms) and reduced LV ejection fraction (≤35%) receiving CRT were retrospectively enrolled. Echocardiography was performed before (baseline) and 6 months after CRT implantation to obtain LV volumes and 18-segment longitudinal strain. A previously developed algorithm was used to generate 45 DTs by personalizing the CircAdapt model to each patient's baseline measurements. From each DT, baseline septal-to-lateral myocardial work difference (MWLW-S,DT) and maximum rate of LV systolic pressure rise (dP/dtmax,DT) were derived. Biventricular pacing was then simulated using patient-specific atrioventricular delay and lead location. Virtual pacing-induced changes ΔMWLW-S,DT and ΔdP/dtmax,DT were correlated with real-world LV end-systolic volume change at 6-month follow-up (ΔLVESV). The DT's baseline MWLW-S,DT and virtual pacing-induced ΔMWLW-S,DT were both significantly associated with the real patient's reverse remodelling ΔLVESV (r = -0.60, P < 0.001 and r = 0.62, P < 0.001, respectively), while correlation between ΔdP/dtmax,DT and ΔLVESV was considerably weaker (r = -0.34, P = 0.02). CONCLUSION: Our results suggest that the reduction of septal-to-lateral work imbalance by virtual pacing in the DT can predict real-world post-CRT LV reverse remodelling. This DT approach could prove to be an additional tool in selecting HF patients for CRT and has the potential to provide valuable insights in optimization of CRT delivery.

Manual vs. automatic assessment of the QT-interval and corrected QT.

AIMS: Sudden cardiac death (SCD) is challenging to predict. Electrocardiogram (ECG)-derived heart rate-corrected QT-interval (QTc) is used for SCD-risk assessment. QTc is preferably determined manually, but vendor-provided automatic results from ECG recorders are convenient. Agreement between manual and automatic assessments is unclear for populations with aberrant QTc. We aimed to systematically assess pairwise agreement of automatic and manual QT-intervals and QTc. METHODS AND RESULTS: A multi-centre cohort enriching aberrant QTc comprised ECGs of healthy controls and long-QT syndrome (LQTS) patients. Manual QT-intervals and QTc were determined by the tangent and threshold methods and compared to automatically generated, vendor-provided values. We assessed agreement globally by intra-class correlation coefficients and pairwise by Bland-Altman analyses and 95% limits of agreement (LoA). Further, manual results were compared to a novel automatic QT-interval algorithm. ECGs of 1263 participants (720 LQTS patients; 543 controls) were available [median age 34 (inter-quartile range 35) years, 55% women]. Comparing cohort means, automatic and manual QT-intervals and QTc were similar. However, pairwise Bland-Altman-based agreement was highly discrepant. For QT-interval, LoAs spanned 95 (tangent) and 92 ms (threshold), respectively. For QTc, the spread was 108 and 105 ms, respectively. LQTS patients exhibited more pronounced differences. For automatic QTc results from 440-540 ms (tangent) and 430-530 ms (threshold), misassessment risk was highest. Novel automatic QT-interval algorithms may narrow this range. CONCLUSION: Pairwise vendor-provided automatic and manual QT-interval and QTc results can be highly discrepant. Novel automatic algorithms may improve agreement. Within the above ranges, automatic QT-interval and QTc results require manual confirmation, particularly if T-wave morphology is challenging.

Synchronization of repolarization after cardiac resynchronization therapy: A combined clinical and modeling study.

INTRODUCTION: The changes in ventricular repolarization after cardiac resynchronization therapy (CRT) are poorly understood. This knowledge gap is addressed using a multimodality approach including electrocardiographic and echocardiographic measurements in patients and using patient-specific computational modeling. METHODS: In 33 patients electrocardiographic and echocardiographic measurements were performed before and at various intervals after CRT, both during CRT-ON and temporary CRT-OFF. T-wave area was calculated from vectorcardiograms, and reconstructed from the 12-lead electrocardiography (ECG). Computer simulations were performed using a patient-specific eikonal model of cardiac activation with spatially varying action potential duration (APD) and repolarization rate, fit to a patient's ECG. RESULTS: During CRT-ON T-wave area diminished within a day and remained stable thereafter, whereas QT-interval did not change significantly. During CRT-OFF T-wave area doubled within 5 days of CRT, while QT-interval and peak-to-end T-wave interval hardly changed. Left ventricular (LV) ejection fraction only increased significantly increased after 1 month of CRT. Computer simulations indicated that the increase in T-wave area during CRT-OFF can be explained by changes in APD following chronic CRT that are opposite to the change in CRT-induced activation time. These APD changes were associated with a reduction in LV dispersion in repolarization during chronic CRT. CONCLUSION: T-wave area during CRT-OFF is a sensitive marker for adaptations in ventricular repolarization during chronic CRT that may include a reduction in LV dispersion of repolarization.

Cardiovascular fetal-to-neonatal transition: an in silico model.

BACKGROUND: Previous models describing the fetal-to-neonatal transition often lack oxygen saturation levels, homeostatic control mechanisms, phasic hemodynamic signals, or describe the heart with a time-varying elastance model. METHODS: We incorporated these elements in the adapted CircAdapt model with the one-fiber model for myocardial contraction, to simulate the hemodynamics of the healthy term human fetal circulation and its transition during the first 24 h after birth. The fetal-to-neonatal model was controlled by a time- and event-based script of changes occurring at birth, such as lung aeration and umbilical cord clamping. Model parameters were based on and validated with human and animal data. RESULTS: The fetal circulation showed low pulmonary blood flow, right ventricular dominance, and inverted mitral and tricuspid flow velocity patterns, as well as high mean ductus venosus flow velocity. The neonatal circulation showed oxygen saturation levels to gradually increase to 98% in the first 15 min after birth as well as temporary left ventricular volume overload. CONCLUSIONS: Hemodynamics of the term fetus and 24-h-old neonate, as well as the events occurring directly after birth and the transition during the first 24 h after birth, were realistically represented, allowing the model to be used for educational purposes and future research. IMPACT: With the addition of oxygen saturation levels, homeostatic pressure-flow control mechanisms, and the one-fiber model for myocardial contraction, a new closed-loop cardiovascular model was constructed to give more insight into the healthy term human fetal circulation and its cardiovascular transition during the first 24 h after birth. Extensive validation confirmed that the hemodynamics of the term fetus and the fetal-to-neonatal transition were realistically represented with the model. This well-validated and versatile model can serve as an education as well as a research platform for in silico investigation of fetal-to-neonatal hemodynamic changes under a wide range of physiological and pathophysiological conditions.

Hemodynamics-driven mathematical model of first and second heart sound generation.

NEW & NOTEWORTHY: To the best of our knowledge, this is the first hemodynamic-based heart sound generation model embedded in a complete real-time computational model of the cardiovascular system. Simulated heart sounds are similar to experimental and clinical measurements, both quantitatively and qualitatively. Our model can be used to investigate the relationships between heart sound acoustic features and hemodynamic factors/anatomical parameters.

Pacing therapy for atrioventricular dromotropathy: a combined computational-experimental-clinical study.

AIMS: Investigate haemodynamic effects, and their mechanisms, of restoring atrioventricular (AV)-coupling using pacemaker therapy in normal and failing hearts in a combined computational-experimental-clinical study. METHODS AND RESULTS: Computer simulations were performed in the CircAdapt model of the normal and failing human heart and circulation. Experiments were performed in a porcine model of AV dromotropathy. In a proof-of-principle clinical study, left ventricular (LV) pressure and volume were measured in 22 heart failure (HF) patients (LV ejection fraction <35%) with prolonged PR interval (>230 ms) and narrow or non-left bundle branch block QRS complex. Computer simulations and animal studies in normal hearts showed that restoring of AV-coupling with unchanged ventricular activation sequence significantly increased LV filling, mean arterial pressure, and cardiac output by 10-15%. In computer simulations of failing hearts and in HF patients, reducing PR interval by biventricular (BiV) pacing (patients: from 300 ± 61 to 137 ± 30 ms) resulted in significant increases in LV stroke volume and stroke work (patients: 34 ± 40% and 26 ± 31%, respectively). However, worsening of ventricular dyssynchrony by using right ventricular (RV) pacing abrogated the benefit of restoring AV-coupling. In model simulations, animals and patients, the increase of LV filling and associated improvement of LV pump function coincided with both larger mitral inflow (E- and A-wave area) and reduction of diastolic mitral regurgitation. CONCLUSION: Restoration of AV-coupling by BiV pacing in normal and failing hearts with prolonged AV conduction leads to considerable haemodynamic improvement. These results indicate that BiV or physiological pacing, but not RV pacing, may improve cardiac function in patients with HF and prolonged PR interval.

Chorioamnionitis induces changes in ovine pulmonary endogenous epithelial stem/progenitor cells in utero.

BACKGROUND: Chorioamnionitis, an intrauterine infection of the placenta and fetal membranes, is a common risk factor for adverse pulmonary outcomes in premature infants including BPD, which is characterized by an arrest in alveolar development. As endogenous epithelial stem/progenitor cells are crucial for organogenesis and tissue repair, we examined whether intrauterine inflammation negatively affects these essential progenitor pools. METHODS: In an ovine chorioamnionitis model, fetuses were intra-amniotically exposed to LPS, 2d or 7d (acute inflammation) before preterm delivery at 125d of gestation, or to intra-amniotic Ureaplasma parvum for 42d (chronic inflammation). Lung function, pulmonary endogenous epithelial stem/progenitor pools, and downstream functional markers were studied. RESULTS: Lung function was improved in the 7d LPS and 42d Ureaplasma groups. However, intrauterine inflammation caused a loss of P63+ basal cells in proximal airways and reduced SOX-9 expression and TTF-1+ Club cells in distal airways. Attenuated type-2 cell numbers were associated with lower proliferation and reduced type-1 cell marker Aqp5 expression, indicative for impaired progenitor function. Chronic Ureaplasma infection only affected distal airways, whereas acute inflammation affected stem/progenitor populations throughout the lungs. CONCLUSIONS: Acute and chronic prenatal inflammation improve lung function at the expense of stem/progenitor alterations that potentially disrupt normal lung development, thereby predisposing to adverse postnatal outcomes. IMPACT: In this study, prenatal inflammation improved lung function at the expense of stem/progenitor alterations that potentially disrupt normal lung development, thereby predisposing to adverse postnatal outcomes. Importantly, we demonstrate that these essential alterations can already be initiated before birth. So far, stem/progenitor dysfunction has only been shown postnatally. This study indicates that clinical protocols to target the consequences of perinatal inflammatory stress for the immature lungs should be initiated as early as possible and ideally in utero. Within this context, our data suggest that interventions, which promote function or repair of endogenous stem cells in the lungs, hold great promise.

Exploring the cause of conduction delays in patients with repaired Tetralogy of Fallot.

AIMS: Cardiac dyssynchrony in patients with repaired Tetralogy of Fallot (rToF) has been attributed to right bundle branch block (RBBB), fibrosis and/or the patches that are inserted during repair surgery. We aimed to investigate the basis of abnormal activation in rToF patients by mapping the electrical activation sequence during sinus rhythm (SR) and right ventricular (RV) pacing. METHODS AND RESULTS: A total of 17 patients were studied [13 with rToF, 2 with left bundle branch block (LBBB), and 2 without RBBB or LBBB (non-BBB)] during medically indicated cardiac surgery. During SR and RV pacing, measurements were performed using 112-electrode RV endocardial balloons (rToF only) and biventricular epicardial sock arrays (four of the rToF and all non-rToF patients). During SR, functional lines of block occurred in five rToF patients, while RV pacing caused functional blocks in four rToF patients. The line of block persisted during both SR and RV pacing in only 2 out of 13 rToF patients. Compared to SR, RV pacing increased dispersion of septal activation, but not dispersion of endocardial and epicardial activation of the RV free wall. During pacing, RV and left ventricular activation dispersion in rToF patients were comparable to that of the non-rToF patients. CONCLUSION: The results of the present study indicate that the delayed activation in the right ventricle of rToF patients is predominantly due to block(s) in the Purkinje system and that conduction in RV tissue is fairly normal.

Electromechanical substrate characterization in arrhythmogenic cardiomyopathy using imaging-based patient-specific computer simulations.

AIMS: Arrhythmogenic cardiomyopathy (AC) is an inherited cardiac disease, characterized by life-threatening ventricular arrhythmias and progressive cardiac dysfunction. The aim of this study is to use computer simulations to non-invasively estimate the individual patient's myocardial tissue substrates underlying regional right ventricular (RV) deformation abnormalities in a cohort of AC mutation carriers. METHODS AND RESULTS: In 68 AC mutation carriers and 20 control subjects, regional longitudinal deformation patterns of the RV free wall (RVfw), interventricular septum (IVS), and left ventricular free wall (LVfw) were obtained using speckle-tracking echocardiography. We developed and used a patient-specific parameter estimation protocol based on the multi-scale CircAdapt cardiovascular system model to create virtual AC subjects. Using the individual's deformation data as model input, this protocol automatically estimated regional RVfw and global IVS and LVfw tissue properties. The computational model was able to reproduce clinically measured regional deformation patterns for all subjects, with highly reproducible parameter estimations. Simulations revealed that regional RVfw heterogeneity of both contractile function and compliance were increased in subjects with clinically advanced disease compared to mutation carriers without clinically established disease (17 ± 13% vs. 8 ± 4%, P = 0.01 and 18 ± 11% vs. 10 ± 7%, P < 0.01, respectively). No significant difference in activation delay was found. CONCLUSION: Regional RV deformation abnormalities in AC mutation carriers were related to reduced regional contractile function and tissue compliance. In clinically advanced disease stages, a characteristic apex-to-base heterogeneity of tissue abnormalities was present in the majority of the subjects, with most pronounced disease in the basal region of the RVfw.

Pulmonary vein isolation in a real-world population does not influence QTc interval.

AIMS: We aimed to examine whether routine pulmonary vein isolation (PVI) induces significant ventricular repolarization changes as suggested earlier. METHODS AND RESULTS: Five-minute electrocardiograms were recorded at hospital's admission (T-1d), 1 day after the PVI-procedure (T+1d) and at 3 months post-procedure (T+3m) from a registry of consecutive atrial fibrillation (AF) patients scheduled for routine PVI with different PVI modalities (radiofrequency, cryo-ablation, and hybrid). Only patients who were in sinus rhythm at all three recordings (n = 117) were included. QT-intervals and QT-dispersion were evaluated with custom-made software and QTc was calculated using Bazett's, Fridericia's, Framingham's, and Hodges' formulas. Both QT- and RR-intervals were significantly shorter at T+1d (399 ± 37 and 870 ± 141 ms) and T+3m (407 ± 36 and 950 ± 140 ms) compared with baseline (417 ± 36 and 1025 ± 164 ms). There was no statistically significant within-subject difference in QTc Fridericia (T-1d 416 ± 28 ms, T+1d 419 ± 33 ms, and T+3m 414 ± 25 ms) and QT-dispersion (T-1d 18 ± 12 ms, T+1d 21 ± 19 ms, and T+3m 17 ± 12 ms) between the recordings. A multiple linear regression model with age, sex, AF type, ablation technique, first/re-do ablation, and AF recurrence to predict the change in QTc at T+3m with respect to QTc at T-1d did not reach significance which indicates that the change in QTc does not differ between all subgroups (age, sex, AF type, ablation technique, first/re-do ablation, and AF recurrence). CONCLUSION: Based on our data a routine PVI does not result in a prolongation of QTc in a real-world population. These findings, therefore, suggest that there is no need to intensify post-PVI QT-interval monitoring.

Incidence, prevalence, and trajectories of repetitive conduction patterns in human atrial fibrillation.

AIMS: Repetitive conduction patterns in atrial fibrillation (AF) may reflect anatomical structures harbouring preferential conduction paths and indicate the presence of stationary sources for AF. Recently, we demonstrated a novel technique to detect repetitive patterns in high-density contact mapping of AF. As a first step towards repetitive pattern mapping to guide AF ablation, we determined the incidence, prevalence, and trajectories of repetitive conduction patterns in epicardial contact mapping of paroxysmal and persistent AF patients. METHODS AND RESULTS: A 256-channel mapping array was used to record epicardial left and right AF electrograms in persistent AF (persAF, n = 9) and paroxysmal AF (pAF, n = 11) patients. Intervals containing repetitive conduction patterns were detected using recurrence plots. Activation movies, preferential conduction direction, and average activation sequence were used to characterize and classify conduction patterns. Repetitive patterns were identified in 33/40 recordings. Repetitive patterns were more prevalent in pAF compared with persAF [pAF: median 59%, inter-quartile range (41-72) vs. persAF: 39% (0-51), P < 0.01], larger [pAF: = 1.54 (1.15-1.96) vs. persAF: 1.16 (0.74-1.56) cm2, P < 0.001), and more stable [normalized preferentiality (0-1) pAF: 0.38 (0.25-0.50) vs. persAF: 0.23 (0-0.33), P < 0.01]. Most repetitive patterns were peripheral waves (87%), often with conduction block (69%), while breakthroughs (9%) and re-entries (2%) occurred less frequently. CONCLUSION: High-density epicardial contact mapping in AF patients reveals frequent repetitive conduction patterns. In persistent AF patients, repetitive patterns were less frequent, smaller, and more variable than in paroxysmal AF patients. Future research should elucidate whether these patterns can help in finding AF ablation targets.

The 'Digital Twin' to enable the vision of precision cardiology.

Providing therapies tailored to each patient is the vision of precision medicine, enabled by the increasing ability to capture extensive data about individual patients. In this position paper, we argue that the second enabling pillar towards this vision is the increasing power of computers and algorithms to learn, reason, and build the 'digital twin' of a patient. Computational models are boosting the capacity to draw diagnosis and prognosis, and future treatments will be tailored not only to current health status and data, but also to an accurate projection of the pathways to restore health by model predictions. The early steps of the digital twin in the area of cardiovascular medicine are reviewed in this article, together with a discussion of the challenges and opportunities ahead. We emphasize the synergies between mechanistic and statistical models in accelerating cardiovascular research and enabling the vision of precision medicine.

The Putative Role of Methylglyoxal in Arterial Stiffening: A Review.

BACKGROUND: Arterial stiffening is a hallmark of vascular ageing and a consequence of many diseases including diabetes mellitus. Methylglyoxal (MGO), a highly reactive α-dicarbonyl mainly formed during glycolysis, has emerged as a potential contributor to the development of arterial stiffness. MGO reacts with arginine and lysine residues in proteins to form stable advanced glycation endproducts (AGEs). AGEs may contribute to arterial stiffening by increased cross-linking of collagen within the extracellular matrix (ECM), by altering the vascular structure, and by triggering inflammatory and oxidative pathways. Although arterial stiffness is mainly determined by ECM and vascular smooth muscle cell function, the effects of MGO and MGO-derived AGEs on these structures have not been thoroughly reviewed to date. METHODS AND RESULTS: We conducted a PubMed search without filtering for publication date which resulted in 16 experimental and 22 clinical studies eligible for inclusion. Remarkably, none of the experimental and only three of the clinical studies specifically mentioned MGO-derived AGEs. Almost all studies reported an association between arterial stiffness and AGE accumulation in the arterial wall or increased plasma AGEs. Other studies report reduced arterial stiffness in experimental models upon administration of AGE-breakers. CONCLUSIONS: No papers published to date directly show an association between MGO or MGO-derived AGEs and arterial stiffening. The relevance of the various underlying mechanisms is not yet clear, which is particularly due to methodological challenges in the detection of MGO and MGO-derived AGEs at the molecular, intra- and pericellular, and structural levels, as well as in challenges in the assessment of intrinsic arterial wall properties at ECM- and tissue levels.

Differentiating the effects of ß-adrenergic stimulation and stretch on calcium and force dynamics using a novel electromechanical cardiomyocyte model.

Cardiac electrophysiology and mechanics are strongly interconnected. Calcium is crucial in this complex interplay through its role in cellular electrophysiology and sarcomere contraction. We aim to differentiate the effects of acute ß-adrenergic stimulation (ß-ARS) and cardiomyocyte stretch (increased sarcomere length) on calcium-transient dynamics and force generation, using a novel computational model of cardiac electromechanics. We implemented a bidirectional coupling between the O'Hara-Rudy model of human ventricular electrophysiology and the MechChem model of sarcomere mechanics through the buffering of calcium by troponin. The coupled model was validated using experimental data from large mammals or human samples. Calcium transient and force were simulated for various degrees of ß-ARS and initial sarcomere lengths. The model reproduced force-frequency, quick-release, and isotonic contraction experiments, validating the bidirectional electromechanical interactions. An increase in ß-ARS increased the amplitudes of force (augmented inotropy) and calcium transient, and shortened both force and calcium-transient duration (lusitropy). An increase in sarcomere length increased force amplitude even more, but decreased calcium-transient amplitude and increased both force and calcium-transient duration. Finally, a gradient in relaxation along the thin filament may explain the nonmonotonic decay in cytosolic calcium observed with high tension. Using a novel coupled human electromechanical model, we identified differential effects of ß-ARS and stretch on calcium and force. Stretch mostly contributed to increased force amplitude and ß-ARS to the reduction of calcium and force duration. We showed that their combination, rather than individual contributions, is key to ensure force generation, rapid relaxation, and low diastolic calcium levels.NEW & NOTEWORTHY This work identifies the contribution of electrical and mechanical alterations to regulation of calcium and force under exercise-like conditions using a novel human electromechanical model integrating ventricular electrophysiology and sarcomere mechanics. By better understanding their individual and combined effects, this can uncover arrhythmogenic mechanisms in exercise-like situations. This publicly available model is a crucial step toward understanding the complex interplay between cardiac electrophysiology and mechanics to improve arrhythmia risk prediction and treatment.

Large vessels as a tree of transmission lines incorporated in the CircAdapt whole-heart model: A computational tool to examine heart-vessel interaction.

We developed a whole-circulation computational model by integrating a transmission line (TL) model describing vascular wave transmission into the established CircAdapt platform of whole-heart mechanics. In the present paper, we verify the numerical framework of our TL model by benchmark comparison to a previously validated pulse wave propagation (PWP) model. Additionally, we showcase the integrated CircAdapt-TL model, which now includes the heart as well as extensive arterial and venous trees with terminal impedances. We present CircAdapt-TL haemodynamics simulations of: 1) a systemic normotensive situation and 2) a systemic hypertensive situation. In the TL-PWP benchmark comparison we found good agreement regarding pressure and flow waveforms (relative errors ≤ 2.9% for pressure, and ≤ 5.6% for flow). CircAdapt-TL simulations reproduced the typically observed haemodynamic changes with hypertension, expressed by increases in mean and pulsatile blood pressures, and increased arterial pulse wave velocity. We observed a change in the timing of pressure augmentation (defined as a late-systolic boost in aortic pressure) from occurring after time of peak systolic pressure in the normotensive situation, to occurring prior to time of peak pressure in the hypertensive situation. The pressure augmentation could not be observed when the systemic circulation was lumped into a (non-linear) three-element windkessel model, instead of using our TL model. Wave intensity analysis at the carotid artery indicated earlier arrival of reflected waves with hypertension as compared to normotension, in good qualitative agreement with findings in patients. In conclusion, we successfully embedded a TL model as a vascular module into the CircAdapt platform. The integrated CircAdapt-TL model allows detailed studies on mechanistic studies on heart-vessel interaction.