From Cycling Between Coupled Reactions to the Cross-Bridge Cycle: Mechanical Power Output as an Integral Part of Energy Metabolism

ATP delivery and its usage are achieved by cycling of respective intermediates through interconnected coupled reactions. At steady state, cycling between coupled reactions always occurs at zero resistance of the whole cycle without dissipation of free energy. The cross-bridge cycle can also be described by a system of coupled reactions: one energising reaction, which energises myosin heads by coupled ATP splitting, and one de-energising reaction, which transduces free energy from myosin heads to coupled actin movement. The whole cycle of myosin heads via cross-bridge formation and dissociation proceeds at zero resistance. Dissipation of free energy from coupled reactions occurs whenever the input potential overcomes the counteracting output potential. In addition, dissipation is produced by uncoupling. This is brought about by a load dependent shortening of the cross-bridge stroke to zero, which allows isometric force generation without mechanical power output. The occurrence of maximal efficiency is caused by uncoupling. Under coupled conditions, Hill’s equation (velocity as a function of load) is fulfilled. In addition, force and shortening velocity both depend on [Ca2+]. Muscular fatigue is triggered when ATP consumption overcomes ATP delivery. As a result, the substrate of the cycle, [MgATP2−], is reduced. This leads to a switch off of cycling and ATP consumption, so that a recovery of [ATP] is possible. In this way a potentially harmful, persistent low energy state of the cell can be avoided.


Introduction
In cellular metabolism, energy transductions are brought about by coupled reactions. The network of energy metabolism is organised in such a way that cycling of respective intermediates, like protons in oxidative phosphorylation (OP), or [Pi], [ADP], and [ATP] in the ATP cycle, is ensured. As was shown previously [1], entropy production during steady state cycling must be zero. This follows from the fact that the line integral taken around a closed path is zero if the integrand is an exact differential. This latter constraint is always fulfilled for potential functions like electro-chemical potentials or affinities. It will be shown that this steady state cycling between coupled reactions is associated with the occurrence of negative resistances.
In a muscle fiber, mechanical power output is coupled to ATP splitting. How this is achieved is not fully understood, although there has been great success in many field endeavours in muscular research such as the structure of the contractile apparatus and its functional correlates [2][3][4][5][6]. Since Huxley's widely accepted sliding filament theory [7,8], the cross-bridge cycle is of central importance, especially in the functional aspects of contraction. This cycle must contain the reactions of free energy transduction from chemical (ATP splitting reaction) to mechanical energy (actin movement against a load force). From the overall reaction, contractile efficiency can be obtained by relating mechanical power output to the dissipation function of ATP splitting, where the mechanical power is given by the product of the force exerted by the load and the shortening velocity. Experimental results show [9][10][11][12], that when efficiency is expressed as a function of v, a curved line with a maximum is obtained. From non-equilibrium thermodynamics (NET, [13]), it is well known that uncoupling is necessary to generate a maximum in efficiency plots (efficiency against reduced force ratio). Thus, to yield such a maximal efficiency, any description of the cross-bridge cycle on a thermodynamic basis must contain an uncoupling mechanism, which uncouples the transduction of free energy from ATP splitting to actin movement.
To describe the cross-bridge cycle in terms of the new flux equations published recently [1], the cross-bridge cycle has to be formulated in relation to this formalism, which combines the basics of NET [13][14][15][16] with Michaelis-Menten-like kinetics of enzyme-catalysed reactions [17]. It will be shown that Hill's equation describing muscular performance [18,19] can be easily deduced by applying the new flux equation.
When compared with other approaches to the energetics of the cross-bridge cycle, the main particularity of the present work may be the fact that this cycle is connected here to energy metabolism of the muscle fiber, i.e., to ATP producing and consuming reactions. The generation of mechanical energy from the free energy of ATP splitting is treated here as one of the parallel reactions of the sarcosol consuming ATP delivered in fast fibers, mainly from glycogenolysis or glycolysis, respectively. This integration into the cell's energy metabolism makes it possible to inspect some variables like ATP and its reaction products and species at high mechanical power output. In addition, concentration changes in metabolites and ions like creatine phosphate, lactate, H + , and Mg 2+ , are of interest under these conditions. This is achieved by formulating, in particular, the ATP splitting reaction according to Alberty [20] as a function of both [H + ] and [Mg 2+ ].
It is the aim of this study to elucidate cycling between coupled reactions, and if such cycling is also involved with the cross-bridge cycle and force generation. In addition, the consequences of uncoupling on power output and efficiency will be shown. For this, a formulation of the cycle in terms of the above mentioned new flux equation had to be derived.
The phenomenon of muscular fatigue at the cellular level occurs when ATP consumption exceeds ATP delivery [21][22][23][24][25]. Under such conditions drastic changes in many metabolite and ion concentrations can be expected. The results of simulations will show to what extent these changes may contribute to fatigue, and if this phenomenon can be explained by such changes alone.

