Isolation of single atrial and ventricular cells from the human heart.

The single isolated heart cell has recently emerged as a model for the study of the structure and function of cardiac cells. Heart muscle cells of adult animals of various species have been successfully isolated by enzymatic digestion of intact cardiac tissue. In this paper a dissociation method that yields living cells from atrial and ventricular tissue of young and adult humans is detailed. The cells retain the morphologic features of cells in intact cardiac tissue, and they generate action potentials and contractions in response to electrical stimulation. The study of isolated human heart cells should make a valuable contribution to knowledge of the normal and diseased heart.     

Dynamic myocardial contractile parameters from left ventricular pressure-volume measurements

A new dynamic model of left ventricular (LV) pressure-volume relationships in beating heart was developed by mathematically linking chamber pressure-volume dynamics with cardiac muscle force-length dynamics. The dynamic LV model accounted for >80% of the measured variation in pressure caused by small-amplitude volume perturbation in an otherwise isovolumically beating, isolated rat heart. The dynamic LV model produced good fits to pressure responses to volume perturbations, but there existed some systematic features in the residual errors of the fits. The issue was whether these residual errors would be damaging to an application where the dynamic LV model was used with LV pressure and volume measurements to estimate myocardial contractile parameters. Good agreement among myocardial parameters responsible for response magnitude was found between those derived by geometric transformations of parameters of the dynamic LV model estimated in beating heart and those found by direct measurement in constantly activated, isolated muscle fibers. Good agreement was also found among myocardial kinetic parameters estimated in each of the two preparations. Thus the small systematic residual errors from fitting the LV model to the dynamic pressure-volume measurements do not interfere with use of the dynamic LV model to estimate contractile parameters of myocardium. Dynamic contractile behavior of cardiac muscle can now be obtained from a beating heart by judicious application of the dynamic LV model to information-rich pressure and volume signals. This provides for the first time a bridge between the dynamics of cardiac muscle function and the dynamics of heart function and allows a beating heart to be used in studies where the relevance of myofilament contractile behavior to cardiovascular system function may be investigated.     

Simple experiments to understand the ionic origins and characteristics of the ventricular cardiac action potential

Electrophysiological experiments are helpful for students to understand the role of electrical activity in heart function. Papillary muscle, which belongs to the ventricle, offers the advantage of being easily studied using glass microelectrodes. In addition, there is commercially available software that simulates ventricular electrical activity and can help overcome some difficulties, such as voltage clamp experiments, which need expensive apparatus when used for studies on living preparations. Here, we present a class practical session that is taken by undergraduate students at our University. In the first part of this class, students record action potentials from papillary muscles with the use of glass microelectrodes, and they change extracellular conditions to study the ionic basis of the action potential. In the second part of the class, students simulate action potentials using the Oxsoft Heart model (v. 4.0) and model their previous experiments on papillary muscle to quantify the effects. In particular, the model is very helpful in promoting understanding of the effect that extracellular potassium has on cardiac action potential by simulating voltage clamp experiments. This twin approach of papillary muscle experiments and computer modeling leads to a good understanding of the functioning of the action potential and can help introduce discussion of some abnormal cardiac functioning.     

Purification and properties of hyaluronidase from human liver. Differences from and similarities to the testicular enzyme.

