Search

Contractile function of cardiac cells is driven by the sliding displacement of myofilaments powered by the cycling myosin crossbridges. Critical to this process is the availability of ATP, which myosin hydrolyzes during the cross-bridge cycle. The diffusion of adenine nucleotides through the myofilament lattice has been shown to be anisotropic, with slower radial diffusion perpendicular to the filament axis relative to parallel, and is attributed to the periodic hexagonal arrangement of the thin (actin) and thick (myosin) filaments. We investigated whether atomistic-resolution details of myofilament proteins can refine coarse-grain estimates of diffusional anisotropy for adenine nucleotides in the cardiac myofibril, using homogenization theory and atomistic thin filament models from the Protein Data Bank. Our results demonstrate considerable anisotropy in ATP and ADP diffusion constants that is consistent with experimental measurements and dependent on lattice spacing and myofilament overlap. A reaction-diffusion model of the half-sarcomere further suggests that diffusional anisotropy may lead to modest adenine nucleotide gradients in the myoplasm under physiological conditions.

Contractile function of cardiac cells is driven by the sliding displacement of myofilaments powered by the cycling myosin crossbridges. Critical to this process is the availability of ATP, which myosin hydrolyzes during the cross-bridge cycle. The diffusion of adenine nucleotides through the myofilament lattice has been shown to be anisotropic, with slower radial diffusion perpendicular to the filament axis relative to parallel, and is attributed to the periodic hexagonal arrangement of the thin (actin) and thick (myosin) filaments. We investigated whether atomistic-resolution details of myofilament proteins can refine coarse-grain estimates of diffusional anisotropy for adenine nucleotides in the cardiac myofibril, using homogenization theory and atomistic thin filament models from the Protein Data Bank. Our results demonstrate considerable anisotropy in ATP and ADP diffusion constants that is consistent with experimental measurements and dependent on lattice spacing and myofilament overlap. A reaction-diffusion model of the half-sarcomere further suggests that diffusional anisotropy may lead to modest adenine nucleotide gradients in the myoplasm under physiological conditions.

The transverse tubular system of rabbit ventricular myocytes consists of cell membrane invaginations (t-tubules) that are essential for efficient cardiac excitation-contraction coupling. In this study, we investigate how t-tubule micro-anatomy, L-type Ca(2+) channel (LCC) clustering, and allosteric activation of Na(+)/Ca(2+) exchanger by L-type Ca(2+) current affects intracellular Ca(2+) dynamics. Our model includes a realistic 3D geometry of a single t-tubule and its surrounding half-sarcomeres for rabbit ventricular myocytes. The effects of spatially distributed membrane ion-transporters (LCC, Na(+)/Ca(2+) exchanger, sarcolemmal Ca(2+) pump, and sarcolemmal Ca(2+) leak), and stationary and mobile Ca(2+) buffers (troponin C, ATP, calmodulin, and Fluo-3) are also considered. We used a coupled reaction-diffusion system to describe the spatio-temporal concentration profiles of free and buffered intracellular Ca(2+). We obtained parameters from voltage-clamp protocols of L-type Ca(2+) current and line-scan recordings of Ca(2+) concentration profiles in rabbit cells, in which the sarcoplasmic reticulum is disabled. Our model results agree with experimental measurements of global Ca(2+) transient in myocytes loaded with 50 μM Fluo-3. We found that local Ca(2+) concentrations within the cytosol and sub-sarcolemma, as well as the local trigger fluxes of Ca(2+) crossing the cell membrane, are sensitive to details of t-tubule micro-structure and membrane Ca(2+) flux distribution. The model additionally predicts that local Ca(2+) trigger fluxes are at least threefold to eightfold higher than the whole-cell Ca(2+) trigger flux. We found also that the activation of allosteric Ca(2+)-binding sites on the Na(+)/Ca(2+) exchanger could provide a mechanism for regulating global and local Ca(2+) trigger fluxes in vivo. Our studies indicate that improved structural and functional models could improve our understanding of the contributions of L-type and Na(+)/Ca(2+) exchanger fluxes to intracellular Ca(2+) dynamics.

