The Virtual Child.

SUMMARY: We are building the world's first Virtual Child-a computer model of normal and cancerous human development at the level of each individual cell. The Virtual Child will "develop cancer" that we will subject to unlimited virtual clinical trials that pinpoint, predict, and prioritize potential new treatments, bringing forward the day when no child dies of cancer, giving each one the opportunity to lead a full and healthy life.

Two-Stage Self-Supervised Contrastive Learning Aided Transformer for Real-Time Medical Image Segmentation.

The availability of large, high-quality annotated datasets in the medical domain poses a substantial challenge in segmentation tasks. To mitigate the reliance on annotated training data, self-supervised pre-training strategies have emerged, particularly employing contrastive learning methods on dense pixel-level representations. In this work, we proposed to capitalize on intrinsic anatomical similarities within medical image data and develop a semantic segmentation framework through a self-supervised fusion network, where the availability of annotated volumes is limited. In a unified training phase, we combine segmentation loss with contrastive loss, enhancing the distinction between significant anatomical regions that adhere to the available annotations. To further improve the segmentation performance, we introduce an efficient parallel transformer module that leverages Multiview multiscale feature fusion and depth-wise features. The proposed transformer architecture, based on multiple encoders, is trained in a self-supervised manner using contrastive loss. Initially, the transformer is trained using an unlabeled dataset. We then fine-tune one encoder using data from the first stage and another encoder using a small set of annotated segmentation masks. These encoder features are subsequently concatenated for the purpose of brain tumor segmentation. The multiencoder-based transformer model yields significantly better outcomes across three medical image segmentation tasks. We validated our proposed solution by fusing images across diverse medical image segmentation challenge datasets, demonstrating its efficacy by outperforming state-of-the-art methodologies.

The effect of ionizing radiation through cardiac stereotactic body radiation therapy on myocardial tissue for refractory ventricular arrhythmias: A review.

Cardiac stereotactic body radiation therapy (cSBRT) is a non-invasive treatment modality that has been recently reported as an effective treatment for ventricular arrhythmias refractory to medical therapy and catheter ablation. The approach leverages tools developed and refined in radiation oncology, where experience has been accumulated in the treatment of a wide variety of malignant conditions. However, important differences exist between rapidly dividing malignant tumor cells and fully differentiated myocytes in pathologically remodeled ventricular myocardium, which represent the respective radiation targets. Despite its initial success, little is known about the radiobiology of the anti-arrhythmic effect cSBRT. Pre-clinical data indicates a late fibrotic effect of that appears between 3 and 4 months following cSBRT, which may result in conduction slowing and block. However, there is clear clinical evidence of an anti-arrhythmic effect of cSBRT that precedes the appearance of radiation induced fibrosis for which the mechanism is unclear. In addition, the data to date suggests that even the late anti-arrhythmic effect of cSBRT is not fully attributable to radiation.-induced fibrosis. Pre-clinical data has identified upregulation of proteins expected to result in both increased cell-to-cell coupling and excitability in the early post cSBRT period and demonstrated an associated increase in myocardial conduction velocity. These observations indicate a complex response to radiotherapy and highlight the lack of clarity regarding the different stages of the anti-arrhythmic mechanism of cSBRT. It may be speculated that in the future cSBRT therapy could be planned to deliver both early and late radiation effects titrated to optimize the combined anti-arrhythmic efficacy of the treatment. In addition to these outstanding mechanistic questions, the optimal patient selection, radiation modality, radiation dose and treatment planning strategy are currently being investigated. In this review, we consider the structural and functional effect of radiation on myocardium and the possible anti-arrhythmic mechanisms of cSBRT. Review of the published data highlights the exciting prospects for the development of knowledge and understanding in this area in which so many outstanding questions exist.

Detection of focal source and arrhythmogenic substrate from body surface potentials to guide atrial fibrillation ablation.

