SUMMARY Organ-specific functions of tissue-resident macrophages in the steady-state heart are unknown. Here we show that cardiac macrophages facilitate electrical conduction through the distal atrioventricular node, where conducting cells densely intersperse with elongated macrophages expressing connexin 43. When coupled to spontaneously beating cardiomyocytes via connexin 43-containing gap junctions, cardiac macrophages have a negative resting membrane potential and depolarize in synchrony with cardiomyocytes. Conversely, macrophages render the resting membrane potential of cardiomyocytes more positive and, according to computational modeling, accelerate their repolarization. Photostimulation of channelrhodopsin 2-expressing macrophages improves atrioventricular conduction, while conditional deletion of connexin 43 in macrophages and congenital lack of macrophages delay atrioventricular conduction. In the Cd11bDTR mouse, macrophage ablation induces progressive atrioventricular block. These observations implicate macrophages in normal and aberrant cardiac conduction.

Mathematical models of the human heart are evolving to become a cornerstone of precision medicine and support clinical decision making by providing a powerful tool to understand the mechanisms underlying pathophysiological conditions. In this study, we present a detailed mathematical description of a fully coupled multi-scale model of the human heart, including electrophysiology, mechanics, and a closed-loop model of circulation. State-of-the-art models based on human physiology are used to describe membrane kinetics, excitation-contraction coupling and active tension generation in the atria and the ventricles. Furthermore, we highlight ways to adapt this framework to patient specific measurements to build digital twins. The validity of the model is demonstrated through simulations on a personalized whole heart geometry based on magnetic resonance imaging data of a healthy volunteer. Additionally, the fully coupled model was employed to evaluate the effects of a typical atrial ablation scar on the cardiovascular system. With this work, we provide an adaptable multi-scale model that allows a comprehensive personalization from ion channels to the organ level enabling digital twin modeling.

Electrophysiological modeling of cardiac tissue is commonly based on functional and structural properties measured in experiments. Our knowledge of these properties is incomplete, in particular their remodeling in disease. Here, we introduce a methodology for quantitative tissue characterization based on fluorescent labeling, 3-D scanning confocal microscopy, image processing and reconstruction of tissue micro-structure at sub-micrometer resolution. We applied this methodology to normal rabbit ventricular tissue and tissue from hearts with myocardial infarction. Our analysis revealed that the volume fraction of fibroblasts increased from 4.83 ± 0.42% (mean ± standard deviation) in normal tissue up to 6.51 ± 0.38% in myocardium from infarcted hearts. The myocyte volume fraction decreased from 76.20 ± 9.89% in normal to 73.48 ± 8.02% adjacent to the infarct. Numerical field calculations on 3-D reconstructions of the extracellular space yielded an extracellular longitudinal conductivity of 0.264 ± 0.082 S/m with an anisotropy ratio of 2.095 ± 1.11 in normal tissue. Adjacent to the infarct, the longitudinal conductivity increased up to 0.400 ± 0.051 S/m, but the anisotropy ratio decreased to 1.295 ± 0.09. Our study indicates an increased density of gap junctions proximal to both fibroblasts and myocytes in infarcted versus normal tissue, supporting previous hypotheses of electrical coupling of fibroblasts and myocytes in infarcted hearts. We suggest that the presented methodology provides an important contribution to modeling normal and diseased tissue. Applications of the methodology include the clinical characterization of disease-associated remodeling.

Aims Macrophages (MΦ), known for immunological roles such as phagocytosis and antigen presentation, have been found to electrotonically couple to cardiomyocytes (CM) of the atrio-ventricular node via Cx43, affecting cardiac conduction in isolated mouse hearts. Here, we characterise passive and active electrophysiological properties of murine cardiac resident MΦ, and model their potential electrophysiological relevance for CM. Methods and Results We combined classic electrophysiological approaches with 3 D florescence imaging, RNA-sequencing, pharmacological interventions and computer simulations. We used Cx3cr1eYFP/+ mice wherein cardiac MΦ were fluorescently labelled. FACS-purified fluorescent MΦ from mouse hearts were studied by whole-cell patch-clamp. MΦ electrophysiological properties include: membrane resistance 2.2 ± 0.1 GΩ (all data mean±SEM), capacitance 18.3 ± 0.1 pF, resting membrane potential -39.6 ± 0.3 mV, and several voltage-activated, outward or inwardly-rectifying potassium currents. Using ion channel blockers (barium, TEA, 4-AP, margatoxin, XEN-D0103, DIDS), flow cytometry, immuno-staining and RNA-sequencing, we identified Kv1.3, Kv1.5 and Kir2.1 as channels contributing to observed ion currents. MΦ displayed four patterns for outward and two for inward-rectifier potassium currents. Additionally, MΦ showed surface expression of Cx43, a prerequisite for homo- and/or heterotypic electrotonic coupling. Experimental results fed into development of an original computational model to describe cardiac MΦ electrophysiology. Computer simulations to quantitatively assess plausible effects of MΦ on electrotonically coupled CM showed that MΦ can depolarise resting CM, shorten early and prolong late action potential duration, with effects depending on coupling strength and individual MΦ electrophysiological properties, in particular resting membrane potential and presence/absence of Kir2.1. Conclusions Our results provide a first electrophysiological characterisation of cardiac resident MΦ, and a computational model to quantitatively explore their relevance in the heterocellular heart. Future work will be focussed at distinguishing electrophysiological effects of MΦ–CM coupling on both cell types during steady-state and in patho-physiological remodelling, when immune cells change their phenotype, proliferate, and/or invade from external sources. Translational Perspective Cardiac tissue contains resident macrophages (MΦ) which, beyond immunological and housekeeping roles, have been found to electrotonically couple via connexins to cardiomyocytes (CM), stabilising atrio-ventricular conduction at high excitation rates. Here, we characterise structure and electrophysiological function of murine cardiac MΦ and provide a computational model to quantitatively probe the potential relevance of MΦ-CM coupling for cardiac electrophysiology. We find that MΦ are unlikely to have major electrophysiological effects in normal tissue, where they would hasten early and slow late CM-repolarisation. Further work will address potential arrhythmogenicity of MΦ in patho-physiologically remodelled tissue containing elevated MΦ-numbers, incl. non-resident recruited cells.