The transverse tubular system of rabbit ventricular myocytes consists of cell membrane invaginations (t-tubules) that are essential for efficient cardiac excitation-contraction coupling. In this study, we investigate how t-tubule micro-anatomy, L-type Ca(2+) channel (LCC) clustering, and allosteric activation of Na(+)/Ca(2+) exchanger by L-type Ca(2+) current affects intracellular Ca(2+) dynamics. Our model includes a realistic 3D geometry of a single t-tubule and its surrounding half-sarcomeres for rabbit ventricular myocytes. The effects of spatially distributed membrane ion-transporters (LCC, Na(+)/Ca(2+) exchanger, sarcolemmal Ca(2+) pump, and sarcolemmal Ca(2+) leak), and stationary and mobile Ca(2+) buffers (troponin C, ATP, calmodulin, and Fluo-3) are also considered. We used a coupled reaction-diffusion system to describe the spatio-temporal concentration profiles of free and buffered intracellular Ca(2+). We obtained parameters from voltage-clamp protocols of L-type Ca(2+) current and line-scan recordings of Ca(2+) concentration profiles in rabbit cells, in which the sarcoplasmic reticulum is disabled. Our model results agree with experimental measurements of global Ca(2+) transient in myocytes loaded with 50 μM Fluo-3. We found that local Ca(2+) concentrations within the cytosol and sub-sarcolemma, as well as the local trigger fluxes of Ca(2+) crossing the cell membrane, are sensitive to details of t-tubule micro-structure and membrane Ca(2+) flux distribution. The model additionally predicts that local Ca(2+) trigger fluxes are at least threefold to eightfold higher than the whole-cell Ca(2+) trigger flux. We found also that the activation of allosteric Ca(2+)-binding sites on the Na(+)/Ca(2+) exchanger could provide a mechanism for regulating global and local Ca(2+) trigger fluxes in vivo. Our studies indicate that improved structural and functional models could improve our understanding of the contributions of L-type and Na(+)/Ca(2+) exchanger fluxes to intracellular Ca(2+) dynamics.

The t-tubules of mammalian ventricular myocytes are invaginations of the cell membrane that occur at each Z-line. These invaginations branch within the cell to form a complex network that allows rapid propagation of the electrical signal, and hence synchronous rise of intracellular calcium (Ca(2+)). To investigate how the t-tubule microanatomy and the distribution of membrane Ca(2+) flux affect cardiac excitation-contraction coupling we developed a 3-D continuum model of Ca(2+) signaling, buffering and diffusion in rat ventricular myocytes. The transverse-axial t-tubule geometry was derived from light microscopy structural data. To solve the nonlinear reaction-diffusion system we extended SMOL software tool (http://mccammon.ucsd.edu/smol/). The analysis suggests that the quantitative understanding of the Ca(2+) signaling requires more accurate knowledge of the t-tubule ultra-structure and Ca(2+) flux distribution along the sarcolemma. The results reveal the important role for mobile and stationary Ca(2+) buffers, including the Ca(2+) indicator dye. In agreement with experiment, in the presence of fluorescence dye and inhibited sarcoplasmic reticulum, the lack of detectible differences in the depolarization-evoked Ca(2+) transients was found when the Ca(2+) flux was heterogeneously distributed along the sarcolemma. In the absence of fluorescence dye, strongly non-uniform Ca(2+) signals are predicted. Even at modest elevation of Ca(2+), reached during Ca(2+) influx, large and steep Ca(2+) gradients are found in the narrow sub-sarcolemmal space. The model predicts that the branched t-tubule structure and changes in the normal Ca(2+) flux density along the cell membrane support initiation and propagation of Ca(2+) waves in rat myocytes.

The t-tubules of mammalian ventricular myocytes are invaginations of the cell membrane that occur at each Z-line. These invaginations branch within the cell to form a complex network that allows rapid propagation of the electrical signal, and hence synchronous rise of intracellular calcium (Ca(2+)). To investigate how the t-tubule microanatomy and the distribution of membrane Ca(2+) flux affect cardiac excitation-contraction coupling we developed a 3-D continuum model of Ca(2+) signaling, buffering and diffusion in rat ventricular myocytes. The transverse-axial t-tubule geometry was derived from light microscopy structural data. To solve the nonlinear reaction-diffusion system we extended SMOL software tool (http://mccammon.ucsd.edu/smol/). The analysis suggests that the quantitative understanding of the Ca(2+) signaling requires more accurate knowledge of the t-tubule ultra-structure and Ca(2+) flux distribution along the sarcolemma. The results reveal the important role for mobile and stationary Ca(2+) buffers, including the Ca(2+) indicator dye. In agreement with experiment, in the presence of fluorescence dye and inhibited sarcoplasmic reticulum, the lack of detectible differences in the depolarization-evoked Ca(2+) transients was found when the Ca(2+) flux was heterogeneously distributed along the sarcolemma. In the absence of fluorescence dye, strongly non-uniform Ca(2+) signals are predicted. Even at modest elevation of Ca(2+), reached during Ca(2+) influx, large and steep Ca(2+) gradients are found in the narrow sub-sarcolemmal space. The model predicts that the branched t-tubule structure and changes in the normal Ca(2+) flux density along the cell membrane support initiation and propagation of Ca(2+) waves in rat myocytes.