Focal sources (FS) are believed to be important triggers and a perpetuation mechanism for paroxysmal atrial fibrillation (AF). Detecting FS and determining AF sustainability in atrial tissue can help guide ablation targeting. We hypothesized that sustained rotors during FS-driven episodes indicate an arrhythmogenic substrate for sustained AF, and that non-invasive electrical recordings, like electrocardiograms (ECGs) or body surface potential maps (BSPMs), could be used to detect FS and AF sustainability. Computer simulations were performed on five bi-atrial geometries. FS were induced by pacing at cycle lengths of 120-270 ms from 32 atrial sites and four pulmonary veins. Self-sustained reentrant activities were also initiated around the same 32 atrial sites with inexcitable cores of radii of 0, 0.5 and 1 cm. FS fired for two seconds and then AF inducibility was tested by whether activation was sustained for another second. ECGs and BSPMs were simulated. Equivalent atrial sources were extracted using second-order blind source separation, and their cycle length, periodicity and contribution, were used as features for random forest classifiers. Longer rotor duration during FS-driven episodes indicates higher AF inducibility (area under ROC curve = 0.83). Our method had accuracy of 90.6±1.0% and 90.6±0.6% in detecting FS presence, and 93.1±0.6% and 94.2±1.2% in identifying AF sustainability, and 80.0±6.6% and 61.0±5.2% in determining the atrium of the focal site, from BSPMs and ECGs of five atria. The detection of FS presence and AF sustainability were insensitive to vest placement (±9.6%). On pre-operative BSPMs of 52 paroxysmal AF patients, patients classified with initiator-type FS on a single atrium resulted in improved two-to-three-year AF-free likelihoods (p-value < 0.01, logrank tests). Detection of FS and arrhythmogenic substrate can be performed from ECGs and BSPMs, enabling non-invasive mapping towards mechanism-targeted AF treatment, and malignant ectopic beat detection with likely AF progression.

Leadless left ventricular endocardial pacing for CRT upgrades in previously failed and high-risk patients in comparison with coronary sinus CRT upgrades.

AIMS: Cardiac resynchronization therapy (CRT) upgrades may be less likely to improve following intervention. Leadless left ventricular (LV) endocardial pacing has been used for patients with previously failed CRT or high-risk upgrades. We compared procedural and long-term outcomes in patients undergoing coronary sinus (CS) CRT upgrades with high-risk and previously failed CRT upgrades undergoing LV endocardial upgrades. METHOD AND RESULTS: Prospective consecutive CS upgrades between 2015 and 2019 were compared with those undergoing WiSE-CRT implantation. Cardiac resynchronization therapy response at 6 months was defined as improvement in clinical composite score (CCS) and a reduction in LV end-systolic volume (LVESV) ≥15%. A total of 225 patients were analysed; 121 CS and 104 endocardial upgrades. Patients receiving WiSE-CRT tended to have more comorbidities and were more likely to have previous cardiac surgery (30.9% vs. 16.5%; P = 0.012), hypertension (59.2% vs. 34.7%; P < 0.001), chronic obstructive airways disease (19.4% vs. 9.9%; P = 0.046), and chronic kidney disease (46.4% vs. 21.5%; P < 0.01) but similar LV ejection fraction (30.0 ± 8.3% vs. 29.5 ± 8.6%; P = 0.678). WiSE-CRT upgrades were successful in 97.1% with procedure-related mortality in 1.9%. Coronary sinus upgrades were successful in 97.5% of cases with a 2.5% rate of CS dissection and 5.6% lead malfunction/displacement. At 6 months, 91 WiSE-CRT upgrades and 107 CS upgrades had similar improvements in CCS (76.3% vs. 68.5%; P = 0.210) and reduction in LVESV ≥15% (54.2% vs. 56.3%; P = 0.835). CONCLUSION: Despite prior failed upgrades and high-risk patients with more comorbidities, WiSE-CRT upgrades had high rates of procedural success and similar improvements in CCS and LV remodelling with CS upgrades.

Clinical effectiveness of a dedicated cardiac resynchronization therapy pre-assessment clinic incorporating cardiac magnetic resonance imaging and cardiopulmonary exercise testing on patient selection and outcomes.

