1. Introduction
Muscles in adult mammals and birds are responsible for most heat generation to keep their organisms thermoregulated [3]. Nevertheless, it has recently been stated that they could also act as a thermosensitive system [4]. Experiments conducted over muscle fibers from cardiac tissue reveal that local thermal perturbation caused using infrared (IR) laser light influences the dynamic response of the fibers. The thermal excitability of real striated myofibrils was addressed before isolating an acting chain and applying a heat pulse [5,6]. Although a skeletal muscle fiber is susceptible to displaying a similar response when stimulated with a spot laser IR, all efforts have been centered on cardiac cells. One example is the increasing number of studies showing how applying laser light pulses makes it possible to pace from groups of cardiomyocytes to an embryonic heart [4,7,8]. Besides muscle systems, local thermal gradients have been recognized to affect intracellular behavior and alter physiological functions [9].
Two explanations to address the thermal influence in cardiomyocytes have been proposed. The first one treats the excitation-contraction process from a global perspective and suggests that the temperature gradient affects ion channel openings; this process is related to Ca2+ releasing from the sarcoplasmic reticulum. Rabbitt's group experiments [10] for photothermal excitation evince transient capacitive currents in a neuromuscular junction. Therefore, we should consider the interplays between myosin, tropomyosin, and actin instead of focusing on single cardiomyocytes or diverting attention away from the sarcomere. In this respect, Ishiwata's group found cardiomyocyte contractions without Ca2+. They adduced contractions to a partial dissociation of tropomyosin from F-actin induced by heating, a self-reliant process from Ca2+ concentration [4]. Additionally, to the mechanism involved in the photothermal contractions, it seems that contractions depend strongly on the magnitude of the temperature increase (
Until now, cardiac tissue has been excited externally using a laser light pulse that punctually heats the water, generating a thermal gradient in its surroundings. Furthermore, one of the primary purposes is to refine this tool to study living muscle tissue with minimum interference, something challenging to achieve through electrical stimulation; or even thinking of possible therapeutic applications. However, considering a muscle fiber, this thermal gradient also could be internally produced by the Ca2+-ATPase sarcoplasmic reticulum, knowing that this pump, under some circumstances, does not follow a strict cycle [11,12]. When this happens, a 'slippage' or skip of the regular steps is observed, traduced in consuming ATP to produce heat, a thermogenic mechanism. In other contexts, for example, under skeletal muscle fatigue, such thermal gradient could be provided by the same intense muscle use.
When working out the skeletal muscle until fatigue, cellular metabolism is under stress, and preserving the energy state becomes a priority to the system [13]. In similar stress-like conditions, Yang et al. [1] have shown that myocytes from rat exhibit glycolytic oscillations and ATP oscillations. Thus, energy-supplying processes can display oscillatory behavior until physiological conditions are favorable again, which, although not confirmed experimentally, could be the case for a skeletal muscle fiber. Considering the internal thermal gradients generated as a consequence of muscle fatigue [14] or thermogenic process, they can be matched advantageously with biochemical oscillations to drive chemical energy production in the muscle tissue at the cellular level.
As presented above, coupling a thermal gradient with glycolytic oscillations is not an unlikely biochemical scenario. Theoretically, all that is needed is a model reaction capable of producing thermal and chemical oscillations. A Sal'nikov-modified model, including a third autocatalytic reaction, fulfills this requirement [2]. In that case, shifting from chemical oscillations (using autocatalytic reaction) to mixed thermal-chemical oscillations is possible.
Furthermore, if this model takes place in each functional unit of a single fiber, i.e., each sarcomere, coupling several Sal'nikov-modified reactors will be equivalent to a single fiber. Each sarcomere is treated as a black box, sensitive to biochemical species and temperature changes. We are not scaling the model to reproduce the exact physical realization in a microscopic fiber scale because the main interest is the formulation of the governing equations to understand the interplay between thermal gradients and biochemical oscillations. With this restriction, we showed that several Sal'nikov-modified units (reactors) coupled through heat exchange let thermal gradients develop within a 1-dimensional macro sarcomere array (an ideal macroscopic analog of a muscle fiber). Moreover, the resulting fiber model is thermally perturbed, and its behavior is analyzed in the biochemical context of a skeletal muscle contraction.
2. Materials and methods
2.1. Modified model Sal'nikov for chemical energy supplied to a muscle sarcomere
From a macroscopic scale, a striated muscle is formed by muscle fibers composed of myofibrils. Zooming a single myofibril reveals that each consists of cylindrical units called sarcomeres [15]. Sarcomeres are the largest ordered multi-protein assemblies; their most remarkable components are thick and thin filaments interlaced, nearly identical in length [16]. The thin filament comprises actin, troponin, tropomyosin, and myosin II motors from the thick filament at this microscopic scale. Each myosin molecule has two cross bridges that attach cyclically to the actin filament and move the actin along the myosin filament. As both filaments overlap, the sarcomere length gradually decreases. The visible result of summing the contribution of all sarcomeres will be the development of macroscopic force and muscle shortening [17,18]. The hydrolysis of ATP provides the energy required to complete the attaching-detaching cycle to ADP and Pi occurring at the myosin catalytic domain. The chemical reaction is coupled to a mechanical response [19].
Sal'nikov's modified model does not pretend to describe sarcomere at these scales. Instead, it focuses on the biochemical aspect of energetics requirement to make muscle contraction possible. In this sense, a sarcomere is comparable to a chemical reactor where a reaction occurs through ATP supply (macro-sarcomere), and the supply rate constrains the reaction time course. In fact, at the molecular level, the myosin cross-bridge is the transducer that interconverts chemical energy (derived from ATP) to mechanical work [20]. Thus, the proposed model intends to understand how thermal gradients can couple to biochemical oscillations of energetic species relevant to muscle function, i.e., ATP, rather than explain details of the sarcomere mechanism.
Sal'nikov’s modified model proposed before [2] was based upon three differential equations (eq. 1-3), including two chemical (Y and Z) species and heat conduction. A parameter γ controls transition from pure thermal instabilities (γ = 0 ) to mixed thermal-chemical instabilities (γ = 0 ).
In this case, thermal dependence is linked to the conversion Zk3 →B (the exothermic reaction with enthalpy Q3 ≠ 0 ), where
and ∇T = 1/Ta - 1/T Table 1 summarizes the parameter values used to simulate system dynamics according to eq. (1)-(4). In the set of governing differential equations,
To couple several of Sal'nikov’s modified systems, the actual model is discretized using the finite differences method. Then eq. (5)-(7) will look like an equation set defined by eq. (8)-(10), keeping in mind that the value for
Physicochemical quantity | Symbol | Value |
---|---|---|
Volume | V | 1.0 dm3 (cylindrical reactor with height h = 0.15 m and r = 0.046 m) |
Surface area |
|
Five dm2 |
Surrounding temperature |
|
400 K |
Activation energy |
|
166 kJ mol-1 |
Heat capacity |
|
0.029 kJ K-1 mol-1 |
Surface heat transfer coefficient |
|
0.055 kJ s-1 dm-2 K-1 |
Enthalpy of exothermic step | Q 3 | 400 kJ mol-1 |
Initial reactant concentration |
|
Three × 10-3 mol dm-3 |
Rate constants |
|
|
Gas constant | 0.008314 | 0.008314 kJ mol-1 K-1 |
Source: Adapted from Barragán and Ochoa-Botache, 2018.
Here, k trans represents the heat transfer coefficient between two adjacent compartments. Note that coupling was made only in Fourier's heat transfer law, and there is no chemical exchange. Doing k trans = 0, in eq. (10), reduces it to eq. (7). It makes sense when thinking about the selectivity of a cell membrane.
2.2. Configuration entropy tool to characterize the system
As mentioned, the present study aims to model chemical energy supply over time to a macro sarcomere unit. After bringing together enough sarcomeres, a myofibril is confirmed. As a result, the chemical energy input of myofibrils is a series of eq. (8)-(10) takes place. A sample of 400 compartments was arranged in a 1-dimensional linear array, and the initial thermal gradient (slope of the linear gradient, m), thermal pulse applying-time (t app ) as well as pulse amplitude (TΔ) were varied to construct a phase diagram (see variables tested in Table 2). The compartment's initial temperature Ti, where i-th is defined according to the equation of a straight line T i = m i + T max where T max = 300 k and the slope m varies between 0.01 and 0.1 K/i-th compartments (Fig. 2). Pulses were applied each 5 s, starting from 40 s to 70 s and affected only the first half of the compartments ( i = 0 to i = 200). A typical 2D contour plot, like that in Figure 1, is drawn when running simulations with the combinations of variables previously described and summarized in Table 2. (Fig. 1 shows the evolution of the 1-dimensional array on the "Compartment" axis projected in the temporal axis and the "Temperature" axis).
Naturally, covering different initial thermal conditions and applying thermal pulses of different intensities or at other times will be reflected in an extensive collection of plots. Each 2D natural contour plot is resized to obtain helpful information from those pictures until receiving sides of the same length L. Thus, raw images are reduced to 5000 × 5000 pixels and displayed in black and white mode, associating each black pixel with one and a white pixel with a 0. This way, the same image is interpreted as a matrix filled with zeroes or ones. All this pre-processing stage prepares a particular image to be analyzed through configuration entropy as described by Andraud et al. [21,22]. The 2D contour plots generated from each simulation by employing configuration entropy are treated as disordered morphologies and characterized with a descriptor that describes their fractal structure [23]. This method has been used to describe fractal landscapes to clay soil samples and even fat, connective tissue [23-25]. This analysis method qualitatively describes the 1-dimensional array response to different conditions. It provides information about the system's behavior only under the restrictions imposed. An analytical stability analysis from eq could achieve a more general analysis. (8)-(10).
The prepared image is analyzed through a sliding square cell of side l. For a particular position of the cell, the number of pixels with value one contained in the cell is counted. Once the whole image is examined, it is computed Nk(l), the number of cells containing k active pixels, and N(l), the total number of visited cells. The probability associated with the state of the cell is thus defined as p k (l) = N k (l)/N(l) and the corresponding entropy is:
which measures the uncertainty associated with the set of possible states the cell of size l may attain. H(l) is an effective descriptor of the likelihood of the random pixel distribution which may occur inside a cell of size
Pulse applying time |
Slope of the thermal gradient |
Pulse amplitude |
---|---|---|
40 45 50 55 60 65 70 |
0.01 0.05 0.08 0.10 |
317 320 323.9 |
Source: Own elaboration
Detailed configuration entropy values are reported in Fig. 7 and 8 when pulse amplitude was fixed to 317 (Fig. 7) and 323.9 K (Fig. 8). Again, the highest bars are contained in the region described above. This zone qualitatively implies a disordered morphological 2D structure associated with ordered structures. As the emergence of black pixels is tied to the number of oscillations exhibited in time, the zone with the highest entropy of Fig. 5 exhibits more oscillations (see 2D plot in Fig. 4a, where the pulse was applied for 40 s and has a higher H*40 =0.4095) than other regions (2D plot in Fig. 4b, where the pulse was applied at 70 s, with lower H*70 = 0.2746).
3. Results and discussion
3.1. Overall effect of a thermal pulse
In these conditions, oscillatory chemical specie Z could be identified as ATP in a real fiber when restoring chemical energy is prioritized. Set of eq. (8)-(10) provides a compartmentalized energy supply model under such stressful conditions. A characteristic course time of temperature developed in a single 1-dimensional array is presented in Fig. 6. Fig. 8(d-f) presents the oscillatory transient in Z specie, which naturally is coupled to the thermal profile.
At this point, asking for the other compartments that complete the ideal muscle fiber outlook becomes logical. Each of them could be characterized with a plot of temperature vs. time or chemical concentration vs. time. However, the image obtained from the 3D plot of the fiber's time history is processed by the constitutive entropy method (a typical 2D contour plot shown in Fig. 1) instead of providing the evolutionary information of the system separately in this way. As a result, each image is numbered from 0 to 1. The advantage of the configuration entropy tool lies in the possibility of discriminating with quantitative criteria an extensive collection of plots. This information is summarized in three 2D phase diagram plots generated after the maximum peak amplitude was set to 317, 320, and 323.9 K (Fig. 5). It can be seen that the zone, with the highest entropy value, is restricted to pulse applying time ranging from tapp = 40 s to 45 s. In the other coordinate, it is contained between m = 0.08 and m = 0.1k/i-compartment.
In quantitative terms, if the thermal Pulse is applied before inherent oscillations start, the system responds with an oscillatory train (this is the case of Fig. 9(a-b)). The number of changes is proportional to the configuration entropy value in the time-compartments space, as shown in Fig. 3-4. However, if perturbation has no effect at the time of application of the Pulse, its amplitude is not comparable with the system's maximum amplitude. Correspondingly, the vibration of Fig. 9c applied at t app = 70 s falls in a valley, making its relative amplitude not identical with the adjacent thermal spikes. Therefore, it does not significantly change the shape of the single-compartment temperature histories and the 2D contour plots are almost unchanged (see Fig. 7-8. Pulse and control applied at 70 s, i.e., Unperturbed systems exhibit similar entropy values). Therefore, it does not significantly change the shape of the single-compartment temperature histories, and the 2D contour plots are almost unchanged (see Figs. 7-8. Pulse and control applied at 70 s, i.e., Unperturbed systems exhibit similar entropy values). Instead, changes disappear once the Pulse is applied (figure not shown).
Besides the aspects concerning the pulse amplitude and when the system is perturbed, initial thermal conditions tend to increase configuration entropy. Linear thermal gradient slope,
3.2. Discussion
When a muscle fiber is exhausted after high exercise loads, depletion of ATP is expected near the Ca2+-ATPases [26]. Oxidative metabolism or pools of energetically representative chemical species like PCr cannot maintain a continuous ATP feed at high rates and will be exhausted sooner or later. This ATP is essential to keep the Ca2+ ions exchange, which plays a fundamental role in the attaching-detaching cycle of myosin-actin filaments. At first sight, an ATP deficit at this level implies stopping all the machinery related to muscle movement. Nevertheless, being quiet is not as bad since an additional chemical impulse could lead to irreversible damage to the tissue [27]. The inconvenience would not be the ceasing of machinery's activity. Still, the time it takes to recover could be reduced if ATP is restored, and not precisely by the oxidative metabolism pathway.
These ideas were developed in a previous paper [2], where a chemical energy inlet was modeled for a single sarcomere. In such cases, it showed the relevance of chemical oscillations coupled with the local thermal gradient in stressing conditions to replenish levels of a chemical specie of interest, i.e., ATP. The set of equations proposed was based on a modified Sal'nikov model and summarized in the eq. (1)-(3). Through numerical analysis, parameters presented in Table 1 were found to reproduce chemical oscillatory behavior in the presence of thermal instabilities, i.e., doing γ = 1. This condition was the starting point for the present work since the main goal was to evaluate a thermal gradient synergistic effect. Either in the form of a heat pulse, as proposed by [4], or a thermal gradient arising from intense muscular work alongside a muscle fiber, non-isothermal conditions modify system dynamics. They could be an emergency mechanism to overwhelm muscle fatigue locally.
Under energetic stressing conditions [1,28], myocytes from rat exhibit glycolytic and sustained ATP oscillations. If, as proposed in the set of eq. (1)-(3), These oscillations are gradually introduced in γ in autocatalysis along an oscillating thermal gradient. A chemical and thermal explosion is observed when the γ = 1 state is reached. It is the portrait of a single reactor unit.
Thus, a single reactor is an oversimplified system hosting some reactions inside a sarcomere with enough detail to capture the relation between the chemical species that promote oscillatory behaviors. A coupling of several identical reactors, connected through heat exchange only, will lead to the idea of a myofibril chemical energy supplier. In other words, the model of eq. (8)-(10) fuels each
Additionally, the 1-dimensional array of reactors (ideal myofibril) is sensitive to thermal changes. Experimentally, the effect of thermal dependence on muscle function is studied with localized thermal stimulus over muscle tissue or single cells (c.f. [7,29,30]) rather than using single myofibrils. This limitation centers the discussion on the qualitative results shown here rather than a numerical comparison against available data.
The way the compartmentalized model distributes chemical energy when no thermal pulse or initial thermal gradient is present is precisely a near identical repetition of the Z profile and thermal profile from the first cell i = 1 to i=400. This ideal case does not show an appreciable change between cells in time despite including thermal exchange from each compartment I with their surrounds. However, when initial thermal conditions of the compartments are forced to behave linearly, following the equation of a straight line, different slopes m reveal a notorious change, especially near the last portion of the compartments. This difference when changing initial thermal conditions is a consequence of Fourier's heat transfer coupling because this law forces two adjacent rooms to exchange heat at the fastest rate possible [5], and being this process slower or more evident between compartments with different temperatures i = 0 and i = 400.
The system is more perturbed with the highest value for m (see Fig. 2). In effect, Fig. 4a shows that, even when a thermal pulse is applied, if m=0.1 K/i-compartment, the last compartments respond with thermal spikes and, conversely, chemical spikes in Z (Fig. 9(a-d)). In general, phase diagrams from Fig. 5 confirm that this value of m is obtained as the highest value for configuration entropy. Even a system unperturbed by a heat pulse has a minimum configuration entropy associated if a thermal profile is imposed since the beginning. Alternatively, in other words, chemical and thermal oscillations arise only as a consequence of the linear initial thermal conditions for each compartment. For this reason, the height of the bars for the undisturbed (control) system in Fig. 7-8 is comparable with bars obtained after a pulse is applied for 70 s (where the Pulse has no visible effect).
Nevertheless, suppose a homogenous thermal gradient would be sufficient to activate all the compartments, i.e., to increase the amount of Z with the thermal burst. In that case, it is more likely that a natural fiber is under a non-homogenous thermal gradient. The linear gradient proposed here is only an ideal case, illustrating how the compartments 'feel' the thermal burst at different times. First compartments are the first affected by thermal explosions, and to the extent we move away to the last boxes, shots are delayed. For this reason, the time course of temperature corresponding to the first compartments seems to be shifted in time concerning the last rooms. Besides, it implies that the chemical supplier model oscillates in time, but the temperature and, conversely, the concentration of Z, varies among compartments at a given time.
In addition to an initial thermal linear gradient, the system is perturbed with a heat pulse, and the observed response varies in function of Pulse applying time and pulse amplitude. Pulse affected only half of the cells, i.e., two hundred compartments. Pulse applying time is a critical parameter if the pulse amplitude is not higher than the maximum amplitude of the oscillation when it is applied. Therefore, between 40 s and 60 s, a pulse starting at 317 K has a significant impact in terms of configuration entropy because, in this range, the amplitude of the Pulse is greater than the temperature of the system (~315 K) (see Fig. 9(a-b)), promoting an oscillatory behavior just as the Pulse is applied. Something similar occurs with initial thermal conditions since for m=0.5 to m=0.1 K/i-compartment, the initial temperature is over the reference temperature of the system, 310 K (see Fig. 2).
Contrary to non-homogenous initial thermal conditions, a thermal pulse impacts homogenous thermal initial conditions similarly. That means synchronizing the fiber's chemical energy production simultaneously and among thermally stimulated compartments. This thermal synchronization has been recognized as a valuable tool, for example, to pace embryonic hearts [7].
4. Conclusion
The thermal excitability of muscle tissue is essential for monitoring the thermoregulatory function and highlighting the thermosensitive function throughout the body. Under fatigue conditions, where there is an ATP deficit, restoring ATP levels could be a good coupling between a biochemical reaction, i.e., glycolysis, and a thermal gradient. Glycolytic oscillations have been reported under stressing energetic conditions, and thermal gradients can be induced externally.
The 1-D ideal fiber governed by Sal’nikov's modified model shows a chemical-thermokinetic coupling. Their response to different scenarios can be studied by varying three parameters (see Table 2): The slope of the initial thermal gradient (comparable to the internal angle developed naturally in a muscle fiber as a result of intense mechanical work); thermal Pulse applying time (external thermal stimulation to the thread); and pulse amplitude. Numerical simulations deliver a collection of 2D plots by combining these parameters. To understand the information quantitatively, we introduce a configuration entropy analysis where a 2D plot is associated with a single H* number, between 0 and 1. 0 means a homogeneous image (no black pixels, no thermal oscillations) and a highly heterogeneous image (black pixels due to thermal oscillations). The initial thermal conditions following a linear profile along the fiber lead to thermal and chemical oscillations (according to the natural coupling implicit in the Sal'nikov modified model for each compartment). Moreover, applying a thermal pulse could develop further chemical changes in the Z species associated with ATP, which could be a workaround to sustain muscle energetics in stressful conditions during short periods.