Modeling transmural heterogeneity of K(ATP) current in rabbit ventricular myocytes

To investigate the mechanisms regulating excitation-metabolic coupling in rabbit epicardial, midmyocardial, and endocardial ventricular myocytes we extended the LabHEART model (Puglisi JL and Bers DM. Am J Physiol Cell Physiol 281: C2049–C2060, 2001). We incorporated equations for Ca²⁺ and Mg²⁺ buffering by ATP and ADP, equations for nucleotide regulation of ATP-sensitive K⁺ channel and L-type Ca²⁺ channel, Na⁺-K⁺-ATPase, and sarcolemmal and sarcoplasmic Ca²⁺-ATPases, and equations describing the basic pathways (creatine and adenylate kinase reactions) known to communicate the flux changes generated by intracellular ATPases. Under normal conditions and during 20 min of ischemia, the three regions were characterized by different I_Na, I_to, I_Kr, I_Ks, and I_Kp channel properties. The results indicate that the ATP-sensitive K⁺ channel is activated by the smallest reduction in ATP in epicardial cells and largest in endocardial cells when cytosolic ADP, AMP, PCr, Cr, P_i, total Mg²⁺, Na⁺, K⁺, Ca²⁺, and pH diastolic levels are normal. The model predicts that only K_ATP ionophore (Kir6.2 subunit) and not the regulatory subunit (SUR2A) might differ from endocardium to epicardium. The analysis suggests that during ischemia, the inhomogeneous accumulation of the metabolites in the tissue sublayers may alter in a very irregular manner the K_ATP channel opening through metabolic interactions with the endogenous PI cascade (PIP₂, PIP) that in turn may cause differential action potential shortening among the ventricular myocyte subtypes. The model predictions are in qualitative agreement with experimental data measured under normal and ischemic conditions in rabbit ventricular myocytes. ATP-sensitive K⁺ channel; creatine and adenylate kinase reactions; phosphatidylinositol phosphates; heart; mathematical model