Теоретическое исследование механизмов функционирования и регуляции цикла Кребса митохондрии и Escherichia coli

Размещено на http://www.allbest.ru//

Размещено на http://www.allbest.ru//

03.00.02 - Биофизика

Автореферат

диссертации на соискание ученой степени кандидата биологических наук

Теоретическое исследование механизмов функционирования и регуляции цикла Кребса митохондрии и Escherichia coli

Могилевская Екатерина Александровна

Красноярск - 2007

Работа выполнена в Институте физико-химической биологии им. А.Н. Белозерского МГУ им. М.В. Ломоносова.

Научный руководитель:

Официальные оппоненты

Ведущая организация: Гематологический научный центр РАМН

Защита диссертации состоится «06» ноября 2007 г. в 1000 ч. мин. на заседании Диссертационного совета Д 003.007.01 в Институте биофизики СО РАН по адресу: 600036, г. Красноярск, Академгородок, д. 50, стр. 50.

С диссертацией можно ознакомиться в библиотеке Института биофизики СО РАН.

Автореферат разослан «15» августа 2007 г.

Ученый секретарь диссертационного совета,

доктор физико-математических наук Н.С. Кудряшева

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность проблемы. В последние годы наблюдается существенный прогресс в области молекулярно-биологических и генетических исследований бактерии Escherichia coli и других организмов. Клеточный метаболизм, включающий взаимопревращения тысяч молекул в виде катализируемых ферментами биохимических реакций, - очень активно изучаемая сейчас система. Для большинства биохимических путей в E. coli известны все молекулы-метаболиты и стехиометрия их взаимопревращений, т.е. имеется так называемая статическая информация о последовательности реакций и о том, какие молекулы и в каких количествах в них участвуют. В данной ситуации может возникнуть иллюзия, что статической информации достаточно для предсказания функционирования биохимических путей при решении фундаментальных и прикладных (например, биоинженерных) задач. Отчасти, это действительно так - в некоторых случаях возможно сделать предсказания на основе только статической информации [Л1]. Однако, как правило, такой подход не оправдывает себя, поскольку внутриклеточные процессы определяются не только последовательностью реакций, но и регуляторными влияниями интермедиатов на ферменты, и генетической регуляцией уровней экспрессии ферментов, что позволяет клетке адаптироваться к изменениям внешней среды. Регуляторные механизмы в клетке ответственны за поддержание гомеостаза и переходы между различными физиологическими состояниями метаболизма. Вот почему крайне важно включать в модели регуляторные механизмы метаболических путей. Для построения таких моделей в данной работе был использован подход кинетического моделирования [Л2], который может быть, в частности, применен для изучения побочных эффектов лекарств. На основе данных по стехиометрии и регуляторным механизмам могут быть реконструированы внутриклеточные процессы, на которые лекарство оказывает негативное действие. Как правило, лекарства имеют множественные эффекты на внутриклеточный метаболизм, например, несколько ферментов могут быть активированы или ингибированы, а также ряд ферментов может быть вовлечен в экскрецию лекарства. Проблема выявления главных и второстепенных механизмов токсического действия лекарств не может быть решена только экспериментально, т. к. анализ различных влияний требует их избирательного выключения, а это, как правило, невозможно сделать в эксперименте. Кинетическое моделирование позволяет исследовать каждый эффект в отдельности и понять, какой из них имеет больший вклад в общий побочный эффект.

Цель данной работы заключалась в выявлении и описании с помощью кинетических моделей особенностей функционирования и регуляции цикла Кребса митохондрии и Escherichia coli.

Для достижения поставленной цели в рамках настоящей работы решались следующие основные задачи:

Разработать кинетические модели цикла Кребса митохондрии, потребляющей глутамат и малат, и Escherichia coli, растущей на ацетате аэробно, основанные на детальном описании функционирования отдельных ферментов. митохондрия цикл кребс

2. Получить с помощью кинетических моделей цикла Кребса митохондрии и Escherichia coli ответы системы на изменение внешних условий - энергетической и биосинтетической нагрузки клетки.

3. Промоделировать совместное и изолированное влияние отдельных механизмов ингибирования салицилатом цикла Кребса для выявления среди них критических и предложить способы восстановления стационарной скорости в цикле Кребса, сниженной салицилатом.

