Human-based computational modelling and simulation are powerful tools to accelerate the mechanistic understanding of cardiac patho-physiology, and to develop and evaluate therapeutic interventions. The aim of this study is to calibrate and evaluate human ventricular electro-mechanical models for investigations on the effect of the electro-mechanical coupling and pharmacological action on human ventricular electrophysiology, calcium dynamics, and active contraction. The most recent models of human ventricular electrophysiology, excitation-contraction coupling, and active contraction were integrated, and the coupled models were calibrated using human experimental data. Simulations were then conducted using the coupled models to quantify the effects of electro-mechanical coupling and drug exposure on electrophysiology and force generation in virtual human ventricular cardiomyocytes and tissue. The resulting calibrated human electro-mechanical models yielded active tension, action potential, and calcium transient metrics that are in agreement with experiments for endocardial, epicardial, and mid-myocardial human samples. Simulation results correctly predicted the inotropic response of different multichannel action reference compounds and demonstrated that the electro-mechanical coupling improves the robustness of repolarisation under drug exposure compared to electrophysiology-only models. They also generated additional evidence to explain the partial mismatch between in-silico and in-vitro experiments on drug-induced electrophysiology changes. The human calibrated and evaluated modelling and simulation framework constructed in this study opens new avenues for future investigations into the complex interplay between the electrical and mechanical cardiac substrates, its modulation by pharmacological action, and its translation to tissue and organ models of cardiac patho-physiology.

In this work, we present a fully coupled fluid-electro-mechanical model of a 50th percentile human heart. The model is implemented on Alya, the BSC multi-physics parallel code, capable of running efficiently in supercomputers. Blood in the cardiac cavities is modeled by the incompressible Navier-Stokes equations and an arbitrary Lagrangian-Eulerian (ALE) scheme. Electrophysiology is modeled with a monodomain scheme and the O'Hara-Rudy cell model. Solid mechanics is modeled with a total Lagrangian formulation for discrete strains using the Holzapfel-Ogden cardiac tissue material model. The three problems are simultaneously and bidirectionally coupled through an electromechanical feedback and a fluid-structure interaction scheme. In this paper, we present the scheme in detail and propose it as a computational cardiac workbench.

No abstract

Aims Develop, calibrate and evaluate with clinical data a human electromechanical modelling and simulation framework for multiscale, mechanistic investigations in healthy and post-myocardial infarction (MI) conditions, from ionic to clinical biomarkers. Methods and results Human healthy and post-MI electromechanical simulations were conducted with a novel biventricular model, calibrated and evaluated with experimental and clinical data, including torso/biventricular anatomy from clinical magnetic resonance, state-of-the-art human-based membrane kinetics, excitation–contraction and active tension models, and orthotropic electromechanical coupling. Electromechanical remodelling of the infarct/ischaemic region and the border zone were simulated for ischaemic, acute, and chronic states in a fully transmural anterior infarct and a subendocardial anterior infarct. The results were compared with clinical electrocardiogram and left ventricular ejection fraction (LVEF) data at similar states. Healthy model simulations show LVEF 63%, with 11% peak systolic wall thickening, QRS duration and QT interval of 100 ms and 330 ms. LVEF in ischaemic, acute, and chronic post-MI states were 56%, 51%, and 52%, respectively. In linking the three post-MI simulations, it was apparent that elevated resting potential due to hyperkalaemia in the infarcted region led to ST-segment elevation, while a large repolarization gradient corresponded to T-wave inversion. Mechanically, the chronic stiffening of the infarct region had the benefit of improving systolic function by reducing infarct bulging at the expense of reducing diastolic function by inhibiting inflation. Conclusion Our human-based multiscale modelling and simulation framework enables mechanistic investigations into patho-physiological electrophysiological and mechanical behaviour and can serve as testbed to guide the optimization of pharmacological and electrical therapies.