Использование комбинированных подходов математического моделирования для описания каталитического цикла fоf1 -атфсинтазы
На правах рукописи
Машковцева Елена Валерьевна ИСПОЛЬЗОВАНИЕ КОМБИНИРОВАННЫХ ПОДХОДОВ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ ДЛЯ ОПИСАНИЯ КАТАЛИТИЧЕСКОГО ЦИКЛА FОF1 -АТФСИНТАЗЫ 03.01.02 – Биофизика
АВТОРЕФЕРАТ
диссертации на соискание ученой степени кандидата физико-математических наук
Москва – 2010
Работа выполнена в Научно-исследовательском институте цитохимии и молекулярной фармакологии
Научный консультант:
кандидат физико-математических наук, доцент Нарциссов Ярослав Рюрикович
Официальные оппоненты:
доктор физико-математических наук, профессор Романовский Юрий Михайлович доктор физико-математических наук Якушевич Людмила Владимировна
Ведущая организация:
Институт математических проблем биологии РАН
Защита состоится 21 октября 2010 года в 14 часов на заседании диссертационного совета Д501.001.96 при Московском Государственном университете им. М.В. Ломоносова по адресу: 119991, Россия, Москва, Ленинские горы, МГУ, биологический факультет, кафедра биофизики, «Новая» аудитория.
С диссертацией можно ознакомиться в библиотеке биологического факультета МГУ им. М.В. Ломоносова.
Автореферат разослан «_» _ 2010 г.
Ученый секретарь диссертационного совета, кандидат биологических наук М.Г. Страховская
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность исследования. Одним из ключевых разделов метаболизма живых клеток является преобразование энергии в ходе окислительно восстановительных реакций и мембранного транспорта. Важнейшую роль в эффективности функционирования энергопроизводящей системы клеток играет АТФсинтаза. Полученная энергия в форме АТФ необходима для протекания всех процессов в живых организмах, от деления до апоптоза. Однако, несмотря на значительные успехи структурной биологии в последние годы, точный механизм функционирования основного энергосинтезирующего фермента клетки до сих пор не детализирован. FoF1-АТФсинтаза относится к АТФазам F типа [Boyer, P.D., 2002] и является уникальным молекулярным мотором, использующим для синтеза АТФ энергию электрохимического градиента, создаваемого за счет переноса электронов при дыхании или фотосинтезе. Этот фермент присутствует в большинстве живых организмов – на внутренней стороне плазматической мембраны бактерий, наружной стороне тилакоидной мембраны хлоропластов, во внутренней мембране митохондрий [Foster et al, 2006;
Тихонов, А.Н., 1997]. Обилие разрозненных структурных и функциональных данных не позволило сформулировать единой теории накопления и преобразования энергии, а также объяснить крайне высокую эффективность молекулярного мотора FoF1-АТФсинтазы.
В последние годы растет число пациентов с генетическими дефектами митохондриальной АТФсинтазы. В настоящее время большинство митохондриальных болезней являются смертельными, что отчасти обусловлено недостаточным пониманием механизма формирования синдромов и неверным выбором терапии заболевания. Поэтому возможность воздействия на активность фермента могла бы помочь скорректировать различные патологические состояния. Одним из возможных подходов может стать моделирование патологически измененного фермента и его гипотетического механизма работы [Das, A.M., 2003].
Высокая степень интеграции фермента с мембраной не позволяет непосредственно исследовать многие параметры его функционирования.
Однако сложность его механизма также осложняет процесс моделирования синтеза/гидролиза АТФ. Накопление и преобразование механической энергии вращения в энергию макроэргической связи АТФ ограничивает применение простого кинетического подхода, при этом описание структурных особенностей фермента не объясняет обратимости катализа. Таким образом, все существующие методы моделируют только часть общего процесса (либо моделирование всего цикла крайне затруднительно по вычислительным ресурсам) и не позволяют представить достаточно подробно весь каталитический цикл АТФсинтазы в целом.
Цель исследования. Построить универсальную модель, описывающую каталитический цикл FоF1-АТФсинтазы, и на ее основе разработать программное обеспечение, имитирующее активность данного белка при различных условиях.
Для достижения поставленной цели было необходимо решить следующие задачи:
1. Получить количественную оценку времени перемещения протонов через мембрану в полуканалах FoF1-АТФсинтазы.
2. На основе уравнений динамики рассчитать основные характеристики поворота (период вращения и моменты сил) для подвижных белковых частей FoF1-АТФсинтазы.
3. Используя вероятностные алгоритмы, получить количественные оценки процессов гидролиза и синтеза АТФ на каталитических сайтах отдельного фермента.
4. На основе полученных результатов моделирования составить алгоритм описания последовательности элементарных событий в ходе полного каталитического цикла FoF1-АТФсинтазы.
5. Формализовать предложенный алгоритм в виде независимого программного обеспечения.
6. Провести сравнение результатов компьютерного моделирования с экспериментально полученными характеристиками синтеза и гидролиза АТФ в различных тканях, происходящего с участием FoF1-АТФсинтазы.
Научная новизна. Впервые разработана физическая модель, позволяющая на основе известных данных о структуре интегрального мембранного белка и предполагаемых принципах его работы описывать его полный каталитический цикл и получать временные зависимости изменения концентрации субстратов и продуктов. На основе предложенной модели впервые был разработан алгоритм, позволяющий симулировать как гидролиз АТФ, так и его синтез АТФсинтазой, а также получать количественные характеристики этих процессов в тканях.
Принципы построения модели и структура алгоритма позволяют адаптировать используемые параметры для различных биологических систем, что значительно расширяет возможности применимости предложенного в работе подхода.
С использованием результатов моделирования разработано оригинальное программное обеспечение (программный пакет PASH ver. 1.0), позволяющее проводить компьютерную симуляцию полного повторяемого каталитического цикла АТФсинтазы в различных условиях. Преимуществом предложенного программного продукта является сочетание детализированного описания элементов структуры мембранного белка и низкие требования к вычислительным мощностям используемых компьютеров.
Практическая значимость работы. Разработанное программное обеспечение на основе предложенного алгоритма описания каталитического цикла FoF1-АТФсинтазы позволяет в рамках исследовательской лаборатории на основе структурных данных и параметров экспериментальной системы получать информацию о количественных характеристиках синтеза и гидролиза АТФ в клетках и тканях. Условия проведения подобной компьютерной симуляции могут быть очень точно адаптированы к неоднородностям оцениваемой системы по измеримым параметрам (вязкость мембран, количество и положение зарядов в белке и т.д.). Подобная адаптация является крайне необходимой для оценки функционирования нативных тканей (органов), что открывает возможности в оценке энергообеспечения клеток организма без проведения дорогостоящих инвазивных процедур. Особенность предложенного алгоритма состоит еще и в том, что осуществление виртуального эксперимента в программном пакете PASH не требует существенных вычислительных мощностей и, следовательно, позволяет значительно снизить временные и финансовые затраты на оценку аналогичных показателей in vivo. В результате применения данного подхода к описанию функционирования фермента могут быть выявлены важные закономерности, последующая экспериментальная проверка которых позволит уточнить известные данные о механизмах образования энергии, что имеет большое значение в исследовании патогенеза и современных методов терапии ряда заболеваний. Упомянутая универсальность и низкая ресурсозатратность разработанного программного обеспечения дает возможность интеграции с другими компьютерными продуктами для одновременного моделирования множества происходящих в живой клетке процессов.
Апробация работы. Материалы диссертации изложены и обсуждены на конференции, посвященной 40-летию медико-биологического факультета РГМУ (г.
Москва, 2003 г.), на 13 Европейской биоэнергетической конференции (г. Пиза, Италия, 2004), на XI Российском национальном конгрессе “Человек и лекарство» (Москва, 2004), на международной конференции «Митохондрии, от молекулярного подхода к физиологии и патологии» (г. Бари, Италия, 2005 г.), на 15 Европейской биоэнергетической конференции (г. Дублин, Ирландия, 2008), на 13 конференции Международной исследовательской группы по системной биологии (г. Эльсинор, Дания, 2008 г.), на 9 международной конференции по системной биологии (г.
Гетеборг, Швеция, 2008 г.), на 14 конференции Международной исследовательской группы по системной биологии (г. Владимир, Россия, 2010 г.).
Публикации. По результатам диссертационной работы опубликовано работ, в том числе 2 статьи в рецензируемых журналах.
Структура и объем диссертации. Диссертационная работа состоит из введения, обзора литературы, и четырех глав, в которых представлены описания использованных методов, результаты исследований и их обсуждение, заключения, выводов и списка использованной литературы. Список литературы включает _ источников, из них _ на русском языке и _ на иностранных языках. Диссертация изложена на страницах машинописного текста, иллюстрирована _ рисунками и _ таблицами.
СОДЕРЖАНИЕ РАБОТЫ
Во введении обоснована актуальность темы, сформулированы цели и задачи исследования, научная новизна и практическая значимость работы.
В первой главе представлен обзор литературных данных. Подробно описаны особенности структуры и молекулярные механизмы каталитического цикла протонной АТФсинтазы F-типа. Кроме структурных данных, освещаются современные представления о каталитическом цикле фермента, проводится сравнение АТФаз, выделенных из различных организмов, анализируется влияние различных химических соединений на активность фермента, а также его функционирование в различных патологических состояниях. Приводится подробное описание и критический анализ существующих моделей FoF1 АТФсинтазы.
Описание функционирования FoF1-АТФсинтазы на основе структурных данных можно проводить с помощью модели ее каталитического цикла.
Поскольку структура множества ферментов одного типа одинакова, возможно применение единого подхода, заключающегося в моделировании работы единичного фермента, которая затем аппроксимируется на все множество работающих белков в данном объеме.
Фермент состоит из внемембранного каталитического домена F1, связанного посредством центральной ножки с мембранным доменом Fo (Рисунок 1).
Глобулярный домен F1 состоит из пяти различных главным образом гидрофильных субъединиц со стехиометрией 33111 [Fillingame et al., 1995;
Senior, 1988]. и субъединицы расположены поочередно вокруг двойной спирали из двух антипараллельных -спиралей -субъединицы. Каталитические участки на -субъединицах расположены на поверхностях между и [Abrahams et al., 1994]. Хвост -субъединицы выступает из комплекса 33 и может связываться с полярными участками с-субъединиц фактора Fo [Watts et al., 1995]. В митохондриях и субъединицы объединены с -субъединицей в центральную ножку [Orriss et al., 1996], как -субъединица бактерий и хлоропластов [Stock et al., 1999].
После отщепления фактора F1 в мембране остается часть FoF1-АТФсинтазы, ответственная за перенос протонов (фактор Fо) [Racker, 1976]. Фактор Fo состоит из большей частью гидрофобных субъединиц а, b и с со стехиометрией a1b2c9-15 [Boyer, 2002].
Рисунок 1. Строение FoF1-АТФсинтазы E.coli.
Весь каталитический цикл FoF1-АТФсинтазы можно условно разделить на несколько процессов, которые в ходе катализа протекают одновременно, но могут происходить и независимо друг от друга: непосредственно катализ синтеза/гидролиза молекул АТФ в факторе F1 на каталитических субъединицах, вращение подвижных белковых частей фермента и перенос протонов через энергосопрягающую мембрану. Все эти процессы необходимы для корректного функционирования фермента, в то же время их временные характеристики различны и могут оцениваться отдельно.
Вторая глава посвящена количественной оценке времени переноса протонов через мембранные полуканалы FoF1-АТФсинтазы.
Определение вероятного расположения протонных полуканалов в факторе Fo FoF1-АТФсинтазы по данным литературы. Были проведены анализ вероятной структуры полуканалов и приблизительная оценка возможного количества молекул воды и протонов, расположенных внутри них. Для этого по данным литературы были выявлены все заряженные и существенные аминокислотные остатки, уточнено их расположение относительно друг друга и предполагаемых полостей. Локализовались отрицательно заряженные атомы карбоксильных групп, которые могли формировать возможные «центры» областей скопления протонов при их движении через пространство внутри белка. Эти остатки выделялись при визуализации pdb-структуры, после чего проводился анализ их взаимного расположения;
по координатам, взятым из структурных данных, вычислялись расстояния между ними. Перед проведением расчетов были также указаны возможные локализации молекул воды.
Рисунок 2. Структура выходного полуканала (слева изображены 4 ТМС а субъединицы, справа – ближайшая с-субъединица). Стрелкой показано предполагаемое направление движения протонов. На спиралях показаны указанные справа аминокислоты, которые являются существенными и/или заряженными.
Возможности оценки среднего времени перемещения протонов в пространстве полуканалов с использованием квантово-механического подхода.
Существует несколько подходов в решении проблемы о перемещении протонов внутри некоторого ограниченного пространства. Особенность данной задачи состоит в том, что, по сути, движение протонов представляет собой не просто стохастическое перемещение частицы под действием взаимодействий с растворителем, но последовательную её локализацию вблизи соответствующих атомов как растворителя (воды), так и аминокислотных остатков белка. В подобных условиях с учетом малых размеров рассматриваемой системы, а также вероятностного характера перемещения частицы можно было бы воспользоваться приближением квантовой механики. Примечательно, что при всей сложности предлагаемого подхода, в сущности для правильной оценки кинетических параметров работы отдельной белковой машины необходимо лишь среднее значение времени перемещения протонов от поверхности мембраны до фиксированных аминокислотных остатков. В дальнейшем этот параметр может быть использован для компьютерного розыгрыша реального времени по методу Монте-Карло. Упомянутый подход позволяет в достаточной степени точно описывать вероятности перемещения частицы в моделируемой системе. Однако с учетом того, что результаты такого моделирования должны быть инкорпорированы в более сложную виртуальную симуляцию каталитического цикла АТФсинтазы, а также должны быть пригодны для постоянного многократного использования в ходе программных циклов, предпочтение должно быть отдано подходу, в котором исследователь сможет избегать множества дополнительных вычислений.
Применение диффузионного подхода к описанию движения протонов по полуканалам FoF1-АТФсинтазы на примере внешнего полуканала. Описание перемещения протонов и получение искомых характеристик движения, пригодных для последующего их включения в программы, происходит с помощью диффузионного приближения. При этом для оценки значения среднего времени диффузионно обусловленного перемещения частицы желательно рассматривать лишь «главные» точки возможного перемещения частицы. В частности, в качестве подобных точек можно выбрать некоторое положение вблизи входа в полуканал и область Asp61 на с-субъединице. В этом случае оценка среднего времени перемещения частицы на расстояние l в l t= % трехмерном пространстве составит: (1), 6 D где ~ - среднее время перемещения, D – коэффициент диффузии. В Таблица t представлены результаты подобных расчетов для случая перемещения протонов по различным «путям» внутри полуканала. Легко видеть, что изменение предполагаемого времени при выборе различных путей не превысит 30% от среднего значения для рассматриваемых примеров.
Предполагаемый путь перемещения D1=1.4·10-10 м/с2 D2=4.1·10-10 м/с протона 6.45·10-9 с 2.20·10-9 с ASP61 SER206 SER202 outside 8.79·10-9 с 3.00·10-9 с ASP61 SER206 outside 11.3·10-9 с 3.86·10-9 с ASP61 outside Таблица 1. Оценка среднего времени (сек) перемещения единичного протона с внешней стороны мембраны на отрицательно заряженные аминокислотные остатки в полуканале. Результаты расчетов в соответствии с выражением (1) представлены для случая двух экспериментально оцененных значений коэффициента диффузии [Albaldawi, N., 1992] и нескольких возможных путей перемещения протона.
Математическое описание протонирования/депротонирования существенных отрицательно заряженных аминокислотных остатков с-субъединиц FoF1 АТФсинтазы в процессе ее каталитического цикла. В соответствии с подходом, предложенным D.A. Cherepanov с соавт., константы скорости протонирования/депротонирования карбоксильного остатка с-субъединицы (kon и koff) в случае равновесия концентрации протонов на поверхности мембраны и внутри протонного полуканала могут быть оценены по эмпирическим 12 pK k on = 5 1012 pH ;
k off = 5 формулам: [Cherepanov, 1999]. С помощью рассчитываемых значений констант нами был предложен метод описания протонированных/депротонированных остатков Asp61. Вероятности нахождения остатка в одном из состояний будут оказывать влияние на совокупное время каталитического цикла.
Предложенный во второй главе подход позволяет количественно оценить время перемещения протонов через полуканалы АТФсинтазы и закладывает основу для моделирования ее полного каталитического цикла.
В третьей главе производится описание вращения подвижных субъединиц FoF1-АТФсинтазы в процессе ее каталитического цикла с помощью уравнений классической механики.
Постановка задачи и описание модели вращения ротора FoF1-АТФсинтазы.
Одной из основных задач данной работы является создание эффективной физической модели, описывающей вращение белкового комплекса АТФсинтазы и позволяющей осуществлять моделирование данного процесса с использованием ограниченных вычислительных мощностей среднего ПК.
Гипотеза о возможном механизме работы фермента предполагает вращение одних частей белкового комплекса относительно других. Существует предположение о том, что вращение может быть вызвано взаимодействием отрицательно заряженного Asp61 второй ТМС с-субъединиц (здесь и далее структурные данные приведены для фермента, выделенного из E.coli) с положительно заряженным остатком Arg210, принадлежащим а-субъединице, входящей в состав статора [Boyer, 2002]. С физической точки зрения вращение ротора можно описать как вращение однородного цилиндра радиусом R и массой m относительно его главной геометрической оси в однородном изотропном диэлектрике с диэлектрической проницаемостью, обусловленное взаимодействием точечного заряда q2 на поверхности цилиндра с фиксированным точечным зарядом ql (расположенным на статоре). Заряды расположены в одной плоскости, параллельной основанию цилиндра, при этом q 1 = 1.6·10 -19 Кл, а q 2 = -1.6·10 -19 Кл.
Общий вид уравнения, описывающего поведение рассматриваемой системы, с учетом силы вязкого трения, силы упругого скручивания, стохастических движений и силы кулоновского взаимодействия следующий:
d 2 (t ) d (t ) = Fq1 q 2 (t ) R sin( (t ) + b) 4 l R I dt 2 dt R 4 E (2) (t ) + Fstoch R cos( stoch ) 4 L (1 + ) Количественное описание динамики вращения ротора FoF1-АТФсинтазы с учетом различных значений параметров модели. Перед тем, как проводить вычисления на основе модели, необходимо провести анализ границ ее применимости. Так как в основе описания вращения цилиндра лежит один из принципов классической механики, границы применения модели не могут выходить за пределы ее применимости. Произведение неопределенностей канонически сопряженных физических величин на четыре порядка больше постоянной Планка. Это дает основание утверждать, что данная модель может быть адекватно описана с использованием представлений классической, а не квантовой механики.
Сначала был рассмотрен предельный случай, описывающий ситуацию, при которой вращение происходит с максимально возможной в рамках предложенной модели угловой скоростью. Такой ситуации соответствует вращение ротора заданной массы в вакууме без упругого натяжения в центральной оси за счет силы взаимодействия двух зарядов. Для данного случая было получено теоретическое решение, в результате которого минимальное время поворота ротора на элементарный угол оценено равным 0,6 нс.
Вязкость митохондриальной внутренней мембраны и модуль Юнга для субъединицы E.coli точно не известны, поэтому в данной работе мы варьировали эти параметры в рамках известных значений для других условий и белков. Был проведен анализ поведения ротора при различных значениях указанных параметров. Полученная динамика вращения ротора представляет собой медленное поступательное движение к положению равновесия. При достаточно малых значениях модуля Юнга этим положением равновесия является положение около угла, около которого затем происходят незначительные колебания ротора за счет стохастической силы. При увеличении модуля Юнга положение равновесия смещается в сторону начального положения, и уже при значениях, на 1-2 порядка меньших известных значений для большинства белков, достигаемый угол значительно меньше необходимого для осуществления каталитического цикла угла.
Изучив зависимость равновесного угла от значения модуля Юнга, можно заключить, что для данного исследуемого белка (E.coli) максимальным значением, при котором может осуществляться каталитический цикл АТФсинтазы при условии соблюдения описанной теории, является значение порядка 5·106 Н/м2. Однако очевидно, что данный результат может быть вызван несовершенством оценки структурных показателей, лежащих в основе модели, например, расстояния между зарядами, уменьшение которого неизбежно приводило бы к увеличению энергии взаимодействия и, таким образом, к увеличению максимально возможных значений показателей упругости. Кроме того, вполне вероятно, что фиксация концов ножки (-субъединицы) не столь жесткая, и позволяет некоторую свободу движения, что может компенсировать получаемую недостаточную упругость центральной ножки, приведенную выше.
Рисунок 3. Зависимости угла поворота ротора от времени при различных значениях вязкости, модуль Юнга -субъединицы Е=3·106 Н/м При допустимой упругости предельный угол поворота ротора, необходимый для протекания каталитического процесса, достигается при любых значениях вязкости. Иными словами, вязкость окружающей ротор мембраны, в отличие от модуля Юнга, не оказывает влияния на амплитуду углового движения в моделируемой системе. В то же время с увеличением вязкости мембраны время достижения Аsp61 области полуканалов, обеспечивающих протонирование/ депротонирование, существенно растет.
Во время поворота белкового ротора происходит изменение величины воздействующих на него моментов сил: параллельно с уменьшением моментов электростатической силы и силы вязкого трения происходит нарастание момента упругого натяжения. Это обусловлено изменением относительного положения зарядов, углового скручивания центральной ножки, а также возрастанием угловой скорости. По мере приближения движущегося отрицательного заряда с-субъединицы к положительно заряженному остатку происходит увеличение вклада стохастической силы, которое выражается в колебаниях значений момента силы вязкого трения.
Результаты оценки времени поворота ротора на один шаг в зависимости от значения вязкости (Е=3·106 Н/м2) представлены в Таблица 2. Легко заметить, что для высоких значений вязкости, например 1 Па·с (вязкость глицерина), время поворота на один шаг увеличивается существенно и составляет 210 мкс, что даже несколько превосходит известные экспериментальные значения [Etzold, 1997].
Таким образом мы можем установить расчетное время полного оборота ротора, в условиях значений Е=3·106 Н/м2 и =0,26 Па·с, оно составляет 1 мс, что близко к экспериментально полученным данным [Junge, W., 2001].
, Па·с 0.03 0.05 0.1 0.25 0.5 0.75 Время, мкс 6 9 20 60 120 170 Таблица 2 Время поворота ротора FoF1-АТФсинтазы на элементарный шаг в зависимости от значения вязкости мембраны.
Применение для описания вращения подвижных белковых субъединиц FoF1 АТФсинтазы уравнения динамики позволяет оценить такие важные характеристики, как время и паттерн вращения. Кроме того, приведенный выше подход позволяет получать характеристики вращения при изменении структуры и свойств фермента, а также параметров среды (в частности, вязкости мембраны).
Рисунок 4. Динамика изменения моментов сил, возникающих в процессе каталитического цикла FoF1-АТФсинтазы за время одного полного оборота фермента. А) Зависимость момента электростатической силы, возникающей между двумя взаимодействующими зарядами, от времени. Б) Зависимость момента силы вязкого трения, возникающей при вращении белковых субъединиц в мембране, от времени В) Зависимость момента силы упругого взаимодействия -субъединицы от времени. Г) Зависимость момента стохастической силы, действующей на ротор при вращении, от времени [Nartsissov, 2006].
Четвертая глава посвящена моделированию каталитического процесса FoF1-АТФсинтазы, происходящего на -субъединицах, с использованием методов Монте-Карло.
Количественные оценки вероятностей состояний -субъединиц фактора F1 в процессе каталитического цикла. Рассматривается процесс синтеза АТФ.
Предполагается, что в каждый момент времени каждая из трех -субъединиц находится в одном из трех конформационных состояний (T, O, L). Весь каталитический цикл разбит на несколько элементарных стадий, т.е. таких событий, которые сами по себе не являются достаточными для функционирования фермента, но без завершения которых катализ невозможен: обратимые процессы связывания с каталитическим сайтом АТФ и АДФ и их отщепление, а также сама реакция образования АТФ из АДФ и неорганического фосфата и наоборот. При описании подобных событий есть два основных параметра, которые характеризуют его протекание: это вероятность самого события и время, за которое данное событие произойдет.
Вероятность для -субъединицы быть в одном из состояний определяется по формулам:
P = [ ADP ] [ ATP ] 1 + ADP + ATP Kd Kd ATP P ATP = ATP K dATP [ ADP ] K (3) 1+ + d [ ATP ] K dADP [ ATP ] ADP P ADP = K ADP K dADP [ ATP ] 1+ d + [ ADP] K dATP [ ADP ] Создание алгоритма для описания последовательности конформационных состояний каталитических -субъединиц фактора F1. После того, как определены значения вероятностей основных событий, следует составить непосредственную схему последовательности того, что будет происходить в каталитических субъединицах параллельно с процессом вращения. Для этого были схематически изображены возможные события на каждой субъединице.
На каждом шаге определялось состояние -субъединицы, наличие и тип связанного нуклеотида, вероятность протекания синтеза или гидролиза АТФ, а также отсоединения нуклеотида. Блок статора был сформирован из трех блоков субъединиц, проходящих взаимопревращения в процессе каталитического цикла с учетом влияния концентрации протонов по обе стороны мембраны.
Рисунок 5. Блок-схема основной части программы по описанию каталитического цикла FoF1-АТФсинтазы.
Подобная вероятностная оценка событий позволяет получить не только количество образованных молекул АТФ и АДФ, но и время, за которое это произошло, т.е. поток нуклеотидов. Кроме того, это первая модель, которая позволяет оценить не только гидролиз АТФ, но и синтез, а также может продемонстрировать состояние каталитического равновесия.
В пятой главе производится описание компьютерного моделирования синтеза и гидролиза АТФ в митохондриях FoF1-АТФсинтазами.
Оценка параметров катализа с помощью компьютерной программы PASH по вероятностной схеме. Сформулированные в работе принципы реализованы в виде аппаратно-независимого программного обеспечения, которое позволяет получать интересующие исследователей параметры. По результатам вычислений можно получить информацию о состоянии каждой -субъединицы на каждом шаге расчета, наличии или отсутствии в них нуклеотидов, реализации процесса синтеза или гидролиза, количестве протонов, перенесенных через мембрану в обоих направлениях, а также количественную характеристику этапа вращения и суммарный угол поворота ротора фермента. Кроме того, для каждого конкретного расчета указываются полученные значения вероятностей, которые корректируются по результатам предшествующих шагов. В конечном итоге в результате работы программы можно оценить макроскопический поток синтезированного или гидролизованного АТФ, а также перенесенных протонов для одного белка (Рисунок 6).
Рисунок 6. Окно программы по расчету вероятностей и описанию состояний -субъединиц. Показаны конформационные состояния субъединиц.
Особый интерес представляет сопоставление результатов моделирования с данными экспериментальных процедур по оценке скорости катализа FoF1 АТФсинтазой синтеза/гидролиза АТФ. Из Таблица 3 очевидно, что с ростом исходной концентрации АТФ скорость гидролиза также увеличивается до некоторого предельного значения. Расчетные значения скорости гидролиза близки к экспериментальным и также демонстрируют увеличение с ростом концентрации субстрата. Однако моделирование проводилось в условиях 100% функционирующих ферментов, что в реальных условиях наблюдается редко, кроме того, может варьировать содержание FoF1-АТФсинтазы в митохондриальном белке (расчет проведен для случая 3,3% [Etzold, 1997]).
Следует отметить, что программа PASH допускает значительно более точную оценку величины скорости гидролиза АТФ АТФсинтазой, однако для ее осуществления необходимо наряду с рассматриваемыми параметрами более точно учитывать количественную характеристику вязкости мембраны.
Экспериментальная Условия Источник Расчетная скорость скорость гидролиза гидролиза АТФ, АТФ, мкмоль/мин/ мкмоль/мин/мг мг белка белка 3,1 Выделение без [Arselin, 1996] 5, фосфолипидов 2,3 [АТФ]=0.1 мМ [Syroeshkin, 1999] 4, 3,5 [АТФ]=0.2 мМ [Syroeshkin, 1999] 4, 5,2 [АТФ]=1 мМ [Syroeshkin, 1999] 5, 6,7 Mg-АТФ, Vmax [Syroeshkin, 1999] 5, Таблица 3. Экспериментальная и расчетная скорости гидролиза АТФ FoF1 АТФсинтазой при различных условиях.
При сопоставлении результатов моделирования с экспериментальной скоростью синтеза АТФ (Таблица 4) очевидно, что в рамках предложенной модели синтез АТФ происходит на несопоставимо меньшем проценте белков, нежели гидролиз. Причина этого феномена может заключаться в том, что для гидролиза необходимо наличие лишь активных сайтов -субъединиц, наличие же потенциала на сопрягающей мембране и отрицательной кооперативности вовсе необязательно (явление унисайтного катализа [Penefsky et al., 1991]). В тоже время в случае синтеза АТФ принципиально наличие нативной сопрягающей мембраны, обладающей достаточной величиной мембранного потенциала. Фактически в случае синтеза наблюдается не просто оценка активности самого фермента, а оценка функциональной активности комплекса ферментов дыхательной цепи и АТФсинтазы. Поскольку выделение препарата митохондрий может приводить к нарушению их нативной структуры, то возникает ситуация, при которой количество реально существующего в пробе белка АТФсинтазы не соответствует количеству реально работающего на синтез фермента.
Вязкость Расчетная скорость Экспериментальная скорость мембраны, синтеза АТФ, синтеза АТФ, мкмоль/мин/мг белка Па·с мкмоль/мин/мг белка 0,03 1, 0,2 [Pitard, 1996] 0,05 0, 0,15-0,2 [Guiseppe Attardi 1995] 0,1 0, 0,1 [Cortes-Hernandez, 2007] 0,25 0, 0,058 [Pepe, 2002] 0,5 0, 0,004 [Massart, 2002] 0,75 0, 0,0029 [Baracca, 2000] 1 0, Таблица 4 Экспериментальная и расчетная скорости синтеза АТФ FoF1-АТФсинтазой при различных условиях. Данные приведены для случая 1% работающих ферментов.
Таким образом, описанная в данном разделе разработанная программа PASH для моделирования каталитического цикла FoF1-АТФсинтазы является законченным программным продуктом, позволяющим исследователям в рамках лабораторий с использованием ПК получать оценку потоков синтеза/гидролиза АТФ ферментом, а также ряд других параметров функционирования описанного белка. С помощью PASH возможно проводить оценку общего количества синтезированного/гидролизованного АТФ в клетке, ткани, органе, во всем организме.
Расчет синтеза/гидролиза АТФ в реальной системе на примере отдельных клеток и тканей. С помощью литературных данных [Безбородкина, 2008] и рассчитанного с помощью PASH количества синтезируемого в единицу времени АТФ, легко оценить общую массу АТФ, синтезируемого любым органом. Так, для печени масса производимого АТФ составляет 1289 г в час (для 5% работающих АТФсинтаз). Известно, что для полноценного развития эмбриона человека необходимо наличие в ооците не менее митохондрий [J.C.St.John, 2004]. Количество энергии, необходимое для данного процесса составляет 1,5 мкг АТФ в сутки, что примерно втрое больше массы самого ооцита (0,4 мкг [Кьюмерле Х.П., 1987]).
Выводы 1. С использованием комбинированных методов математического моделирования построена универсальная модель, описывающая совокупность элементарных процессов, происходящих на белковых субъединицах АТФсинтазы.
2. На основе предложенной физической модели, базирующейся на данных рентгеноструктурного анализа, сформулирован алгоритм, представляющий полный каталитический цикл исследуемого белка.
3. Разработано независимое программное обеспечение, позволяющее проводить компьютерные эксперименты, описывающие синтез и гидролиз АТФ моделируемым белком при различных условиях.
4. Показано, что реализация предложенного подхода позволяет без дополнительных условий с помощью компьютерной симуляции получать экспериментально измеримые характеристики исследуемой системы в условиях обратимости процессов, происходящих на АТФсинтазе.
5. На основе сопоставления результатов моделирования с экспериментальными данными продемонстрирована возможность использования разработанного программного продукта для получения количественных характеристик энергетического обеспечения тканей живого организма.
Список работ, опубликованных по теме диссертации 1. Машковцева Е.В. Теоретическое описание механо-конформационных изменений подвижных белковых составляющих FoF1-АТФсинтазы. // Вестник РГМУ – 2003. №3 (34).- С. 2. Машковцева Е.В., Нарциссов Я.Р. Возможности применения теоретического моделирования энергосинтезирующих систем клеток для изучения патологических состояний различной этиологии. // Тезисы докладов ХI Российского Национального Конгресса "Человек и Лекарство".-2004.- С.813.
3. Нарциссов Я.Р., Машковцева Е.В. Математическая модель механо конформационных изменений в интегральных мембранных белковых комплексах FОF1-АТФсинтазы. // Биофизика.- 2005.- Т.50.- № 6.- С.1048 1054.
4. Nartsissov Ya.R., Mashkovtseva E.V. A novel mechanical model of ATP synthesis and hydrolysis in FOF1-ATP synthase. // Biochim. Biophys. Acta. 2004.- EBEC Short report.-V.13.- P.117.
5. Mashkovtseva E.V., Nartsissov Ya.R. Theoretical description of rotation within FOF1-ATP synthase using solid-state mechanics approach. // Abstracts of Bari International Conference on Mitochondria, from Molecular insight to Physiology and Pathology.- 2005.- Bary - P.36.
6. Nartsissov Ya. R., Mashkovtseva E.V. Application of rigid body mechanics to theoretical description of rotation within FOF1-ATP synthase. // J. Theor. Biol. 2006.- V.242.- №2.- P.300-308.
7. Mashkovtseva E.V., Nartsissov Ya.R. Quantum-mechanical approach to the description of proton diffusion in the membrane protein pore of ATP-synthase. // BBA-Bioenergetics. – 2008. V.1777, suppl. 1. – P.S18.
8. Mashkovtseva E.V., Boronovsky S.E., Nartsissov Ya.R. Pulsing charge model of the proton movement through a membrane protein channel of ATP synthase. // In book: Will bottom up meet top-down? - 2008. – P.15-20.
9. Mashkovtseva E.V., Boronovsky S.E., Nartsissov Ya.R. Quantum-mechanical description of step-by-step protons motion through membrane protein channels.// The 9th International Conference on System Biology.
Abstract
book. – 2008. – P.160.
10. Mashkovtseva E.V., Boronovsky S.E., Nartsissov Ya.R. Mechano-probability model of FoF1-ATPsynthase Catalytic Cycle: a novel approach to ATP synthesis and hydrolysis description. // In book: Modern trends in Systems biology. Virtual modeling and regulation. - 2010. – P.153-155.