4. Показать, возможно ли, объединяя в модели in vitro экспериментальные данные из различных источников, описать измеренное in vivo распределение потоков в цикле Кребса E. coli, потребляющей ацетат как источник углерода.

Научная новизна. Впервые разработана кинетическая модель сегмента цикла Кребса митохондрии, который функционирует в состоянии повышенной активности митохондрий, на основе информации об отдельных ферментах. Значения кинетических параметров, входящих в уравнения скорости, оценены по in vitro литературным экспериментальным данным. Концентрации ферментов определены из экспериментальных данных по дыханию суспензии митохондрий на глутамате и малате. На модели изучены эффекты салицилата на энергетический метаболизм митохондрии. Было показано, что ингибирование сукцинатдегидрогеназы и ?-кетоглутаратдегидрогеназы вносит существенный вклад в общее ингибирующее действие салицилата, тогда как разобщение окислительного фосфорилирования и потребление коэнзима А в реакциях трансформации салицилата незначительно влияют на скорость окисления субстратов в цикле Кребса. Модель позволяет предсказать, что заингибированный салицилатом поток в цикле Кребса может быть увеличен путем перераспределения потоков в цикле увеличением концентраций внемитохондриальных глутамата и малата и снижением концентраций внемитохондриального??-кетоглутарата и внутримитохондриального глицина. Разработана также детальная кинетическая модель цикла Кребса Escherichia coli, растущей на ацетате аэробно. На основе in vitro данных дано подробное описание функционирования и регуляции ферментов полного цикла Кребса и глиоксилатного шунта, учтена регуляция изоцитратдегидрогеназы путем фосфорилирования. С помощью модели показано, как будет изменяться распределение потоков между циклом Кребса и глиоксилатным шунтом при изменении энергетической и биосинтетической активностей клетки.

Практическое значение. С помощью построенной модели сегмента цикла Кребса митохондрии, функционирующей при повышенном энергопотреблении, изучено ингибирующее влияние салицилата на стационарный поток по циклу Кребса и предложены возможные способы предотвращения его уменьшения. На модели цикла Кребса Escherichia coli показана возможность объединения in vitro экспериментальных данных по отдельным ферментам для описания поведения системы in vivo. Предложенный подход кинетического моделирования позволяет решать практические фармакологические (показать, как можно уменьшить токсические эффекты уже существующих лекарств, а также предсказать возможные побочные эффекты новых лекарственных веществ, которые находятся в стадии разработки) и биоинженерные задачи.

Апробация работы. Результаты диссертации докладывались на 2 съезде токсикологов России (Москва, 2003), на 3 европейской конференции по вычислительной биологии (Глазго, 2004), на 11 и 12 международных конференциях по биотермокинетике (Оксфорд, 2004, Тракай, 2006), на 12 международной конференции «Математика. Компьютер. Образование» (Пущино, 2005), на конференции «Российская биоэнергетика: от молекул к клетке» (Москва, 2005).

Структура и объем диссертации. Работа состоит из введения, пяти глав исследований, заключения и списка литературы. Работа представляет собой рукопись на 166 страницах, включая 57 рисунков и 5 таблиц.

СОДЕРЖАНИЕ РАБОТЫ

Глава 1. Обзор литературы

В первой главе представлен обзор литературы, посвященной структуре цикла Кребса, гепатотоксическому действию салицилатов и математическим моделям цикла Кребса митохондрии и E. coli.

Цикл трикарбоновых кислот - это заключительный общий путь окисления молекул-источников энергии: углеводов, аминокислот и жирных кислот, обеспечивающий генерирование ATP. Кольцевая стехиометрия цикла трикарбоновых кислот была установлена Х.А. Кребсом и У.А. Джонсоном в 1937 году [Л3]. К настоящему времени накоплен большой массив данных, описывающих отдельные стадии цикла. Для всех ферментов имеются экспериментально полученные зависимости скорости реакций от концентраций субстратов, продуктов, изучена регуляция отдельных ферментов интермедиатами цикла. Для многих ферментов предложены возможные механизмы их функционирования, описывающие экспериментальные данные.