BACKGROUND: Pre-procedural assessment of patients undergoing cardiac resynchronization therapy (CRT) is heterogenous and patients implanted with unfavorable characteristics may account for non-response. A dedicated CRT pre-assessment clinic (CRT PAC) was developed to standardize the review process and undertake structured pre-procedural evaluation. The aim of this analysis was to determine the effectiveness on patient selection and outcomes. METHODS: A prospective database of consecutive patients attending the CRT PAC between 2013 and 2018 was analyzed. Pre-operative assessment included cardiac magnetic resonance (CMR) and cardiopulmonary exercise testing (CPET). Patients were considered CRT responders based on improvement in clinical composite score (CCS) and/or reduction in left ventricular end-systolic volume (LVESV) ≥ 15% at 6-months follow-up. RESULTS: Of 252 patients reviewed in the CRT PAC during the analysis period, 192 fulfilled consensus guidelines for implantation. Of the patients receiving CRT, 82% showed improvement in their CCS and 57% had a reduction in LVESV ≥ 15%. The presence of subendocardial scar on CMR and a peak VO2 ≤ 12 ml/kg/min on CPET predicted CRT non-response. Two patients were unsuitable for CRT as they had end-stage heart failure and died during follow-up. The majority of patients initially deemed unsuitable for CRT did not suffer from unexpected hospitalization for decompensated heart failure or died from cardiovascular disease; only 8 patients (13%) received CRT devices during follow-up because of symptomatic left ventricular systolic impairment. CONCLUSION: A dedicated CRT PAC is able to appropriately select patients for CRT. Pre-procedural investigation/imaging can identify patients unlikely to respond to, or may not yet be suitable for CRT.

Standardised computed tomographic assessment of left atrial morphology and tissue thickness in humans.

AIMS: Left atrial (LA) remodelling is a common feature of many cardiovascular pathologies and is a sensitive marker of adverse cardiovascular outcomes. The aim of this study was to establish normal ranges for LA parameters derived from coronary computed tomographic angiography (CCTA) imaging using a standardised image processing pipeline to establish normal ranges in a previously described cohort. METHODS: CCTA imaging from 193 subjects recruited to the Budapest GLOBAL twin study was analysed. Indexed LA cavity volume (LACVi), LA surface area (LASAi), wall thickness and LA tissue volume (LATVi) were calculated. Wall thickness maps were combined into an atlas. Indexed LA parameters were compared with clinical variables to identify early markers of pathological remodelling. RESULTS: LACVi is similar between sexes (31 ml/m2 v 30 ml/m2) and increased in hypertension (33 ml/m2 v 29 ml/m2, p = 0.009). LASAi is greater in females than males (47.8 ml/m2 v 45.8 ml/m2 male, p = 0.031). Median LAWT was 1.45 mm. LAWT was lowest at the inferior portion of the posterior LA wall (1.14 mm) and greatest in the septum (median = 2.0 mm) (p < 0.001). Conditions known to predispose to the development of AF were not associated with differences in tissue thickness. CONCLUSIONS: The reported LACVi, LASAi, LATVi and tissue thickness derived from CCTA may serve as reference values for this age group and clinical characteristics for future studies. Increased LASAi in females in the absence of differences in LACVi or LATVi may indicate differential LA shape changes between the sexes. AF predisposing conditions, other than sex, were not associated with detectable changes in LAWT.Clinical trial registration:http://www.ClinicalTrials.gov/NCT01738828.

Hypokalemia Promotes Arrhythmia by Distinct Mechanisms in Atrial and Ventricular Myocytes.

RATIONALE: Hypokalemia occurs in up to 20% of hospitalized patients and is associated with increased incidence of ventricular and atrial fibrillation. It is unclear whether these differing types of arrhythmia result from direct and perhaps distinct effects of hypokalemia on cardiomyocytes. OBJECTIVE: To investigate proarrhythmic mechanisms of hypokalemia in ventricular and atrial myocytes. METHODS AND RESULTS: Experiments were performed in isolated rat myocytes exposed to simulated hypokalemia conditions (reduction of extracellular [K+] from 5.0 to 2.7 mmol/L) and supported by mathematical modeling studies. Ventricular cells subjected to hypokalemia exhibited Ca2+ overload and increased generation of both spontaneous Ca2+ waves and delayed afterdepolarizations. However, similar Ca2+-dependent spontaneous activity during hypokalemia was only observed in a minority of atrial cells that were observed to contain t-tubules. This effect was attributed to close functional pairing of the Na+-K+ ATPase and Na+-Ca2+ exchanger proteins within these structures, as reduction in Na+ pump activity locally inhibited Ca2+ extrusion. Ventricular myocytes and tubulated atrial myocytes additionally exhibited early afterdepolarizations during hypokalemia, associated with Ca2+ overload. However, early afterdepolarizations also occurred in untubulated atrial cells, despite Ca2+ quiescence. These phase-3 early afterdepolarizations were rather linked to reactivation of nonequilibrium Na+ current, as they were rapidly blocked by tetrodotoxin. Na+ current-driven early afterdepolarizations in untubulated atrial cells were enabled by membrane hyperpolarization during hypokalemia and short action potential configurations. Brief action potentials were in turn maintained by ultra-rapid K+ current (IKur); a current which was found to be absent in tubulated atrial myocytes and ventricular myocytes. CONCLUSIONS: Distinct mechanisms underlie hypokalemia-induced arrhythmia in the ventricle and atrium but also vary between atrial myocytes depending on subcellular structure and electrophysiology.

