Введение
Водоросли рода Porphyridium известны тем, что продуцируют фикобилипротеины, полиненасыщенные жирные кислоты и экзополисахариды [1]. Фикобилипротеины это группа окрашенных, водорастворимых белков, составляющих основной комплекс светособирающих пигментов цианобактерий, красных водорослей, глаукоцистофитов и криптофитов. Фикоэритрин является основным пигментом красных морских водорослей, составляющим 6080 % от общего количества растворимых белков [2]. Также пигменты представлены хлорофиллом а и каротиноидами (преимущественно бета-каротин и зеоксантин). Хлорофилл а основной фотосинтетический пигмент микроводорослей, а такие каротиноиды, как бета-каротин и зеоксантин, являются фотопротекторными пигментами и передают энергию на хлорофилл при фотоситезе [3].
Количественное содержание пигментов в культуре микроводорослей представляет ценную информацию о физиологическом состоянии, однако их определение достаточно сложная и трудоемкая задача [4]. Среди множества доступных современных методов наиболее используемыми являются микроволновая экстракция, ультразвуковая экстракция, гомогенизация под высоким давлением, экстракция сверхкритическими жидкостями и ускоренная экстракция растворителем (также называемая экстракцией жидкостью растворителем под давлением или субкритической экстракцией растворителем) [5; 6]. При работе с интенсивными культурами микроводорослей достаточно часто возникает необходимость определения пигментов, не прибегая к экстракции и не вмешиваясь в процессы роста культуры. Определение пигментов in vivo это быстрый метод без разрушения клеток. Одним из способов вычисления концентрации пигментов в нативном состоянии является анализ спектров поглощения культуры микроводорослей с помощью математического моделирования [7].
Таким образом, целью данной статьи является разработка математической модели спектра поглощения культуры Porphyridium purpureum, которая позволяла бы определять концентрации основных фотосинтетических пигментов: хлорофилла а, В-фикоэритрина и суммарных каротиноидов.
1. Материалы и методы исследования
Экспериментальные работы выполняли на базе кафедры "Физика" СевГУ. В работе использовалась красная морская водоросль Porphyridium purpureum (Bory) Drew et Ross (Rhodophyta). Культуру P. purpureum выращивали в накопительном режиме, используя среду для красных морских водорослей [8], в плоскопараллельном фотобиореакторе объемом 1,2 л, рабочая толщина культуры 2 см [9]. Нижняя грань фотобиореактора расположена под углом с целью улучшения перемешивания суспензий микроводорослей и цианобактерий. Сверху культиватор закрывался пластиковой крышкой, в которой было выполнено отверстие для подачи воздуха, а также был оснащен системой охлаждения ("водяной рубашкой"), обеспечивающей поддержание температуры.
В качестве источника освещения использовали горизонтальную световую решетку из светодиодных ламп LCD Feron LB-213 мощностью 10 Вт. В эксперименте освещенность составляла 5 клк или 17 Вт·м . Температуру стабилизировали на уровне 25 1 C. Барботаж культуры осуществляли аквариумным компрессором Hailea ACO-308 воздухом через аквариумный распылитель, представляющий собой пластиковую трубку длиной 5 см, диаметром 5 мм, у которой диаметр пор не превышает 0,1 мм. Скорость подачи воздуха составляла 1 л·л культуры в минуту. Дополнительного введения углекислого газа не производилось.
Отбор проб для определения оптической плотности проводили из разных точек внутри фотобиореактора: отбирали по 5 мл суспензии клеток водорослей, получая таким образом "среднюю пробу" объемом 30 мл. В средней пробе после перемешивания определяли коэффициент пропускания. Оптическую плотность рассчитывали по формуле: D = lg(T), где Т величина пропускания, определяемая на фотометре КФК-2 при длине волны 750 нм, погрешность измерения величины пропускания не превышала 1 %. Кюветы располагали максимально близко к фотоприемнику, что позволяло снизить ошибку измерения оптической плотности культуры, связанную со светорассеянием. Пробы c оптической плотностью выше единицы предварительно разбавляли свежей питательной средой, подбирая коэффициент разбавления таким образом, чтобы показания КФК-2 попадали в диапазон наименьшей погрешности (0,2 0,6 единиц).
Для определения сухого веса 510 мл суспензии центрифугировали в течение 10 минут при 3000 об·мин , сливали надосадочную жидкость, осадок промывали дистиллированной водой, повторно центрифугировали и сушили в течение суток при 55 C. В результате получена линейная зависимость сухого веса от оптической плотности D , а коэффициент пропорциональности составил г СВ·л ·ед. опт. пл .
Пробы для определения содержания пигментов отбирали ежедневно на различных фазах роста накопительной культуры после тщательного перемешивания. Концентрации хлорофилла а, В-фикоэритрина и суммарных каротиноидов (далее хл а, В-ФЭ, кр) определяли согласно стандартным методикам [1012]. После отбора проб из фотобиореактора проводили центрифугирование при 3000 об·мин , надосадочную жидкость сливали. Далее проводили несколько циклов "замораживание-оттаивание" в морозильной камере при -18 C. Затем В-ФЭ экстрагировали водным буферным раствором с рН 77,5. После изменения цвета осадка (из красного на зеленый), с помощью ацетона выделяли хл а. Спектры экстрактов пигментов записывали на спектрофотометре Unico в кварцевых 1 см кюветах в диапазоне длин волн 400 800 нм с шагом 0,5 нм. Регистрировали оптическую плотность полученных экстрактов в области характеристических максимумов поглощения В-ФЭ (545 нм), R-фикоцианина (615 нм), аллофикоцианина (650 нм), хл а (663 нм) и кр (480 нм). Концентрацию В-ФЭ определяли по [10], хл а [11], кр [12] с учетом объемов экстрактов и массы навесок:
(1)
(2)
(3)
где D , D , D , D , D значения оптических плотностей в соответствующей длине волны; экстинкция В-фикоэритрина 10 (л·г ·см ); экстинкция R-фикоцианина 7 (л·г ·см ); экстинкция аллофикоцианина 5,8 (л·г ·см ); экстинкция хлорофилла а 88,15 (л·г ·см ); экстинкция суммарных каротиноидов 2500 (л·г ·см ).
Спектры поглощения культуры регистрировались в диапазоне от 400 до 800 нм с шагом 0,5 нм на двухлучевом спектрофотометре Lambda 365 Double Beam UV-Visible (производитель: Perkin Elmer, Индия), который оснащен интегрирующей сферой (ИС) диаметром 60 мм (внешнее покрытие BaSO ). Для определения концентрации пигментов методом кривых Гаусса необходим истинный спектр поглощения, компенсированный на рассеяние, который определяли по формуле [13]:
(4)
где D( ;r) оптическая плотность образца, расположенного на некотором расстоянии r от ИС; L(r;0) поправочный коэффициент ослабления света;D( ;0) оптическая плотность образца при стандартном положении кюветы.
2. Результаты и их обсуждение
Экспериментальные результаты
На рис. 1, A представлена накопительная кривая роста культуры P. purpureum, которая характеризуется наличием ярко выраженных экспоненциальной и линейной фаз роста.
Рис. 1. Накопительная кривая роста культуры P. purpureum (А). Сплошные линии аппроксимации экспоненциальной (0 3 сутки) и линейной (3 8 сутки) фаз. Истинные спектры поглощения культуры P. purpureum, скомпенсированные на рассеяние (В)
Fig. 1. Batch growth curves of P. purpureum culture (А). Solid lines are approximations of exponential (0 3 days) and linear (3 8 days) phases. True absorption spectra of P. purpureum culture compensated for scattering (В)
Биомасса культуры за две недели эксперимента увеличилась почти в 50 раз и достигла 1,7 г СВ·л , что с учетом толщины фотобиореактора (1 см) составляет 17 г СВ·м . Экспоненциальная и линейная фазы роста культуры описывались выражениями:
(5)
где В начальная плотность культуры,г СВ·л ; максимальная удельная скорость роста, которая составила 0,7 сут ; В плотность культуры на момент начала линейного роста t ; P максимальная продуктивность, значение которой составило 0,21 г СВ·л ·сут .
Проведенные расчеты позволяют предположить, что рост порфиридиума начиная с третьих суток эксперимента был ограничен потоком углекислого газа. В работе [14] показано, что максимальная продуктивность культуры P. purpureum достигала 0,42 г СВ·л ·сут при такой же освещенности и рабочем слое культуры.
На рис. 1, B представлена динамика истинных спектров поглощения, компенсированных на рассеяния, определенные по формуле (4). Получим математическую модель, позволяющую описать такие данные и определить концентрации основных фотосинтетических пигментов.
Модель спектра поглощения
При анализе нормированного в точке максимума (678 нм) спектра поглощения культуры P. purpureum, для достоверного выявления пиков была вычислена производная второго порядка. Контур второй производной похож на спектр поглощения, но его полосы имеют более тонкую организацию [15]. Производная второго порядка (отрицательные области) позволяет точно выявить скрытые дополнительные максимумы в общем спектре поглощения (рис. 2, A). Можно предположить, что выявленные явные максимумы связаны с определенными полосами поглощения основных пигментов красной микроводоросли P. purpureum. Полученные данные послужили начальными точками максимумов пиков при аппроксимации спектра кривыми Гаусса. Каждый отдельный пик был описан формулой (6):
(6)
где D( ) оптическая плотность, отн. ед; D амплитуда пика, отн. ед; длина волны, нм; положение максимума пика, нм; полуширина пика, нм.
Рис. 2. Спектр поглощения культуры микроводоросли P. purpureum, нормированный по значению D и вторая производная спектра (умноженная на (-100)) (A). Математическая модель спектра хл а (B), кр (C), В-ФЭ (D), R-фикоцианина (E) и аллофикоцианина (F)
Fig. 2. Absorption spectrum of microalgae P. purpureum normalized by D value and the second derivative of the spectrum (multiplied by (-100))(multiplied by (-100)) (A). Mathematical model of the spectrum of chl a (B), car (C), B-PE (D), R-ficocyanin (E), and allophicocyanin (F)
Для описания спектра поглощения нативной формы хлорофилла а в работах [7; 16] приводится 4 кривых Гаусса с максимумами 415, 435, 623, 675 нм. Пик хл а в области 680 нм необходимо разделить на 2 кривых (673 и 683 нм). Такое разделение основано на дифференциальном спектре появление двух максимумов в области 680 нм. Аппроксимация спектра поглощения культуры P. purpureum подтвердила, что в области 680 нм пик хл а необходимо разделять на две составляющие, при этом коэффициент детерминации приближается к 1 (рис. 2, B). Нативные каротиноиды у [16] описываются двумя кривыми Гаусса с максимумами 489, 532 нм. В нашем случае (рис. 2, C) в процессе приближения функции к нативному спектру поглощения получены 2 кривые Гаусса с максимумами 465 и 490 нм. У фикобилиновых пигментов форма и положение максимумов практически не изменяются при переходе от экстрактов к нативным формам [10; 17]. На основе полученных данных о спектрах поглощения пигментов на рис. 2, (B F) приведены модели отдельных пигментов.
Таким образом, модель спектра поглощения культуры P. purpureum представляет собой сумму спектров хл а и В-ФЭ, R-фикоцианина, аллофикоцианина и кр:
(7)
где D( ) общая оптическая плотность; D ( ) оптическая плотность В-фикоэритрина; D ( ) оптическая плотность R-фикоцианина; D ( ) оптическая плотность аллофикоцианина; D ( ) оптическая плотность хлорофилла а; D ( ) оптическая плотность суммарных каротиноидов; C концентрация В-фикоэритрина (г·л ); C концентрация R-фикоцианина (г·л ); C концентрация аллофикоцианина (г·л ); C концентрация хлорофилла а (г·л ); C концентрация суммарных каротиноидов (г·л ); экстинкция В-фикоэритрина 10 л·г ·см ; экстинкция R-фикоцианина 7 л·г ·см ; экстинкция аллофикоцианина 5,8 л·г ·см ; экстинкция хлорофилла а 88,15 л·г ·см ; экстинкция суммарных каротиноидов 2500 л·г ·см .
Рис. 3. Пример разделения спектра поглощения культуры микроводоросли P. purpureum на отдельные кривые Гаусса методом аппроксимации. Черными точками показан спектр поглощения культуры; зеленая линия спектр хлорофилла а, розовая В-фикоэритирина; фиолетовая R-фикоцианина; синяя аллофикоцианина; желтая суммарных каротиноидов; красная суммарная модель спектра
Fig. 3. Example of separation of the absorption spectrum of the microalgae P. purpureum into separate Gaussian curves by approximation. Black dots show the absorption spectrum of the culture; green line chlorophyll a spectrum, pink B-phycoerythyrin; purple R-phycocyanin; blue allophycocyanin; yellow total carotenoids; red total spectrum model
На рис. 3 представлен пример разделения спектра поглощения культуры микроводоросли P. purpureum на отдельные кривые Гаусса методом аппроксимации моделью (7). Предлагаемый способ декомпозиции спектра может быть использован как экспресс-метод определения концентрации пигментов без их выделения в чистом виде, что важно при работе с накопительными культурами без вмешательства в процессы их роста.
Верификация модели
Для верификации полученных результатов была проведена серия параллельных измерений концентрации хл а, кр и В-ФЭ по стандартным методикам и предлагаемой модели (7). Сравнение результатов представлено на рис. 4.
Рис. 4. Концентрация В-фикоэритрина (A), хлорофилла а (B) и суммарных каротиноидов (C), определенных по стандартной методике (13) и при помощи модели (6). Значения углов наклона прямой во всех вариантах близки к 1 (R >0,95)
Fig. 4. Concentrations of B-phycoerythrin (A), chlorophyll a (B) and total carotenoids (C) determined using the standard technique (13) and model (5). The values of the angles of slope of the straight line in all variants are close to 1 (R >0,95)
Согласно рис. 4, значения концентраций В-ФЭ, хл а и кр, полученные с помощью стандартного метода, а также концентрации, определенные методом разложения истинных спектров поглощения на кривые Гаусса, коррелируют между собой с высокой точностью (R >0,95). Для определения концентрации хл а в обоих методах использовался коэффициент экстинкции в ацетоновом экстракте 88,15 л·г ·см [18]. Известно, что нативный коэффициент экстинкции не совпадает с эталонным значением для ацетонового экстракта. В работе [19] коэффициент экстинкции хлорофилла а в нативной форме равен 76 л·г ·см , а по данным [20] нативный коэффициент экстинкции хлорофилла а составляет 67,5 л·г ·см в максимуме поглощения 678 нм. Одним из объяснений различия между нативным коэффициентом экстинкции и коэффициента экстинкции в ацетоновом растворе является эффект "упаковки" пигментов в клетках. Результаты расчетов показывают, что при экстинкции хл а 88,15 л·г ·см коэффициент корреляции (тангенс угла наклона) равен 0,9 (рис. 4 Б). Если при расчете концентрации хл а с помощью модели (7) использовать экстинкции 100 л·г ·см , то коэффициент корреляции будет равен 1. Коэффициенты экстинкции В-ФЭ и кр не изменяются при переходе от экстрактов к нативным формам, поэтому в модели их значения составили 10 л·г ·см и 2500 л·г ·см соответственно. Следует отметить, что для определения концентрации фикобилипротеинов и суммарных каротиноидов необходимо аппроксимировать спектр поглощения культуры моделью (7). Для расчета концентрации хл а достаточно использовать значение оптической плотности при длине волны 680 нм, которое необходимо разделить на экстинкцию, так как поглощение других пигментов не влияет на оптическую плотность в этой области.
Заключение
В работе на основе анализа методов определения концентрации основных фотосинтетических пигментов: В-фикоэритрина, хлорофилла а и суммарных каротиноидов разработана модель спектра поглощения культуры P. purpureum. Установлено, что метод кривых Гаусса с высокой точностью коррелирует со стандартными биохимическими методиками и может быть использован как экспресс-метод определения концентрации основных фотосинтетических пигментов P. purpureum без вмешательства в ход эксперимента. Однако предложенная модель не позволяет провести разделение каротиноидов (бета-каротина и зеоксантина). Отметим, что при работе с моделью (7) необходимо использовать только истинный спектр поглощения культуры, который, в свою очередь, должен быть записан на спектрофотометре, оснащенном интегрирующей сферой.
Работа выполнена в рамках госзадания ФИЦ "Институт биологии южных морей имени А.О. Ковалевского РАН" , "Комплексное исследование механизмов функционирования биотехнологических комплексов с целью получения активных веществ из гидробионтов" ( гос. регистрации 1023032700554-2-1.6.16)