Anomalous ion diffusion within skeletal muscle transverse tubule networks
© Shorten and Soboleva; licensee BioMed Central Ltd. 2007
Received: 08 March 2007
Accepted: 17 May 2007
Published: 17 May 2007
Skeletal muscle fibres contain transverse tubular (t-tubule) networks that allow electrical signals to rapidly propagate into the fibre. These electrical signals are generated by the transport of ions across the t-tubule membranes and this can result in significant changes in ion concentrations within the t-tubules during muscle excitation. During periods of repeated high-frequency activation of skeletal muscle the t-tubule K+ concentration is believed to increase significantly and diffusive K+ transport from the t-tubules into the interstitial space provides a mechanism for alleviating muscle membrane depolarization. However, the tortuous nature of the highly branched space-filling t-tubule network impedes the diffusion of material through the network. The effective diffusion coefficient for ions in the t-tubules has been measured to be approximately five times lower than in free solution, which is significantly different from existing theoretical values of the effective diffusion coefficient that range from 2–3 times lower than in free solution. To resolve this discrepancy, in this paper we study the process of diffusion within electron microscope scanned sections of the skeletal muscle t-tubule network using mathematical modelling and computer simulation techniques. Our model includes t-tubule geometry, tautness, hydrodynamic and non-planar network factors.
Using our model we found that the t-tubule network geometry reduced the K+ diffusion coefficient to 19–27% of its value in free solution, which is consistent with the experimentally observed value of 21% and is significantly smaller than existing theoretical values that range from 32–50%. We also found that diffusion in the t-tubules is anomalous for skeletal muscle fibres with a diameter of less than approximately 10–20 μm as a result of obstructed diffusion. We also observed that the [K+] within the interior of the t-tubule network during high-frequency activation is greater for fibres with a larger diameter. Smaller skeletal muscle fibres are therefore more resistant to membrane depolarization. Because the t-tubule network is anisotropic and inhomogeneous, we also found that the [K+] distribution generated within the network was irregular for fibres of small diameter.
Our model explains the measured effective diffusion coefficient for ions in skeletal muscle t-tubules.
Skeletal muscle fibres contain transverse tubular (t-tubule) networks that provide for rapid propagation of electrical signals into the fibre and ensures near simultaneous contraction in the constituent myofibrils. These t-tubule networks are highly branched space-filling networks that are located near sarcomere Z-lines in amphibians and the A-I junction in mammals. The t-tubule network therefore largely lies in a plane perpendicular to the axis of the muscle fibre. The electrical signals that propagate along the planar t-tubules are generated by the transport of ions across the t-tubule membranes and this can result in significant changes in ion concentrations within the t-tubules during muscle excitation. In particular, K+ ions accumulate within the t-tubule network as a result of voltage-gated K+ channels that repolarize the membrane during action potentials. This K+ accumulation can result in membrane depolarization and reduced membrane excitability as a result of Na+ channel inactivation and consequently reduced muscle force output . The K+ ionic gradients across the t-tubule membranes are re-established by the Na+-K+ exchanger following an action potential . During periods of repeated high-frequency activation of skeletal muscle the t-tubule K+ concentration ([K+] t ) is believed to increase significantly [3, 4] and diffusive K+ transport from the t-tubules into the interstitial space provides a further mechanism for alleviating muscle membrane depolarization. In this paper we study the process of K+ diffusion within the skeletal muscle t-tubule network using mathematical modelling and computer simulation techniques.
The structure of biological tissue impedes the diffusion and transport of material through the tissue. In particular it has been calculated that the effective diffusion coefficient for ions in the t-tubules is approximately five times lower than in free solution [5, 6]. In this paper we analyze diffusion on the skeletal muscle t-tubule networks imaged by Peachey and Eisenberg  and Hayashi et al. . We have found that the t-tubule network restricts the diffusion of material through the network as a result of geometric, hydrodynamic and network tautness factors. We also find that diffusion in the t-tubules over short distances is anomalous (i.e. does not obey Brownian motion).
Diffusion within the t-tubule network
where D is the temperature-dependent solute diffusion coefficient in free solution, c is the solute concentration in the tissue, t is time and D app is the apparent diffusion coefficient. The tortuosity factor then determines the difference between the apparent diffusion coefficient in the tissue and the diffusion coefficient in free solution. This equation then describes the macroscopic transport of the solute through the tissue. The mean squared displacement (MSD) of the solute is⟨r2(t)⟩ = 2dτDt
where d is the spatial dimensionality. Several regular two-dimensional networks of one-dimensional t-tubules have been used to model the t-tubule geometry. Tortuosity factors have been calculated for these networks using flux density analysis and homogenization methods and it has been found that τ = 0.5 for regular hexagonal, square and trigonal network geometries and τ = 0.4 for staggered square network geometries [13, 17–19]. Furthermore, using theoretical methods, Mathias et al.  have estimated that τ = 0.32 from the morphometric parameters of frog skeletal muscle t-tubule networks. The irregular structure of the t-tubule network therefore significantly reduces diffusion on the network.
Experimental measurements of the t-tubule tortuosity factor have been made for a number of different ionic species. Tortuosity measurements range from τ = 0.11 for Ca2+ to τ = 0.21 for K+, which are significantly lower than the tortuosity factors calculated on theoretical grounds (0.32 ≤ τ ≤ 0.5). For example, Fujishiro and Kawata  investigated Na+ diffusion into the t-tubule network when frog skeletal muscle fibres were moved from a Na+ free solution into normal Ringer solution with [Na+] = 111 mM. The authors employed a diffusion equation and found that the diffusion of Na+ along the transverse tubules was consistent with an apparent diffusion coefficient of = 3.11 × 10-6 cm2/s, which is significantly less than that for an aqueous solution at 25°C (14.8 × 10-6 cm2/s). Using Markov chain Monte Carlo we found that the SEM associated with this calculated apparent diffusion coefficient is 0.3 cm2/s. This data is consistent with a t-tubule network total tortuosity factor of τ = 0.21 ± 0.02. Almers  investigated K+ diffusion in frog sartorius t-tubules using a voltage-clamp technique at 22°C. Using a computer model they found that their data was consistent with an apparent t-tubule K+ diffusion coefficient of = 3.8 × 10-6 cm2/s. We found using Markov chain Monte Carlo that the SEM associated with this calculated apparent diffusion coefficient is 0.1 cm2/s. The K+ diffusion coefficient in free solution at 25°C is 18.3 × 10-6 cm2/s . This data is consistent with a t-tubule network total tortuosity factor of τ = 0.21 ± 0.005. Almers et al.  investigated Ca2+ depletion in frog muscle tubules using a voltage-clamp technique at 20–24°C. They found that their data was consistent with an apparent t-tubule Ca2+ diffusion coefficient of = 0.87 × 10-6 cm2/s with a SEM of 0.18 cm2/s. The Ca2+ diffusion coefficient in free solution at 25°C is 7.7 × 10-6 cm2/s . This data is therefore consistent with a t-tubule network total tortuosity factor of τ = 0.11 ± 0.02. Ca2+ binding may possibly explain the lower tortuosity factor for Ca2+. The experimentally observed t-tubule network total tortuosity factors (0.11 ≤ τ ≤ 0.21) are therefore significantly smaller than existing estimates of the tortuosity factor based on theoretical considerations (0.32 ≤ τ ≤ 0.5).
Anomalous and obstructed K+ diffusion in the t-tubule network
⟨r2(t)⟩ = K α t α , α ≠ 1
where Γ(α) is the Gamma function.
On sections of the t-tubule network with a diameter greater than approximately 11–15 μm diffusion can be described by the standard homogenised diffusion equation (Eq. 1). The random walk simulations can therefore be used to calculate the tissue tortuosity factors. T-tubule network tortuosity factors of τ = 0.25 and τ = 0.35 were calculated for the human and frog skeletal muscle t-tubule geometries respectively. The geometric tortuosity factors for these networks calculated by random walk simulation are consistent with the previous findings of Mathias et al.  who used flux density analysis to estimate that τ = 0.32 from the morphometric parameters of frog skeletal muscle t-tubule networks. These theoretical estimates of the t-tubule network tortuosity factor (τ = 0.25, 0.35, 0.32) are nonetheless significantly larger than experimentally observed total tortuosity factors (0.11 ≤ τ ≤ 0.21). We believe that this inconsistency is due to non-planar and hydrodynamic effects. If we include these effects into our theoretical model then we obtain a tortuosity factor that is consistent with the experimental data.
Skeletal muscle t-tubule networks are generally conceptualised as planar networks perpendicular to the fibre direction that are located at the boundary between the A-band and I-band regions in vertebrate fibres. Nevertheless, electron microscope measurements of longitudinal cross-sections of human skeletal muscle t-tubule networks indicate that the network is not planar and lies within a region about 300–400 nm in width . T-tubule branches are therefore greater in length than they appear in two-dimensional transverse images of t-tubule networks. We estimate from images of longitudinal cross-sections of t-tubule networks by Hayashi et al.  that the t-tubule branch lengths are approximately 10% greater than depicted in Fig. 2. From the relationship between the diffusion coefficient, space and time defined by Eq. (2), it follows that that the effective diffusion coefficient on this non-planar t-tubule network is reduced by the factor 1/1.12 = 0.83.
and . K+ diffusion in the t-tubules is therefore 8% slower than in free solution. Hydrodynamic wall effects within the small diameter t-tubules therefore have a small but significant effect on K+ diffusion within the t-tubules.
Our modelling analysis therefore predicts a total tortuosity factor for K+ within the t-tubule network of 0.30 × 0.92 × 0.83 = 0.23 (geometric, hydrodynamic and non-planar factors), which agrees well with the experimental observations of Almers et al.  who found a t-tubule network total tortuosity factor for K+ of τ = 0.21. Although there is variability in the t-tubule branch diameters within the t-tubule network [10, 29], and therefore variability in the level of hydrodynamic wall effects at different locations in the network, we found that this has a negligible effect on the t-tubule network total tortuosity factor for K+. We use a standard deviation in the diameter of the t-tubule branches of 1–3 nm [10, 29] and assume that each t-tubule branch contains 3 segments of different diameter . T-tubule network geometries are also different in different muscle fibres. For example, the t-tubule network from slow-twitch guinea pig skeletal muscle is more non-planar than that from fast-twitch skeletal muscle . T-tubule branch lengths are 3% greater in length than in transverse images of t-tubules in fast-twitch guinea pig skeletal muscle and 30% greater in length than in transverse images of t-tubules in slow twitch fibres . From the relationship between the diffusion coefficient, space and time defined by Eq. (2), it follows that the effective diffusion coefficient on a non-planar t-tubule network is reduced by the factor 1/1.032 = 0.94 in fast-twitch guinea pig skeletal muscle and by 1/1.32= 0.59 in slow-twitch fibres. This suggests that the effective diffusion coefficient on the t-tubule network in slow-twitch fibres is 37% smaller than in fast-twitch fibres.
The effect of the t-tubule network geometry on the [K+] distribution within t-tubules and membrane depolarization
where Ω denotes the region within the t-tubule network, Ω' denotes the t-tubule membrane, Γ denotes the boundary between the t-tubule network and the interstitial space, c0 = 4 mM is the [K+] in the interstitial space, is the outward normal at the t-tubule membrane, B = 0.002 μm-1 represents the net rate of K+ release across the t-tubule membrane, A = 60 mM represents the reduction in K+ release from voltage-gated K+ channels due to t-tubule K+ accumulation and D is the K+ diffusion coefficient in free solution. The parameters A and B were chosen to reproduce qualitatively the macroscopic t-tubule [K+] gradient generated in the model by Wallinga et al. . Eq. (11) was solved using the Matlab software package (The MathWorks) with the method of lines .
Experiments show that the skeletal muscle t-tubule network significantly impedes the diffusive transport of material through the t-tubules. Tortuosity measurements range from τ = 0.11 for Ca2+ to τ = 0.21 for K+ and Na+ [5, 6, 21], which are significantly lower than the tortuosity factors calculated on theoretical grounds (0.32 ≤ τ ≤ 0.5). Our work demonstrated why previous theoretical studies overestimate the t-tubule tortuosity factor. We calculated total tortuosity factors of τ = 0.19 and τ = 0.27 respectively for K+ transport in human  and frog  skeletal muscle t-tubule geometries. Our theoretical estimate of the tortuosity factor, which includes geometric, hydrodynamic and non-planar effects, is consistent with experimental observations. A contributing reason for the difference in the tortuosity factor τ between the human and frog networks is that the latter network is a traced taut reconstruction of the network (i.e. the t-tubule branches do not contain wiggles). However, it is evident from scanned transverse sections in Fig. 2 that the t-tubule branches do contain wiggles and are therefore not taut. These t-tubule network wiggles were predicted to exist by Almers  based on the observation that muscle fibres can significantly shorten without damage (if muscle contracts under constant volume the t-tubule network must stretch). The network parameters obtained for the human network are therefore potentially more representative of skeletal muscle tissue.
We also observed that diffusion in the t-tubules is anomalous over short distances (i.e. does not obey Brownian motion) as a result of obstructed diffusion. The level of anomalous diffusion is described by an anomalous exponent (α = 1 represents Brownian diffusion) and we found anomalous exponents of α = 0.86 and α = 0.91 for the human and frog networks respectively. Obstructed diffusion is a localised phenomenon for an anomalous exponent (α) greater than the percolation threshold 0.697 . For this reason the anomalous exponent α is also a function of the local t-tubule network neighbourhood defined by R*. The anomalous exponents reported in this paper are averages over the network geometry and for the frog network we found that 0.88 ≤ α ≤ 0.95 and 5 ≤ R* ≤11 μm, which is a measure of network inhomogeneity. For small muscle fibres (diameter less than 10–20 μm) diffusion on the t-tubule network can be described by the fractional diffusion equation (Eq. 7), whereas for large muscle fibres (diameter greater than 20 μm) the homogenised diffusion equation (Eq. 1) provides an accurate approximation to diffusion on the t-tubule network.
A significant [K+] gradient is believed to be generated within the t-tubule network of skeletal muscle fibres during repeated high-frequency activation . We found that the [K+] within the interior of the t-tubule network during high-frequency activation is greater for fibres with a larger diameter and therefore smaller skeletal muscle fibres are more resistant to membrane depolarization. We also found that fast-twitch muscle fibres have larger t-tubule tortuosity factors than slow-twitch fibres and are therefore better able to remove K+ from the t-tubules via diffusive transport. Fast-twitch fibres are more susceptible to membrane depolarization and faster K+ diffusion in the t-tubules therefore allows fast-twitch fibres to better control membrane depolarization. Because the t-tubule network is anisotropic and inhomogeneous, we also observed that the [K+] distribution generated within the network was irregular. Although the level of irregularity in the [K+] distribution within the network was small for larger sized networks (> 20 μm), it was significant for fibres of small diameter (< 20 μm). T-tubule networks within large muscle fibres can therefore be considered as isotropic irregular networks and a homogenized description of the network can be safely employed to model the t-tubule system [4, 32].
- Sejersted OM, Sjogaard G: Dynamics and consequences of potassium shifts in skeletal muscle and heart during exercise. Physiol Rev. 2000, 80 (4): 1411-1481.PubMedGoogle Scholar
- Neilson OB, Ortenblad N, Lamb GD, Stephenson DG: Excitability of the t-tubular system in rat skeletal muscle: roles of K+ and Na+ gradients and Na+-K+ pump activity. J Physiol. 2004, 557 (1): 133-146. 10.1113/jphysiol.2003.059014.View ArticleGoogle Scholar
- Westerblad H, Lee JA, Lamb AG, Bolsover SR, Allen DG: Spatial gradients of intracellular calcium in skeletal muscle during fatigue. Pflug Arch Eur J Phy. 1990, 415: 734-740. 10.1007/BF02584013.View ArticleGoogle Scholar
- Wallinga W, Meijer SL, Alberink MJ, Vliek M, Wienk ED, Ypey DL: Modelling action potentials and membrane currents of mammalian skeletal muscle fibres in coherence with potassium concentration changes in the T-tubular system. Eur Biophys J. 1999, 28: 317-329. 10.1007/s002490050214.View ArticlePubMedGoogle Scholar
- Almers W: Potassium concentration changes in the transverse tubules of vertebrate skeletal muscle. Fed Proc. 1980, 39 (5): 1527-1532.PubMedGoogle Scholar
- Fujishiro N, Kawata H: Tetanus responses under rapid bath solution change: Electrotonic depolarization of transverse tubules may release Ca2+ from sarcoplasmic reticulum of Rana japonica skeletal muscle. Comp Biochem Phys A. 1992, 103 (4): 661-666. 10.1016/0300-9629(92)90163-K.View ArticleGoogle Scholar
- Peachey LD, Eisenberg BR: Helicoids in the t system and striations of frog skeletal muscle fibres seen by high voltage electron microscopy. Biophys J. 1978, 22: 145-154.PubMed CentralView ArticlePubMedGoogle Scholar
- Hayashi K, Miller RG, Brownell AKW: Three-dimensional architecture of sarcoplasmic reticulum and t-system in human skeletal muscle. Anat Rec. 1987, 218: 275-283. 10.1002/ar.1092180308.View ArticlePubMedGoogle Scholar
- Soeller C, Cannell MB: Examination of the transverse tubular system in living cardiac rat myocytes by 2-photon microscopy and digital image-processing techniques. Circ Res. 1999, 84: 266-275.View ArticlePubMedGoogle Scholar
- Ploug T, van Deurs B, Ai H, Cushman SW, Ralston E: Analysis of GLUT4 distribution in whole skeletal muscle fibers: Identification of distinct storage compartments that are recruited by insulin and muscle contractions. J Cell Biol. 1998, 142 (6): 1429-1446. 10.1083/jcb.142.6.1429.PubMed CentralView ArticlePubMedGoogle Scholar
- Eisenberg BR, Kuda AM, Peter JB: Stereological analysis of mammalian skeletal muscle I. Soleus Muscle of the Adult Guinea Pig. J Cell Biol. 1974, 60: 732-754. 10.1083/jcb.60.3.732.PubMed CentralView ArticlePubMedGoogle Scholar
- Franzini-Armstrong C, Ferguson DG, Champ C: Discrimination between fast- and slow-twitch fibres of guinea pig skeletal muscle using relative surface density of junctional transverse tubule membrane. J Muscle Res Cell M. 1988, 9: 403-414. 10.1007/BF01774067.View ArticleGoogle Scholar
- Mathias RT, Eisenberg RS, Valdiosera R: Electrical properties of frog skeletal muscle fibers interpreted with a mesh model of the tubular system. Biophys J. 1977, 17: 57-93.PubMed CentralView ArticlePubMedGoogle Scholar
- Mathias RT: Effect of tortuous extracellular pathways on resistance measurements. Biophys J. 1983, 42: 55-59.PubMed CentralView ArticlePubMedGoogle Scholar
- El-Kareh AW, Braunstein SL, Secomb TW: Effect of cell arrangement and interstitial volume fraction on the diffusivity of monoclonal antibodies in tissue. Biophys J. 1993, 64: 1638-1646.PubMed CentralView ArticlePubMedGoogle Scholar
- Friedrich O, Ehmer T, Uttenweiler D, Vogel M, Barry PH: Numerical analysis of Ca2+ depletion in the transverse tubular system of mammalian muscle. Biophys J. 2001, 80: 2046-2055.PubMed CentralView ArticlePubMedGoogle Scholar
- Adrian RH, Chandler WK, Hodgkin AL: The kinetics of mechanical activation in frog muscle. J Physiol. 1969, 204: 207-230.PubMed CentralView ArticlePubMedGoogle Scholar
- Schneider MF: Linear electrical properties of the transverse tubules and surface membrane of skeletal muscle fibres. J Gen Physiol. 1970, 56: 640-671. 10.1085/jgp.56.5.640.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen KC, Nicholson C: Changes in brain cell shape create residual extracellular space volume and explain tortuosity behavior during osmotic challenge. P Natl Acad Sci USA. 2000, 97: 8306-8311. 10.1073/pnas.150338197.View ArticleGoogle Scholar
- Koneshan S, Rasaiah JC, Lynden-Bell RM, Lee SH: Solvent structure, dynamics, and ion mobility in aqueous solutions at 25°C. J Phys Chem B. 1998, 102: 4193-4204. 10.1021/jp980642x.View ArticleGoogle Scholar
- Almers W, Fink R, Palade PT: Calcium depletion in frog muscle tubules: the decline of calcium current under maintained depolarization. J Physiol. 1981, 312: 177-207.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang JH: Tracer-diffusion in liquids. IV. Self-diffusion of calcium ion and chloride ion in aqueous calcium chloride solutions. J Am Chem Soc. 1953, 75: 1769-1770. 10.1021/ja01103a539.View ArticleGoogle Scholar
- Metzler R, Klafter J: The random walks guide to anomalous diffusion: A fractional dynamics approach. Phys Rep. 2000, 339: 1-77. 10.1016/S0370-1573(00)00070-3.View ArticleGoogle Scholar
- Saxton MJ: Anomalous diffusion due to obstacles: A Monte Carlo study. Biophys J. 1994, 66: 394-401.PubMed CentralView ArticlePubMedGoogle Scholar
- Saxton MJ: Anomalous subdiffusion in fluorescence photobleaching recovery: A Monte Carlo study. Biophys J. 2001, 81: 2226-2240.PubMed CentralView ArticlePubMedGoogle Scholar
- Olveczky BP, Verkman AS: Monte Carlo analysis of obstructed diffusion in three dimensions: Application to molecular diffusion in organelles. Biophys J. 1998, 74: 2722-2730.PubMed CentralView ArticlePubMedGoogle Scholar
- Havlin S, Ben-Avraham D: Diffusion in disordered media. Adv Phys. 1987, 36 (6): 695-798. 10.1080/00018738700101072.View ArticleGoogle Scholar
- Nitsche JM, Balgi G: Hindered Brownian diffusion of spherical solutes within circular cylindrical pores. Ind Eng Chem Res. 1994, 33: 2242-2247. 10.1021/ie00033a030.View ArticleGoogle Scholar
- Fujimaki N, Peachey LD, Murakami T, Ishikawa H: Three-dimensional visualization of the T-system in fixed and embedded frog skeletal muscle fibers by confocal laser scanning reflection microscopy. Bioimages. 1993, 1: 167-174.Google Scholar
- Lee HS, Matthews CJ, Braddock RD, Sander GC, Gandola F: A MATLAB method of lines template for transport equations. Environ Model Software. 2004, 19: 603-614. 10.1016/j.envsoft.2003.08.017.View ArticleGoogle Scholar
- Sutherland H, Salmons S, Ramnarine IR, Capoccia M, Walsh AA, Jarvis JC: Adaptive conditioning of skeletal muscle in a large animal model (Sus domesticus). J Anat. 2005, 209: 165-177. 10.1111/j.1469-7580.2006.00598.x.View ArticleGoogle Scholar
- Henneberg K, Roberge FA: Simulation of propagation along an isolated skeletal muscle fiber in an isotropic volume conductor. Ann Biomed Eng. 1997, 25: 5-28. 10.1007/BF02647346.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.