The impact of wall thickness and curvature on wall stress in patient-specific electromechanical models of the left atrium.

The left atrium (LA) has a complex anatomy with heterogeneous wall thickness and curvature. The anatomy plays an important role in determining local wall stress; however, the relative contribution of wall thickness and curvature in determining wall stress in the LA is unknown. We have developed electromechanical finite element (FE) models of the LA using patient-specific anatomical FE meshes with rule-based myofiber directions. The models of the LA were passively inflated to 10mmHg followed by simulation of the contraction phase of the atrial cardiac cycle. The FE models predicted maximum LA volumes of 156.5 mL, 99.3 mL and 83.4 mL and ejection fractions of 36.9%, 32.0% and 25.2%. The median wall thickness in the 3 cases was calculated as [Formula: see text] mm, [Formula: see text] mm, and [Formula: see text] mm. The median curvature was determined as [Formula: see text] [Formula: see text], [Formula: see text], and [Formula: see text]. Following passive inflation, the correlation of wall stress with the inverse of wall thickness and curvature was 0.55-0.62 and 0.20-0.25, respectively. At peak contraction, the correlation of wall stress with the inverse of wall thickness and curvature was 0.38-0.44 and 0.16-0.34, respectively. In the LA, the 1st principal Cauchy stress is more dependent on wall thickness than curvature during passive inflation and both correlations decrease during active contraction. This emphasizes the importance of including the heterogeneous wall thickness in electromechanical FE simulations of the LA. Overall, simulation results and sensitivity analyses show that in complex atrial anatomy it is unlikely that a simple anatomical-based law can be used to estimate local wall stress, demonstrating the importance of FE analyses.

Evaluation of a real-time magnetic resonance imaging-guided electrophysiology system for structural and electrophysiological ventricular tachycardia substrate assessment.

AIMS: Potential advantages of real-time magnetic resonance imaging (MRI)-guided electrophysiology (MR-EP) include contemporaneous three-dimensional substrate assessment at the time of intervention, improved procedural guidance, and ablation lesion assessment. We evaluated a novel real-time MR-EP system to perform endocardial voltage mapping and assessment of delayed conduction in a porcine ischaemia-reperfusion model. METHODS AND RESULTS: Sites of low voltage and slow conduction identified using the system were registered and compared to regions of late gadolinium enhancement (LGE) on MRI. The Sorensen-Dice similarity coefficient (DSC) between LGE scar maps and voltage maps was computed on a nodal basis. A total of 445 electrograms were recorded in sinus rhythm (range: 30-186) using the MR-EP system including 138 electrograms from LGE regions. Pacing captured at 103 sites; 47 (45.6%) sites had a stimulus-to-QRS (S-QRS) delay of ≥40 ms. Using conventional (0.5-1.5 mV) bipolar voltage thresholds, the sensitivity and specificity of voltage mapping using the MR-EP system to identify MR-derived LGE was 57% and 96%, respectively. Voltage mapping had a better predictive ability in detecting LGE compared to S-QRS measurements using this system (area under curve: 0.907 vs. 0.840). Using an electrical threshold of 1.5 mV to define abnormal myocardium, the total DSC, scar DSC, and normal myocardium DSC between voltage maps and LGE scar maps was 79.0 ± 6.0%, 35.0 ± 10.1%, and 90.4 ± 8.6%, respectively. CONCLUSION: Low-voltage zones and regions of delayed conduction determined using a real-time MR-EP system are moderately associated with LGE areas identified on MRI.

Changes in contractility determine coronary haemodynamics in dyssynchronous left ventricular heart failure, not vice versa.