How Negative Conductances Are Generated
In a previous article [1] it was shown that at steady state all affinities and dissipation functions of closed pathways associated with coupled in series reactions must vanish. In the following, it will be demonstrated that the overall resistance (=1/conductance) of such cycles must also be zero. As a consequence, the existence of negative conductances (or resistances) has to be called for.
According to [1] a coupled two-flux-system can be described as: , and (1) because A 1 usually is negative, 1  must also be negative. Expanding the right hand terms yields: , and (2d) The term: represents that partial conductance of L c , which is associated with A 1 , while belongs to A 2 . They relate to the usual different forms of energy being processed through the coupling reaction. Obviously, when A 1 is negative, L c1 must also be negative to yield a negative 1  .
The same result can also be derived by starting from flux equations, yielding: , and So, to yield a positive J 1 , L c1 has to be negative for a negative A 1 . A 1 and A 2 are in series, hence, L c can be regarded as the equivalent conductance of both in series conductances L c1 and L c2 , yielding: These theoretical results are confirmed by simulations.

Conductances in Cycles between Coupled Reactions
In a reaction sequence in which two coupled reactions in series are involved, the output force of the A negative). Because both reactions are coupled, the conductance (resistance) of the whole cycling process is brought about by both partial (in series) conductances associated with the input and load affinity, respectively. In oxidative phosphorylation (OP), as described in detail in [1], a proton cycle is generated over the inner mitochondrial membrane. At steady state, coupled outward proton pumping by redox (NAD red and FAD red ) reactions of the respiratory chain (J NA and J FA ) equals the back flow of a given fraction formulated. Opposite equality of partial conductances is also fulfilled for this cycle (the above results were obtained by using the simulation SIM GlOx from reference [1]). For a further illustration, an analytically solvable example is given in the Appendix section. Simple electric circuits consisting of one battery connected to an outer conductance, or of two batteries in series, are analysed. These examples show very clearly the behavior of coupled in series reactions.
Further evidence of such an equality of conductances comes from the known fact that for a coupled reaction with an attached load, conductance matching (L Ld = L c ) is needed to achieve a maximal power output [1]. At total coupling, the output power is given by: (6) The maximal P out is found by differentiation with respect to the variable A 1 , while A 2 remains constant, and by setting the derivative equal to zero: In addition to cycling at the inner mitochondrial membrane, other types of cycles occur in metabolism. Especially in skeletal muscle cells, the phosphofructokinase (PFK) reaction in conjunction with the fructose-1,6-biphosphatase (FBPase) operating anti-parallel represent a substrate cycle, which may control the pathway of glycolysis (GLY) more sensitively than would be possible by PFK alone.
In this cycle, fructose-1,6-biphosphate (FBP), which is produced by ATP-coupled formation from fructose-6-phosphate (F6P), is cycled back via FBPase to F6P. However, usually both fluxes are not equal. Also to demonstrate the opposite equality of partial conductances for this kind of cycle, only equal fluxes can be used for this purpose.
As a further example, the phosphocreatine shuttle will be considered. The creatine kinase (CK) reaction can also be regarded as a coupled reaction. Here, ATP splitting powers phosphocreatine (PCr) formation from creatine (Cr), which may proceed near equilibrium. As described in detail in reference [1], ATP is shuttled between locations of ATP formation (for instance in the inter-membrane space in mitochondria) and locations of high ATP demand like myofibrils. By analogy to an electric circuit built by two in series batteries with an outer circuit conductance (see Appendix (A4)), the output affinity of PCr formation in the inter-membrane space of mitochondria corresponds to To ensure diffusional flow of PCr and Cr between both locations, an additional driving force (corresponding to U e ; see (A4)) with associated conductance must be present. Under such conditions partial conductances do not match. Only when the additional conductance corresponding to the diffusional process (L e ) is added to I c1 L does this sum become opposite and equal to II c2 L , as is shown in (A4). L e depends greatly on structural features. So, to achieve a high diffusional conductance, diffusional paths must be as short as possible, which in turn requires a high grade of structural organization [26][27][28].
It seems worth mentioning that coupled systems like pump and leak cycles are often not in a steady state. For instance, steady state cycling through sarco/endoplasmatic reticulum Ca 2+ ATPases (SERCA) and Ca 2+ release channels of the sarcoplasmatic reticulum (SR) breaks off during activation of contraction. There is an enormous Ca 2+ efflux through release channels; meanwhile the pumping rate of SERCAs may be low. Under these conditions, respective conductances may greatly differ; however, when a new steady state cycling is reached, the partial conductance of SERCA must be opposite and equal to the conductance of the Ca 2+ release channels. The opposite has to be expected, when release channels close again, and the Ca 2+ pumping rate exceeds the release rate.