The effect of T3 (3,3',5-tri-iodothyronine) on protein turnover in skeletal and cardiac muscle was measured in intact rats by means of a 6 h [14C]tyrosine-infusion technique. Treatment with 25-30 micrograms of T3/100 g body wt. daily for 4-7 days increased the fractional rate of protein synthesis in skeletal muscle. Since the fractional growth rate of the muscle was decreased or unchanged, T3 treatment increased the rate of muscle protein breakdown. These findings suggest that increased protein degradation is an important factor in decreasing skeletal-muscle mass in hyperthyroidism. In contrast with skeletal muscle, T3 treatment for 7 days caused an equivalent increase in the rate of cardiac muscle growth and protein synthesis. This suggests that hyperthyroidism does not increase protein breakdown in heart muscle as it does in skeletal muscle. The failure of T3 to increase proteolysis in heart muscle may be due to a different action on the cardiac myocyte or to systemic effects of T3 which increase cardiac work.     

Fluid Dynamics of Heart Development

The morphology, muscle mechanics, fluid dynamics, conduction properties, and molecular biology of the developing embryonic heart have received much attention in recent years due to the importance of both fluid and elastic forces in shaping the heart as well as the striking relationship between the heart’s evolution and development. Although few studies have directly addressed the connection between fluid dynamics and heart development, a number of studies suggest that fluids may play a key role in morphogenic signaling. For example, fluid shear stress may trigger biochemical cascades within the endothelial cells of the developing heart that regulate chamber and valve morphogenesis. Myocardial activity generates forces on the intracardiac blood, creating pressure gradients across the cardiac wall. These pressures may also serve as epigenetic signals. In this article, the fluid dynamics of the early stages of heart development is reviewed. The relevant work in cardiac morphology, muscle mechanics, regulatory networks, and electrophysiology is also reviewed in the context of intracardial fluid dynamics.     

Theoretical approach to blood ejection from the human left ventricle

An ejection dynamics mathematical model of human left ventricle (LV) based on physiological data of human heart is proposed in this study. The mathematical equations were expressed in terms of vorticity-stream function equations in a prolate spheroidal coordinate system. These equations combined with specified boundary conditions were numerically solved by using an alternating-direction-implicit (ADI) algorithm with second order accuracy. The unsteady aspects of the ejection process were subsequently introduced into the numerical simulation. The numerical results have shown that the present ellipsoidal model could be available to simulate the ejection process of the human LV. Such a model combined with cardiac muscle mechanics could be studied further to determine altered left ventricular function in cardiac diseases.     

Models of cardiac excitation-contraction coupling in ventricular myocytes

Mathematical and computational modeling of cardiac excitation-contraction coupling has produced considerable insights into how the heart muscle contracts. With the increase in biophysical and physiological data available, the modeling has become more sophisticated with investigations spanning in scale from molecular components to whole cells. These modeling efforts have provided insight into cardiac excitation-contraction coupling that advanced and complemented experimental studies. One goal is to extend these detailed cellular models to model the whole heart. While this has been done with mechanical and electophysiological models, the complexity and fast time course of calcium dynamics have made inclusion of detailed calcium dynamics in whole heart models impractical. Novel methods such as the probability density approach and moment closure technique which increase computational efficiency might make this tractable.     

Regulation of Cardiac Muscle Contractility

The heart's physiological performance, unlike that of skeletal muscle, is regulated primarily by variations in the contractile force developed by the individual myocardial fibers. In an attempt to identify the basis for the characteristic properties of myocardial contraction, the individual cardiac contractile proteins and their behavior in contractile models in vitro have been examined. The low shortening velocity of heart muscle appears to reflect the weak ATPase activity of cardiac myosin, but this enzymatic activity probably does not determine active state intensity. Quantification of the effects of Ca⁺⁺ upon cardiac actomyosin supports the view that myocardial contractility can be modified by changes in the amount of calcium released during excitation-contraction coupling. Exchange of intracellular K⁺ with Na⁺ derived from the extracellular space also could enhance myocardial contractility directly, as highly purified cardiac actomyosin is stimulated when K⁺ is replaced by an equimolar amount of Na⁺. On the other hand, cardiac glycosides and catecholamines, agents which greatly increase the contractility of the intact heart, were found to be without significant actions upon highly purified reconstituted cardiac actomyosin.     

Modelling the mechanical properties of cardiac muscle

A model of passive and active cardiac muscle mechanics is presented, suitable for use in continuum mechanics models of the whole heart. The model is based on an extensive review of experimental data from a variety of preparations (intact trabeculae, skinned fibres and myofibrils) and species (mainly rat and ferret) at temperatures from 20 to 27°C. Experimental tests include isometric tension development, isotonic loading, quick-release/restretch, length step and sinusoidal perturbations. We show that all of these experiments can be interpreted with a four state variable model which includes (i) the passive elasticity of myocardial tissue, (ii) the rapid binding of Ca²⁺ to troponin C and its slower tension-dependent release, (iii) the kinetics of tropomyosin movement and availability of crossbridge binding sites and the length dependence of this process and (iv) the kinetics of crossbridge tension development under perturbations of myofilament length.     

Regulation of oxidative phosphorylation in intact mammalian heart in vivo

A dynamic computer model of oxidative phosphorylation in intact heart was developed by modifying the model of oxidative phosphorylation in intact skeletal muscle published previously. Next, this model was used for theoretical studies on the regulation of oxidative phosphorylation in intact heart in vivo during transition between different work intensities. It is shown that neither a direct activation of ATP usage alone nor a direct activation of both ATP usage and substrate dehydrogenation, including the calcium-activated tricarboxylate acid cycle dehydrogenases, can account for the constancy of [ADP], [PCr], [P(i)] and [NADH] during a significant increase in oxygen consumption and ATP turnover encountered in intact heart in vivo. Only a direct activation of all oxidative phosphorylation complexes in parallel with a stimulation of ATP usage and substrate dehydrogenation enabled to reproduce the experimental data concerning the constancy of metabolite concentrations. The molecular background of the differences between heart and skeletal muscle in the kinetic behaviour of the oxidative phosphorylation system is also discussed.     

4D subject-specific inverse modeling of the chick embryonic heart outflow tract hemodynamics

Blood flow plays a critical role in regulating embryonic cardiac growth and development, with altered flow leading to congenital heart disease. Progress in the field, however, is hindered by a lack of quantification of hemodynamic conditions in the developing heart. In this study, we present a methodology to quantify blood flow dynamics in the embryonic heart using subject-specific computational fluid dynamics (CFD) models. While the methodology is general, we focused on a model of the chick embryonic heart outflow tract (OFT), which distally connects the heart to the arterial system, and is the region of origin of many congenital cardiac defects. Using structural and Doppler velocity data collected from optical coherence tomography, we generated 4D (\(\hbox {3D}\,+\,\hbox {time}\)) embryo-specific CFD models of the heart OFT. To replicate the blood flow dynamics over time during the cardiac cycle, we developed an iterative inverse-method optimization algorithm, which determines the CFD model boundary conditions such that differences between computed velocities and measured velocities at one point within the OFT lumen are minimized. Results from our developed CFD model agree with previously measured hemodynamics in the OFT. Further, computed velocities and measured velocities differ by \(<\)15 % at locations that were not used in the optimization, validating the model. The presented methodology can be used in quantifications of embryonic cardiac hemodynamics under normal and altered blood flow conditions, enabling an in-depth quantitative study of how blood flow influences cardiac development.     

Isolation of the ryanodine receptor from cardiac sarcoplasmic reticulum and identity with the feet structures

Ryanodine, a highly toxic alkaloid, reacts specifically with the Ca2+ release channels which are localized in the terminal cisternae of sarcoplasmic reticulum (SR). In this study, the ryanodine receptor from cardiac SR has been purified, characterized, and compared with that of skeletal muscle SR. The ryanodine receptor was solubilized with 3-[(3-cholamidopropyl)dimethylammonio]-1-propanesulfonate (CHAPS) in the presence of phospholipids. Purification was performed by sequential affinity chromatography followed by gel permeation chromatography in the presence of CHAPS and phospholipids. The enrichment of the receptor from cardiac microsomes was about 110-fold. The purified receptor contained a major polypeptide band of Mr 340,000 with a minor band of Mr 300,000 (absorbance ratio 100/8) on sodium dodecyl sulfate-polyacrylamide gel electrophoresis. Electron microscopy of the purified receptor from heart showed square structures of 222 +/- 21 A/side, which is the unique characteristic of feet structures of junctional face membrane of terminal cisternae of SR. Recently, we isolated the ryanodine receptor from skeletal muscle (Inui, M., Saito, A., and Fleischer, S. (1987) J. Biol. Chem. 262, 1740-1747). The ryanodine receptors from heart and skeletal muscle have similar characteristics in terms of protein composition, morphology, chromatographic behavior, and Ca2+, salt, and phospholipid dependence of ryanodine binding. However, there are distinct differences: 1) the Mr of the receptor is slightly larger for skeletal muscle (Mr approximately 360,000); 2) the purified receptor from heart contains two different affinities for ryanodine binding with Kd values in the nanomolar and micromolar ranges, contrasting with that of skeletal muscle SR which shows only the high affinity binding; 3) the affinity of the purified cardiac receptor for ryanodine was 4-5-fold higher than that of skeletal muscle, measured under identical conditions. The greater sensitivity in ryanodine in intact heart can be directly explained by the tighter binding of the ryanodine receptor from heart. The present study suggests that basically similar machinery (the ryanodine receptor and foot structure) is involved in triggering Ca2+ release from cardiac and skeletal muscle SR, albeit there are distinct differences in the sensitivity to ryanodine and other ligands in heart versus skeletal muscle.     