Функционирование цикла Кребса в целом также изучалось экспериментально и теоретически. Было установлено, что субстратами цикла Кребса могут быть различные субстраты - сукцинат, глутамат, малат, ?-кетоглутарат, а также пируват, жирные кислоты и аминокислоты - источники молекулы ацетил-коэнзима A. Известно, что в состоянии покоя и активности ткани окисление субстратов в митохондриях может идти не по идентичным маршрутам. Один из субстратов ЦТК - янтарная кислота (сукцинат) - по интенсивности окисления и аккумуляции энергии резко превышает остальные. Кинетическое преимущество окисления сукцината особенно важно при энергетических затратах. М.Н. Кондрашовой было показано [Л4], что при активной деятельности происходит переключение на преимущественное использование сукцината. Активация обмена в организме (мышечной нагрузкой, стрессом, введением катехоламинов, гипоксическим или холодовым воздействием, при пробуждении от спячки) приводит к избирательному усилению окисления и образования сукцината, в то время как окисление НАД-зависимых субстратов не меняется. Показано, что эффективность использования сукцината, образующегося в митохондриях из ?-кетоглутарата, приблизительно в 100 раз выше, чем добавляемого в среду инкубации [Л5]. В 1989 г. М.Н. Кондрашовой [Л4] был рассмотрен путь притока высокоэффективного сукцината через сукцинил-КоА из ?-кетоглутарата, образующегося не через участок лимонных кислот ЦТК, а путем переаминирования глутамата и оксалоацетата, осуществляемого глутамат-оксалоацетат трансаминазой. Этот путь является важным путем интенсивного образования сукцината при активации функций митохондрии, в частности при напряжении, связанном с патологическими процессами. Модель именно такого варианта функционирования цикла Кребса была разработана в данной работе для изучения побочных эффектов салицилата. Существование такого укороченного цикла Кребса было показано в митохондриях гепатомы Морриса 3924A [Л6]. В ряде экспериментальных работ [Л7] было также установлено, что митохондриальный цикл Кребса может функционировать в отличном от традиционного варианте.

Функционирование цикла Кребса в бактериях и, в частности, в Escherichia coli, также изучено достаточно полно. При аэробном росте E. coli на ацетате наблюдается экспрессия всех ферментов цикла Кребса, а также ферментов глиоксилатного шунта [Л8]. В этом случае два атома углерода молекулы ацетата не потребляются полностью в декарбоксилирующих реакциях цикла Кребса, а частично направляются в глиоксилатный шунт для использования в биосинтетических процессах [Л9]. После формирования изоцитрата поток углерода разделяется между изоцитратдегидогеназой и изоцитратлиазой. Распределение изоцитрата между этими двумя реакциями регулируется обратимым фосфорилированием изоцитратдегидрогеназы ее киназой/фосфатазой. Этот белок, в свою очередь, подвержен влиянию ряда центральных метаболитов, чьи концентрации служат сигналами энергетической и биосинтетической потребностей клетки.

В настоящей работе кинетическая модель сегмента цикла Кребса митохондрии, функционирующего при потреблении глутамата и малата, применялась для изучения гепатотоксических эффектов салицилата, который относится к группе нестероидных противовоспалительных средств (НПВС). В то время как механизмы возникновения язвообразования при его применении достаточно хорошо изучены, изменения в клетке, развивающиеся при токсических воздействиях, в частности, на печень, остаются не до конца выясненными. Из литературы известно, что большие дозы аспирина могут приводить к некрозу гепатоцитов, развитию истинной печеночной недостаточности, к последующей эндогенной интоксикации, печеночной коме и смерти [Л10]. Поражение печени как побочный эффект может проявляться и при применении аспирина и других салицилатов в терапевтических дозах [Л11]. Одной из причин этих тяжелых последствий являются нарушения в функционировании и регуляции энергетического метаболизма гепатоцитов, в частности, цикла Кребса. В процессе метаболизма в печени аспирин и его производные могут ингибировать ?-окисление жирных кислот [Л12], уменьшать пул кофермента А [Л13], ингибировать ?-кетоглутаратдегидрогеназу и сукцинатдегидрогеназу [Л14]. Известно также, что салицилаты повышают проницаемость внутренней митохондриальной мембраны для протонов, уменьшая тем самым трансмембранный потенциал [Л15]. Все эти негативные воздействия могут приводить к нарушениям энергетического метаболизма.

Ранее был построен ряд математических моделей разной степени детализации для изучения функционирования цикла Кребса. В работах Е.Е. Селькова и соавт. [Л16], а также в работе R. Ramakrishna et al. [Л17], исследовалась лишь стехиометрия системы. В работах группы В.В. Дынника [Л18] изучались различные регуляторные механизмы в