From Chemical Potentials to Mechanical Force Generation
In striated muscle cells like ventricular muscle cells (VMs) or skeletal muscle fibers (SMFs), force generation as well as shortening is brought about by the cyclic action of cross bridges. It is a known fact that this process is powered by ATP splitting. The underlying mechanism of the energy transduction process, however, is not completely understood. Here, a thermodynamic description of the cycle is derived using a formalism recently published [1]. It takes into account the basic energetics of enzyme-catalysed reactions, which states that the overall affinity of the catalysed and non-catalysed processes must be equal. For an enzyme-catalysed reaction like: (S = substrate, E = enzyme, ES = enzyme-substrate complex, P = product), this means that at steady state the sum of the affinities of substrate binding, transition, and product release must yield the affinity of the non-catalysed reaction, which is given by the reaction affinity of all involved compounds in the bulk solution: (7a) or, after contraction of the first two terms: and R K are equilibrium constants of the binding, transition, and release reaction, respectively, whereas r K  denotes that of the non-catalysed reaction).
An analogous reaction sequence is used here to describe the cross-bridge cycle. The following cycle is given in chemical notation, i.e., the charges of involved species are taken into account. The cycle begins with the splitting reaction of the de-energised actomyosin complex (A-M) by MgATP 2in the diffusional space of myofibrils: This first reaction yields dissociated actomyosin with MgATP 2bound to myosin (the bold point denotes binding to myosin). Two negative charges develop on the dissociated actin, which are neutralised by potassium ions, K  , stemming from free MgATP 2− , which is now bound to myosin heads. On the dissociated myosin heads, it neutralises both emerging positive charges. This first actomyosin dissociation and binding of MgATP 2− to myosin is followed by ATP splitting on the myosin heads. This transition reaction is described by It is coupled to the formation of energised myosin ( 2+ M  ), which is characterised by a tilting of the myosin head from a more bent arms position by an angle of about 60° towards the respective Z disc, so that now the myosin head builds a right angle with the opposing actin filament. 2+ M  contains free energy from reaction R 2 as conformational energy. The force generating stroke of the myosin head is triggered by the association reaction to form the energised actomyosin complex (cross-bridge): Because of uncompensated charges, the resulting intermediate in curly brackets lacks firmness. A more stable conformation is obtained by the last reaction of the cycle, which restores the electroneutrality of the cross-bridges: The existence of a less stable interaction of myosin with actin has been shown previously [29][30][31].
In a coupled reaction, de-energised actomyosin is restored by dissociating first MgADP  and 2 4 H PO  from cross-bridges and then by releasing the stored conformational energy. During this reaction the cross-bridge tilts back by 60° towards the sarcomere centre, whereby free energy is transferred to the actin filament as mechanical energy. If the constraint which has to be overcome during stroking. The quantity Because ionic species are involved, the reaction sequence of the cycle should be markedly enforced by electrostatic interactions. So MgATP 2− binding can proceed only if actomyosin dissociates, whereas release of products becomes possible only when at the same time cross-bridge formation occurs. Moreover, the conformational change in the myosin head forces it into a new position, which favours an interaction with actin at a new actin binding site displaced a certain distance towards the Z-disc. During stroking, binding of a new MgATP 2− molecule and detaching of cross-bridges may preferentially occur at the end of the power stroke, when cross-bridges form an angle of about 60° with the actin filament (see below for uncoupling by stroke shortening).
The contractile performance of whole muscle and of SMFs is exceptionally well reproduced by Hill's equation [19]. This equation relates the shortening velocity v to the mechanical load force Ld F which has to be overcome during shortening. (9a) The above function represents a hyperbola, which fits remarkably well with experimental data obtained under isotonic conditions.
To obtain an equivalent expression from the flux equation Km l N  (see below for a mechanistic interpretation of uncoupling and  values). Conversion to mechanical units can then be done in the same way as shown above.
In Figure 1  , where the coupled flux must be zero and only uncoupled fluxes are possible.
In the following, an attempt has been made to interpret the above results, which were gained from a phenomenological approach, mechanistically by relating coupled and uncoupled fluxes to possible cross-bridge actions.
, coupled reactions with associated actin filament movement come to a halt, because the driving force has vanished. As already mentioned above, now only uncoupled fluxes can occur. Such a situation may also be realised with isometric contraction, which is known to be associated with ATP splitting and heat production, but without power output. That is, a mechanism has to be found which explains the identity of the isometric force F 0 with F P , which was merely formally derived from the input affinity P Str A by a conversion factor. This is achieved by defining the uncoupling mechanism by a shortening of the stroke length Str l of the power stroke. Under totally coupled conditions, the input flux is given by: Uncoupling by stroke shortening dissipates free energy, which can be expressed by a leak dissipation function: The leak conductance P StrL L can be replaced by Str L , because this latter conductance may depend mainly on the formation mechanism of the actomyosin bond. The stroke reaction associated with conformational changes of the myosin head is assumed to proceed at a high conductance, since the energising reaction ( En J ), which is coupled to the same conformational change in the reverse direction, also proceeds at a very high conductance. So an increase by stroke shortening of a high conductance The input flux then is given by: The output flux is reduced by stroke shortening as if it were uncoupled.    , Ld Str J formally could be negative, which would mean that actin filaments were moving in the direction of stretching. This is, however, impossible, because actomyosin bonds would have to be broken by a load force, which is smaller than F 0 . Therefore, in this region of loads, Ld Str J cannot be negative; it must remain zero.