Induction of myofibrillogenesis in cardiac lethal mutant axolotl hearts rescued by RNA derived from normal endoderm

A strain of axolotl, Ambystoma mexicanum, that carries the cardiac lethal or c gene presents an excellent model system in which to study inductive interactions during heart development. Embryos homozygous for gene c contain hearts that fail to beat and do not form sarcomeric myofibrils even though muscle proteins are present. Although they can survive for approximately three weeks, mutant embryos inevitably die due to lack of circulation. Embryonic axolotl hearts can be maintained easily in organ culture using only Holtfreter's solution as a culture medium. Mutant hearts can be induced to differentiate in vitro into functional cardiac muscle containing sarcomeric myofibrils by coculturing the mutant heart tube with anterior endoderm from a normal embryo. The induction of muscle differentiation can also be mediated through organ culture of mutant heart tubes in medium 'conditioned' by normal anterior endoderm. Ribonuclease was shown to abolish the ability of endoderm-conditioned medium to induce cardiac muscle differentiation. The addition of RNA extracted from normal early embryonic anterior endoderm to organ cultures of mutant hearts stimulated the differentiation of these tissues into contractile cardiac muscle containing well-organized sarcomeric myofibrils, while RNA extracted from early embryonic liver or neural tube did not induce either muscular contraction or myofibrillogenesis. Thus, RNA from anterior endoderm of normal embryos induces myofibrillogenesis and the development of contractile activity in mutant hearts, thereby correcting the genetic defect.     