цикле Кребса. Также были разработаны детальные модели, включающие описание кинетики ферментов [Л19]. Спектр моделей цикла Кребса E. coli не так широк [Л20]. С помощью этих моделей были выявлены особенности работы и регуляции цикла Кребса, однако к их недостаткам можно отнести произвольную запись уравнений скорости, которые не отражают механизмов каталитических циклов ферментов, а также отсутствие учета действия некоторых эффекторов на ферменты. Также в перечисленных работах значения параметров, входящих в уравнения скорости, которые не могут быть измерены в эксперименте, выбираются произвольным образом. Представленный в данной работе подход основан на выводе уравнений скорости ферментов согласно механизмам их работы, учете известных регуляторных связей в цикле Кребса, а также определении неизвестных параметров из литературных экспериментальных данных, что позволяет делать более достоверные выводы и предсказания.

Глава 2. Материалы и методы

2.1 Кинетическая модель сегмента цикла Кребса митохондрии, потребляющей глутамат и малат.

Кинетическая модель представляет собой систему обыкновенных дифференциальных уравнений, которая в каждый момент времени определяет состояние рассматриваемой системы химических реакций, т.е. задает концентрации метаболитов этой совокупности реакций как функции времени [Л22]:

(2.1)

Здесь X - концентрация метаболита, и - суммарные скорости его продукции и потребления.

Построенная математическая модель сегмента цикла Кребса митохондрии описывает потребление глутамата и малата как субстратов (см. схему цикла на рис. 1). Переменными модели являются: Gluin, Aspin, OAA, KGin, SucCoA, СoA, Suc, Fum, Malin, SDH, SDH-OAA. Концентрации Gluout, Aspout, Hout, Hin, Malout, KGout, P, ATP, ADP, Ca2+, GTP, GDP, Q, QH2, NAD,

NADH не изменяются со временем, т.е. являются параметрами модели. Модель сегмента цикла Кребса митохондрии описывается системой обыкновенных дифференциальных уравнений (2.1) в соответствии со схемой (Рис. 1). Например, изменение концентрации внутримитохондриального -кетоглутарата во времени определяется следующим образом:

(2.2)

В системе существует четыре уравнения детального баланса:

(сохранение азота аминогрупп);

(сохранение коэнзима А);

(сохранение четырехуглеродного скелета);

(сохранение сукцинатдегидрогеназы).

2.2 Кинетическая модель цикла Кребса E.coli.

Схема функционирования цикла Кребса E. coli при аэробном росте на

ацетате представлена на Рис. 2. Поскольку при росте E. coli в таких

условиях экспрессируются ферменты глиоксилатного шунта - изоцитратлиаза (aceA) и малатсинтаза (aceB) - появляется разветвление на уровне изоцитрата, который потребляется либо изоцитратдегидрогеназой (icd, IDH), либо изоцитратлиазой. Распределение изоцитрата между этими двумя реакциями регулируется киназой/фосфатазой IDH (aceK). Если большая часть IDH фосфорилирована, т.е. является неактивной (IDHP), изоцитрат в основном направляется в глиоксилатный шунт для пополнения сукцината (Suc) и малата (Mal), расходуемых на биосинтезы. С другой стороны, если IDH находится в нефосфорилированной (активной) форме, изоцитрат поступает в нижний сегмент цикла Кребса и дважды декарбоксилируется ферментами IDH и KGDH. В последнем случае цикл Кребса реализует свою энергетическую функцию, обеспечивая восстановительными эквивалентами дыхательную цепь. Функционирование aceK как киназы/фосфатазы, в свою очередь, регулируется уровнями таких центральных метаболитов как AMP, пируват (Pyr), 3-фосфоглицерат (PG) [Л23]. Также в модель входят оттоки из малата - малик-фермент (ME), оксалоацетата - фосфоенолпируваткарбоксилаза (PEPCL), 2-кетоглутарата (KG) - глутаматдегидрогеназа (GDH). Учтены также реакции окисления NADH дыхательной цепью, потребление NADPH на биосинтезы, синтез ATP ATP-синтазой и потребление ATP на клеточные нужды. Переменными модели являются AcP, AcCoA, CoA, OAA, Cit, Aco, iCit, KG, NADP, NADPH, IDH, IDHP, сукцинил-КоА (SucCoA), NAD, NADH, ADP, ATP, Suc, фумарат (Fum), Mal, глиоксилат (Glx). Концентрации протонов (H), P, Q, QH2, PEP, PG, AMP, Pyr не изменяются со временем, т.е. являются постоянными параметрами модели. Модель цикла Кребса E. coli описывается системой обыкновенных дифференциальных уравнений (2.1), в которой изменение во времени каждой переменной определяется в соответствии со схемой (Рис. 2). Например, уравнение, определяющее концентрацию оксалоацетата, записывается следующим образом:

