Математические модели процессов термической инактивации микроорганизмов.
Прежде всего остановимся на количественных характеристиках эффективности термической стерилизации, которые наиболее часто употребляются в литературе (табл. 6).
Таблица 6. Наиболее употребляемые характеристики эффективности термической стерилизации
| Характеристика | Физический смысл | Условное обозначение |
| Смертельное время | Время, необходимое для инактивации всех микроорганизмов данной популяции | ТС |
| Время гибели 90% популяции | Время, необходимое для снижения концентрации микроорганизмов в 10 раз | Ио |
| Зависимость отмирания микроорганизмов от температуры | Разность температур, необходимая для изменения т1о в 10 раз | Z |
| Удельная скорость гибели | Количество микроорганизмов, отмирающих в единицу времени | k |
| •Фактор инактивации | Логарифм отношения начального числа спор (Ng) к конечному числу спор (У), выживших после нагревания | п |
| Критерий стерилизации | Натуральный логарифм отношения количества жизнеспособных спор Уо до стерилизации к количеству N жизнеспособных спор в момент времени т после начала стерилизации | V |
Значения некоторых характеристик приведены в табл.
7—11.Таблица 7. Значения т10 (в с) для суспензий спор некоторых аэробных бактерий [148]
Температура, °С
| Бактерия | 100 | 120 | 130 | 140 | 150 | 160 |
| Вас. cereus | 12,1 | 4,2 | 2,6 | 1,3 | 1,0 | 0,7 |
| Вас. subtil is | 27,8 | 4,5 | 3,1 | 2,1 | 1,1 | 0,5 |
| Вас- stearothermophilus | 2857,0 | 38,6 | 8,8 | 3,9 | 2,4 | 1,4 |
| Микроорганизм и штамм | Суспендирующая среда | Температура, °С | Т19, мин | ||||
| Вас. cereus | |||||||
| NCTC 5893 | Фосфатный буфер (pH 7,0)' | 100 | 5,5 | ||||
| То же | То же | 121,1 | 0,0065 | ||||
| Вас. licheniformis | |||||||
| NCTC 7589 'Вас. pantothenticus | Не указана | 100 | 13,5 | ||||
| NCTC 8162 | 100 | 7,3 | |||||
| Вас. megatherium | |||||||
| АТСС 8245 | Вода | 100 | 0,83 | ||||
| NCTC 7581 | » | 100 | 1,0 | ||||
| То же | » | 115 | 0,025 | ||||
| Бас. coagulans | Фосфатный буфер (pH 7,0) | 121,1 | 0,04 | ||||
| Sp. 33 | Не указана | 106 | 6,3—7,3 | ||||
| То же | » | 108 | 4,2—4,6 | ||||
| » | 110 | 2,2—2,5 | |||||
| » | » | 112 | 1,2—1,8 | ||||
| Фосфатный буфер (pH 7,0) | 121,1 | 3,0 | |||||
| Бас. coagulans var. therrnoacidurans | |||||||
| АТСС 8038 | То же | 89 | 32,05 | ||||
| То же | 93 | 18,8 | |||||
| » | 96 | 8,31 | |||||
| » Бас. stearothermophilus | » | 86 | 71,14 | ||||
| NCIB 8919 | Вода | 115 | 25 . | ||||
| То же | 115 | 24 | |||||
| » | Дрожжи | 115 | 34 | ||||
| NCIB 8157 | Вода | 121 | 4 | ||||
| То же | » | 115 | 11 | ||||
| NCA 1518 | Не указана | 115,6 | 7,03 | ||||
| То же | Снятое молоко | 127 | 0,13 | ||||
| NCIB 8224 | Фосфатный буфер (pH 7,0) | 121,1 | 14,0 | ||||
| Вас. subtil is | |||||||
| NCTC 8236 | Вода | 100 | 11,3 | ||||
| То же | То же | 105 | 2—2.5 | ||||
| Bacillus | |||||||
| 1-11 | Фосфатный буфер (pH 6,8) | 107,5 | 14 | ||||
| То же | То же | ПО | 6,9 | ||||
| » | 112,5 | 4,3 | |||||
| 115 | 2,2 | ||||||
| Микроорганизм и штамм | Суспендирующая среда | Температура, СС | Т,о» мин | ||||
| То же | То же | 117,5 | 1,3 | ||||
| » | 121 | 0,57 | |||||
| 1-12 | 107,5 | 16,6 | |||||
| То же | 110 | 9,1 | |||||
| » | 112,5 | 4,9 | |||||
| 115 | 2,75 | ||||||
| » | 117,5 | 1,55 | |||||
| » | 121 | 0,71 | |||||
| 5230 (15U) | Фосфатный буфер (pH 7,0) | 121,1 | 0,48 | ||||
| То же | То же | 121,1 | 0,08 | ||||
| 121 | 0.35 | ||||||
| Cl. perfringens | |||||||
| АТСС 3624 | Не указана | 100 | 0,31—0,39 | ||||
| Ш-1168 | То же | 100 | 1,0—1,5 | ||||
| Т-65 | » | 100 | 5,2—6,6 | ||||
| NCTC-8238 | 100 | 13,0—17,6 | |||||
| С1. sporogenes | |||||||
| РА 3679 | Фосфатный буфер (pH 7,0) | 107 | 13,4—16,1 | ||||
| То же | То же | 111 | 5,84—8,04 | ||||
| » | » | 115 | 2,77—3,63 | ||||
| » | 90 | 1000 | |||||
| » | » | 121 | 0,84 | ||||
| NCTC 532 | Вода | 70 | 1170 | ||||
| То же | То же | 80 | 245,8 | ||||
| Cl. botulinum | » | 90 | 34,2 | ||||
| Е | Не указана | 77 | 1,8 | ||||
| 62 А | Буфер | 112,8 | 1,23 | ||||
| 12885 А | То же | 112,8 | 1,09 | ||||
| 213 В | 112,8 | 1,32 | |||||
| 32 В | 112,8 | 0.15 | |||||
| Cl. butyricum | |||||||
| 555 | Фосфатный буфер (pH 6,5) | 85 | 23 | ||||
| То же | То же (pH 7,0) | 85 | 23 | ||||
| » | Молоко | 85 | 14 | ||||
| Cl. butyricum | |||||||
| 556 | Фосфатный буфер (pH 6,5) | 85 | 18 | ||||
| То же | То же (PH 7,0) | 85 | 12 | ||||
| Молоко | 85 | 21 | |||||
| Суспендирующая | Темпера- | ||
| Микроорганизм и штамм О. tyrobutyricum 500 То же » 514 То же » 518 То же » » » » €1. aureofoetidum NCTC 505 То же Cl. histolyticum NCIB 503 То же 3» | среда Фосфатный буфер (pH 6,5) То же (pH 7,0) (рН7,5) Молоко Фосфатный буфер (pH 7,0) Молоко Фосфатный буфер (pH 6,5) То же (pH 7,0) (рН»7,5) Молоко Фосфатный буфер (pH 7,0) Молоко Фосфатный буфер (pH 6,5) То же (pH 7,0) (РН7.5) Молоко Фосфатный буфер (pH 7,0) Молоко Вода То же » » » | тура, °С 90 90 90 90 121 121 90 90 90 90 121 121 90 90 90 90 121 121 80 90 95 70 80 85 90 | Т10, мин 6.5 13 13 9 0,012 0,006 15 20,9 18 14.1 0,012 0,008 11 18 21 10 0,003 0,007 4205 139.1 27,2 1097 243.1 40.5 11.5 |
Более полные сведения о значении тю по набору микроорганизмов даны в работе [80], однако представленный температурный диапазон не охватывает всех значений, характерных для режимов термической стерилизации в современных микробиологических производствах. Тем не менее эти данные позволяют заклю- 23
чить, что с повышением температуры значения тю уменьшаются; наиболее резистентными к теплу являются споры микроорганизма Вас. stearothermophilus.
Терморезистентность различных типов микроорганизма (например, Cl. botulinum) варьирует в достаточно широком диапазоне. Наконец, на значения тк> существенно влияет суспендирующая среда, в которой находятся микробные клетки.
Даже приведенные весьма неполные данные о показателях термической устойчивости микроорганизмов указывают на достаточно большой объем информации, полученной в этой области. Заметим, что значения тщ
Таблица 9. Значения 7. для спор микроорганизмов, экспонированных к влажному теплу [80]
| Микроорганизм | Температура, °С | Микроорганизм | Температура, °С |
| Вас. stearothermophilus | 7 | Cl. tetani | 9,0 |
| Термофилы | 10 | Cl. botulinum В-364 | 6,2 |
| Вас. subtil is | 13 | Cl. botulinum | 12 |
| То же | 8,2 | в фосфатном буфере | 8,4—10 |
| 9,4 | Cl. tyrobutiricum | ||
| Вас. cereus | 7,4* | Cl. perfringens | 10 |
| Вас. licheniformis | 6,0 | Cl. aureofoetidum | 6,8 |
| Вас. megatherium | 8,8 | Cl. histolyticum | 10 |
| Вас. pantothenticus | 10,6 | Cl. sporogenes Cl. sporogenes PA 3679 | 13 10,6 |
Таблица 10. Значение удельной скорости гибели спор Вас. stearothermophilus штамм 1618 при различных температурах [42, 145]
| Температура, CC | K, хМИН'~1 | Температура, °С | к, мин-1 | Температура, °С | к, мин--1 | Температура, °С | к, мин"1 |
| 100 | 0,013 | 111 | 0,193 | 122 | 2,44 | 133 | 23,4 |
| 101 | 0,017 | 112 | 0,234 | 123 | 3,07 | 134 | 29,5 |
| 102 | 0,023 | 113 | 0,302 | 124 | 3,76 | 135 | 37,2 |
| 103 | 0,030 | 114 | 0,412 | 125 | 4,57 | 136 | 46,8 |
| 104 | 0,036 | 115 | 0,540 | 126 | 5,90 | 137 | 59,0 |
| 105 | 0,048 | 116 | 0,653 | 127 | 7,40 | 138 | 74,2 |
| 106 | 0,062 | 117 | 0,810 | 128 | 9,35 | 139 | 93,5 |
| 107 | 0,083 | 118 | 1,002 | 129 | 11,70 | 140 | 104,8 |
| 108 | 0,109 | 119 | 1,210 | 130 | 14,80 | 141 | 118,0 |
| 109 | 0,135 | 120 | 1,480 | 131 | 16,60 | 142 | 148,0 |
| 110 | 0,163 | 121 | 1,830 | 132 | 18,60 | 143 | 186,0 |
Таблица 11. Значения критерия стерилизации для спор Вас. stearothermophilus штамм 1518 при различной температуре (42, 145]
| Температура, °С | V | Температура, °С | V | Температура, °С | V | Температура, °С | V |
| 100 | 0,013 | 111 | 0,91 | 122 | 11,8 | 133 | 130,9 |
| 101 | 0,030 | 112 | 1,15 | 123 | 14,9 | 134 | 160,4 |
| 102 | 0,053 | ИЗ | 1,41 | 124 | 18,3 | 135 | 197,6 |
| 103 | 0,083 | 114 | 1,86 | 125 | 23,2 | 136 | 244,4 |
| 104 | 0,120 | 115 | 2,40 | 126 | 29,1 | 137 | 298,4 |
| 105 | 0,170 | 116 | 3,05 | 127 | 36,5 | 138 | 377,6 |
| 106 | 0,230 | 117 | 3,86 | 128 | 45,8 | 139 | 471,1 |
| 107 | 0,310 | 118 | 4,86 | 129 | 57,5 | 140 | 576,0 |
| 108 | 0,420 | 119 | 6,07 | 130 | 72,3 | 141 | 694,0 |
| 109 | 0,560 | 120 | 7,55 | 131 | 88,9 | 142 | 842,0 |
| 110 | 0,720 | 121 | 9,36 | 132 | 107,5 | 143 | 1028‘0 |
получены в основном при температурах, характерных для режимов стерилизации различных продуктов в пищевой промышленности. Они не охватывают температурного диапазона, используемого в микробиологической промышленности при стерилизации оборудования и различных жидкостей. Однако, исходя уже из этих данных, можно заметить, что даже микроорганизмы одного вида существенно отличаются по термоустойчивости в зависимости от штамма и характеристик среды экспонирования.
Одним из основных в математических моделях, основанных на понятии смертельного времени, является соотношение [97]
где тс — время отмирания всех микроорганизмов в популяции при температуре (с; т — то же, но при температуре t; Z — см. табл. 6. Очевидно, что
Расчет эффективности режимов стерилизации основан на следующих предпосылках [72, 84]. Убыль количества жизнеспособных бактериальных спор во времени при повышенной температуре подчиняется логарифмической зависимости. Величина тю не зависит от кон- 25
центрации спор. Наибольшее количество выживших во> время стерилизации спор соответствует наименее прогреваемой точке. Конечный эффект процесса стерилизации зависит от концентрации наиболее термоустойчивых бактериальных спор.
Иногда в практических расчетах используют следующее соотношение:
где Гс — приведенный стерилизующий эффект (время термической гибели); / — температура в наименее прогреваемой точке; т — продолжительность поддержания температуры /; kF — переводной коэффициент, называемый иногда коэффициентом летальности.
Часто в расчетах, особенно в технологии пищевых производств, пользуются экспериментально определенными характеристиками термической устойчивости спор Cl. botulinum, для которых Z=10°C. Значения коэффициента летальности для некоторых спор при температурах, характерных для режимов стерилизации, приведены в табл. 12.
Таблица 12. Значения коэффициента летальности
kp для спор [32]
Стерилизующий эффект определяют как сумму значений kF при температуре, достигаемой в наименее прогреваемой точке в периоды нагревания, выдержки и охлаждения стерилизуемого продукта во время замеров, которые повторяют через равные промежутки времени (ти). Затем полученное значение ^kF умножают на ти и получают фактическое значение стерилизующего эффекта [84, 97]. Такие расчеты полу- 26
чили широкое распространение в пищевой промышленности, особенно при стерилизации консервов.
Во многих расчетах, особенно в микробиологической и медицинской промышленности, пользуются моделями, в которых исходным является уравнение
dN/dx = — kN,
(4)
где ft-—удельная скорость гибели микроорганизмов; N'—количество жизнеспособных микроорганизмов в момент времени т.
Интегрирование уравнения (4) в пределах от 0 до т и от No до N приводит к решению:
т
V = In N„/N = J kdx, о
а при k = const (Z = const) N = No exp (—kx), (5)
где jV0 — количество жизнеспособных спор в популяции до стерилизации (т=0).
Исходя из приведенного выше определения Тю,
Wo= (/Ш=ехр(-£т10), (б)
откуда
tw = 2,303/6 (7)
В некоторых случаях [31—33] термоустойчивость
микроорганизмов выражают соотношением
Дс = «т10, (8)
n=lgN0/N (см. табл. 6). (9)
Следовательно, фактор инактивации n (V = 2,3n) характеризует долю инактивированных при стерилизации микроорганизмов. Например, при п — 8 после стерилизации останется 1•10~8 часть от исходного числа жизнеспособных микроорганизмов данной популяции.
В пищевой промышленности требуемое значение п определяют из выражения
где %о — концентрация спор в продукте, например в питательной среде, до стерилизации; V—объем среды; Ь — допускаемый брак по нестерильности (в процентах).
Все приведенные соотношения для расчета эффективности режимов стерилизации предполагают, что процесс отмирания микроорганизмов во времени под влиянием повышенной температуры описывается экспоненциальной зависимостью. Однако накопленные экспери- 27
ментальные данные свидетельствуют о том, что кривые отмирания часто имеют S-образный вид [42, 72, 134„ 161]. Биологический смысл таких кривых еще не выяснен [134], однако при рассмотрении различных задач, в том числе и при построении математических моделей процесса инактивации, нужно учитывать отклонение от экспоненты, требующее ужесточения режима стерилизации.
В каких же условиях получают в основном экспоненциальную зависимость между количеством оставшихся жизнеспособных микроорганизмов и продолжительностью воздействия повышенной температуры? Как правило, это лабораторные данные, полученные в опытах на разбавленных гомогенных суспензиях микроорганизмов. Более того, в лабораториях часто экстраполируют полученные результаты, не доводя опыт до полного отмирания микроорганизмов данной популяции, а ограничиваясь лишь существенным снижением количества жизнеспособных спор. В то же время именно наличие весьма термостойких спор на практике часто приводит к нестерильным операциям, а отработанные в лабораториях режимы нуждаются в существенной корректировке в сторону увеличения температуры и (или) продолжительности ее воздействия. Это положение подтверждается уже полученными результатами [72]. На их основе предлагается оценивать стерилизующий эффект по уточненной формуле
где хп — поправка, учитывающая отклонение кривой гибели от экспоненты.
Условия, при которых учитывается выживаемость наиболее термостойких спор, следующие:
где ЛГТ ■— общее число спор в каждой серии образцов на одно соотношение температуры и времени; т — число образцов (например, капилляров); с — число спор в одном образце (капилляре).
В случаях, не отвечающих требованиям неравенств, величина Fc не будет учитывать выживаемость наиболее термостойких спор [72]. Таким образом, при отработке режимов стерилизации вопрос о представительности исследуемых выборок является исключительно важным.
В пищевой промышленности иногда для расчета, смертельного времени пользуются следующим соотношением [35]:
где аир — величины, характеризующие термическую устойчивость микроорганизмов и определяемые экспериментально.
В литературе еще довольно часто встречаются расчеты, основанные на характеристике
Qi0:
где kt И Т( — удельная скорость гибели спор и время гибели при температуре t соответственно; &і+ю и Tf+io — то же при температуре Г+10.
Величина Qio связана с величиной Z следующим соотношением [64]:
Принципиальный недостаток «метода Qi0» заключается в том, что величина Qio не является константой, а зависит от температуры:
In Qlo = 10£/(ЯЛ + ю^Г), (16),
где Е — энергия активации; R — газовая постоянная; Т — температура стерилизации, 7’=/+273.
Необходимо показать связь величины Z с величинами Qio и Fc:
где t2 и б — конечная и начальная температуры материала; вовремя десятикратного уменьшения популяции; t10=10tjo; Fe— время термической гибели; Fi = 10 F2.
Большое число опубликованных в литературе моделей процессов инактивации базируется на знании удельной скорости гибели спор и энергии их активации. При этом учитывается экспериментально установленный факт, свидетельствующий о том, что изменение, удельной скорости гибели спор в зависимости от температуры описывается известным уравнением Аррениуса [106, 107, 123].
где А — стерический фактор (коэффициент).
Для двух отличающихся значений температуры: 
В этом виде уравнение Аррениуса можно использовать непосредственно для того, чтобы рассчитать энергию активации реакции гибели клеток.
Следует подчеркнуть, что уравнение Аррениуса является эмпирическим. Более точная зависимость скорости реакции гибели от термодинамических параметров определяется уравнением Эйринга, лежащим в основе теории абсолютных скоростей химических реакций:
В приведенных уравнениях:
■къ— константа Больцмана; ha —постоянная Планка; kF*— изменение свободной чнрпгии активации; Д5* — изменение энтропии активации реакции; — затраты энергии активации.
Подставляя уравнение (21) в (20), получаем:
Сравнение уравнения (22) с уравнением (18) показывает, что член А непосредственно связан с энтропией активации реакции:
а член
— с энергией активации
Попытки использовать теорию скоростей химических реакций применительно к процессам инактивации микробных популяций предпринимались неоднократно {106, 107, 127]. При этом полагают, что в клетке происходит тепловая инактивация какой-то одной молекулы. Именно такой подход позволяет рассматривать процесс тепловой.гибели как реакцию 1-го порядка. Для этого имеются следующие основания. Прямая активация молекул под воздействием тепловой энергии приводит к разрушению молекулярных связей, а по теории химических реакций подобный процесс протекает как реакция 1-го порядка. В условиях избытка кислорода при тепловой стерилизации реакцию окисления в микробных клетках можно рассматривать как реакцию 1-го порядка. Реакция с водой или парами воды в микробной клетке при значительном избытке 30
одного компонента (в данном случае воды) протекает не как реакция 2-го порядка, а подобно реакции окисления, как реакция 1-го порядка. В подобном случае в уравнении скорости снижения концентрации непрореагировавших молекул См можно вместо См подставить N (количество сохранившихся жизнеспособных микробных клеток), о чем уже сказано выше.
Хорошо известен подход к изучению гибели популяции в тепловом поле, основанный на теории мишеней [127]. По этой теории предполагается существование фактора гибели в виде частиц, подобных радиоактивным. В процессе стерилизации эти частицы сталкиваются с мишенями в структуре чувствительных К НИМ' молекул ДНК, РНК или белка, вызывают глубокие изменения или разрушения, которыми и обусловливается гибель микробных клеток. Таким фактором гибели могут быть, например, молекулы воды в микробных клетках, возбужденные тепловой энергией.
Если р — вероятность того, что в результате столкновений частиц фактора гибели с мишенями в какой- то промежуток времени т в каждой из мишеней произойдут изменения, в результате которых будут утрачены нормальные функции, то, взяв вместо числа клеток число мишеней, получим:
где Лф —число мишеней в момент времени т=0, т. е. до начала стерилизации; М' — то же, но в момент времени т>0.
Если принять, что в каждойклетке существует только одна мишень, т. е,
В этом случае p = k рассматривают как временную плотность вероятности возникновения условий гибели.
Соответствующие модели гибели клеток можно получить и на основе наиболее вероятного распределения клеток по теплоустойчивости [127].
Допустим, что теплоустойчивость клеток, измеряемая в данном случае продолжительностью их выживания при определенной температуре, варьирует от ті до Tj. Следовательно, по теплоустойчивости все клетки популяции можно разделить на і-групп и количество клеток, относящихся к каждой группе, обозначить как 3£
Таким образом, применение самых различных подходов, основанных на различных представлениях о механизме гибели клеток под действием теплового ПОЛЯ, позволяет получать экспоненциальную зависимость между количеством жизнеспособных клеток и продолжительностью их нагревания. Этот факт и многочисленные экспериментальные данные послужили основанием считать, что гибель клеток в тепловом поле подчиняется логарифмическому закону [35, 85, 123, 145, 148]. Однако не менее многочисленные данные указывают на существенное отклонение кривых гибели от экспоненты, особенно если речь идет о значительных объемах популяций микроорганизмов в производственных условиях. Существует множество объяснений S-образного характера кривых отмирания, предполагающих наличие периодов задержки (фаза термоактивации) и отмирания наиболее термостойких спор [42, 72, 99, 134]. Остановимся лишь на наиболее существенных из них, а также на соответствующих моделях.
Известно, что для микробных популяций характерно явление гетерархии, т. е. гетерогенности по многим признакам (продолжительность жизни, прочность обо- 32
лочек, размеры клеток и др.), например гетерогенность по признаку терморезистентности, которая в наибольшей степени проявляется у культур, полученных при периодическом культивировании. Таким образом, в условиях промышленного производства при выборе конкретных условий термической инактивации той или иной микробной популяции мы встречаемся с необходимостью уничтожения спор, значительно отличающихся по своей устойчивости к тепловому воздействию. Из результатов, полученных некоторыми авторами [72], следует, что продолжительность отмирания наиболее термостойких спор в несколько раз больше, чем основной массы спор, отмирающих по экспоненте. Наибольшая разница наблюдается при температуре, характерной для работы установок непрерывной стерилизации (£>130° С), что выдвигает повышенные требования к выбору режимов, предусматривающих сочетание температуры и времени ее воздействия на жидкость.
К причинам, приводящим к отмиранию популяций по S-образной кривой, относится и наличие не единичных клеток, а их конгломератов. На этот факт редко обращают внимание в условиях лабораторий, однако в промышленном производстве практически не встречаются популяции, состоящие из единичных микроорганизмов. Учет конгломератов совершенно необходим при обосновании режимов стерилизации.
Конгломераты клеток в реальных популяциях, как правило, покрыты пленками поверхностно-активных веществ (ПАВ), непосредственно влияющих на достижение требуемой температуры во всех точках стерилизуемого объема. На устойчивость микроорганизмов к действию тепла влияет не только их вид, но и возраст, а также условия культивирования, в отдельных случаях—и условия временного хранения и транспортировки.
В последние годы появились сообщения о том, что погибающие клетки выделяют в окружающую среду- специфические защитные вещества [134].
По-видимому, отмеченные причины, приводящие к отклонению от экспоненциального закона отмирания популяций, не являются исчерпывающими. Однако они свидетельствуют о разном характере и многообразии факторов, в той или иной мере влияющих на выживае- 2—1950 33
мость микробных популяций и вызывающих значительные трудности при создании математических моделей.
Так называемая полимолекулярная модель позволяет анализировать два возможных типа изменений в процессе гибели. Предполагают, что гибель клеток происходит в результате разрушения какого-то одного вида молекул в этих клетках. Причем, если в клетке содержится т' подобных молекул, то клетка погибает при разрушении п' из этих т' молекул. Если р — вероятность того, что в какое-то время т будет разрушено какое-то число из т' молекул, a q — вероятность того, что они останутся неразрушенными, процесс можно описать следующим уравнением:
чения константы скорости реакции разрушения молекул. Вероятность разрушения молекул:
а вероятность того, что молекулы не
будут разрушены
В соответствии с этими уравнениями для системы (т', п'), в которой гибель клетки обусловлена разрушением п' из т' имеющихся молекул, можно записать:
где Пі'Сі — коэффициент.
Из уравнения (32) выведено несколько уравнений для частных случаев:
Применение полимолекулярной модели существенна усложняется, когда гибель микроорганизмов связана с. разрушением более чем одного типа молекул [127].
Достаточно известной является так называемая модель активированных спор [123]. Ее автор, анализируя 34
кривую гибели микроорганизмов, полученную при исследовании процесса тепловой гибели спор, связывает характер кривой с явлением замедленной гибели в начале нагревания. В самом деле, если k\ — константа скорости перехода спор (N) в состояние прорастающих спор (N*), a k2—константа скорости гибели спор, то
При более детальном исследовании этой модели полагают [151], что начальная популяция спор является смесью активированных (способных к образованию колоний на подходящей среде) и неактивированных спор No, которые при нагревании переходят (причем обратимо) в активированное состояние (Ао). Активированные споры под действием тепла погибают (удельная скорость гибели k), а неактивированные переходят в активное состояние (константа скорости тепловой активации а) и только после этого погибают (удельная скорость гибели k). Наглядно это представлено на схеме 1.
Схема 1. Упрощенная модель инактивации с учетом активированных бактериальных спор.
личество неактивированных спор в момент времени т, изменяющееся в соответствии с уравнением
Изменение и Lx во времени можно представить
При сравнении результатов опытов с Вас. stearo- X thermophilus и расчетных данных форма кривых оказалась сходной, однако наблюдались и некоторые различия [151]. Для устранения последних авторы попытались отказаться от допущения, что тепловая активация
протекает по типу реакции 1-го порядка.
Известны модели, учитывающие влияние предыдущего нагревания (история нагревания) в ходе стерилизации на тепловую чувствительность клеток. Стерилизуемую систему можно разделить на две следующие системы [127]:
— без предыстории нагревания — предшествующий нагрев в процессе стерилизации не влияет на тепловую чувствительность жизнеспособных клеток;
— с предысторией нагревания — предшествующий нагрев в процессе стерилизации влияет на тепловую чувствительность жизнеспособных клеток. Схема процесса представляется следующей:
Отмечено, что наиболее простой моделью системы с предысторией нагревания является модель, соответствующая реакции 1-го порядка. Теплочувствительность 36
молекул при нагревании изменяется от более низкой (2?т) к более высокой (5т) или наоборот, но так или иначе молекулы оказываются с измененной тепловой чувствительностью (рт). На приведенной выше схеме k — константа скорости изменений в направлении DT, X — константа скорости изменений в направлении
или R.
Данная модель позволяет описать процесс следующим обпазом:
Такую модель использовали для анализа тепловой гибели спор Вас. subtilis var. niger и Вас. pumilus. Заслуживает внимания тот факт, что предположительная величина энергии активации соответствует реакции на этапе X. Для спор Вас. pumilus этап %s представляет собой реакцию защитного (в отношении нагревания) типа *
~ при этом энергия активации
составляет около 105 кДж/моль. Для спор Вас. subtilis var. niger этап hR,наоборот, является реакцией тепловой активации
при этом энергия актива
ции составляет 264,6 кДж/моль. Данная величина соответствует энергии активации тепловой денатурации белков (в отличие от энергии активации этапа Z у Вас. pumilus).
В приведенном выше рассуждении ys и уд принимались равными единице, однако подобный анализ возможен и при меньших значениях этих показателей, т. е. в случае, когда в исходной популяции спор одновременно присутствуют S- и Д-споры. Некоторые авторы исходят из допущения, что в процессе термической 37
стерилизации среди активированных спор содержатся споры, уже инактивированные теплом [127, 151].
Если индекс і указывает на активированные споры, то общее количество жизнеспособных спор во время т определяется как общая сумма Np
На схеме 2 показана выведенная из этой общей мо-
Схема 2. Упрощенная модель интерпретации процесса тепловой гибели спор Вас. subtil is.
дели упрощенная модель, которая оказалась пригодной для интерпретации процесса тепловой гибели спор Вас. subtilis. В этой упрощенной модели Ni и N2—исходные группы спор, находящихся соответственно в активном и менее активном состоянии. Это согласуется с предположением о том, что в целом гибель жизнеспособных спор Вас. subtilis (Лг = Аі + А2 + Лг3) описывается логарифмической зависимостью, а неповрежденных (неинактивированных) спор (N* = Ni+N2)—S-образной кривой. Величины N и N* во время т определяются по следующим уравнениям:

В некоторых работах рассмотрены выводы рациональных, по мнению авторов, моделей, основанных на знании кинетики химических реакций и служащих для определения условий выживания микробов при термическом воздействии [44, 100, 149, 150].
Учитывая сложность процессов, протекающих в стерилизуемой жидкости, большинство исследователей принимают некоторые допущения, упрощающие задачу создания моделей. Часто постулируют, что причиной гибели клеток являются необратимые химические реакции, среди которых особое место занимают реакции денатурации некоторых белковых компонентов. При этом признают, что распределение точек, в которых протекают реакции денатурации, носит статистический характер [44, 113]. Затем рассматривают идеализированный вариант: летальные условия однородны во всех частях стерилизуемой жидкости, а размножение клеток в созданных условиях подавляется мгновенно и полностью. Допускают, что популяция клеток гомогенна, а статистические особенности, имеющие отношение к гибели клеток, одинаковы во всех точках стерилизуемого объема жидкости. Весьма рационально предположение р том, что гибель каждой отдельной клетки является независимым событием. Следовательно, на вероятность гибели данной клетки в определенный момент времени не влияют те процессы, которые происходят в других клетках в данный момент или уже произошли до рассматриваемого момента времени.
Учитывая изложенные выше допущения, для описания процесса термической инактивации микроорганизмов можно использовать закон Пуассона [44]. Этот закон в теории вероятностей является предельным случаем биномиального распределения случайной величины и обеспечивает высокую точность, при которой число испытаний бывает относительно велико, а вероятность появления ожидаемого события в одном испытании весьма мала.
Применительно к процессу отмирания _ микробной популяции под действием повышенной температуры вероятность Pd возникновения d реакций денатурации в некотором объеме. V (при среднем числе п центров реакций в единице объема) может выражаться формулой
где d— целое число (d—X, 2, 3...).
Если для инактивации клеток необходимо, чтобы в объеме возникло не менее d реакций, то все клетки, в которых возникнет реакций меньше (с?—1, d—2, d—3,...), останутся неповрежденными. Очевидно, что при d=i неповрежденными или выжившими останутся только те клетки, в которых не возникнет ни одной реакции денатурации. Вероятность наличия таких клеток равна: і
При d—2 останутся неповрежденными клетки, в которых возникло не более одной реакции денатурации, а вероятность наличия таких клеток равна:
Аналогичные рассуждения можно повторить для d = 3, 4, 5... Клетки, которые поражаются при d = l, условно называют одноударными, при d — 2— двухударными и т. д. При «многоударном» порядке отмирания популяции микроорганизмов (d = m) суммарная вероятность наличия неповрежденных клеток равна:
Число погибших клеток равно произведению вероятностей на общее число объектов (Ао) в объеме жидкости, подвергаемой тепловой обработке. Очевидно, что число выживших клеток равно произведению No на Лї-і:

Первый член в уравнении (55) определяет число выживших одноударных клеток, у которых в объеме V не произошли реакции денатурации:
т. е. описывает кинетику отмирания микроорганизмов по экспоненте. Однако уравнение (57) не отражает зависимости количества выживших клеток от времени, поэтому его нельзя использовать в практических расчетах. Для выяснения этой связи целесообразно записать уравнение (57) в следующем виде [441:
Величина, стоящая в левой части уравнения (58), представляет собой безразмерный критерий стерилизации V, характеризующий убыль количества жизнеспособных микроорганизмов во времени т в зависимости от температуры и констант термической устойчивости микроорганизмов:
— для переменной температуры стерилизации
о
— для постоянной температуры стерилизации
С учетом приведенных соотношений изменение количества жизнеспособных клеток в объеме жидкости в зависимости от режимов тепловой обработки (сочетание Т и т) и термической устойчивости микроорганизмов записывается следующим образом:
— для случая с переменной температурой стерилизации:
— для случая с постоянной температурой стерилизации:
Уравнение (62) можно использовать для расчета количества жизнеспособных клеток в момент времени т в заданном объеме жидкости, подвергаемой стерилизации в изотермических условиях. Оно помогает уточ- 41
нить тип (одноударный или многоударный) обрабатываемой популяции микроорганизмов. Это можно сделать путем сравнения экспериментальных данных, отражающих зависимость N = f(x), с данными, рассчитанными по уравнению (62) для одноударных, двухударных и т. д. систем. Очевидно, что при 4=1 кривая гибели клеток представляет собой экспоненту, а при дальнейшем увеличении d до 2, 3, 4 и т. д. она все больше будет отклоняться от экспоненциальной и приближаться к S-образной кривой.
Необходимо отметить, что значения удельной скорости гибели спор зависят только от температуры, а при изотермических условиях стерилизации могут рассматриваться в качестве постоянных величин, характерных для популяции данного вида. Это является одним из преимуществ рассмотренной модели.
В настоящее время получены данные об изменении в зависимости от температуры значений удельной скорости гибели спор некоторых микроорганизмов, наиболее часто встречающихся в микробиологических производствах и отличающихся повышенной термической устойчивостью (Вас. megatherium, Вас. mycoides, Вас. stearothermophilus и др. [123, 145]).
Существование в реальных производствах S-образных кривых отмирания микроорганизмов в тепловом поле свидетельствует о необходимости тщательной оценки пригодности режимов стерилизации, разработанных и оказавшихся высокоэффективными в лабораторных условиях. Часто лабораторные эксперименты, проводимые с разбавленными гомогенными суспензиями микроорганизмов, указывают на экспоненциальную кривую отмирания, в то время как на производстве отмечается S-образная кривая. На практике это приводит к необходимости «ужесточения» лабораторных режимов, увеличению времени выдержки (при t = const) в 2—3 раза или к соответствующему повышению температуры (TB = const) (см. [44, 50, 51]).
Приведенный далеко не полный перечень моделей, описывающих процесс инактивации микроорганизмов под действием тепла, позволяет не только представить их многообразие и разнохарактерность (различие в исходных посылках, используемом математическом аппарате и др.), но и оценить некоторые возможности прак- 42
тического использования таких моделей, особенно в практике микробиологических производств. Следует отметить, что тип применяемой модели, как правило, не только зависит от существа решаемой задачи и ее особенностей, но и имеет черты установившихся традиций. Примером служат соотношения, широко используемые в практике пищевых производств и базирующиеся на понятиях смертельного времени, стерилизующего эффекта и др. Такие соотношения основаны преимущественно на предположении о том, что гибель клеток в тепловом поле подчиняется экспоненциальному закону, хотя в последние годы накоплено много фактов, указывающих на весьма частые отклонения кривых отмирания от экспоненты. Попытки учесть эти данные в традиционно существующих моделях заслуживают сами по себе по* ложительной оценки, однако не устраняют другой существенный недостаток, заключающийся в том, что. для характеристики термоустойчивости микроорганизмов и эффективности режимов стерилизации используются показатели (Z, тю), зависящие от температуры.
В то же время очевидно, что переход к достаточно строгим количественным зависимостям предъявляет повышенные требования к характеристикам термоустойчивости. В этом плане более предпочтительны методы, основанные на изучении термодинамических характеристик микроорганизмов, например энергии активации. Следует, однако, подчеркнуть, что уравнения, заимствованные из кинетики химических реакций и примененные к процессу инактивации микроорганизмов, теряют ясный физический смысл, а соответствующие константы термоустойчивости можно рассматривать лишь как формализованные характеристики, полезные для анализа процессов термической стерилизации.
3.2.