BACKGROUND: Biventricular pacing has been shown to increase both cardiac contractility and coronary flow acutely but the causal relationship is unclear. We hypothesised that changes in coronary flow are secondary to changes in cardiac contractility. We sought to examine this relationship by modulating coronary flow and cardiac contractility. METHODS: Contractility and lusitropy were altered by varying the location of pacing in 8 patients. Coronary autoregulation was transiently disabled with intracoronary adenosine. Simultaneous coronary flow velocity, coronary pressure and left ventricular pressure data were measured in the different pacing settings with and without hyperaemia and wave intensity analysis performed. RESULTS: Multisite pacing was effective at altering left ventricular contractility and lusitropy (pos. dp/dtmax -13% to +10% and neg. dp/dtmax -15% to +17% compared to baseline). Intracoronary adenosine decreased microvascular resistance (362.5 mm Hg/s/m to 156.7 mm Hg/s/m, p < 0.001) and increased LAD flow velocity (22 cm/s vs 45 cm/s, p < 0.001) but did not acutely change contractility or lusitropy. The magnitude of the dominant accelerating wave, the Backward Expansion Wave, was proportional to the degree of contractility as well as lusitropy (r = 0.47, p < 0.01 and r = -0.50, p < 0.01). Perfusion efficiency (the proportion of accelerating waves) increased at hyperaemia (76% rest vs 81% hyperaemia, p = 0.04). Perfusion efficiency correlated with contractility and lusitropy at rest (r = 0.43 & -0.50 respectively, p = 0.01) and hyperaemia (r = 0.59 & -0.6, p < 0.01). CONCLUSIONS: Acutely increasing coronary flow with adenosine in patients with systolic heart failure does not increase contractility. Changes in coronary flow with biventricular pacing are likely to be a consequence of enhanced cardiac contractility from resynchronization and not vice versa.

Bringing in vitro analysis closer to in vivo: Studying doxorubicin toxicity and associated mechanisms in 3D human microtissues with PBPK-based dose modelling.

Doxorubicin (DOX) is a chemotherapeutic agent of which the medical use is limited due to cardiotoxicity. While acute cardiotoxicity is reversible, chronic cardiotoxicity is persistent or progressive, dose-dependent and irreversible. While DOX mechanisms of action are not fully understood yet, 3 toxicity processes are known to occur in vivo: cardiomyocyte dysfunction, mitochondrial dysfunction and cell death. We present an in vitro experimental design aimed at detecting DOX-induced cardiotoxicity by obtaining a global view of the induced molecular mechanisms through RNA-sequencing. To better reflect the in vivo situation, human 3D cardiac microtissues were exposed to physiologically-based pharmacokinetic (PBPK) relevant doses of DOX for 2 weeks. We analysed a therapeutic and a toxic dosing profile. Transcriptomics analysis revealed significant gene expression changes in pathways related to "striated muscle contraction" and "respiratory electron transport", thus suggesting mitochondrial dysfunction as an underlying mechanism for cardiotoxicity. Furthermore, expression changes in mitochondrial processes differed significantly between the doses. Therapeutic dose reflects processes resembling the phenotype of delayed chronic cardiotoxicity, while toxic doses resembled acute cardiotoxicity. Overall, these results demonstrate the capability of our innovative in vitro approach to detect the three known mechanisms of DOX leading to toxicity, thus suggesting its potential relevance for reflecting the patient situation. Our study also demonstrated the importance of applying physiologically relevant doses during toxicological research, since mechanisms of acute and chronic toxicity differ.

Transcatheter mitral valve replacement in mitral annulus calcification - "The art of computer simulation".

There is considerable interest in transcatheter prosthetic valve treatment for mitral valve disease in high-risk individuals. Although the presence of mitral annular calcium (MAC) may provide an anchoring zone for such devices, results to date have been modest with reported technical failure rates approaching 30% in specialist centers. This in part relates to the risk of left ventricular outflow tract obstruction and device dislodgment but also to the lack of specific imaging guidelines to plan for such procedures. We present the use of finite element analysis and computer simulation based on cardiac CT in three patients with severe MAC in whom transcatheter devices were considered. In the first two cases, the computer simulations were performed after the clinical procedure and were concordant with the clinical outcome. For the third case, computer simulation was performed prior to the clinical procedure. This indicated unsuitability for transcatheter device deployment and a subsequent medical management was adopted. Overall, our initial results suggest that computer simulation may have the potential to improve patient selection for transcatheter mitral valve replacement in the presence of significant MAC.

Normoxic cells remotely regulate the acid-base balance of cells at the hypoxic core of connexin-coupled tumor growths.