(2.3)

В системе содержатся пять уравнений детального баланса:

(сохранение адениновых нуклеотидов);

(сохранение пиридиновых нуклеотидов);

(сохранение кофермента А);

(сохранение фосфорилированных пиридиновых нуклеотидов);

(сохранение изоцитратдегидрогеназы).

2.3 Методы описания ферментов цикла Кребса в моделях.

Для описания механизмов функционирования отдельных ферментов использовались литературные экспериментальные данные по изучению очищенных ферментов in vitro. Вывод уравнений скорости реакций ферментов основывался на принципах квазистационарного состояния и/или быстрого равновесия [Л24] и состоял из нескольких этапов:

Построение каталитического цикла фермента;

Вывод уравнения стационарной скорости реакции, катализируемой ферментом, с использованием параметров каталитического цикла (констант скоростей и констант диссоциации/равновесия отдельных элементарных реакций каталитического цикла);

Вывод соотношений, связывающих кинетические параметры ферментативной реакции (константы Михаэлиса, константы ингибирования, каталитические константы и др.) с параметрами каталитического цикла;

Вывод уравнения скорости с использованием классических кинетических параметров ферментативной реакции.

После того как уравнения скорости были выведены, определялись значения фигурирующих в них кинетических параметров. Этот процесс также состоял из нескольких этапов:

В литературе были найдены все доступные экспериментальные данные по исследованию кинетических свойств данного фермента in vitro. Эти данные представляют собой либо зависимости начальной скорости работы фермента от концентраций субстратов, продуктов, эффекторов, либо зависимости изменения концентраций субстратов и/или продуктов рассматриваемой ферментативной реакции от времени, а также данные по изотопному обмену, происходящему в присутствии фермента;

Для количественного описания экспериментов строилась кинетическая минимодель с использованием выведенного уравнения скорости работы фермента. Например, для описания экспериментально полученных зависимостей начальных скоростей от концентраций субстратов, продуктов и эффекторов минимоделью являлось выведенное ранее уравнение скорости, а для описания экспериментальных зависимостей субстратов (продуктов) ферментативной реакции от времени в качестве минимодели использовалась система обыкновенных дифференциальных уравнений, в правую часть которой входило уравнение скорости рассматриваемого фермента;

Значения параметров уравнения скорости определялись из условия наилучшего совпадения экспериментальных данных с соответствующими им результатами численного решения минимоделей. Для этого использовалась программа DBSolve7.0 [Л25], в которой реализована идентификация параметров по алгоритму Hook-Jeeves [Л26].

2.4 Методы исследования поведения моделей цикла Кребса.

Модели цикла Кребса исследовались также с помощью программы DBSolve 7.0 [Л25], которая, используя методы численного интегрирования, по заданной системе дифференциальных уравнений, описывающей модель, и заданным начальным условиям позволяет получать зависимости переменных от времени. Изучалась также зависимость стационарной скорости цикла от различных параметров.

Глава 3. Описание кинетики ферментов цикла Кребса

В этой главе дается подробное описание кинетики ферментов цикла Кребса митохондрии и E. coli, вывод уравнений скорости работы ферментов и определения параметров, входящих в уравнение, согласно алгоритму, приведенному в Главе 2. Учтено влияние ингибиторов и активаторов на ферменты. Для тех ферментов цикла Кребса E. coli, по которым в литературе имелись данные о зависимости их активности от pH, мы вводили в уравнение скорости pH-зависимость, что позволило объединять разнородные экспериментальные данные по исследованию фермента in vitro при разных pH. Например, функционирование 2-кетоглутаратдегидрогеназы E. coli описывалось в соответствии с необратимым механизмом Ping Pong. Субстратное ингибирование фермента 2-кетоглутаратом [Л27] было описано с помощью введения в модель двух центров связывания 2-кетоглутарата, между которыми существует кооперативное взаимодействие. Схема каталитического цикла приведена на Рис. 3. Кроме того, мы описали