Power Output and Efficiency
In experiments, mechanical power output is often represented in relation to shortening velocity. In also, En1 Str2 L L   is fulfilled, and therefore, cross-bridge cycling at zero resistance.
In addition, the equality of , which describes the conductance of the whole cycle including coupled inputs and outputs is nearly exactly obeyed. The overall efficiency of the cross-bridge cycle is obeyed: (17) as is the overall dissipation function given by:   On the one hand, the driving force is changed by the load potential (see Figure 1,     It seems plausible to suggest that during shortening it is not always the same group of cross-bridges that is active, but that, e.g., at 1.08 µM [Ca 2+ ], four different groups may alternately be involved with contraction. The cycling frequency of an individual cross-bridge would then be much lower than the frequency of ATP splitting, which might be advantageous, especially at high velocities. Furthermore, an alternating involvement of groups may be absolutely necessary for a smooth shortening. How this might be accomplished is so far not known. An involvement of special filaments of the sarcomere [CB] [µM] cytoskeleton [34][35][36], which may be responsible for a subtle sensing of load forces and an undisturbed takeover of a given load by a new fraction of cross-bridges during synchronous stroking, seems indispensable. The values of maximal tension (=force/unit area in N/m 2 = Pa, Pa = Pasqual) obtained from SIM GLYgen (A16) in the present study are comparable to experimental values. For instance, a value of 372 kPa (from F 0 = 7.756 × 10 −4 N, [Ca 2+ ] = 1.06 µM, 37°C) found here, seems to be in reasonable agreement with about 320 kPa resulting from measurements with a fast-twitch mouse fiber at 25°C [37].  [12].
All these parameters of contractile performance may, however, be reduced to a certain extent by dissipative frictional processes associated with v which are not addressed in the present simulation. Such dissipation during fiber shortening may be produced mainly by viscous deformations of membranes and the filament lattice. Figure 5A shows the time courses of rates of [H + ] changes. Interestingly, [H + ] production by ATP splitting is practically compensated by [H + ] consumption by the CK reaction. The contribution by the AK reaction is negligible. A similar behavior is found for Mg 2+ ( Figure 5B). A concentration increase in this ion is mainly brought about by acidification. A second source of protons is given by the disturbance of lactate production by glycogenolysis (or glycolysis) and lactate efflux via lactate/H symport at the sarcolemma. Especially when lactate and H + accumulate in the glycocalyx (the outer aspect of the sarcolemma), the concentrations of these compounds also increase drastically in the sarcosol. This seems to be the main mechanism of sarcosolic acidification. Muscular fatigue at the cellular level can be defined as a phase of markedly reduced contractile performance, which largely recovers after a period of rest [38]. Because metabolites like creatine, ADP, Pi, H + , and lactate accumulate during conditions of fatigue in a similar way as can be observed during ischemia or hypoxia, which are known to be the result of impaired ATP production, it seems justified to suggest that the preconditioning for fatigue may also be initiated by a deterioration of the energy metabolism of the muscle fibers. Whenever ATP delivery does not match ATP consumption, such a situation may arise.
These effects can be easily demonstrated with a simulation of glycogenolytic or glycolytic ATP production in the absence of mitochondrial metabolism (SIM GLYgen, see (A16)), which is related to the energy metabolism of fast muscle fibers. However, when a back pressure on glycogenolysis (or glycolysis) is produced by accumulated extracellular [Lac] e and [H + ] e , the flux through this pathway may become reduced. In addition, efficiency has been reduced by switching from glycogenolysis to glycolysis. The power output of ATP production is markedly reduced by these combined effects. As a result, the power of ATP production begins to fall, so that ATP consumption may overcome ATP production. Steady state cycling through ATP consuming and producing pathways can now no longer be maintained. Figure 6 shows   Under these conditions all myosin heads form cross-bridges, which however are unable to perform the power stroke, since the input force is equal to the opposed load force. In such a situation myosin heads may be bound to actin and may have dissociated H 2 PO 4 − and MgADP − similar to an isometric contraction, but in contrast to those latter conditions, equilibrium of forces is now brought about at a much lower load force ( P Str A = − Ld Str A = 0.375 × 10 4 J/mol at 1.08 µM [Ca 2+ ]). A load-dependent actomyosin splitting by MgATP 2− at the beginning of the stroke, that is uncoupling, is impossible under these conditions. So cross-bridge cycling with concomitant ATP consumption may be completely prevented.
[ATP], therefore, can recover rapidly, even if the conditions leading to fatigue first remain unchanged. By this mechanism the fatigued skeletal muscle fiber is capable of protecting itself from the dangerous risk of irreversible cell damage. This seems to be necessary, since this cell type is voluntarily controlled without any protecting mechanism against an unbridled consumption of ATP, as is known to occur with other ATP coupled reactions such as, for instance, ion pumps. These are controlled mainly by the ion concentration that they are transporting. For example, when [Na + ] in the sarcosol is lowered by the Na/K pump to values below 10.0 mM, the reaction rate of this transport process is increasingly deactivated by the decreasing [Na + ], so that ATP consumption also is reduced.
Such a protective mechanism is not known, however, for the cross-bridge cycle. Contractions with concomitant ATP splitting would be incessantly initiated, as long as firing of nervous impulses persisted. The voluntary muscle fiber would obey this parent command up to exhaustion or even up to cell death, if the fatigue producing mechanism were absent. Obviously it represents that special control mechanism which is necessary to protect voluntary muscle from dangerous ATP depletion during phases of high energetic demands. The above results are obtained from a simulation, in which ATP production is confined solely to glycolysis. A whole muscle, however, is constructed from many types of functionally different fibers, with slow fibers having densely packed mitochondria, and fast fibers in which mitochondrial density can be very low. It is a known fact that especially fast fibers with a very low content of mitochondria and, therefore, a mainly anaerobic ATP production, are much more liable to be affected by overpowering than slow fibers. This may be brought about primarily by the preconditioning effects especially associated with this fiber type. At very high energetic demands, fast fibers produce much more lactate and protons through glycolysis than slow fibers, which can oxidise pyruvate by mitochondria. That is, the glycolytic ATP production rate of fast fibers may be decelerated by a back pressure, which may be generated by accumulating lactate and protons during high power output. In slow fibers with a high rate of oxidative glucose metabolism, such a back pressure cannot be produced as easily. Therefore, the metabolic changes leading to fatigue are simulated here with respect to fast fibers with a low resistance to fatigue. This weakness may be best demonstrated with an extreme fiber type, which can produce ATP solely by GLY. However, fibers completely devoid of mitochondria may not exist in vertebrate muscle. The results of this fatigue model, therefore, can only be taken as an approximation of real fast muscle fibers. Muscular contraction on the basis of NET has been treated theoretically by Caplan and Essig [13]. These authors explained the curvature of Hill's equation by an uncoupling. They defined the degree of coupling by 0.2 and 0.3). This latter equality is also fulfilled by the present model (  = 0.309). However, in contrast to the approach of the above authors, the hyperbolic form of Hill's equation is produced here by introducing a Michaelis-Menten-like inhibition factor into the respective conductance, as already mentioned above. Moreover, uncoupling in the present model is produced through a load-dependent stroke shortening, which generates the maximum obtained with power plots.

Methods
Here the energy metabolism of a fast twitch muscle fiber is treated. That is, ATP production by this fiber type is solely brought about by metabolism of glycogen and/or glucose. Mitochondria are absent. Glycolytically produced [NAD red ] has to be reoxidised by the lactate dehydrogenase (LDH) reaction, and the lactate plus proton formed thereby is released to the extracellular space via Lac/H symport. Electrophysiological reactions at the cell membrane (sarcolemma) are omitted. Also, most reactions of the sarcoplasmatic reticulum (SR) are not addressed. Only Ca 2+ pumping by the sarco/endoplasmatic reticulum Ca 2+ ATPase (SERCA) as an ATP consuming reaction is included in simulations besides several other reactions of ADP production (see SIM GLYgen (A16)) taken over from reference [1]. Therefore, [Ca 2+ ] is treated as an adjustable constant.
To determine the fractional fiber volume V Cell , a cylindrical geometry of the muscle cell is assumed. With radius R Cell = 25.76 µm, and a length L = 10 3 µm (fraction of whole fiber length), V Cell = 2.0847 × 10 6 µm 3 or 2.0847 nL, and A Cell = 2.0847 × 10 3 µm 2 . From data of Aliev et al. [39] for the heart, the volume of the sarcosol, V c , can be determined by adding the mitochondrial to the fibrillar volume, yielding V c /V Cell = (321 + 195 + 55)/758.5 = 0.7528 or V c ≈ 1.57 nL. Then c  can be obtained using c  = 10 −12 /(F×V c ) = 6.6024×10 −9 µM/C (F = Faraday's constant, C = Coulomb). That is, to yield the corresponding flux in µM/ms from an electric current entering the sarcosol, this current in fA (= pS×mV = 10 −18 C/ms; pS = pico Siemens) has to be multiplied by c  .
For calculation of force and velocity, the dimensions of the force generating cross-sectional area and the number of half-sarcomeres (HS) of the fibrils must be known. The contractile machinery of skeletal muscle fibers is organised in fibrils having diameters between 1.0 and 2.0 µm, which are built up from in series sarcomeres connected by Z-discs over the whole length of a fibril, i.e., from end to end of the fiber. The functional unit is given by the HS. The principal filaments of an HS are the thick myosin and the thin actin filaments. In cross-sections, myosin filaments show hexagonal geometry. From this symmetry the fibrillar volume V Fibr can be obtained. One hexagon is composed of six equilateral triangles of side length l Tri = 41.0 nm [12] and equal angles of 60°. The area of a hexagonal fibril (or HS) of radius R Fibr = 18.0 × 41.0 = 738 nm is given by: The total volume of fibrils is given by 0.866×V Cell (see reference [39]). The number N Fibr is then given by: For the determination of the total number of myosin heads of an HS, the number of myosin filaments of an HS must be known. The above hexagonal area of an HS can be constructed from equilateral triangles of l Tri = 41.0 nm. A large triangle of the hexagonal HS with n-fold side length contains: F CB is the force of one single cross-bridge. It is obtained here from the stroke input potential A P (see Results) under these conditions, by assuming that the stroke length l Str = 12.0 nm. Then F CB is given by:

Conclusions
Cycling between coupled reactions occurs, especially in energy metabolism. It is shown that the overall resistance of such cycles must vanish, and that the resistance or conductance associated with the negative output affinity of a coupled reaction also has to be negative. The following may be illuminating: The entropy change of a spontaneously proceeding reaction ( r i S  ) is always positive.
When a reaction is forced against spontaneity, r i S  must become negative. All reaction parameters associated with entropy changes, like affinities and conductances, must inevitably follow a sign change of r i S  whenever such a change occurs. This is not a contradiction to Ohm's law, but a consequence of the phenomenological definition of a conductance through L J A    . It can be concluded that the occurrence of negative conductances is realised not only with biochemical reactions in living cells, but represents a fundamental concept of coupled processes. This concept is realised here for the cross-bridge cycle. The reactions of the cycle are described on a thermodynamic basis using the kinetic approach of enzyme-catalysed reactions. Hill's equation for muscular performance can be derived on this basis. However, uncoupling has to be introduced to yield a maximal efficiency of power output. Here the uncoupling mechanism is not an accidental process during energy transduction, but a necessary interference during force generation, which ultimately produces an isometric contraction.
Although mechanical acceleration may also be possible on a cellular basis by changes in sarcosolic [Ca 2+ ], it seems highly unlikely, however, that this may be sufficient to allow normal locomotion of a subject. Only the control by the nervous system can bring about coordinated actions of several muscle fibers, groups of fibers, or even several different muscles. In this way, accelerated and decelerated motion becomes possible. To achieve this, the number of force generating cross-bridges is varied by a change in cross-sectional area, that is, by altering the number of fibers recruited for contraction. Thereby the locomotion at high efficiency or maximal power output can be controlled by will. Also, isometric contractions are indispensible for coordinated actions. They are produced by reducing the cross-sectional area to such an extent that a load dependent uncoupling is initiated to stop fiber shortening. In many species nervous control of muscles is not a capability which is present from birth on. To reach a certain level of adroitness an individual has to learn-often during a long lasting phase of exercise-to control muscle action by coordination.

Negative Resistances in Simple Electric Circuits
Reactions occurring in a common car battery can be considered as coupled. A redox reaction is started by introducing a catalyst (the electrodes), which couples the affinity A R of the redox reaction to the formation of an electrical potential difference  L c1 and L c2 are given by -5/2 and 5/3 Ω −1 , respectively, or R i1 = −0.4, and R i2 = 0.6 Ω, which fulfils R i1 + R i2 = R i = 0.2 Ω. Setting R e = 4.0 Ω yields R i1 = −4.0, R i2 = 4.2 Ω, and again R i = 0.2 Ω. Moreover, in the electric circuit the overall resistance R i1 + R e vanishes. It should be noticed that partial resistances are not constant, although R i is a constant. They depend both on R e and L e , respectively.
The total dissipation function of reactions in the battery and of the outer circuit is given by: (240 J/s or 0.24 kW), or (A3a) This latter result means that the total entropy production of the circuit is given by that of the redox reaction. This is also valid with respect to heat production.
Power output as a function of R e is given by: In simulations, instead of complete d[H + ]/dt, only those fluxes producing or consuming protons are considered, because changes of [H + ] depend mainly on these fluxes (see Figure 5A).
[Mg 2+ ] is introduced as a variable only in those simulations that deal with muscular fatigue. Because changes of [Mg 2+ ] depend mainly on acidification, and pH does not change markedly even under conditions of high power output, this variable is set constant to 800 µM for all other simulations.
In the above equations, methods of calculus are used so formulas can be held compact. In simulations, however, these equations must be incorporated in an explicit form, which often results in very voluminous expressions.

Simulation of Glycogenolysis and Glycolysis
Most flux equations of glycogenolysis are congruent with those of a simulation of glycolysis given in [1]; they are taken over from that article. Glucose-6-phosphate (G6P) formation by hexokinase (HK) and glycogen phosphorylase is now included. The new flux equations used here are as follows.