ATP fuels the removal of metabolic end-products, including H+ ions that profoundly modulate biological activities. Energetic resources in hypoxic tumor regions are constrained by low-yielding glycolysis, and any means of reducing the cost of acid extrusion, without compromising pH homeostasis, would therefore be advantageous for cancer cells. Some cancers express connexin channels that allow solute exchange between cells, and we propose that, via this route, normoxic cells supply hypoxic neighbors with acid-neutralizing HCO3- ions. This hypothesis was tested by imaging cytoplasmic pH in spheroidal tissue growths of connexin43-positive pancreatic cancer Colo357 cells during light-controlled H+ uncaging at the hypoxic core. Cytoplasmic acid retention at the core was halved in the presence of CO2/HCO3-, but this process requires a restorative HCO3- flux. The effect of CO2/HCO3- was ablated by connexin43 inhibition or knockdown. In connexin-decoupled spheroids, 4,4'-diisothiocyano-2,2'-stilbenedisulfonic acid (DIDS), an inhibitor of HCO3- uptake, had no effect on cytoplasmic [H+] in the H+-uncaging region, indicating that DIDS-sensitive transport is not an adequate pH-regulatory strategy therein. With intact connexin-coupling, acid retention at the core increased upon DIDS treatment, indicating that HCO3- ions are taken up actively by peripheral cells and then transmitted passively to cells at the hypoxic core. Thus, the energetic burden of pH regulation is offloaded from hypoxic cells onto metabolically altruistic normoxic neighbors.-Dovmark, T. H., Hulikova, A., Niederer, S. A., Vaughan-Jones, R. D., Swietach, P. Normoxic cells remotely regulate the acid-base balance of cells at the hypoxic core of connexin-coupled tumor growths.

A Biophysical Systems Approach to Identifying the Pathways of Acute and Chronic Doxorubicin Mitochondrial Cardiotoxicity.

The clinical use of the anthracycline doxorubicin is limited by its cardiotoxicity which is associated with mitochondrial dysfunction. Redox cycling, mitochondrial DNA damage and electron transport chain inhibition have been identified as potential mechanisms of toxicity. However, the relative roles of each of these proposed mechanisms are still not fully understood. The purpose of this study is to identify which of these pathways independently or in combination are responsible for doxorubicin toxicity. A state of the art mathematical model of the mitochondria including the citric acid cycle, electron transport chain and ROS production and scavenging systems was extended by incorporating a novel representation for mitochondrial DNA damage and repair. In silico experiments were performed to quantify the contributions of each of the toxicity mechanisms to mitochondrial dysfunction during the acute and chronic stages of toxicity. Simulations predict that redox cycling has a minor role in doxorubicin cardiotoxicity. Electron transport chain inhibition is the main pathway for acute toxicity for supratherapeutic doses, being lethal at mitochondrial concentrations higher than 200µM. Direct mitochondrial DNA damage is the principal pathway of chronic cardiotoxicity for therapeutic doses, leading to a progressive and irreversible long term mitochondrial dysfunction.

Species-dependent adaptation of the cardiac Na+/K+ pump kinetics to the intracellular Na+ concentration.

The Na(+)/K(+) ATPase (NKA) plays a critical role in maintaining ionic homeostasis and dynamic function in cardiac myocytes, within both the in vivo cell and in silico models. Physiological conditions differ significantly between mammalian species. However, most existing formulations of NKA used to simulate cardiac function in computational models are derived from a broad range of experimental sources spanning many animal species. The resultant inability of these models to discern species-specific features is a significant obstacle to achieving a detailed quantitative and comparative understanding of physiological behaviour in different biological contexts. Here we present a framework for characterising the steady-state NKA current using a biophysical mechanistic model specifically designed to provide a mechanistic explanation of the NKA flux supported by self-consistent species-specific data. We thus compared NKA kinetics specific to guinea- pig and rat ventricular myocytes. We observe that the apparent binding affinity for sodium in the rat is significantly lower, whereas the overall pump cycle rate is doubled, in comparison to the guinea pig. This sensitivity of NKA to its regulatory substrates compensates for the differences in Na(+) concentrations between the cell types. NKA is thereby maintained within its dynamic range over a wide range of pacing frequencies in these two species, despite significant disparities in sodium concentration. Hence, by replacing a conventional generic NKA model with our rat-specific NKA formula into a whole-cell simulation, we have, for the first time, been able to accurately reproduce the action potential duration and the steady-state sodium concentration as functions of pacing frequency.