зависимость активности фермента от pH, исходя из классического предположения [Л24], что фермент может протонироваться в активном центре, причем активной является единожды протонированная форма, а депротонированная и дважды протонированные формы являются неактивными. На Рис. 3 показано протонирование свободной формы фермента Е1; таким же образом происходит протонирование всех других форм фермента. Уравнение скорости фермента было выведено согласно каталитическому циклу (Рис. 3) в следующем виде:

(3.1)

Константы отщепления ионов водорода от их комплекса с ферментом были определены по экспериментальным данным [Л28] (см. Рис. 4).

Ряд неизвестных параметров определялся из экспериментальных данных [Л27] (см. Рис. 5).

Глава 4. Исследование функционирования цикла Кребса при повышенной активности митохондрий и влияния на него салицилата с помощью модели

После того как модель сегмента цикла Кребса митохондрии, в котором потребляются глутамат и малат, была построена, исследовались ответы системы на изменение внешних условий. Было показано, как зависит скорость работы такого укороченного цикла Кребса от ATPазной нагрузки и степени восстановленности пиридиновых нуклеотидов.

Для описания эффекта салицилата на энергетический метаболизм использовались литературные экспериментальные данные. Были найдены данные о том, как изменяется трансмембранный потенциал и pH в митохондрии при воздействии салицилата. Константы ингибирования салицилатом двух ферментов цикла Кребса (?-кетоглутаратдегидрогеназы и сукцинатдегидрогеназы) определялись из описания полной моделью экспериментальных данных по влиянию салицилата на дыхание суспензии митохондрий (см. Рис. 6). Для того чтобы оценить вклад каждого из

вышеперечисленных механизмов ингибирования салицилатом потока по описываемому сегменту цикла Кребса, рассматривалось влияние как отдельных механизмов ингибирования, так и их совместное действие на стационарный поток. На рис. 7,а показано, как скорость потребления глутамата зависит от его внемитохондриальной концентрации при отсутствии салицилата (кривая 1), при учете отдельных механизмов влияния салицилатов (кривые 2 - 5) и при совместном учете всех механизмов (кривая 6). Анализируя полученные результаты (рис. 7,а), можно заключить, что присутствие реакций, потребляющих CoA, - салицил-CoA лигазы и ацил-CoA-глицин ацилтрансферазы, практически не изменяет потока в цикле Кребса (см. Рис. 7,а, кривая 2). Разобщение окислительного фосфорилирования салицилатами (Рис. 7,а, кривая 3) оказывает незначительный эффект. Тогда как включение по отдельности остальных механизмов значительно снижает поток в цикле. В самом деле, принимая во внимание ингибирование либо кетоглутаратдегидрогеназы (Рис. 7,а, кривая 4), либо сукцинатдегидрогеназы (Рис. 7,а, кривая 5), получаем снижение скорости окисления глутамата практически в 20 раз. Включение всех возможных ингибирующих влияний салицилатов снижает поток почти до нулевого уровня (Рис. 7,а, кривая 6). Таким образом,

сравнение различных механизмов ингибирующего эффекта салицилата позволяет заключить, что главным образом этот эффект создается путем ингибирования сукцинатдегидрогеназы и кетоглутаратдегидрогеназы. На рис. 7,б показано, как меняется стационарное распределение потоков в рассматриваемом сегменте цикла Кребса при добавлении салицилата. Видно, что в контроле при нулевой концентрации салицилата практически весь ??кетоглутарат поглощается кетоглутаратдегидрогеназой (см. схему на Рис. 1), тогда как при увеличении концентрации салицилата поток через кетоглутаратдегидрогеназу (VKGDH) падает, а поток через -кетоглутаратмалатный переносчик (VKMC) увеличивается. Это происходит из-за ингибирования салицилатом двух ферментов нижнего сегмента цикла Кребса (см. Рис. 1), что приводит к перенаправлению потока через -кетоглутаратмалатный шунт.

Предсказание возможных способов предотвращения снижения скорости потребления глутамата в цикле Кребса при воздействии салицилата.

