CC BY
26
10. Ivanova A.V., Bobkina E.M., Il'ina V.N. K flore pamjatnika prirody «Gora Krasnaja» Krasnojarskogo rajona Samarskoj oblasti // Samarskaja Luka: problemy regional'noj i global'noj jekologii. 2011. T. 20. № 3. S. 88-105.
11. Il'ina N.S., Il'ina V.N., Mitroshenkova A.E. Prirodnyj kompleks «Verhov'ja reki Binaradki»: sovremennoe sostojanie i ohrana (Krasnojarskij rajon, Samarskaja oblast') // Vestnik Volzhskogo universiteta im. V.N. Tatishheva. 2011. Vyp. 12. S. 35-41.
12. Ustinova A.A., Matveev V.I., Il'ina N.S., Solov'eva V.V., Mitroshenkova A.E., Rodionova G.N., Shishova T.K., Il'ina V.N. Ohranjaemye prirodnye territorii Samarskoj oblasti: vydelenie, monitoring, rastitel'nyj pokrov // Izvestija Samarskogo nauchnogo centra Rossijskoj akademii nauk. 2011. T. 13. № 6. S. 1523-1528.
13. Il'ina V.N. Jekologicheskaja plastichnost' vidov flory urochishha «Verhov'ja reki Binaradki» // Reproduktivnaja biologija, geografija i jekologija rastenij i soobshhestv Srednego Povolzh'ja: Materialy Vserossijskoj konferencii (27-29 nojabrja 2012 g.). Ul'janovsk: UlGPU, 2012. S. 107-109.
14. Il'ina V.N. Jekologicheskaja plastichnost' flory urochishha «Lysaja gora» (Studenyj ovrag, Krasnoglinskij rajon g. o. Samara) // Nauchnyj dialog. 2013. № 3 (15). S. 43-56.
15. Il'ina V.N., Mitroshenkova A.E., Ustinova A.A. Organizacija i monitoring osobo ohranjaemyh prirodnyh territorij v Samarskoj oblasti // Samarskii nauchnyj vestnik. 2013. № 3 (4). S. 41-44.
16. Mitroshenkova A.E., Il'ina V.N., Ustinova A.A. Prirodnyj kompleks «Igonev dol»: sovremennoe sostojanie i ohrana (Kinel'skij rajon, Samarskaja oblast') // Izvestija Samarskogo nauchnogo centra Rossijskoj akademii nauk. 2013. T. 15. № 3-2. S. 852-855.
17. Il'ina V.N. Opredelenie prirodoohrannogo statusa redkih vidov rastenij Krasnoj knigi Samarskoj oblasti (vtoroe izdanie) na osnove osobennostej ih ontogeneza i populjacionnoj struktury // Fitoraznoobrazie Vostochnoj Evropy. 2014. T. VIII. № 4. S. 98-113.
18. Il'ina V.N. Osobennosti struktury cenoticheskih populjacij ostrolodochnika kolosistogo Oxytropis spicata (Pall.) O. et B. Fedtsch. (Fabaceae) v Samarskoj oblasti // Izvestija Samarskogo nauchnogo centra Rossijskoj akademii nauk. 2014. T. 16. № 5-5. S. 1637-1643.
19. Il'ina V.N., Mitroshenkova A.E. Sohranenie fitoraznoobrazija na osobo ohranjaemyh prirodnyh territorijah Samarskoj oblasti // Problemy sovremennoj biologii. 2014. № XII. S. 20-26.
20. Il'ina V.N., Mitroshenkova A.E. Rol' pamjatnikov prirody regional'nogo znachenija v sohranenii fitoraznoobrazija v Samarskoj oblasti // Izvestija Samarskogo nauchnogo centra Rossijskoj akademii nauk. 2014. T. 16. № 1-4. S. 1205-1208.
21. Il'ina V.N., Savchenko A.A. Soderzhanie razlichnyh veshhestv v pochvah i po-verhnostnyh vodah na territorii nekotoryh pamjatnikov prirody Kinel'-Cherkasskogo rajona Samarskoj oblasti // Karel'skij nauchnyj zhurnal. 2014. № 1(6). S. 119-121.
22. Mitroshenkova A.E., Il'ina V.N. Botanicheskoe kraevedenie Samarskoj oblasti: aktual'nye problemy i perspektivy razvitija // Samarskij nauchnyj vestnik. 2014. № 2 (7). C. 71-74.
23. Il'ina V.N. Vedenie Krasnoj knigi Samarskoj oblasti: k opredeleniju pri-rodoohrannogo statusa redkih vidov rastenij // Strukturno-funkcional'naja organizacija i dinamika rastitel'nogo pokrova: materialy II Vserossijskoj nauchno-prakticheskoj konferencii s mezhdunarodnym uchastiem, posv. 80-letiju so dnja rozhdenija d.b.n., prof. V.I. Matveeva, 30-31 janvarja 2015 goda, Samara. Samara: PGSGA, 2015. S. 131-137.
24. Aladinskaja A.R., Anopchenko T.Ju., Afonina I.A., Ahmedenov K.M., Domashenko Ju.E., Drogobuzhskaja S.V., Ivanova T.K., Il'ina V.N., Karaeva Ju.V., Kirsanov S.A., Korostiev A.V., Koshim A.G., Kravcova M.V., Krapchin I.P., Kremeneckaja I.P., Kucherov V.S., Lashhuk V.V., Mitroshenkova A.E., Murzin A.D., Murzina S.M., Chernyshev M.A. Ohrana okruzhajushhej sredy ot negativnogo vozdejstvija hozjajstvennoj dejatel'nosti: nauchnaja monografija; pod red. D.V. Eliseeva. Novosibirsk, 2015. 260 s.
Кокшенев В.Б.
Кандидат физико-математических наук, профессор физического факультета Федерального Университета Минас Жерайс,
г. Бело Горизонти, Штат Минас Жерайс, Бразилия
НОВЫЙ ПОДХОД В ТЕОРИИ ПОПЕРЕЧНЫХ МОСТИКОВ СТАЦИОНАРНО СОКРАЩАЮЩИХСЯ
И РАСТЯГИВАЮЩИХСЯ МЫШЦ.
Аннотация
Несмотря на достаточно широкое применение в биомеханике, модель скользящих нитей, предложенная Хаксли, не позволяет описать поведение активированных скелетных мышц в условиях стационарного растяжения. В настоящей статье впервые предложено новое самосогласованное решение фундаментального кинетического уравнения Хаксли для распределения актино-миозиновых ниточных связей (тянуще-тормозящих поперечных мостиков) в рамках обобщенной термо- динамической теории флуктуаций. Тем самым модифицирован традиционный подход, представляющий собой упрощенный расчет частот образования и разрыва мостиков на основе теории переходных состояний. Предложено обобщенное выражение для мышечного усилия, пригодное как для состояния сжатия, так и для растяжения мышцы. Единый механизм сокращения мышцы интерпретируется взаимодействием между головками миозина, связанными c нитями актина, находящимися в термодинамически равновесных состояниях, и вращающимися неуравновновешенными головками миозина, находящимися в ассиметрично распределенных, локально равновесных состояниях. Переход от стационарного режима активного сокращения к активному удлинению мышцы объясняется реконструкцией поперечных мостиков и изменением частоты их образования и разрыва при сохранении основных упругих характеристик мостиков. Предлагается обобщенное уравнение состояния «сила-скорость» с использованием лишь одного подгоночного параметра в каждом из режимов генерации силы, в то время как в подходе Хаксли (1957) и его последующих модификациях требуется четыре параметра для описания экпериментальных данных для каждого режима актуации.
Ключевые слова: мышечное сокращение, уравнения состояния, cтатистическая механика, теория скользящих нитей, соотношения «сила-скорость».
Kokshenev V.B.
PhD, Professor, Physics Department of the Federal University of Minas Gerais, Belo Horizonte, Brazil,
NEW INSIGHTS INTO THE THEORY OF CROSS-BRIDGE TENSION OF STEADILY SHORTENING
AND STRETCHING MUSCLES
Abstract
Despite the great success, seminal Huxley's sliding filament model broadly fails to explain the steady state behavior of stretched activated skeletal muscles. Here, a new self-consistent solution to the fundamental kinetic equation for the distribution of actin-myosin linkages (cross-bridges) is proposed in light of the generalized thermodynamic theory of fluctuations, thus substituting the transition state theory traditionally used for the transition rates. The unified description for mechanism of the force output in both shortening and stretching regimes is attributed to the interplay between the uniformly-state distributed, thermodynamically equilibrated myosin heads attached to actin filaments and asymmetrically-state distributed, mechanically equilibrated rotating myosin heads. The crossover between two steady regimes is associated with a reconstruction of the cross-bridge domains and change their attach-detach rates, leaving unchanged basic mechanical characteristics. Theory suggests a unified generic force-velocity equation, using only a single combination of the cross-bridge parameters for each regime, while Huxley's approach (1957) as well as its subsequent modifications, requires four adjustable parameters to fit the same data.
Keywords: muscle contraction, equations of state, statistical mechanics, sliding filament theory, force-velocity relations.
27
Introduction
The early studies of muscle fibers under the light microscope revealed cross-striations running normal to the fiber axis during concentric and eccentric muscle contractions and two researcher groups [13, 14] laid the foundations for the sliding filament theory. They simultaneously suggested that muscle contractions occur due to the relative sliding of the thick myosin filaments past the thin rigid actin filaments, mediated by the ATP-dependent elastic actin-myosin linkages (cross-bridges, CBs) [13], working as independent force generators [14].
Based on the idea of the existence of mechanical equilibrium of elastic myosin molecules at the regular sites of relatively rigid actin filaments, A.F. Huxley explained the mechanism of generation of elastic muscle force via strain variation of myosin heads bound to actin molecules [15], resulting in lever-arm picture [16]. However, in spite of the great success in illuminating the force generation and power liberation during muscle shortening [15], Huxley's seminal approach generally failed to describe the ascending branch of the force-velocity curve established phenomenologically for the steady muscle lengthening. Considerable efforts have been therefore made to find an exact numerical solution to the partial differential equation [15] via ad hoc modifications of the rate functions [3, 9, 36, 37], including an arbitrary rate-velocity dependence provided the "exact" steady state solutions [3]. All such modifications explored the fact that the sliding filament model [15] leaves a free choice of the load-dependent attachment and detachment rate functions with respect of CB strains. Many researchers simulated numerically a range of muscle properties by fitting empirical data on the muscle tension in lengthening by adopting either linear-strain functions or constants [21, 22, 35], or their appropriate combinations [4, 15, 38]. Moreover, the linear-strain exponential rate functions, also suggested by Huxley [15], were either extended [23, 26] or simply substituted [29, 30] by the bilinear-strain exponential functions.
An extension of the number of adjustable parameters in sliding filament models via arbitrary chosen rate functions does not improve our knowledge on the nature of muscle force generation, thus many physical inconsistencies of modified Huxley's models were reported. For example, the requirement to Huxley's CBs to be attached at long distances during steadily muscle stretching was found as unreasonable [9], when studied within the scope of the distribution moment model by Zahalak [38], employed recently in [35]. Thereby, the central question was raised on the existence of Huxley's CBs, favoring myosin head attachment events at large displacements with the increase in contraction velocity. Likewise [9], controversies surrounding Huxley's approach were brought forth in careful studies of the force exposed by a single CB during lengthening [24].
Treating the sliding filament model by Huxley as an incomplete one, further generalizations over many-step mechanochemical myosin cycle, relating the kinetic rate constants with the standard Gibbs energies of actin-myosin states, were proposed in [6, 8, 11, 12, 16, 22, 2931, 36]. Given that the employed transition state theory treats i) the entire system to be in a global thermodynamic equilibrium and ii) all CB degrees of freedom as being in ergodic states, such a simplified combination of the equilibrated chemical thermodynamics and mechanically equilibrated muscle dynamics seems to be rather dubious.
The approach developed in this study employs a general tendency towards establishing the thermodynamic equilibrium between all CB states, getting some of states out of the equilibrium. The treatment of statistically random strain-dependent CB states on the basis of the general theory of fluctuations removes the aforesaid uncertainties of the transition rate functions. A new analytical solution to the standard kinetic equation of the sliding filament model by Huxley [15] is proposed in the form unifying mechanisms of the force generation during both regimes of muscle contractions.
THEORY
Huxley’s approach revisited
At a fixed muscle contraction velocity V, the occupation number nf of CB states, combining in time t myosin filament with actin
filament of the total number of sites N, obeys the common "balance" kinetic equation
d dNV
-Nv ( x, t) - — ■
fiN dx
■~NL ■ ~r = f (x)(N - Nv ) - g( x)Nv ox dt
(1)
Here f (x) and g(x) are attachment and detachment rates of the bound state located at a displacement (strain) x estimated from the
actin filament site x=0 equilibrated by CB elastic energy. The steady state Nv (X, ro) , achieved on practice by late times t » f 1 , g 1 provided qnv(x,t)/dt = 0 in equation (1) reduces it to the fundamental equation (4) in [15], namely
- Vdnv(x) = f (x)(1 - nv) - g(x)nv, (2)
2 dx
suggested for CB steady-state occupation probability n^ (x) = NF (x, <x>)/ N called by CB proportion [15]. According to Huxley [15], the CB force output is generated when x decreases at a positive velocity of sliding of the actin filament va and a negative velocity of
myosin filament VM , i.e. VA = -VM = -dx / dt > 0. The contraction velocity per one-half sarcomere V/2 determines the contraction velocity V of the muscle as a whole, when modeled by V = у — VM = 2УА ■ As can be derived from Huxley's equation (6) in [15] with in part preservation of his notations, the overall generated force
F(tota‘) = lim —
L^ro l
I
dx
F ( x)ny ( x)—
21 Л
(3)
l.
was evaluated via the force F(x) per one myosin site, as one actin site is carried past it. Here 5 is the sarcomere length, A is the trial distance between the nearest sites in the actin filament evaluated in equation (15) in [15], and N' is the number of sites in the thick filament in the overlapping zone of length L.
The solution to the steady equation (2) for the shortening regime (hereafter distinguished by index 1) was proposed as a combination of the localized (short-range) and delocalized (large-range) bound sub-states located respectively in positive and negative domains. The corresponding CB proportions, namely
n{y°c)( x) = Щ-
1 - exp
У11 x_ -1
V I h2
and
(deloc) ( x) = n nV (x) = no;
1 - exp I- У
exp | 2xg;
0 < x < h
- ro < x < 0,
are reproduced exactly from (7) and (8) in [15] and parameterized here by
„ -Л
f1 + g;
and V;= h(f; + g;).
(4)
(5)
(6)
In turn, this description of the two CB sub-states follows from the rates postulated by linear functions, namely
01
28
(7)
x x
f (x) = /1- and g(x) = glT, for 0< x <xt: h h
h,
and two piecewise constants f(x)=0, g\x)=g1', for -<»<x<0. The muscle concentric steady force pftn) reduced to the isometric force F0 found on the basis of (3)-(7), namely
VV ^ „ (8)
F ( short)
F
Vi L
1 - exP I - V
1+-
2h2 gi2
for V > 0,
was fitted in [15] by the four CB kinetic parameters: f = 43.3 s *, gt = 10.0 s *, ft = 0 , and gt' =209 s f In addition, two more adjustable parametersh » 15nm and V = V(exp) /4, where V(exp) is the empirical maximum shortening velocity, were directly and indirectly
employed when tested by Hill's empirical equation [10] (for details, see chapter IV in [15]). It is noteworthy that the nearest-site distance in the actin filament treated as a free parameter was estimated as l a h, i.e. close to the known nearest-molecular distance in the myosin filament ft = 14.5 nm [5]. However, the ratio p(stref) / ft = (f + g )/g = 5.33 reported by Huxley [15] for the muscle lengthening regime,
contrasts to the observed ratios falling between 1.8 and 2.0 (see e.g. [9, 14, 26]).
A new solution to the fundamental kinetic equation
Beyond any specific suggestions, the formal solution to the steady-state equation (2):
nr(x) = n0(x)+Anv(x) = n0 +(1 -n0)cr expf -- f \f(x’) +g(x’)\dx’\x = ^- = +^-, ^
\ xJ0 ) at 2
is valid for any contraction shortening velocity V (= - 2 X >0) and lengthening velocity V (= 2x < 0), leaving arbitrary a choice of the rate functionsf x) and g(x). The first order differential equation possesses, as common, has only one free constant (denoted in (9) by ft ), whereas another one, namely
»0(x) = -/i^, ('0)
f (x)+g( x)
controlled by muscle e intrinsic rates, straightforwardly follows from (2) taken at V = 0. One can see that in Huxley's version of (9) the constant = -n (1 -^1)-1eXp(-V1 /V) results from the periodic boundary condition ft(h) = 0, providing the non-Gaussian CB
proportion shown in (4).
Besides the boundary conditions re-considered below, the property of periodicity in the overlapping part of the actin filament of total length Nd, should be taken into consideration. Moreover, as noticed by Hill et al. ([12], p. 346]), the quantity n (x)dx has a meaning of the
late-time probability to find one of the two myosin heads in the bound state within the domain [x , x + dx]. Consequently, the total force output in the finite-range overlapped zone is
, , (• + Nd dx' ?d dx tl It
FSzone) = r p(x')«„ (x') — = N f F(x)«„ (x)—, where x'= xN. ( )
V pNd - - -- - 2d F- - - -v 2d’
Here F (x) is the active force per one actin cell, substituting that in (3) estimated per one myosin site. Such a consideration suggests the statistical equivalence of all the occupied cells in the actin filament treated as a one-dimensional crystal of the lattice constant d. If one takes into account that the sum of strain-dependent occupation numbers NK (x, да) of the one-molecule steady states x is limited by the total
number N of available actin filament sites, the following normalization condition, namely
£>>2*=j>)f = l>(x)f =L (12)
constraining the CB occupation probability, should be applied. Such a normalization is consistent with (3) and (11).
It was proposed by T. L. Hill [11] to extend Huxley's idea on the minimum of CB elastic energy at x = 0 over the minimum of Gibbs energy [8, 12]. More generally, one deals with a thermodynamic system associated with the electrically activated, loaded muscle and mechanically equilibrated system steadily contracting a fixed temperature.
Having the configurational entropy ft (x) , a given CB contracting at the velocity V tends towards the global equilibrium characterized
by the maximum ft(0) . The mean of a random strain x of the bound myosin molecule is treated as the thermodynamic quantity, linear with the generated force (see e.g. equations (3) and (16) below), where the x-state probability ft(x) is proportional to exp(SK(x)) . The latter is due to Boltzmann's principle stating that the logarithm of the probability of fluctuations is proportional to its entropy, introduced by Einstein in 1910 as the methodological foundation for the fluctuation theory [20, 25]. As common, the fluctuation part of the entropy ASV(x) = ft (x) - ft (0) can be related to the negative minimum work required to bring the given CB from a certain non-equilibrated state, associated with x strain, to the "true equilibrium state" at x = 0 under interaction with the rest of the whole system. Taking into account the maximum of entropy at x = 0, the first derivative of ft(x) is zero and the second derivative is negative that yields the Gaussian form for ft (x) , centered at x = 0.
One can see that the requirement on ft (x) to be Gaussian can be satisfied in (9) solely by the linear parameterization of the rate functions, namely
f (x) = L—, g(x) = gm—, for xm = x+ ^ x ^ 0 or xm x x_
- x < x < 0.
(13)
The shown positive and negative domains [0 , xm] arise from the velocity-displacement-sign agreement xx > 0, also required by the Gaussian form. Finally, the normalization constant in equation (9)
ki herf
. , V 2x
with и = — =-----------
Vm (f m + gm ) xm
>0,
(14)
readily follows from the normalization condition discussed in (12). In (14), the standard error function erf (y) = (2/>h)Jexp(-t2)dt,
0
acting within the domain [ erf (0) = 0 , erf (да) = 1 ], is employed.
-1
c =
Boundary conditions
29
The periodic boundary conditions for CB states should be additionally applied to the solution of (2), thereby introducing additional correlations and reducing the number of independent parameters. As seen in (9), a new solution generic for both steady regimes is composed of the static state occupied with n0(x) = fm / (fm + gm) , following from (10) and (13), and the dynamic state with Anv(x), disappearing in the static limit V^ ± 0 at any finite x located either in positive (d > x+ > x > 0) or negative (-d > - x > x > 0) strain domains described in (13). The static CB states are generally permitted to exist in a positive domain 0 < x < x0 or a negative domain - x0 < x < 0 determined by the new CB structure parameter x0 (< d) . The dynamic-state domain limitations arise from the mechanical constraints discussed below.
Bearing in mind that the positivity of the velocity U(‘'Acri) = VA - VM = -2x , in shortening regime, and the negativity of the stretching velocity V(M) =VA — VM = 2x
(see text below equation (9)) are constrained by the requirement xx > 0 , one infers that the only negative domains are mechanically available for the CB dynamic states. The provided analysis summarized in the trial (9) yields
f x
exp I----2
x0 d 2________v Vxm
d x_^VV J 1
nv (x) = n0®0(x) + (1 - n0 -гУ
erf
-Jv
®v ( x)
(15)
obtained with the help of (14). The discussed above boundary conditions now are incorporated into the novel solution to Huxley's model via the "static-state" ©0(x) and "dynamic-state" ®v(x) auxiliary functions, specifying the CB domains via the standard Heaviside (step)
function ®(y). Specifically, ®(y) is one for y > 0 and zero for y <0; ®0 (x) = ®(x) -®(x - x0) , for the positive-strain CB domain; ©0 (x) = ®(x) - ®(x + x) and ®v(x) = ®(-x) - ®(x + x_) , for negative domains.
In Fig.1, possible CB domains, simultaneously allowed by kinetics and thermodynamics for the dynamic-state fluctuations near the equilibrated static states, are shown along with the self-consistent mechanical constraints imposed on directions of the force output and sliding filament velocities.
RESULTS
Muscle force
The mean force produced by one CB occupied in a periodic cell of the actin filament
F'
1 Т/
(cell)
i »2
= F +AF, = kf^mv(x)dx =F - kx-(1 - n -г)Ф(v), F0 = + --0n0
(16)
r+d dx _ _ x
-d ^ '2d "~2" 41 d^ ~ 2d
is calculated using (11) and (15) on the basis of the CB stiffness k = F(x) / x [15, 16], velocity-independent intrinsic muscle quantity, e.g. [14]. The static part F0 is found to be positive during shortening regime and negative during stretching regime, as demonstrated in Fig. 1. The negative-domain dynamic part AF is linear with the universal dynamic function
Fig. 1 - A mechanical scheme of the steady force generation by the bound actin-myosin linkages resulting from (15) and (16). Each of the two equivalent myosin heads of the effective CB may be bound with the actin filament stationary or dynamically, forming respectively the equilibrated attached state A (shown by the open circle) within domains [0 , ± xo] or the non-equilibrium steady dynamic rotated state R (closed circle), within the negative domains [0 , -xt_] and [0 , -x2_], shown respectively for the muscle shortening and stretching. The arrows indicate directions of the sliding velocity of the actin filament VA and the myosin filament VM .
Ф(и) = 2.-
£
1 - exp
erf
1
4v
with Ф(0) = 0, Ф(1) = 0.486, and Ф(ю) = 1,
where the scaling positive parameter V is defined in (14). In this way, a unique equation
F0 n0x0
(17)
(18)
30
substituting Huxley's model (8), is proposed for the CB force output reduced to CB isometric force F0, unifying both concentric (the case of negative sign in equation (18)) and eccentric (positive sign in equation (18)) muscle contractions driven respectively by the positive and negative steady velocities V = uVm . Below, the one-parameter fitting analysis of the proposed theory is conducted on the basis of known experimental data on muscle force.
Fig. 2 - Analysis of the force-velocity branch predicted for the steady muscle shortening. The points and dashed-point curve are Hill's empirical data [10] shown in (22) and those fitted by Huxley's curve discussed in (8). The solid line is (19) that fits Hill's data at X = 0.85. Inset: The reduced CB occupation numbers modeled in the Appendix for different shortening velocities reduced to the maximum velocity V (shown by solid lines) and the corresponding attach and detach rates (dashed lines). The dotted line schematically drawn for V / V^ = 0.1 indicates a correction due to the possible dynamic-static (R-A) transient state discussed in the Discussion.
In Fig. 2, the tension-velocity data on steady shortening muscle are described by the equivalent equation following from (18), namely
F
(short)
ry(short)
= 1
Ф(^У / Vnaxl)
, 0 < V < Vnaxl, Vnaxl = ^, V»1 = (fl— + El- )■
(19)
F Po Ф№
Here, the well known property of the zero force at the maximum velocity Vmax is employed to relate V^“p) with the characteristic CB velocity Vml specified in (19), through a new auxiliary parameter X. The two adjustable parameters in the equivalent (18) and (19) are connected via the universal dynamic function ст(П°а) = Ф(Я)-l.
Fig. 3 - The steady force output by one cross-bridge occupying one cell of the actin filament versus the stretching velocity. The open circles are the mean datapoints of the forces (re-scaled by F = —l .95 pN) measured by Lombardi and Piazzesi in frog muscle fibers maintained at 50°C and negative velocities, lying between 75 and 1030 nm/s (see Fig. 7 in [14]), scaled here by d = 36 nm. The closed square
31
indicates the force per one CB reported in [24] for unspecified velocity. The theoretical curve is drawn via (20) taken at V£m,0i> = -8d / 5 . Inset: The CB reduced occupation numbers (solid, dotted and dashed lines) and attach-detach rates (dashed lines) predicted for three distinct velocities reduced to the characteristic velocity Vmod). For details, see equation (25) and Appendix. The dotted line has the meaning of that
in the inset in Fig. 2.
Similarly, the data on the reduced force output are fitted in Fig. 3 by
— (ste f V ^
= 1 + ф|— , -»< V < 0, Vm2 =-x2- (f2-+ g2-) < 0, (20)
- { Vm2 )
using the characteristic negative stretching velocity Vm2 as a free parameter. For simplicity, the model parameter (т(т2ф in (18) is taken to be one in order to fit approximately the experimental data on <T(“p) generally lying between 0.8 and 1.0, that follows from the mentioned above maximal force F_„ . Combining results of the fitting analyses in figs. 2 and 3, one has
tf-ir > 1 ^>. (21)
The consequences of the constrains in (21) imposed on kinetic and structure CB characteristics are modeled in the Appendix and shown in the insets in Figs. 2 and 3.
DISCUSSION
Model assumptions and relation to other theories
A description of the muscle action broadly requires the synthesis of mechanical, structural and biochemical concepts. The sliding filament model determined here in the simplest form includes one mechanical parameter k, related to CB stiffness in the fully activated muscle, and four structural parameters: a number N of actin sites available, a period d of the thin filament lattice, associated with the CB domain, and two sub-domain distances x0 and xm. The two kinetic parameters fm and gm characterize the slopes of the attach and detach rate functions. The CB parameters are additionally distinguished by indexes 1 and 2, specifying the shortening and stretching muscle regimes, as well as by positive and negative signs, indicating positive-strain and negative-strain domains. As common, the longitudinal CB stiffness is adopted to be state independent, i.e. k0±=k1_=k2_=k.
Similar to a number of studies [6, 15, 16, 29, 30, 35, 36], it is also assumed that the myosin head can bind to any point on the thin rigid filament. The found single strain-dependent actin-myosin bound state, characteristic of muscle cycling periodic steady-velocity process, emerges as a complex CB state. It effectively accounts for all possible translation-rotation CB random states, presented here by the minimal set of thermodynamic and mechanical characteristics governed by a single variable strain x. Since the occupation-number nV(x) has the meaning of the probability of CB located within the periodic domain of actin filament it is additionally (to other researches) normalized. Moreover, the proposed a strain-state description allows to relate the occupation-state probability n^x) to corresponding Boltzmann's entropy SVX), thereby synthesizing CB thermodynamics and mechanics.
The net muscle tension would be zero, if the thermal equilibrium was possible between all CBs degrees of freedom and the surrounding thermal bath [28]. Unlike the thermal equilibrium postulated by transition state theory [6, 8, 11, 12, 16, 22, 29, 30, 36], the found actin-myosin bound state is favored solely by a trend towards the global thermal equilibrium. The attach-detach rates found on the basis theory of thermodynamic fluctuations appear to be limited by linear-strain functions, remarkably suggested by Huxley's [15]. Consequently, the complex CB state appears to be combined from the equilibrated static and non-equilibrated dynamic bound sub-states distributed respectively via the uniform n0 and strain-dependent An^x) occupation probabilities, as shown in (9).
In order to reveal a link between the complex CB state and more realistic detail CB states, let us address the sliding filament (3G) model by Smith and Geeves [29, 30] dealing with four actin-myosin bound states: the two attached states (A-MDP and A-MT) and two rotated states (A.MDP and A.MT) called respectively by the attached A and rotated R states.
Analysis of the occupation probabilities in the shortening regime at moderate and high velocities indicates (see, respectively, Fig. 5 in [29] and Fig. 6B in [30]) that the attached and rotated states exist respectively in the negative and positive domains, in accord with our analysis shown in the inset in Fig. 2. Moreover, the overall four-state domain spreads over approximately 36 nm, the distance associated here, as well as in [12, 27, 33], with the actin filament period d (though, a distance of 38.5 nm was adopted in [29, 30] for d). One therefore arrives at more "mechanistic" treatment of the complex CB state simplified by the two-state kinetics, on the one hand, and more specific thermodynamic treatment of the A and R states oversimplified by the transition state theory, on the other hand.
Thus, mutually correlated (via occupation numbers) the static and dynamic sub-states are attributed, respectively, to the head-myosin attached, thermally (globally) equilibrated bound states and the head-myosin rotated, mechanically (locally) equilibrated bound states.
It was also suggested in [29] to use strong-binding and weak-binding states as synonymous with attached and rotated states. Moreover, the static-attached and dynamic-rotated states were independently revealed by Rayment et al. [27], via the observations "of catalytic domains of myosin being initially weakly attached to actin, and the following structural changes resulting in tight binding of actin-myosin linkages". Thereby, structural studies of force generators in contracting muscles strongly collaborate with main properties of CBs of this study obtained through the application of the second law of thermodynamics to the new static and dynamic solutions of the two-state differential equation pioneering by Huxley.
The force generation mechanism
As shown in (16), the CB force output FV is due to the interplay between the velocity-independent force F0 and contractive force AFV produced by, respectively, static A and dynamic R sub-states. During concentric muscle contractions with a positive velocity V, the so-called isometric force F0 is also positive, whereas the contractive force AFV is negative as found in (19) and that collaborates with the ATP hydrolysis results. During eccentric contractions, commonly associated with the negative direction of the velocity V, both the forces are also negative (see Fig. 1).
A question arises on the role of A and R states in the muscle force generation revealed by force-velocity steady behavior by other researches. Although detail analysis on the basis of the 3G model was made neither for steady muscle shortening [29] nor stretching [30], it was mentioned in [30] that A and R states are respectively responsible for the low-force and high-force output. As for Duke's model [6], it was only reported on the fitting of the data by Edman [7] on muscle steady shortening by four dimensionless combinations of model parameters, but no explicit force-velocity relation was provided or discussed.
The seminal comparative analysis of A. F. Huxley's and A. V. Hill's data on muscle concentric tension (see [15], p. 287) is revisited in Fig. 2, where the fitted experimental data are presented by the famous phenomenological relation [10], namely
p(exp> p(exp> 1 у / тл(ехр>
-P— = = a----------, (22)
-0(cxp) -0(cxp) a + V / Vj;“p)
where a = 0.25 corresponds to the force-velocity data —(exp) on the titanic stimulated shortening frog muscle at 0°C. The analysis of the
high-velocity wing of the tension-velocity data provided with the help of (19) indicates that the attached A states, having a positive domain, control the curve above V/Vmax = 0.5. This conclusion is due to the occupation probabilities of the non-equilibrated rotated states decrease with sliding velocity. Such a high-velocity behavior, shown in the inset in Fig. 2, is remarkably in agreement with the negatively located R
32
states shown in Fig. 6B in [30]. It is noteworthy that Huxley's delocalized state determined in (5), also fits the high-velocity data by A. V. Hill [10] and Edman [7]. However, the low-velocity Hill's data within 0.2 < V/Vmax < 0.5 reveal the inconsistency with Huxley's CBs (4), located in the positive domain and exerted negative force AFr. As shown by current study, the CB positive-strain rotations contradict to the overall trend of achieving of maximum entropy. Consequently, instead of the retarded detachment g < f1) of Huxley's CBs discussed in (4) and (7), the negatively dislocated rotations of myosin heads require faster detachment than attachment (gm > fm1) during muscle steady shortening, and slower detachment (gm2 < fm2 ) during muscle steady lengthening, as illustrated in the insets in Figs. 2 and 3.
As for the regular deviation of the solid line in Fig. 2 from the data at very low velocities , it can be understood by an incomplete description of the almost-equilibrated slow rotations lying in the vicinity of x ~ 0 where R and A states overlap within the distance S of ~ 5.5 nm, as seen in Fig. 5 in [29] and also in Fig. 6B in [30]. Our fit analysis of the low-velocity data could be therefore improved when the CB occupation probabilities will be locally symmetrized within the range -S < x <S, as shown by the dotted line for the particular case V/Vmax =
0.1 in the inset in Fig. 2. As for the discrepancy between the data by A. V. Hill [27] and Edman [7], the likely distinct experimental conditions underlying the preparation of equilibrated and non-equilibrated CB states during steady muscle shortening need further analysis.
In Fig. 3, the upper branch of the generic force-velocity curve (18) fits well the data on CB force during the stretch [29], when taken under experimental condition F)stret* = 2F0 and drawn via (20) using a single adjustable parameter Vf2xp). Likewise the analysis in Fig. 2, a deviation of the theoretical curve from the data at low velocities can be improved when a certain R-A overlapping bound state is introduced, similarly to that suggested by Mehta and Herzog [24] in their Fig. 3.
Stretching against shortening
In view of the additivity of the CB force output, the provided study on the steady-state force-velocity equation (18) can be equally applied to the reduced tension due to a single CB, sarcomere, fiber, or muscle as a whole during its steady shortening or lengthening. As shown in (25) in the Appendix, the model conditions (21) of observation of the data via the universal dynamic function Ф(и) are reduced to the inequalities ni(tre,) > 1/2 > n(shhrt), limiting occupation numbers of the attached CB states. Since n0 = fm / (fm + gm), this finding is reformulated into the aforesaid inequalities f(shhrt) < g((^hort) and f(ss'el) > g{(l'e,). In contrast to Huxley-like models, see e.g. [9], this finding is supported by recent data on the cycle duty ratio в ( = 1 - n0, the characteristic time of attachment fm1 related to total CB cycling time fm~1+ gm_1) [24] generally indicating that pstret < pshort (see table 1 in [24]). Therefore, one can extend the above inequalities to
(short) Ashort) Astret) (stret)
gm > Jm > Jm > gm ’ (23)
in accord with the insets in Figs. 2 and 3. Such a hierarchy of the rates established for the new CBs corroborates the working hypothesis in [24] on "that a stretched cross-bridge might remain attached longer than a cross-bridge that had been shortened while attached", i.e.
(.tfret)-1 (short )-1
gm gm .
As may be expected, this deduction conflicts with the conclusion that the "reattachment (in steady lengthening is)... faster than attachment in the isometric condition or during shortening" [14], one of the major results obtained in [14] on the basis of Huxley's (4) and (5). A crossover between the stretching and shortening muscle regimes is shown to be due to a reconstruction of the bound sub-states, including structural domains and exchange kinetic parameters, leaving unchanged the mechanical CB characteristics.
Implications of the approach
In spite of revealed by different researches inconsistencies of Huxley's CBs with the experimental data on both muscle steady shortening and stretching, one can see that the major mechanical concept on the CB compliance [15] remains to be common to all modern sliding filament models. The incorporation of conformational changes of myosin molecules into the two-state model [15] improved kinetic cycle by multiple CB states [6, 14, 21, 22, 29, 30, 35], characteristic of real striated and smooth muscles. Further generalization of the rigid sliding filament model [15] over tropomyosin deformations resulted in strong correlations between originally "independent force generators" [14]. The corresponding incorporation of the filament compliance resulted in the effect of cooperatively in CB tension production, now understood at both molecular and macroscopic levels of consideration, i.e. via the filament crystal-like periodicity, e.g. [21], and through the elastic continues medium, e.g. [31].
In the same context, a simplified consideration for the complex, attached-and-rotated CB state, propagating along periodic actin filament lattice of period d = 36 nm [10, 27, 33] is considered here in the continues-strain approximation, likewise the studies [6, 8, 11, 12, 14-16, 26, 29, 30, 36, 38], naturally excluding the so-called "target zones", favorably oriented sites on the helices of the actin filament separated by distances of about 5.5 nm [21, 33], 11 nm [21, 34], or 36 nm [21, 33]. One can see that unlike the target-zone domains, associated with certain load-dependent bound sub-states, ultimately determining the preferable step distance [21], the load-independent distance d plays the role of an external parameter, that may scale structural, mechanical and dynamic characteristics, as shown in (18). Given that the force-velocity steady muscle behavior is not affected directly by variations of the structural or energetic parameters, a new additional domain activated during stretching in the myosin filament [1] can be readily included.
The utility of the proposed improvement of the simplest sliding filament model may be also clarified by the fact that any temporal nonsteady muscular behavior, including residual force effects [35], is underlaid by the steady-state solutions to kinetic equations [3, 14, 26, 35]. The principal suggestion, differing the current study from other analytical approaches explaining the force production during muscle shortening [6, 12, 22, 24, 30] and stretching [29, 30], is that instead of the thermally equilibrated CB states postulated by the transition-state theory, it is adopted here that not all of them are allowed to be in thermal equilibrium. According to the generalized fluctuation theory [17, 20, 25], directly exploring the second law of thermodynamics, the configurational part of the entropy drives the muscle loaded system towards the global equilibrium throughout the locally equilibrated CB rotated states. As the result, the mechanically-equilibrated rotations of 51 myosins, attaching and detaching via linear-strain-dependent rates, are predicted to be distributed via the symmetric Gaussian form. The existence of such dynamic rotations located in the negative-strain CB domains is revealed here via both branches of the generic force-velocity curve. The extension of the proposed approach on the specific-state transient CB kinetics will not affect the overall-state domain and therefore the general picture of force generation.
Another prediction, challenging further experimental verification, is that the entropy production during the steady stretching regime is lower than that during shortening regime. This finding arises from the constraint imposed on the mean numbers for the rotated states, i.e. &hplrel) < bhy"l), opposing the case of the attached static states. Thereby, a reduction of the energy consumption during lengthening can be understood within the same framework, excluding a speculation about the existence of two distinct mechanisms for muscle steady shortening and steady lengthening, as initially suggested by Huxley ([15], p. 292) and then modeled by Lombardi and Piazzesi [14, 26].
Alternative approach in the framework of the sliding filament has been recently proposed Shilko and co-workers in [2] and developed in [28]. On the basis of the application of Huxley’s kinetic equations to pulling and braking cross-bridges, the stationary solutions for the mean statistical occupation numbers of bridges was obtained [2], leading to the phenomenological Hill equation and thus revealing the origin of equation (22) and its free parameters at a molecular level. As a next step, it is interesting to test the proposed model for experimental data for the muscles during steady elongation and thereby to generalize the same Hill's equation for the phase of muscular stretching.
CONCLUSION
The provided study clearly demonstrates how the force response to external loads, emerging in the steadily shortening or stretching muscle, yields universal correlations between mechanical, structural and kinetic model muscle characteristics, when rationalized in light of
33
the fundamental law of thermodynamics. Likewise, the correlations caused by accommodation of the generated force to distinct patterns of external loads during efficient eccentric, isometric, or concentric contractions were revealed in terms of the primary locomotor muscle functions via the most general laws of the theory of elasticity [18,19].
Acknowledgements
The author is grateful Scott Medler for helpful comments.
APPENDIX
Modeling CB domains
The inequalities shown in (21) following from the fitting analysis provided in figs. 2 and 3 are now specified by
x,d . x~d „„
% < ——11----------- and — >-----—----------- (24)
x0i(x01 + xi- ) x02- (x02- ^ x2- )
with the help of (18), where negative signs indicate negative CB-strain domains. Then, additional geometrical constraints Xi_< d / 2 and x2- < d, corresponding to the fact that the tail of the myosin molecules is longer than heads (see Fig. 1), generally results in x2_ > xt_, thereby providing и02 > «oi in (24). If a simplified model requirement of periodicity of CB states x(m°d) + x<m°d) = d is also adopted, along with a condition of the coexistence of the sub-states (e.g. x(m°d) = x(m°d) = d, one has the requirement
«<m°d) > 1/2>«<m°d), (25)
For the case of muscle shortening, the chosen CB structure parameters are x,(m°d) = 2d /3 and xf™*® = d /3 , which in combination with the fitting parameter <r(le1xp) = 1.22 shown in (18) specify the numbers n(m°d) = 0.44 and n(m°dl(x) for the static (attached A) and dynamic (rotated R) states illustrated in the inset in Fig. 2. Moreover, when the maximum velocity ¥(£?! = 3200 nm/s for the experimental data analyzed in Fig. 2 is employed, the model attachment and detachment rates y'l(m°d) = 138s_1 and gj™*® = 176 s_1 are straightforwardly follow from the relationsf_ = Vmx1 «oi/Xxi_ and gi_ = fi_(1 - «oi) / «oi provided by (9), (10), and (19).
For the case of activated muscle stretching, the fitting analysis in Fig. 3 provides Vjje4’) = - 288 nm/s). In the inset in Fig. 3, the used
estimate n(m°d) = 0 75 for the static-state occupation data is obtained with the help of the relation «02_ = f2_X2_ / |Vm2|, following from (10) and (20), where (/2_)i is associated with the characteristic dwell time for the CB attachment found about 0.167 s [24]. Such a simplified modeling provides the model data /2(m°d) = 6 s“i and g(m°d) = 2 s“i estimated via g2_=f2_(1 - «02_) /«02_ and illustrated in the inset in Fig. 3.
References
1. Brunello E., Reconditi M., Elangovan R., Linari M., Sun Y.B., Narayanan T., Panine P., Piazzesi G., Irving M., Lombardi V. Skeletal muscle resists stretch by rapid binding of the second motor domain of myosin to actin // Proc. Natl. Acad. Sci. USA. - 2007. - No 50. - P. 20114-20119.
2. Chernous D. A., Shilko S. V. Modelling of contractive activity of the muscle tissue // Russian Journal of Biomechanics. - 2006. -No. 3. - P. 51-60.
3. Chin L., Yue P., Feng J. J., Seow C.Y. Mathematical simulation of muscle cross-bridge cycle and force-velocity relationship // Biophys. J. - 2006. - No 10. - P. 3653 -3663.
4. Cole G.K., Bogert A.J., Herzog W., Gerritsen K.G. M. Modelling of force production in skeletal muscle undergoing stretch // J. Biomech. - 1996. - No 8. - P. 1091-1104.
5. Craig R., Woodhead J. L. Structure and function of myosin filaments // Curr. Opin. Struc. Biol. -2006.- No 2. - P. 204-212.
6. Duke T. A. J. Molecular model of muscle contraction // Proc. Natl. Acad. Sci. UsA - 1999. - No 6. - P. 2770-2775.
7. Edman K. A. P. Double-hyperbolic force-velocity relation in frog muscle fibres // J. Physiol. - 1988.- No 10. - P. 301-321.
8. Eisenberg E, Hill T. L., Chen Y. D. Cross-bridge model of muscle contraction. Quantitative analysis // Biophys. J. - 1980.- No 2. -P. 195 -227.
9. Harry J. D., Ward A. W., Heglund N. C., Morgan D. L., McMahon T. A. Cross-bridge cycling theories cannot explain high-velocity lengthening behavior in frog muscle // Biophys. J. - 1990.- No2. - P. 201-208.
10. Hill A.V. The heat of shortening and the dynamic constants of muscle // Proc. Royal Soc. London - 1938.- No 843. - P. 136-195.
11. Hill T. L. Theoretical formalism for the sliding filament model of contraction of striated muscle, part I. // Prog. Biophys. Molec. Biol. -1974.-P. 267-340.
12. Hill T. L., Eisenberg E., Chen Y., Podolsky R. J. Some self-consistent two-state sliding filament models of muscle contraction // Biophys. J. - 1975. - P. 335-372.
13. Huxley H. E., Hanson J. Changes in the cross-striations of muscle during contraction and stretch and their structural interpretation // Nature - 1954.- No 4412. - P. 973-976.
14. Huxley A. F., Niedergerke R. Interefrenec microscopy of living muscle fibres // Nature - 1954.- No 4412. - P. 971-973.
15. Huxley A. F. Muscle structure and theories of contraction // Prog. Biophys. Biophys. Chem. - 1957.- No 7, - P. 255-318.
16. Huxley A.F., Simmons R.M. Proposed mechanism of force generation in striated muscle Nature - 1971.- No 5321. - P. 533-538.
17. Kestin J., Dorfman R. A course in statistical thermodynamics . - New York, San Francisco: Academic Press, 1971.
18. Kokshenev V. B. A force-similarity model of the activated muscle is able to predict primary locomotor functions // J. Biomech. -2008.- No 4. - P. 912-915.
19. Kokshenev V. B., Modeling mechanical patterns for striated muscles // J. Biomedical Science and Engineering, - 2014.- No 8. - P. 473-496.
20. Landau L. D., Lifshitz E. M. Statistical Physics, third ed. - London: Pergamon Press, 1989.
21. Lan G., Sun S. X. Flexible light-chain and helical structure of F-actin explain the movement and step size of myosin-VI // Biophys. J. - 2006.- No 11. - P. 4002-4013.
22. Lan G., Sun S.X. Dynamics of myosin-drive skeletal muscle contraction: I. steady-state force generation // Biophys. J. - 2005-No 6. - P. 4107-4117.
23. Lombardi V., Piazzesi G. The contractile response during steady lengthening of stimulated frog muscle fibres J. Physiol. - 1990-No 12. - P. 141-171.
24. Mehta A., Herzog W. Cross-bridge induced force enhancement? // J. Biomech. - 2008.- No 7. - P. 1611-1615.
25. Onsager L., Machlup S. Fluctuations and irreversible processes // Phys. Rev. - 1953.- No 6. - P. 1505-1512.
26. Piazzesi G., Francini F., Linary M., Lombardi V. Tension transients during lengthening of tetanized muscle fibre of the frog // J. Physiol. - 1992.- No 1. - P. 659-771.
27. Rayment I., Holden H. M., Whittaker M., Yohn C. B., Lorenz M., Holmes K. C., Milligan R. A. Structure of the actin-myosin complex and its implications for muscle contraction // Science - 1993.- No 5117. - P. 58-65.
34
28. Shilko S. V., Chemous D. A., Pleskachevsky Yu. M. Nonlinear deformation of skeletal muscles in a passive state and in isotonic contraction // Mechanics of Composite Materials. - 2012. - No 3. - P. 331-342.
29. Smith D. A., Geeves M. A. Strain-dependent cross-bridge cycle for muscle // Biophys. J. - 1995.- No 2. - P. 524-537.
30. Smith D. A., Geeves M. A. Strain-dependent cross-bridge cycle for muscle I. Steady-state behavior // Biophys. J. - 1995.- No 2. -P. 538-552.
31. Smith D. A., Maytum R., Geeves M. A. Cooperative regulation of myosin-activation interactions by a continuous flexible chain I: Actin-tropomyosin systems // Biophys. J. - 2003.- No 5. - P. 3155-3167.
32. Smith D. A., Geeves M. A. Cooperative regulation of myosin-activation regulations by a continuous actin-tropomyosin-troponin and regulation by calcium // Biophys. J. - 2003.- No 5. - P. 3168-3180.
33. Steffen W., Smith D., Simmons R., Sleep J. Mapping the actin filament with myosin // Proc. Natl. Acad. Sci. USA - 2001.- No 26. -P. 14949-14954.
34. Veigel C., Molloy J. E., Schmitz S., Kendrick-Jones J. Load-dependent kinetics of force production by smooth muscle myosin measured with optical tweezers // Nature Cell Biol. - 2003.- No 11. - P. 980-986.
35. Walcott S., Herzog W. Modeling residual force enhancement with generic cross-bridge models // Math. Biosci. - 2008.- No 2. - P. 172-186.
36. Walcott S., Sun S. X. Hysteresis in cross-bridge models of muscle // Phys. Chem. Chem. Phys. - 2009.- No 24. - P. 4871-4881.
37. Wu J. Z., Herzog W., Cole G. K. Modeling dynamic contraction of muscle using the cross-bridge theory // Math. Biosci. - 1997-No 1. - P. 69-78.
38. Zahalak G. I. A distribution-moment approximation for kinetic theories of muscular contraction // Math. Biosci. - 1981.- No 1-2. -P. 89-114.
Сафонова В.Ю.
Доктор биологических наук, доцент, Оренбургский педагогический университет НЕКОТОРЫЕ ПОДХОДЫ К БЕЗОПАСНОСТИ ПРОДУКТОВ ПИТАНИЯ В ЧРЕЗВЫЧАЙНЫХ СИТУАЦИЯХ
Аннотация
В статье рассматриваются вопросы безопасности продуктов питания в чрезвычайных ситуациях, связанных с загрязнением продукции сельскохозяйственного производства радионуклидами. Эффективные способы консервирования мяса в условиях отсутствия холодильных установок.
Ключевые слова: безопасность, чрезвычайные ситуации, консервирование мяса.
Safonova V.Ju.
PhD in Biology, associate professor, Orenburg State Pedagogical University SOME APPROACHES TO FOOD SAFETY IN EMERGENCIES
Abstract
In the article the questions of safety offoodstuffs are examined in the emergencies related to contamination ofproducts of agricultural production by радионуклидами. Effective methods of canning of meat in the conditions of absence of refrigeration units.
Keywords: safety, emergencies, canning of meat.
Если рассматривать продовольственную безопасность страны в целом, то она складывается из физической, экономической доступности продовольствия и безопасности питания.
В первом случае она предусматривает наличие продуктов питания в нужном объеме и ассортименте на всей территории страны при бесперебойном их поступлении в сферу потребления. Это условие возможно за счет обеспечения государственного контроля, как за внешними, так и за внутренними поставками продовольствия, включая имеющиеся запасы продуктов питания.
Во втором - предусматривает экономическую доступность продовольствия. Это то, что каждый гражданин страны должен иметь такой уровень доходов, который бы позволял ему приобрести минимальный набор продуктов питания. Данное условие обеспечивается, как за счет достатка доходов населения, так и оптимального само обеспечения населения продовольствием за счет подсобных хозяйств и дачных участков.
В третьем случае продовольственная безопасность предусматривает безопасность питания. При этом большое значение имеет качество сырья и продуктов питания. Они должны соответствовать установленным требованиям и гарантировать безопасное их потребление. В этом плане интерес представляет уровень загрязнения продукции сельскохозяйственного производства радионуклидами, появление которых обусловлено аварийными ситуациями на предприятиях ядерного топливного цикла [1; 3].
В чрезвычайных ситуациях при масштабном радиационном поражении сельскохозяйственных животных при отсутствии холодильных камер и консервирующего оборудования наиболее эффективным и экономически выгодным методом является посол [2]. Традиционный метод консервирования посолом предусматривает температуру воздуха при хранении не выше +40С. При более высоких температурах он не обеспечивает получения доброкачественной солонины. В нашей работе использовалось мясо, полученное от внешне облученных ионизирующим излучением сельскохозяйственных животных и подвергнутое нетрадиционным способам консервирования и хранения при температуре +15...+25 0С. Учитывая особенности температурных условий консервирования и хранения мяса (говядина) были использованы две рецептуры консерванта.
Первая рецептура состоит из поваренной соли, соляной кислоты и бензоата натрия в соответствующих пропорциях в питьевой водопроводной воде.
Вторая - состоит из поваренной соли, молочной кислоты, бензоата натрия, воды питьевой водопроводной в соответствующих пропорциях. Тот и другой консервант включает в обязательном порядке поваренную соль, кислоту, либо соляную, либо молочную, бензоат натрия и воду. Очень простой и доступный состав.
В работе оценивали ветеринарно-санитарные показатели говядины, полученной от внешне облученных животных и подвергнутых нетрадиционным способам консервирования. После созревания туш облученных и контрольных животных брались куски мяса весом 1 -2 кг, которые подвергались последующему консервированию и последующему хранению. В качестве консервантов использовались две выше представленные рецептуры, обеспечивающие сохранение доброкачественности говядины, свинины и баранины по существующим ГОСТам в условиях консервирования и хранения мяса при температуре воздуха +15-+250 С. При этом мясо облученных животных помещали в емкости с полиэтиленовыми прокладками и заливали консервантом с жидкостным коэффициентом 1:3 и хранили в течение 60 суток при температуре комнатной температуре. Исследования качества продукта проводили через 30 и 60 суток консервации.
В результате исследований установлено, что через 30 суток хранения консервированное мясо имело характерный для свежей солонины из говядины запах. Микроскопия мазков-отпечатков не выявила наличия микроорганизмов. Следы распада мышечной ткани отсутствовали. При постановке реакции с сернокислой медью бульон оставался прозрачным. Содержание аминоаммиачного азота в опытных группах составило 1,02±0,09 мг и 1,07±0,04 мг соответственно, что соответствовало контрольным величинам. Реакция на пероксидазу во всех группах была положительной. Микробная обсемененность мышечной ткани в опытных образцах не
35