Uniform sarcomere shortening behavior in isolated cardiac muscle cells

We have observed the dynamics of sarcomere shortening and the diffracting action of single, functionally intact, unattached cardiac muscle cells enzymatically isolated from the ventricular tissue of adult rats. Sarcomere length was measured either (a) continuously by a light diffraction method or (b) by direct inspection of the cell's striated image as recorded on videotape or by cinemicroscopy (120--400 frames/s). At physiological levels of added CaCl2 (0.5--2.0 mM), many cells were quiescent (i.e., they did not beat spontaneously) and contracted in response to electrical stimulation (less than or equal to 1.0-ms pulse width). Sarcomere length in the quiescent, unstimulated cells (1.93 +/- 0.10 [SD] micrometers), at peak shortening (1.57 +/- 0.13 micrometers, n = 49), and the maximum velocity of sarcomere shortening and relengthening were comparable to previous observations in intact heart muscle preparations. The dispersion of light diffracted by the cell remained narrow, and individual striations remained distinct and laterally well registered throughout the shortening- relengthening cycle. In contrast, appreciable nonuniformity and internal buckling were seen at sarcomere lengths < 1.8 micrometers when the resting cell, embedded in gelatin, was longitudinally compressed These results indicate (a) that shortening and relengthening is characterized by uniform activation between myofibrils within the cardiac cell and (b) that physiologically significant relengthening forces in living heart muscle originate at the level of the cell rather than in extracellular connections. First-order diffracted light intensity, extremely variable during sarcomere shortening, was always greatest during midrelaxation preceding the onset of a very slow and uniform phase of sarcomere relengthening.     

Peroxynitrite irreversibly decreases diastolic and systolic function in cardiac muscle

Much of the damaging action of nitric oxide in heart may be due to its diffusion-limited reaction with superoxide to form peroxynitrite. Direct infusion of peroxynitrite into isolated perfused hearts fails to model the effects of in situ formation because the bulk of peroxynitrite decomposes before reaching the myocytes. To examine the direct effects of peroxynitrite on the contractile apparatus of the heart, we exposed intact and skinned rat papillary muscles to a steady state concentration of 4-microM peroxynitrite for 5 min, followed by a 30-min recovery period to monitor irreversible effects. In intact muscles developed force fell immediately to 26% of initial force, recovering to 43% by 30 min. Resting tension increased by 600% immediately, and was still elevated 500% by 30 min. Nitrotyrosine immunochemistry showed that peroxynitrite can induce tyrosine nitration at low concentrations and is capable of penetrating 200-380 microm into the papillary muscle after a 5-min infusion. Decomposed peroxynitrite had no effect on either intact or skinned muscle developed force or resting tension. Our results show that peroxynitrite directly damages both developed force and resting tension of isolated heart muscle, which can be extrapolated to systolic and diastolic injury in intact hearts.     