В этом разделе ставится вопрос: возможно ли найти такие изменения концентраций внемитохондриальных субстратов, которые бы позволили увеличить потребление глутамата в цикле Кребса и компенсировать ингибирующее действие салицилата? Для решения этой задачи все влияния салицилатов были включены в полную модель (Рис. 1). На рис. 8 показано, что одновременное увеличение концентрации внешнего малата (с 0.495 мМ до 10 мМ), снижение концентрации внешнего ??кетоглутарата (с 0.54 мМ до нуля) и снижение концентрации внутримитохондриального глицина (с 1 мМ до 1e-4 ?М) приводит к значительному восстановлению стационарной скорости потребления глутамата аспартатглутаматным переносчиком. Этот результат может быть объяснен следующим образом. Как было показано выше, все механизмы влияния салицилата, кроме разобщающего, воздействуют на нижнюю часть описываемого сегмента цикла Кребса (реакции KGDH, STK, SDH, FUM на Рис. 1) и не влияют на верхнюю часть (реакции AspAT, MDH и KMC). Это означает, что перенаправление потока в цикле Кребса из нижней его части в реакцию-шунт, осуществляемую кетоглутаратмалатным переносчиком, может повлечь повышение потока в цикле Кребса, заингибированном салицилатом. На рис. 8,б показано, что повышение суммарного потока аспартатглутаматного переносчика осуществляется именно за счет перераспределения потоков в пользу кетоглутаратмалатного шунта.

Глава 5. Описание с помощью модели цикла Кребса E.coli in vivo данных и предсказания модели

Значения тех свободных параметров модели цикла Кребса E. coli, функционирующего на ацетате в аэробных условиях, которые не могли быть определены из имеющихся экспериментальных данных, а именно константы равновесия для фосфоенолпируваткарбоксилазы, малик-фермента, киназы/фосфатазы изоцитратдегидрогеназы и константы скорости для суммарной биосинтетической нагрузки, были оценены из in vivo данных по распределению стационарных потоков в цикле Кребса [Л21] как такие значения перечисленных параметров, которые позволили воспроизвести на модели экспериментальные данные. Для этого также пришлось увеличить концентрации изоцитратлиазы, фосфоенолпируваткарбоксилазы и глутаматдегидрогеназы, которые ранее были вычислены по специфической активности клеточного экстракта E. coli. На рис. 9 представлено сравнение распределения потоков, измеренного in vivo (Рис. 9,а) и рассчитанного с помощью модели (Рис. 9,б). Видно, что модель хорошо воспроизводит экспериментально измеренное распределение потоков в цикле Кребса. Также с помощью модели было изучено поведение системы в зависимости от ATPазной и биосинтетической нагрузок и от интенсивности работы дыхательной цепи.

ВЫВОДЫ

1. Построены кинетические модели сегмента цикла Кребса митохондрии, потребляющей глутамат и малат, и цикла Кребса Escherichia coli, растущей на ацетате аэробно, основанные на детальном описании функционирования отдельных ферментов.

2. С помощью модели сегмента цикла Кребса митохондрии получены ответы системы на изменение внешних условий - энергетической и биосинтетической нагрузки клетки. На модели цикла Кребса E.coli получены зависимости стационарных концентраций метаболитов и стационарных потоков от концентрации субстрата ацетата, уровня биосинтетической и ATPазной нагрузки и от уровня потребления NADH комплексом I дыхательной цепи.

3. Показано, что при воздействии салицилата подавление функции митохондрий обусловлено главным образом ингибированием сукцинатдегидрогеназы и кетоглутаратдегидрогеназы. Модель позволяет предсказать, что одновременное увеличение концентрации малата и глутамата и снижение концентрации кетоглутарата в цитозоле, а также снижение концентрации внутримитохондриального глицина приводит к значительному восстановлению стационарной скорости потребления глутамата в результате перераспределения потоков в сегменте цикла Кребса.

4. Модель цикла Кребса E.coli, построенная на основе in vitro экспериментальных данных из различных источников, позволяет описать экспериментальное распределение потоков в цикле при условии увеличения концентраций ферментов изоцитратлиазы, фосфоенолпируваткарбоксилазы и глутаматдегидрогеназы относительно концентраций этих ферментов, рассчитанных из соответствующих специфических активностей экстракта клеток E.coli, выращенных на ацетате.

СПИСОК РАБОТ, ОПУБЛИКОВАННЫХ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Могилевская (Зобова) Е.А. Исследование гепатотоксичности салицилатов и поиск способов ее предотвращения с помощью кинетической модели цикла трикарбоновых кислот // Открытый всероссийский конкурс на лучшую научную студенческую работу в разделе медицинские науки. - М. - 2002. - С. 59.

