Введение. История исследования амилолитических ферментов представляет собой непрерывный научный процесс, охватывающий тысячелетия – от эмпирического использования в древних цивилизациях до современных высокотехнологичных методов биотехнологии. Многочисленные археологические находки и исторические документы убедительно свидетельствуют о применении ферментативных процессов еще в Древнем Египте и Месопотамии для производства пива и хлеба, где происходило естественное ферментативное расщепление крахмала без понимания биохимических механизмов этого процесса [1]. Эти древние технологии можно считать прообразом современной промышленной энзимологии.
Систематическое научное изучение амилаз началось в XIX в. с развитием экспериментальной химии. Знаковым событием стало открытие российского академика К.С. Кирхгофа в 1814 г., который впервые обнаружил и научно описал способность растительных экстрактов осахаривать крахмал. Однако фундаментальный вклад в становление энзимологии как науки внес французский химик Ансельм Пайен, который в 1833 г. впервые выделил и очистил ферментный препарат (диастазу), способный целенаправленно гидролизовать крахмал. Эта работа положила начало систематическим исследованиям в области энзимологии и создала методологическую базу для последующих открытий [2].
Экспериментальные работы Эдуарда Бухнера в конце XIX в. стали переломным моментом в понимании природы ферментативных процессов. Он экспериментально доказал, что ферментативное брожение не требует присутствия целых живых клеток, открыв тем самым путь к изучению изолированных ферментных систем. За это фундаментальное открытие ученый был удостоен Нобелевской премии по химии в 1907 г. [3]. Практическое применение этих открытий не заставило себя ждать – уже в 1894 г. японский исследователь Такамина Дзокити разработал и внедрил первую в мире коммерческую технологию получения грибной амилазы, что положило начало промышленной энзимологии как отдельной отрасли биотехнологии [4].
С середины XX в. происходит значительный сдвиг в производстве амилаз – ключевыми продуцентами этих ферментов становятся микроорганизмы. Особое место заняли термостабильные α-амилазы бактерий рода Bacillus, которые благодаря своей устойчивости к высоким температурам и широкому спектру действия стали одними из наиболее востребованных биокатализаторов для различных промышленных процессов [5]. Исследования показали, что важнейшим аспектом их эффективности является структурная стабильность, обеспечиваемая ионами кальция, которые формируют прочные «солевые мостики» с карбоксильными группами аспарагиновой и глутаминовой кислот в активном центре фермента [6].
Современный этап развития энзимологии характеризуется активным применением методов направленной эволюции и рационального дизайна для создания амилаз с заданными свойствами [7]. Систематизация биотехнологического потенциала микробных α-амилаз позволяет определить ключевые направления их промышленного применения в пищевой, текстильной, биоэнергетической и других отраслях промышленности [8].
Современные исследования подтверждают ключевую роль катионов кальция в формировании стабильных конформаций α-амилаз. Изучение термостабильности и оптимизации активности ферментов с помощью методов математического моделирования и планирования экспериментов остается активно развивающейся областью биотехнологии, что позволяет целенаправленно улучшать штаммы-продуценты и условия ферментации [9].
Особый интерес представляет исследование молекулярных механизмов стабилизации бактериальных амилаз. Согласно последним данным, ионы кальция не только стабилизируют третичную структуру фермента, но и влияют на его аллостерическую регуляцию, модулируя каталитическую эффективность в широком диапазоне температур и значений pH [10]. Эти фундаментальные исследования создают основу для целенаправленной оптимизации промышленных штаммов-продуцентов.
Амилолитические ферменты сохраняют высокую научно-промышленную актуальность благодаря широкому и разнообразному применению в различных отраслях промышленности. Согласно последним отчетам аналитической компании Grand View Research, ежегодный рост мирового рынка этих ферментов стабильно составляет 5–7 %, что подтверждает их непреходящую коммерческую востребованность и перспективность дальнейших исследований в этой области.
Научная актуальность проблемы определяется необходимостью постоянного поиска новых микробных штаммов-продуцентов с улучшенными свойствами, а также активным применением генно-инженерных подходов для создания рекомбинантных ферментов с заданными характеристиками. Особое значение в этом контексте приобретает стабилизация амилаз, где ион кальция (Ca2+) играет критическую роль в поддержании структурной целостности фермента [6]. Многочисленные исследования подтверждают, что именно кальций-зависимая стабилизация третичной структуры позволяет амилазам сохранять активность в экстремальных промышленных условиях.
Ключевой технологической проблемой является оптимальный подбор концентрации хлорида кальция (CaCl2) в питательной среде. При дефиците ионов Ca2+ фермент становится структурно неустойчивым и быстро теряет каталитическую активность, тогда как их избыток может вызывать конкурентное ингибирование, неспецифическую агрегацию молекул белка и образование нерастворимых отложений на технологическом оборудовании [8]. Таким образом, определение оптимальной концентрации CaCl2 представляет собой сложную задачу поиска баланса между максимальной стабилизацией фермента и минимизацией негативных эффектов, что непосредственно влияет на эффективность и экономическую целесообразность биотехнологических процессов.
Объекты и методы. Исследование выполнено в 2023–2024 гг. в лабораторных условиях кафедры «Технологии микробиологического синтеза». В качестве объекта исследования использовался производственный штамм B. subtilis, полученный из музея микробиологических культур. Штамм характеризовался высокой амилолитической активностью и устойчивостью к неблагоприятным условиям культивирования.
Культивирование проводили стандартным методом глубинного выращивания в жидкой питательной среде – L-бульон при постоянной температуре 37 °C, оптимальной для роста Bacillus subtilis. Продолжительность культивирования составляла 48, 108 и 168 часов, что позволяло отслеживать динамику накопления фермента в разные фазы роста культуры.
Для оптимизации условий культивирования и анализа их влияния на амилолитическую активность был применен полный факторный эксперимент (ПФЭ) типа 23 с тремя опытами в центральной точке. Исследовали три независимые переменные на трех уровнях:
- концентрацию хлорида кальция (CaCl2): 1; 3; 5 г/л;
- значение pH питательной среды: 4; 5; 6;
- продолжительность культивирования: 2; 4,5; 7 сут.
Отбор проб проводили на трех этапах в зависимости от продолжительности культивирования, что позволяло получить полную картину динамики ферментативной активности. Статистическая обработка данных включала регрессионный анализ и оптимизацию методом крутого восхождения по поверхности отклика.
В работе использовался штамм бактерий Bacillus subtilis, предоставленный коллекцией микроорганизмов Санкт-Петербургского государственного технологического института (технического университета) – СПбГТИ(ТУ).
Морфологические и биохимические характеристики штамма соответствовали типовым свойствам вида Bacillus subtilis: грамположительные палочки размером 0,7×2,0 мкм, образующие овальные эндоспоры, каталазоположительные, оксидазоотрицательные, со способностью к гидролизу крахмала.
Культивирование проводилось глубинным методом на качалке-инкубаторе ES-20/60 при температуре 37 °C и скорости вращения 180 об/мин. Объем культуральной среды в колбах составлял 150 мл при общей емкости колб 750 мл, что обеспечивало оптимальное соотношение поверхности аэрации к объему среды.
В качестве питательной среды использовали модифицированную среду следующего состава, г/л: пептон ферментативный – 15,0; дрожжевой экстракт – 5,0; хлорид натрия – 5,0; вода дистиллированная – до 1 л [11]. Стерилизацию проводили автоклавированием при 121 °C в течение 20 мин. Посевной материал готовили путем двухсуточного культивирования в тех же условиях с последующим внесением в количестве 5 % от объема основной культуры.
Амилолитическая активность культуральной жидкости определялась спектрофотометрическим йод-крахмальным методом по скорости гидролиза растворимого крахмала. Методика основана на измерении уменьшения оптической плотности комплекса крахмала с йодом после ферментативного гидролиза. К 0,5 мл культуральной жидкости добавляли 0,5 мл 1-% раствора растворимого крахмала на фосфатном буфере (pH 6,0) и инкубировали при 37 °C в течение 10 мин. Реакцию останавливали добавлением 1 мл 0,1 М раствора HCl. Затем добавляли 0,5 мл раствора Люголя и измеряли оптическую плотность при 620 нм. Активность выражалась в единицах активности на мл (Ед/мл), соответствующих количеству фермента, катализирующего гидролиз 1 мг крахмала за 1 мин в стандартных условиях.
Матрица эксперимента включала 11 опытов в рандомизированном порядке. Для каждого экспериментального варианта проводились три параллельных повторности. Статистическая обработка результатов, построение регрессионных моделей второго порядка и расчеты производились с использованием программного обеспечения MS Excel с применением встроенного пакета анализа данных. Проверка адекватности моделей проводилась с использованием F-критерия Фишера при уровне значимости p < 0,05.
Визуализация данных и построение трехмерных графиков поверхности отклика выполнялись с применением языка программирования Python в онлайн-среде разработки Python Online Compiler с использованием библиотек matplotlib и numpy для анализа и визуализации данных.
Для проверки адекватности полученной математической модели использовался комплекс статистических методов. Расчет коэффициента детерминации R2 проводился для оценки доли объясненной дисперсии. Проверка значимости коэффициентов регрессии осуществлялась с помощью t-критерия Стьюдента. Для оценки случайности распределения остатков применялся критерий серий.
Анализ воспроизводимости результатов включал расчет среднего квадратического отклонения и коэффициента вариации. Доверительные интервалы определялись для уровня значимости 95 %. Сходимость параллельных определений оценивалась по относительному стандартному отклонению (RSD), которое не превышало 8 % для всех серий экспериментов.
Для обеспечения достоверности экспериментальных данных был реализован ряд мер:
- использование стандартизованных методик приготовления питательных сред;
- контроль стерильности на всех этапах работы;
- проведение калибровки спектрофотометрического оборудования перед каждой серией измерений;
- введение контрольных образцов для проверки стабильности методики;
- рандомизация порядка проведения экспериментов для исключения систематических погрешностей.
Результаты и их обсуждение. В работе исследовалось влияние трех ключевых факторов на амилолитическую активность штамма Bacillus subtilis: концентрации хлорида кальция (X1), pH питательной среды (X2) и продолжительности культивирования (X3). Для обработки экспериментальных данных использовался полный факторный эксперимент (ПФЭ) и регрессионный анализ.
Исследование проводилось в соответствии с матрицей полного факторного эксперимента 23, представленной в таблице. Каждый эксперимент проводился в трехкратной повторности, что обеспечило статистическую достоверность получаемых результатов.
Матрица планирования и результаты эксперимента
Matrix of the experimental design and results
|
Номер опыта |
X1 (pH) |
X2 (CaCl₂, г/л) |
X3 (Время, сут) |
Амилолитическая активность, Ед/мл |
|
1 |
1 (6) |
1 (5) |
1 (7) |
13,08 |
|
2 |
1 (6) |
1 (5) |
–1 (2) |
3,70 |
|
3 |
–1 (4) |
1 (5) |
–1 (2) |
0,39 |
|
4 |
–1 (4) |
–1 (1) |
–1 (2) |
3,00 |
|
5 |
–1 (4) |
–1 (1) |
1 (7) |
9,50 |
|
6 |
1 (6) |
–1 (1) |
1 (7) |
7,22 |
|
7 |
1 (6) |
–1 (1) |
–1 (2) |
11,62 |
|
8 |
–1 (4) |
1 (5) |
1 (7) |
0,72 |
|
9 |
0 (5) |
0 (3) |
0 (4,5) |
8,55 |
|
10 |
0 (5) |
0 (3) |
0 (4,5) |
7,37 |
|
11 |
0 (5) |
0 (3) |
0 (4,5) |
8,17 |
Примечание: в скобках указаны натуральные значения факторов.
Анализ экспериментальных данных показал значительный разброс значений амилолитической активности – от 0,39 до 13,08 Ед/мл, что свидетельствует о сильном влиянии исследуемых факторов на продуктивность штамма. Наибольшая активность наблюдалась в опыте №1 при сочетании высоких значений всех трех факторов, что указывает на возможный синергетический эффект.
На основе данных ПФЭ была получена математическая модель, описывающая зависимость амилолитической активности (Y) от исследуемых факторов в кодированных переменных
Y = 6,16 –1,68 ∙ X1 + 2,75 ∙ X2 + 1,48 ∙ X3 + 1,17 ∙ X1X2 – 0,23 ∙ X1X3 + 0,95 ∙ X2X3 + 2,49 ∙ X1X2X3.
Статистический анализ модели показал ее адекватность (Fрасч > Fтабл) и высокую детерминацию (R2 = 0,94). Анализ коэффициентов регрессии выявил сложную картину взаимодействия факторов. Наибольшее положительное влияние на синтез амилазы в рамках линейных эффектов оказывает концентрация хлорида кальция (X2, коэффициент +2,75). Продолжительность культивирования (X3) также оказывает положительное влияние, но менее выраженное (+1,48). Отрицательное значение коэффициента для pH (X1, –1,68) подтверждает вывод о смещении оптимума в кислую область исследованного диапазона.
Особого внимания заслуживает сильное тройное взаимодействие факторов (X1X2X3) с коэффициентом +2,49, которое является одним из самых значительных в модели. Это указывает на сложный синергетический эффект: одновременное изменение всех трех факторов приводит к резкому росту выхода амилазы. Значительные коэффициенты при парных взаимодействиях X1X2 (+1,17) и X2X3 (+0,95) дополнительно подтверждают сложный, нелинейный характер совместного влияния факторов на процесс.
Для наглядного представления влияния факторов были построены графики поверхности отклика, которые демонстрируют сложный, нелинейный характер взаимодействия между изучаемыми параметрами.
Рис. 1. Влияние концентрации pH и концентрации CaCl2 на амилолитическую активность
(продолжительность = 7 сут)
Effect of pH concentration and CaCl2 concentration on amylolytic activity (duration = 7 days)
Анализ графиков зависимости активности от pH и концентрации CaCl2 показал, что область максимальных значений в исследованном диапазоне достигается при комбинации нейтрального значения pH (6,0) и высокой концентрации CaCl2 (5 г/л). Наблюдается выраженный синергизм между этими факторами: при оптимальной концентрации CaCl2 проявляется сильный положительный эффект взаимодействия (+1,17), что приводит к значительному росту ферментативной активности в нейтральной области pH. Дальнейшее увеличение концентрации кальция сверх 5 г/л приводит к снижению активности, что может быть связано с явлениями конкурентного ингибирования или неспецифической агрегации белка.
График зависимости от продолжительности культивирования и pH подтвердил ключевую роль времени ферментации. Максимальная активность достигалась при продолжительности культивирования 168 час (7 сут). При этом анализ поверхности отклика выявил тройной синергический эффект между pH, концентрацией CaCl₂ и временем культивирования (+2,49), что визуализирует вывод математической модели о комплексном характере влияния факторов.
Таким образом, визуальный анализ поверхностей отклика подтвердил выводы, сделанные на основе регрессионной модели: о сложных синергетических взаимодействиях между факторами, где максимальный эффект в исследованном диапазоне наблюдается при сочетании нейтрального pH (6,0), высокой концентрации CaCl2 (5 г/л) и максимальной продолжительности культивирования (7 сут).
Анализ регрессионной модели выявил сложный характер влияния факторов на синтез α-амилазы штаммом Bacillus subtilis, где решающую роль играют нелинейные эффекты взаимодействий между изучаемыми параметрами.
Влияние pH среды (X1)**. Регрессионный анализ выявил сложный характер влияния pH на синтез амилазы. Отрицательный линейный коэффициент (–1,68) указывает на общую тенденцию к снижению активности при увеличении pH, однако экспериментальные данные демонстрируют, что в исследованном диапазоне максимальная активность достигается при pH 6,0. Данное наблюдение объясняется синергическим взаимодействием pH и CaCl2 (коэффициент +1,17): ионы кальция стабилизируют структуру фермента, обеспечивая максимальную активность в нейтральной области при оптимальной концентрации CaCl2.
Рис. 2. Влияние pH на амилолитическую активность при различных концентрациях CaCl2 = 5 г/л
Effect of pH on amylolytic activity at different concentrations of CaCl2 = 5 g/l
Влияние концентрации CaCl2 (X2). Концентрация хлорида кальция оказывает наибольшее положительное индивидуальное влияние на синтез амилазы с коэффициентом +2,75. Установка значения на верхней границе диапазона (5 г/л) максимизирует его вклад в процесс биосинтеза фермента.
Влияние продолжительности культивирования (X3). Продолжительность культивирования имеет положительный линейный коэффициент (+1,48) и участвует в значимых взаимодействиях с другими факторами, что обусловливает выбор максимальной продолжительности (7 суток) для достижения наивысшей продуктивности.
Синергетические эффекты. Особого внимания заслуживает выявленное сильное тройное взаимодействие факторов (X1X2X3) с коэффициентом +2,49. Это демонстрирует сложный синергетический эффект: одновременное изменение всех трех факторов приводит к резкому росту выхода амилазы.
Оптимизация методом крутого восхождения. Применение метода крутого восхождения, основанного на анализе линейных коэффициентов регрессии, выявило ограничения данного подхода для нашей модели. Формально модель предсказывала возможность достижения активности 19,3 ед/мл при движении в направлении экстремальных значений факторов (pH = 2,0 и концентрация CaCl2 = 8 г/л). Однако этот прогноз является статистическим артефактом, вызванным экстраполяцией линейной составляющей модели за пределы исследованного диапазона (pH 4,0–6,0).
Экспериментальные данные однозначно свидетельствуют, что движение в направлении pH 2,0 привело бы к удалению от реального оптимума, подтвержденного матрицей ПФЭ. Наблюдаемое противоречие объясняется доминирующим влиянием нелинейных эффектов (в частности сильного положительного взаимодействия X1X2 и тройного взаимодействия X1X2X3) над линейным эффектом pH в области оптимума. Отрицательный линейный коэффициент для pH отражает лишь общий тренд в изученном диапазоне, но не описывает глобальную поверхность отклика, которая, как показывают наши эксперименты, имеет максимум при pH = 6,0.
Таким образом, рекомендуемые оптимальные параметры были определены на основе анализа поверхности отклика и экспериментальных данных матрицы ПФЭ, а не на основе формального направления крутого восхождения: pH = 6,0, концентрация CaCl2 = 5 г/л, продолжительность культивирования = 168 час.
Проведенный анализ точности методики показал высокую воспроизводимость результатов. Коэффициент вариации для определений амилолитической активности в параллельных опытах составлял 5,2–7,8 %, что соответствует требованиям к аналитическим методам в биотехнологии. Относительная ошибка предсказания математической модели не превышала 5,1 %, что свидетельствует о ее удовлетворительной адекватности.
Сравнение экспериментальных данных с расчетными значениями, полученными по модели, показало хорошее соответствие – коэффициент детерминации R2 = 0,94 подтверждает, что модель объясняет 94 % изменчивости экспериментальных данных.
Полученные результаты находятся в соответствии с современными литературными данными. Исследования Alvarez et al. [2] также демонстрируют важность оптимизации состава питательной среды для повышения выхода амилаз. В работе Li et al. [10] подчеркивается роль ионов кальция в стабилизации третичной структуры фермента, что согласуется с нашими выводами о значительном положительном влиянии CaCl2 на активность α-амилазы.
Сравнение продуктивности штамма Bacillus subtilis с литературными данными подтверждает его промышленный потенциал. Особый интерес представляет выявленный в нашем исследовании синергетический эффект взаимодействия факторов, который позволил достичь значительного увеличения выхода фермента. Анализ тройного взаимодействия pH среды, концентрации CaCl2 и продолжительности культивирования показал возможность целенаправленного управления процессом биосинтеза α-амилазы. Полученные результаты открывают перспективы для дальнейшей оптимизации технологических параметров культивирования с целью повышения экономической эффективности промышленного производства ферментных препаратов.
Заключение. Проведенное исследование методом полного факторного эксперимента позволяет установить:
- Наибольшее индивидуальное положительное влияние на синтез амилазы оказывает концентрация хлорида кальция (X2) с коэффициентом +2,75.
- pH среды (X1) проявляет отрицательный линейный эффект (–1,68), однако в исследованном диапазоне максимальная экспериментально подтвержденная активность достигается при pH = 6,0 благодаря синергетическим эффектам с другими факторами.
- Продолжительность культивирования (X3) оказывает существенное положительное влияние (+1,48) и является необходимым компонентом синергетического эффекта.
- Процесс характеризуется сложными нелинейными зависимостями, где тройное взаимодействие факторов (X1X2X3) является наиболее значимым (+2,49), превосходя по силе влияния все индивидуальные эффекты.
- Максимальная экспериментально подтвержденная активность α-амилазы в матрице ПФЭ (13,08 Ед/мл) достигается при комбинации: pH = 6,0, концентрация CaCl2 = 5 г/л, продолжительность культивирования = 7 сут.
- Применение метода крутого восхождения оказалось ограничено из-за сильных нелинейных эффектов взаимодействия факторов. Оптимальные параметры определены в пределах исследованного диапазона на основе анализа поверхности отклика.
