Biophysical models are increasingly used to gain mechanistic insights by fitting and reproducing experimental and clinical data. The inherent variability in the recorded datasets, however, presents a key challenge. In this study, we present a novel approach, which integrates mechanistic modeling and machine learning to analyze in vitro cardiac mechanics data and solve the inverse problem of model parameter inference. We designed a novel generative adversarial network (GAN) and employed it to construct virtual populations of cardiac ventricular myocyte models in order to study the action of Omecamtiv Mecarbil (OM), a positive cardiac inotrope. Populations of models were calibrated from mechanically unloaded myocyte shortening recordings obtained in experiments on rat myocytes in the presence and absence of OM. The GAN was able to infer model parameters while incorporating prior information about which model parameters OM targets. The generated populations of models reproduced variations in myocyte contraction recorded during in vitro experiments and provided improved understanding of OM’s mechanism of action. Inverse mapping of the experimental data using our approach suggests a novel action of OM, whereby it modifies interactions between myosin and tropomyosin proteins. To validate our approach, the inferred model parameters were used to replicate other in vitro experimental protocols, such as skinned preparations demonstrating an increase in calcium sensitivity and a decrease in the Hill coefficient of the force–calcium (F–Ca) curve under OM action. Our approach thereby facilitated the identification of the mechanistic underpinnings of experimental observations and the exploration of different hypotheses regarding variability in this complex biological system.

Single cardiomyocytes are widely used for investigations of the cellular and molecular mechanisms of regulation and modulation of cardiac performance. Intact cardiomyocytes allow one to study in detail cell function avoiding the effects of extracellular matrix and neighboring cells. The most established protocols of cardiomyocyte isolation are based on the isolated heart perfusion using a Langendorff-apparatus or on intraventricular perfusion using a syringe. However, the yield of single cardiomyocytes obtained by these methods may be low due to the cell injury following non-uniform enzyme digestion of connective tissue in different heart chambers. Moreover, isolation of atrial cardiomyocytes is challenging because of their small size and complex geometric shape. Here we present a new protocol for simultaneous isolation of high quality cardiomyocytes from the atria, ventricular free walls and interventricular septum. The protocol is based on the combination of the Langendorff perfusion method with the intraventricular and intra-atrial injection technique taking into account the collagen content variation between the different heart chambers. Obtained cells demonstrate rod-shaped morphology, a clear and regular sarcomere striation pattern and rat-specific frequency-dependence of contraction and calcium transient parameters. Our protocol provides gentle cell isolation that increases the yield of single cardiomyocytes suitable for biophysical researches .

Acetylcholine (ACh) is the neurotransmitter of the parasympathetic nervous system that modulates cardiac function, and its high concentrations may induce atrial fibrillation. We compared the ACh action on the mechanical function of single cardiomyocytes from the left atria (LA) and the right atria (RA). We exposed single rat LA and RA cardiomyocytes to 1, 10, and 100 µM ACh for 10–15 min and measured the parameters of sarcomere shortening–relengthening and cytosolic calcium ([Ca²⁺]_i) transients during cell contractions. We also studied the effects of ACh on cardiac myosin function using an in vitro motility assay and analyzed the phosphorylation level of sarcomeric proteins. In LA cardiomyocytes, ACh decreased the time to peak sarcomere shortening, time to 50% relengthening, and time to peak [Ca²⁺]_i transients. In RA cardiomyocytes, ACh affected the time of shortening and relengthening only at 10 µM. In the in vitro motility assay, ACh reduced to a greater extent the sliding velocity of F-actin over myosin from LA cardiomyocytes, which was accompanied by a more pronounced decrease in phosphorylation of the myosin regulatory light chain (RLC) in LA cardiomyocytes than in RA cardiomyocytes. Our findings indicate that ACh plays an important role in modulating the contractile function of LA and RA, provoking more pronounced changes in the time course of sarcomere shortening–relengthening and the kinetics of actin–myosin interaction in LA cardiomyocytes.

17β-estradiol (E2) is the most active metabolite of estrogen with a wide range of physiological action on cardiac muscle. Previous studies have reported E2 effects predominantly for the ventricles, while the E2 impact on the atria has been less examined. In this study, we focused on the direct E2 effects on atrial and ventricular contractility at the cellular and molecular levels. Single atrial and ventricular cardiomyocytes (CM) from adult (24 weeks-old) female Wistar rats were incubated with 10 nM E2 for 15 min. Sarcomere length and cytosolic [Ca²⁺]_i transients were measured in mechanically non-loaded CM, and the tension–length relationship was studied in CM mechanically loaded by carbon fibers. The actin–myosin interaction and sarcomeric protein phosphorylation were analyzed using an in vitro motility assay and gel electrophoresis with Pro-Q Diamond phosphoprotein stain. E2 had chamber-specific effects on the contractile function of CM with a pronounced influence on ventricular CM. The characteristics of [Ca²⁺]_i transients did not change in both atrial and ventricular CM. However, in ventricular CM, E2 reduced the amplitude and maximum velocity of sarcomere shortening and decreased the slope of the passive tension–length relationship that was associated with increased TnI and cMyBP-C phosphorylation. E2 treatment accelerated the cross-bridge cycle of both atrial and ventricular myosin that was associated with increased phosphorylation of the myosin essential light chain. This study shows that E2 impairs the mechanical function of the ventricular myocardium while atrial contractility remains mostly preserved. Hormonal replacement therapy (HRT) with estrogen is by far the most effective therapy for treating climacteric symptoms experienced during menopause. Here we found a chamber specificity of myocardial contractile function to E2 that should be taken into account for the potential side effects of HRT.

Atrial fibrillation (AF) is the most common cardiaс arrhythmia accounting for significant morbidity and mortality. Cardiac alternans, defined at the cellular level as cyclic beat-to-beat alternations in cell contraction, action potential morphology, and intracellular Ca²⁺ concentration, have been associated with sustained AF. Although alternans of action potential and intracellular Ca²⁺ transients in atrial cardiomyocytes in AF have been elucidated, little is known about the mechanical alternans, cyclic alternations in contraction force, and cardiomyocyte shortening. Moreover, the susceptibility of different regions of atria to alternans formation is unclear. The aim of this study is to analyze and compare the parameters of alternans in sarcomere shortening-relengthening in single cardiomyocytes isolated from the left and right atria (LA, RA) of rats with acetylcholine-CaCl₂ induced AF. We found that autonomic AF provoked mechanical alternans in single LA and RA cardiomyocytes. We also showed that the degree and magnitude of mechanical alternans were greater in LA cardiomyocytes than in RA cardiomyocytes. Our results propose that LA cardiomyocytes have a greater susceptibility to AF-induced beat-to-beat alternans in sarcomere shortening amplitudes than RA myocytes that may cause LA mechanical dysfunction.