2. Могилевская (Зобова) Е.А. Исследование механизмов гепатотоксичности салицилатов с помощью кинетической модели цикла трикарбоновых кислот / Могилевская (Зобова) Е.А., Демин О.В. // Тезисы докладов 2-го съезда токсикологов России. - М. - 2003. - С. 453.

3. Могилевская Е.А. Кинетическая модель функционирования 2-кетоглутаратдегидрогеназы Escherichia coli / Могилевская Е.А., Лебедева Г.В., Демин О.В. // Российский Биомедицинский журнал Medline.ru - Электрон. журнал. - 2006. - Т. 7, Ст. 44. - С. 442-449.

4. Могилевская Е.А. Кинетическая модель функционирования и регуляции изоцитратдегидрогеназы Escherichia coli / Могилевская Е.А., Лебедева Г.В., Горянин И.И., Демин О.В. // Биофизика. - 2007. - Т. 52, В. 1. - С. 47-56.

5. Могилевская Е.А. Кинетическая модель функционирования цитратсинтазы E. coli / Могилевская Е.А., Лебедева Г.В., Демин О.В. // Математика, Компьютер, Образование. Труды XII международной конференции. - Т. 3. - Пущино, 2005. - С. 934-944.

6. Могилевская (Зобова) Е.А. Кинетическая модель цикла Кребса E.coli / Могилевская (Зобова) Е.А., Лебедева Г.В., Демин О.В. // Математика, Компьютер, Образование. Тезисы XII международной конференции. - Пущино, 2005. - С. 187.

7. Могилевская (Зобова) Е.А., Демин О.В. Кинетическая модель цикла Кребса митохондрии / Могилевская (Зобова) Е.А., Демин О.В. // Российская биоэнергетика: от молекул к клетке. Тезисы конференции. - Москва, 2005. - С. 68.

8. Mogilevskaya E.A. Application of mitochondrial Krebs cycle kinetic modeling to investigate salicylate hepatotoxic effect // Systems Biology: redefining BioThermoKinetics. Trakai. - 2006. - P. 57.

9. Mogilevskaya E.A. Cellular kinetic modeling of the microbial metabolism / Goryanin I.I., Lebedeva G.V., Mogilevskaya E.A., Metelkin E.A., Demin O.V. // Methods Biochem. Anal. - 2006. - V. 49- P. 437-488.

10. Mogilevskaya (Zobova) E.A. Kinetic Model of E.coli Krebs Cycle / Mogilevskaya (Zobova) E.A., Lebedeva G.V., Demin O.V. // Developing Concepts for Systems Biology. 11th Workshop of the BioThermoKinetics Study Group. - Oxford, 2004. - P.55.

11. Mogilevskaya (Zobova) E.A. Kinetic Model of E.coli Krebs Cycle / Mogilevskaya (Zobova) E.A., Lebedeva G.V., Demin O.V. // 12th International Conference on Intelligent Systems for Molecular Biology, 3rd European Conference on Computational Biology. - Glasgow, 2004. - P.217.

12. Mogilevskaya E.A. Kinetic Model of Mitochondrial Krebs Cycle: Unraveling the Mechanism of Salicylate Hepatotoxic Effects / Mogilevskaya E.A., Demin O.V., Goryanin I. // Journal of Biological Physics. - 2006. - V. 32. - P. 245-271.

13. Mogilevskaya (Zobova) E.A. Kinetic Modelling as a Modern Technology to Explore and Modify Living Cells / Demin O.V., Lebedeva G.V., Kolupaev A.G., Mogilevskaya (Zobova) E.A., Plyusnina T.Yu., Lavrova A.I., Dubinsky A., Goryacheva E.A., Tobin F., Goryanin I.I. // Modelling in Molecular Biology. Natural Computing Series. - Springer, 2004. - P. 59-103.

14. Mogilevskaya (Zobova) E.A. Kinetic modelling of the E. coli metabolism / Demin O.V., Plyusnina T.Y., Lebedeva G.V., Mogilevskaya (Zobova) E.A., Metelkin E.A., Kolupaev A.G., Goryanin I.I., Tobin F. // Topics in Current Genetics. - Springer, 2005. - P. 31-67.

Размещено на Allbest.ru