In silico identification of potential calcium dynamics and sarcomere targets for recovering left ventricular function in rat heart failure with preserved ejection fraction

Heart failure with preserved ejection fraction (HFpEF) is a complex disease associated with multiple co-morbidities, where impaired cardiac mechanics are often the end effect. At the cellular level, cardiac mechanics can be pharmacologically manipulated by altering calcium signalling and the sarcomere. However, the link between cellular level modulations and whole organ pump function is incompletely understood. Our goal is to develop and use a multi-scale computational cardiac mechanics model of the obese ZSF1 HFpEF rat to identify important biomechanical mechanisms that underpin impaired cardiac function and to predict how whole-heart mechanical function can be recovered through altering cellular calcium dynamics and/or cellular contraction. The rat heart was modelled using a 3D biventricular biomechanics model. Biomechanics were described by 16 parameters, corresponding to intracellular calcium transient, sarcomere dynamics, cardiac tissue and hemodynamics properties. The model simulated left ventricular (LV) pressure-volume loops that were described by 14 scalar features. We trained a Gaussian process emulator to map the 16 input parameters to each of the 14 outputs. A global sensitivity analysis was performed, and identified calcium dynamics and thin and thick filament kinetics as key determinants of the organ scale pump function. We employed Bayesian history matching to build a model of the ZSF1 rat heart. Next, we recovered the LV function, described by ejection fraction, peak pressure, maximum rate of pressure rise and isovolumetric relaxation time constant. We found that by manipulating calcium, thin and thick filament properties we can recover 34%, 28% and 24% of the LV function in the ZSF1 rat heart, respectively, and 39% if we manipulate all of them together. We demonstrated how a combination of biophysically based models and their derived emulators can be used to identify potential pharmacological targets. We predicted that cardiac function can be best recovered in ZSF1 rats by desensitising the myofilament and reducing the affinity to intracellular calcium concentration and overall prolonging the sarcomere staying in the active force generating state.     

Viscoelastic material properties of the myocardium and cardiac jelly in the looping chick heart

Accurate material properties of developing embryonic tissues are a crucial factor in studies of the mechanics of morphogenesis. In the present work, we characterize the viscoelastic material properties of the looping heart tube in the chick embryo through nonlinear finite element modeling and microindentation experiments. Both hysteresis and ramp-hold experiments were performed on the intact heart and isolated cardiac jelly (extracellular matrix). An inverse computational method was used to determine the constitutive relations for the myocardium and cardiac jelly. With both layers assumed to be quasilinear viscoelastic, material coefficients for an Ogden type strain-energy density function combined with Prony series of two terms or less were determined by fitting numerical results from a simplified model of a heart segment to experimental data. The experimental and modeling techniques can be applied generally for determining viscoelastic material properties of embryonic tissues.     

Effect of Global Cardiac Ischemia on Human Ventricular Fibrillation: Insights from a Multi-scale Mechanistic Model of the Human Heart

Acute regional ischemia in the heart can lead to cardiac arrhythmias such as ventricular fibrillation (VF), which in turn compromise cardiac output and result in secondary global cardiac ischemia. The secondary ischemia may influence the underlying arrhythmia mechanism. A recent clinical study documents the effect of global cardiac ischaemia on the mechanisms of VF. During 150 seconds of global ischemia the dominant frequency of activation decreased, while after reperfusion it increased rapidly. At the same time the complexity of epicardial excitation, measured as the number of epicardical phase singularity points, remained approximately constant during ischemia. Here we perform numerical studies based on these clinical data and propose explanations for the observed dynamics of the period and complexity of activation patterns. In particular, we study the effects on ischemia in pseudo-1D and 2D cardiac tissue models as well as in an anatomically accurate model of human heart ventricles. We demonstrate that the fall of dominant frequency in VF during secondary ischemia can be explained by an increase in extracellular potassium, while the increase during reperfusion is consistent with washout of potassium and continued activation of the ATP-dependent potassium channels. We also suggest that memory effects are responsible for the observed complexity dynamics. In addition, we present unpublished clinical results of individual patient recordings and propose a way of estimating extracellular potassium and activation of ATP-dependent potassium channels from these measurements.