На правах рукописи МОГИЛЕВСКАЯ Екатерина Александровна ТЕОРЕТИЧЕСКОЕ ИССЛЕДОВАНИЕ МЕХАНИЗМОВ ФУНКЦИОНИРОВАНИЯ И РЕГУЛЯЦИИ ЦИКЛА КРЕБСА МИТОХОНДРИИ И ESCHERICHIA COLI 03.00.02 – Биофизика АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата биологических наук Красноярск – 2007 2 Работа выполнена в Институте физико-химической биологии им. А.Н. Белозерского МГУ им. М.В. Ломоносова. Научный руководитель: кандидат физико-математических наук, Олег Владимирович Демин старший научный сотрудник Института физико-химической биологии им. А.Н. Белозерского Официальные оппоненты: доктор биологических наук, Мария Николаевна Кондрашова профессор Института теоретической и экспериментальной биофизики РАН доктор физико-математических наук, Юрий Леонидович Гуревич ведущий научный сотрудник Института биофизики СО РАН Ведущая организация: Гематологический научный центр РАМН Защита диссертации состоится «06» ноября 2007 г. в 1000 ч. мин. на заседании Диссертационного совета Д 003.007.01 в Институте биофизики СО РАН по адресу: 600036, г. Красноярск, Академгородок, д. 50, стр. 50. С диссертацией можно ознакомиться в библиотеке Института биофизики СО РАН. Автореферат разослан «15» августа 2007 г. Ученый секретарь диссертационного совета, доктор физико-математических наук Н.С. Кудряшева 3 ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ Актуальность проблемы. В последние годы наблюдается существенный прогресс в области молекулярно-биологических и генетических исследований бактерии Escherichia coli и других организмов. Клеточный метаболизм, включающий взаимопревращения тысяч молекул в виде катализируемых ферментами биохимических реакций, - очень активно изучаемая сейчас система. Для большинства биохимических путей в E. coli известны все молекулы-метаболиты и стехиометрия их взаимопревращений, т.е. имеется так называемая статическая информация о последовательности реакций и о том, какие молекулы и в каких количествах в них участвуют. В данной ситуации может возникнуть иллюзия, что статической информации достаточно для предсказания функционирования биохимических путей при решении фундаментальных и прикладных (например, биоинженерных) задач. Отчасти, это действительно так – в некоторых случаях возможно сделать предсказания на основе только статической информации [Л1]. Однако, как правило, такой подход не оправдывает себя, поскольку внутриклеточные процессы определяются не только последовательностью реакций, но и регуляторными влияниями интермедиатов на ферменты, и генетической регуляцией уровней экспрессии ферментов, что позволяет клетке адаптироваться к изменениям внешней среды. Регуляторные механизмы в клетке ответственны за поддержание гомеостаза и переходы между различными физиологическими состояниями метаболизма. Вот почему крайне важно включать в модели регуляторные механизмы метаболических путей. Для построения таких моделей в данной работе был использован подход кинетического моделирования [Л2], который может быть, в частности, применен для изучения побочных эффектов лекарств. На основе данных по стехиометрии и регуляторным механизмам могут быть реконструированы внутриклеточные процессы, на которые 4 лекарство оказывает негативное действие. Как правило, лекарства имеют множественные эффекты на внутриклеточный метаболизм, например, несколько ферментов могут быть активированы или ингибированы, а также ряд ферментов может быть вовлечен в экскрецию лекарства. Проблема выявления главных и второстепенных механизмов токсического действия лекарств не может быть решена только экспериментально, т. к. анализ различных влияний требует их избирательного выключения, а это, как правило, невозможно сделать в эксперименте. Кинетическое моделирование позволяет исследовать каждый эффект в отдельности и понять, какой из них имеет больший вклад в общий побочный эффект. Цель данной работы заключалась в выявлении и описании с помощью кинетических моделей особенностей функционирования и регуляции цикла Кребса митохондрии и Escherichia coli. Для достижения поставленной цели в рамках настоящей работы решались следующие основные задачи: 1. Разработать кинетические модели цикла Кребса митохондрии, потребляющей глутамат и малат, и Escherichia coli, растущей на ацетате аэробно, основанные на детальном описании функционирования отдельных ферментов. 2. Получить с помощью кинетических моделей цикла Кребса митохондрии и Escherichia coli ответы системы на изменение внешних условий энергетической и биосинтетической нагрузки клетки. 3. Промоделировать совместное и изолированное влияние отдельных механизмов ингибирования салицилатом цикла Кребса для выявления среди них критических и предложить способы восстановления стационарной скорости в цикле Кребса, сниженной салицилатом. 4. Показать, возможно ли, объединяя в модели in vitro экспериментальные данные из различных источников, описать измеренное in vivo распределение потоков в цикле Кребса E. coli, потребляющей ацетат как источник углерода. 5 Научная новизна. Впервые разработана кинетическая модель сегмента цикла Кребса митохондрии, который функционирует в состоянии повышенной активности митохондрий, на основе информации об отдельных ферментах. Значения кинетических параметров, входящих в уравнения скорости, оценены по in vitro литературным экспериментальным данным. Концентрации ферментов определены из экспериментальных данных по дыханию суспензии митохондрий на глутамате и малате. На модели изучены эффекты салицилата на энергетический метаболизм митохондрии. Было показано, что ингибирование сукцинатдегидрогеназы и -кетоглутаратдегидрогеназы вносит существенный вклад в общее ингибирующее действие салицилата, тогда как разобщение окислительного фосфорилирования и потребление коэнзима А в реакциях трансформации салицилата незначительно влияют на скорость окисления субстратов в цикле Кребса. Модель позволяет предсказать, что заингибированный салицилатом поток в цикле Кребса может быть увеличен путем перераспределения потоков в цикле увеличением концентраций внемитохондриальных глутамата и малата и снижением концентраций внутримитохондриального внемитохондриального-кетоглутарата глицина. Разработана также и детальная кинетическая модель цикла Кребса Escherichia coli, растущей на ацетате аэробно. На основе in vitro данных дано подробное описание функционирования и регуляции ферментов полного цикла Кребса и глиоксилатного шунта, учтена регуляция изоцитратдегидрогеназы путем фосфорилирования. С помощью модели показано, как будет изменяться распределение потоков между циклом Кребса и глиоксилатным шунтом при изменении энергетической и биосинтетической активностей клетки. Практическое значение. С помощью построенной модели сегмента цикла Кребса митохондрии, функционирующей при повышенном энергопотреблении, изучено ингибирующее влияние салицилата на стационарный поток по циклу Кребса и предложены возможные способы 6 предотвращения его уменьшения. На модели цикла Кребса Escherichia coli показана возможность объединения in vitro экспериментальных данных по отдельным ферментам для описания поведения системы in vivo. Предложенный подход кинетического моделирования позволяет решать практические фармакологические (показать, как можно уменьшить токсические эффекты уже существующих лекарств, а также предсказать возможные побочные эффекты новых лекарственных веществ, которые находятся в стадии разработки) и биоинженерные задачи. Апробация работы. Результаты диссертации докладывались на 2 съезде токсикологов России (Москва, 2003), на 3 европейской конференции по вычислительной биологии (Глазго, 2004), на 11 и 12 международных конференциях по биотермокинетике (Оксфорд, 2004, Тракай, 2006), на 12 международной конференции «Математика. Компьютер. Образование» (Пущино, 2005), на конференции «Российская биоэнергетика: от молекул к клетке» (Москва, 2005). Структура и объем диссертации. Работа состоит из введения, пяти глав исследований, заключения и списка литературы. Работа представляет собой рукопись на 166 страницах, включая 57 рисунков и 5 таблиц. 7 СОДЕРЖАНИЕ РАБОТЫ Глава 1. Обзор литературы В первой главе представлен обзор литературы, посвященной структуре цикла Кребса, гепатотоксическому действию салицилатов и математическим моделям цикла Кребса митохондрии и E. coli. Цикл трикарбоновых кислот – это заключительный общий путь окисления молекул-источников энергии: углеводов, аминокислот и жирных кислот, стехиометрия обеспечивающий цикла генерирование трикарбоновых кислот ATP. была Кольцевая установлена Х.А. Кребсом и У.А. Джонсоном в 1937 году [Л3]. К настоящему времени накоплен большой массив данных, описывающих отдельные стадии цикла. Для всех ферментов имеются экспериментально полученные зависимости скорости реакций от концентраций субстратов, продуктов, изучена регуляция отдельных ферментов интермедиатами цикла. Для многих ферментов предложены возможные механизмы их функционирования, описывающие экспериментальные данные. Функционирование цикла Кребса в целом также изучалось экспериментально и теоретически. Было установлено, что субстратами цикла Кребса могут быть различные субстраты - сукцинат, глутамат, малат, -кетоглутарат, а также пируват, жирные кислоты и аминокислоты источники молекулы ацетил-коэнзима A. Известно, что в состоянии покоя и активности ткани окисление субстратов в митохондриях может идти не по идентичным маршрутам. Один из субстратов ЦТК – янтарная кислота (сукцинат) – по интенсивности окисления и аккумуляции энергии резко превышает остальные. Кинетическое преимущество окисления сукцината особенно важно при энергетических затратах. М.Н. Кондрашовой было показано [Л4], что при активной деятельности происходит переключение на преимущественное использование сукцината. Активация обмена в организме (мышечной нагрузкой, стрессом, введением катехоламинов, гипоксическим или холодовым воздействием, при пробуждении от спячки) 8 приводит к избирательному усилению окисления и образования сукцината, в то время как окисление НАД-зависимых субстратов не меняется. Показано, что эффективность использования сукцината, образующегося в митохондриях из -кетоглутарата, приблизительно в 100 раз выше, чем добавляемого в среду инкубации [Л5]. В 1989 г. М.Н. Кондрашовой [Л4] был рассмотрен путь притока высокоэффективного сукцината через сукцинил-КоА из -кетоглутарата, образующегося не через участок лимонных кислот ЦТК, а путем переаминирования глутамата и оксалоацетата, осуществляемого глутамат-оксалоацетат трансаминазой. Этот путь является важным путем интенсивного образования сукцината при активации функций митохондрии, в частности при напряжении, связанном с патологическими процессами. Модель именно такого варианта функционирования цикла Кребса была разработана в данной работе для изучения побочных эффектов салицилата. Существование такого укороченного цикла Кребса было показано в митохондриях гепатомы Морриса 3924A [Л6]. В ряде экспериментальных работ [Л7] было также установлено, что митохондриальный цикл Кребса может функционировать в отличном от традиционного варианте. Функционирование цикла Кребса в бактериях и, в частности, в Escherichia coli, также изучено достаточно полно. При аэробном росте E. coli на ацетате наблюдается экспрессия всех ферментов цикла Кребса, а также ферментов глиоксилатного шунта [Л8]. В этом случае два атома углерода молекулы ацетата не потребляются полностью в декарбоксилирующих реакциях цикла Кребса, а частично направляются в глиоксилатный шунт для использования в биосинтетических процессах [Л9]. После формирования изоцитрата поток углерода разделяется между изоцитратдегидогеназой и изоцитратлиазой. Распределение изоцитрата между этими двумя реакциями регулируется обратимым фосфорилированием изоцитратдегидрогеназы ее киназой/фосфатазой. Этот белок, в свою очередь, подвержен влиянию ряда центральных 9 метаболитов, чьи концентрации служат сигналами энергетической и биосинтетической потребностей клетки. В настоящей работе кинетическая модель сегмента цикла Кребса митохондрии, функционирующего при потреблении глутамата и малата, применялась для изучения гепатотоксических эффектов салицилата, который относится к группе нестероидных противовоспалительных средств (НПВС). В то время как механизмы возникновения язвообразования при его применении достаточно хорошо изучены, изменения в клетке, развивающиеся при токсических воздействиях, в частности, на печень, остаются не до конца выясненными. Из литературы известно, что большие дозы аспирина могут приводить к некрозу гепатоцитов, развитию истинной печеночной недостаточности, к последующей эндогенной интоксикации, печеночной коме и смерти [Л10]. Поражение печени как побочный эффект может проявляться и при применении аспирина и других салицилатов в терапевтических дозах [Л11]. Одной из причин этих тяжелых последствий являются нарушения в функционировании и регуляции энергетического метаболизма гепатоцитов, в частности, цикла Кребса. В процессе метаболизма в печени аспирин и его производные могут ингибировать -окисление жирных кислот [Л12], уменьшать пул кофермента А [Л13], ингибировать кетоглутаратдегидрогеназу и сукцинатдегидрогеназу [Л14]. Известно также, что салицилаты повышают проницаемость внутренней митохондриальной мембраны для протонов, уменьшая тем самым трансмембранный потенциал [Л15]. Все эти негативные воздействия могут приводить к нарушениям энергетического метаболизма. Ранее был построен ряд математических моделей разной степени детализации для изучения функционирования цикла Кребса. В работах Е.Е. Селькова и соавт. [Л16], а также в работе R. Ramakrishna et al. [Л17], исследовалась лишь стехиометрия системы. В работах группы В.В. Дынника [Л18] изучались различные регуляторные механизмы в 10 цикле Кребса. Также были разработаны детальные модели, включающие описание кинетики ферментов [Л19]. Спектр моделей цикла Кребса E. coli не так широк [Л20]. С помощью этих моделей были выявлены особенности работы и регуляции цикла Кребса, однако к их недостаткам можно отнести произвольную запись уравнений скорости, которые не отражают механизмов каталитических циклов ферментов, а также отсутствие учета действия некоторых эффекторов на ферменты. Также в перечисленных работах значения параметров, входящих в уравнения скорости, которые не могут быть измерены в эксперименте, выбираются произвольным образом. Представленный в данной работе подход основан на выводе уравнений скорости ферментов согласно механизмам их работы, учете известных регуляторных связей в цикле Кребса, а также определении неизвестных параметров из литературных экспериментальных данных, что позволяет делать более достоверные выводы и предсказания. Глава 2. Материалы и методы 2.1 Кинетическая модель сегмента цикла Кребса митохондрии, потребляющей глутамат и малат. Кинетическая модель представляет собой систему обыкновенных дифференциальных уравнений, которая в каждый момент времени определяет состояние рассматриваемой системы химических реакций, т.е. задает концентрации метаболитов этой совокупности реакций как функции времени [Л22]: dX VпродукцииХ VпотребленияХ dt (2.1) Здесь X – концентрация метаболита, V продукцииХ и V потребл енияХ - суммарные скорости его продукции и потребления. Построенная математическая модель сегмента цикла Кребса митохондрии описывает потребление глутамата и малата как субстратов (см. схему цикла на рис. 1). Переменными модели являются: Gluin, Aspin, 11 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, Обозначения: AGC-аспартат-глутаматный переносчик; AspAT-аспартатаминотрансфераза; KGDH- кетоглутаратдегидрогеназа; STK-сукцинаттиокиназа; SDH-сукцинатдегидрогеназа; FUM-фумараза; MDHмалатдегидрогеназа; KMC-кетоглутарат-малатный переносчик; ISDH-реакция связывания оксалоацетата с сукцинатдегидрогеназой; SLсалицил-CoA лигаза; SGT – салицил-CoA-глицин ацилтрансфераза; Gluin – внутримитохондриальный глутамат; Aspin – внутримитохондриальный аспартат; OAA оксалоацетат, KGin внутримитохондриальныйкетоглутарат, SucCoA сукцинилКоА, CoA – коэнзим А, Suc - сукцинат, Fum фумарат, Malin – внутримитохондриальный малат. Пунктиром обозначены ингибирующие воздействия салицилата (Sal). Рис 1. Схема сегмента цикла Кребса, потребляющего глутамат и малат при активном функционировании митохондрий, с учетом влияния салицилата. NADH не изменяются со временем, т.е. являются параметрами модели. Модель сегмента цикла Кребса митохондрии описывается системой обыкновенных дифференциальных уравнений (2.1) в соответствии со схемой (Рис. 1). внутримитохондриального Например, изменение -кетоглутарата во времени определяется следующим образом: dKGin V AspAT VKGDH VKMC dt концентрации (2.2) В системе существует четыре уравнения детального баланса: 12 Glu in Aspin N tot (сохранение азота аминогрупп); CoA SucCoA CoAtot (сохранение коэнзима А); OAA KGin SucCoA Suc Fum Mal in SDH _ OAA Ctot (сохранение четырехуглеродного скелета); SDH SDH _ OAA SDH tot (сохранение сукцинатдегидрогеназы). 2.2 Кинетическая модель цикла Кребса E.coli. Схема функционирования цикла Кребса E. coli при аэробном росте на ацетате представлена на Рис. 2. Поскольку при росте E. coli в таких Рис. 2. Схема функционирования цикла Кребса E.coli при аэробном росте на ацетате. Обозначения: ack - Ацетаткиназа; pta - Фосфотрансацетилаза; gltA - Цитратсинтаза; acn Аконитаза; icd,IDH - Изоцитратдегидрогеназа; aceK - Киназа/Фосфатаза IDH; sucAB,lpd 2-кетоглутаратдегидрогеназа; sucCD – СукцинилКоА-лигаза; sdhABCD Сукцинатдегидрогеназа; fumA - Фумараза; mdh - Малатдегидрогеназа; aceA Изоцитратлиаза; aceB - Малатсинтаза; PEPCL - Фосфоенолпируваткарбоксилаза; ME – Малик-фермент; GDH - Глутаматдегидрогеназа; cI – комплекс I; Bs – NADPH потребление на биосинтезы; As – ATP-синтаза; Al – ATP нагрузка. Пунктирные стрелки с острыми концами обозначают активирующие влияния, а пунктирные стрелки с тупыми концами – ингибирующие влияния метаболитов на ферменты. 13 условиях экспрессируются ферменты глиоксилатного шунта – изоцитратлиаза (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). Например, уравнение, определяющее концентрацию оксалоацетата, записывается следующим образом: 14 dOAA VMDH VCS VPEPCK dt (2.3) В системе содержатся пять уравнений детального баланса: ATP ADP Atot (сохранение адениновых нуклеотидов); NAD NADH NADtot (сохранение пиридиновых нуклеотидов); AcCoA CoA SucCoA CoAtot NADP NADPH NADPtot (сохранение кофермента А); (сохранение фосфорилированных пиридиновых нуклеотидов); IDH IDHP IDH tot (сохранение изоцитратдегидрогеназы). 2.3 Методы описания ферментов цикла Кребса в моделях. Для описания механизмов функционирования отдельных ферментов использовались литературные экспериментальные данные по изучению очищенных ферментов in vitro. Вывод уравнений скорости реакций ферментов основывался на принципах квазистационарного состояния и/или быстрого равновесия [Л24] и состоял из нескольких этапов: 1) Построение каталитического цикла фермента; 2) Вывод уравнения стационарной скорости реакции, катализируемой ферментом, с использованием параметров каталитического цикла (констант скоростей и констант диссоциации/равновесия отдельных элементарных реакций каталитического цикла); 3) Вывод соотношений, ферментативной ингибирования, реакции связывающих кинетические параметры (константы Михаэлиса, константы каталитические константы и др.) с параметрами использованием классических каталитического цикла; 4) Вывод уравнения скорости с кинетических параметров ферментативной реакции. После того как уравнения скорости были выведены, определялись значения фигурирующих в них кинетических параметров. Этот процесс также состоял из нескольких этапов: 15 В литературе были найдены все доступные экспериментальные 1) данные по исследованию кинетических свойств данного фермента in vitro. Эти данные представляют собой либо зависимости начальной скорости работы фермента от концентраций субстратов, продуктов, эффекторов, либо зависимости изменения концентраций субстратов и/или продуктов рассматриваемой ферментативной реакции от времени, а также данные по изотопному обмену, происходящему в присутствии фермента; Для 2) количественного описания экспериментов строилась кинетическая минимодель с использованием выведенного уравнения скорости работы фермента. Например, для описания экспериментально полученных зависимостей начальных скоростей от концентраций субстратов, продуктов и эффекторов минимоделью являлось выведенное ранее уравнение скорости, а для описания экспериментальных зависимостей субстратов (продуктов) ферментативной реакции от времени в качестве минимодели использовалась система обыкновенных дифференциальных уравнений, в правую часть которой входило уравнение скорости рассматриваемого фермента; 3) Значения параметров уравнения скорости определялись из условия наилучшего совпадения экспериментальных данных с соответствующими им результатами численного решения минимоделей. Для этого использовалась программа DBSolve7.0 [Л25], в которой реализована идентификация параметров по алгоритму Hook–Jeeves [Л26]. 2.4 Методы исследования поведения моделей цикла Кребса. Модели цикла Кребса исследовались также с помощью программы DBSolve 7.0 [Л25], которая, используя методы численного интегрирования, по заданной системе дифференциальных уравнений, описывающей модель, и заданным начальным условиям позволяет получать зависимости переменных от времени. Изучалась также зависимость стационарной скорости цикла от различных параметров. 16 Глава 3. Описание кинетики ферментов цикла Кребса В этой главе дается подробное описание кинетики ферментов цикла Кребса митохондрии и E. coli, вывод уравнений скорости работы ферментов и определения параметров, входящих в уравнение, согласно алгоритму, приведенному в Главе 2. Учтено влияние ингибиторов и активаторов на ферменты. Для тех ферментов цикла Кребса E. coli, по которым в литературе имелись данные о зависимости их активности от pH, мы вводили в уравнение скорости pH-зависимость, что позволило объединять разнородные экспериментальные данные по исследованию фермента in vitro при разных pH. Например, функционирование 2кетоглутаратдегидрогеназы необратимым механизмом E. coli Ping описывалось Pong. в соответствии Субстратное с ингибирование фермента 2-кетоглутаратом [Л27] было описано с помощью введения в модель двух центров связывания 2-кетоглутарата, между которыми существует кооперативное взаимодействие. Схема каталитического цикла приведена на Рис. 3. Кроме того, мы описали Обозначения: E1, E2, E3 – состояния фермента; E1,0депротонированная форма фермента E1; E1,2- дважды протонированная форма фермента E1. Рис. 3. Схема каталитического цикла 2-кетоглутаратдегидрогеназы E. coli. зависимость активности фермента от pH, исходя из классического предположения [Л24], что фермент может протонироваться в активном центре, причем активной является единожды протонированная форма, а депротонированная и дважды протонированные формы являются неактивными. На Рис. 3 показано протонирование свободной формы фермента Е1; таким же образом происходит протонирование всех других 17 форм фермента. Уравнение скорости фермента было выведено согласно каталитическому циклу (Рис. 3) в следующем виде: KG CoA NAD KG (k1 k 4 KG ) KG CoA NAD K d1 K d K d Kd2 CoA NAD NAD ( CoA NAD (k 2 k 3 k 4 k 3 k 2 k 4 ) k 4 k 3 NAD Kd Kd Kd k 2 k3 V KGDH KGDH * K 1H KG 2 H (1 H ) ( KG KG K d1 K d 2 H K2 CoA KG CoA NAD k 2 k 4 CoA ) KG ( CoA NAD (k 2 k 3 k1 k 3 k 2 k1 ) Kd K d1 K d K d k1 k 3 NAD CoA CoA NAD k 2 k1 CoA ) CoA NAD k 2 k 3 NAD Kd Kd Kd Kd (1 CoA KG KG 2 )( k k SucCoA * k k SucCoA * (3.1) 2 1 2 4 KG K dCoA K dKG K dKG 1 1 Kd2 k 2 k 3 SucCoA * NADH ) (1 NAD CoA )k 2 k 3 NADH * CoA ) NAD Kd Kd Константы отщепления ионов водорода от их комплекса с ферментом были определены по экспериментальным данным [Л28] (см. Рис. 4). Рис. 4. Зависимость максимальной активности 2-кетоглутаратдегидрогеназы от pH, представленная экспериментальными точками [Л28], и описываемая кривой согласно уравнению скорости (3.1) при t=25oC. Ряд неизвестных параметров определялся из экспериментальных данных [Л27] (см. Рис. 5). 18 CoA=0.5 мМ; NAD=3 мМ; KGDH=5.3 нМ; pH7.0; t=4oC Рис. 5. Зависимость начальной скорости 2-кетоглутаратдегидрогеназной реакции от концентрации 2-кетоглутарата, представленная экспериментальными точками [Л27] и описываемая теоретической кривой согласно уравнению скорости (3.1). Глава 4. Исследование функционирования цикла Кребса при повышенной активности митохондрий и влияния на него салицилата с помощью модели После того как модель сегмента цикла Кребса митохондрии, в котором потребляются глутамат и малат, была построена, исследовались ответы системы на изменение внешних условий. Было показано, как зависит скорость работы такого укороченного цикла Кребса от ATPазной нагрузки и степени восстановленности пиридиновых нуклеотидов. Для описания эффекта салицилата на энергетический метаболизм использовались литературные экспериментальные данные. Были найдены данные о том, как изменяется трансмембранный потенциал и pH в митохондрии при воздействии салицилата. Константы ингибирования салицилатом двух ферментов цикла Кребса (-кетоглутаратдегидрогеназы и сукцинатдегидрогеназы) определялись из описания полной моделью экспериментальных данных по влиянию салицилата на дыхание суспензии митохондрий (см. Рис. 6). Для того чтобы оценить вклад каждого из 19 1) KGout=10 мМ, pHout =7.4, T=30oC; (белые квадраты) 2) KGout=10 мМ, pHout =7.4, T=30oC, Sal=6,7 мМ (черные квадраты). Рис. 6. Ингибирование салицилатом скорости дыхания митохондрий на кетоглутарате, описываемая моделью и экспериментальными точками [Л14]. вышеперечисленных механизмов ингибирования салицилатом потока по описываемому сегменту цикла Кребса, рассматривалось влияние как отдельных механизмов ингибирования, так и их совместное действие на стационарный поток. На рис. 7,а показано, как скорость потребления глутамата зависит от его внемитохондриальной концентрации при отсутствии салицилата (кривая 1), при учете отдельных механизмов влияния салицилатов (кривые 2 - 5) и при совместном учете всех механизмов (кривая 6). Анализируя полученные результаты (рис. 7,а), можно заключить, что присутствие реакций, потребляющих CoA, – салицил-CoA лигазы и ацил-CoA-глицин ацилтрансферазы, практически не изменяет потока в цикле Кребса (см. Рис. 7,а, кривая 2). Разобщение окислительного фосфорилирования салицилатами (Рис. 7,а, кривая 3) оказывает незначительный эффект. Тогда как включение по отдельности остальных механизмов значительно снижает поток в цикле. В самом деле, принимая во внимание ингибирование либо кетоглутаратдегидрогеназы (Рис. 7,а, кривая 4), либо сукцинатдегидрогеназы (Рис. 7,а, кривая 5), получаем снижение скорости окисления глутамата практически в 20 раз. Включение всех возможных ингибирующих влияний салицилатов снижает поток почти до нулевого уровня (Рис. 7,а, кривая 6). Таким образом, 20 а) 1) влияние салицилатов не учитывается; 2) потребление CoA в процессе биотрансформации салицилатов; 3) разобщающее действие салицилатов; 4) ингибирование салицилатом -кетоглутаратдегидрогеназы; 5) ингибирование салицилатом сукцинатдегидрогеназы; 6) учтены все влияния салицилатов (Sal = 5 мМ) б) Gluout=20 мМ; KGout=0; Malout=0.5 мМ; Aspout=0 Рис. 7. а) Значимость различных механизмов ингибирования салицилатами цикла Кребса митохондрии, функционирующего на глутамате и малате. Зависимости скорости потребления глутамата от его внемитохондриальной концентрации соответствуют моделям, в которых учтены различные механизмы влияния салицилатов; б) Зависимость величины потоков через аспартат-глутаматный переносчик (AGC), кетоглутаратдегидрогеназу (-KGDH) и -кетоглутаратмалатный переносчик (-KMC) от концентрации салицилата. сравнение различных механизмов ингибирующего эффекта салицилата позволяет заключить, что главным образом этот эффект создается путем ингибирования сукцинатдегидрогеназы и кетоглутаратдегидрогеназы. На рис. 7,б показано, как меняется стационарное распределение потоков в рассматриваемом сегменте цикла Кребса при добавлении салицилата. Видно, что в контроле при нулевой концентрации салицилата практически весь кетоглутарат поглощается кетоглутаратдегидрогеназой (см. схему на Рис. 1), тогда как при увеличении концентрации салицилата 21 поток через кетоглутаратдегидрогеназу (VKGDH) падает, а поток через кетоглутаратмалатный переносчик (VKMC) увеличивается. Это происходит из-за ингибирования салицилатом двух ферментов нижнего сегмента цикла Кребса (см. Рис. 1), что приводит к перенаправлению потока через -кетоглутаратмалатный шунт. Предсказание возможных способов предотвращения снижения скорости потребления глутамата в цикле Кребса при воздействии салицилата. В этом разделе ставится вопрос: возможно ли найти такие изменения концентраций внемитохондриальных субстратов, которые бы позволили увеличить потребление глутамата в цикле Кребса и компенсировать ингибирующее действие салицилата? Для решения этой задачи все влияния салицилатов были включены в полную модель (Рис. 1). На рис. 8 показано, что одновременное увеличение концентрации внешнего малата (с 0.495 мМ до 10 мМ), снижение концентрации внешнего кетоглутарата (с 0.54 мМ до нуля) и снижение концентрации внутримитохондриального глицина (с 1 мМ до 1e-4 М) приводит к значительному восстановлению стационарной скорости потребления глутамата аспартатглутаматным переносчиком. Этот результат может быть объяснен следующим образом. Как было показано выше, все механизмы влияния салицилата, кроме разобщающего, воздействуют на нижнюю часть описываемого сегмента цикла Кребса (реакции KGDH, STK, SDH, FUM на Рис. 1) и не влияют на верхнюю часть (реакции AspAT, MDH и KMC). Это означает, что перенаправление потока в цикле Кребса из нижней его части в реакциюшунт, осуществляемую кетоглутаратмалатным переносчиком, может повлечь повышение потока в цикле Кребса, заингибированном салицилатом. На рис. 8,б показано, что повышение суммарного потока аспартатглутаматного переносчика осуществляется именно за счет перераспределения потоков в пользу кетоглутаратмалатного шунта. 22 а) 1 – Sal=5 мМ; KGout=0.54 мМ; Malout=0.495 мМ; Aspout=0; Gly=1 мМ; 2 – Sal=5 мМ; KGout=0; Malout=10 мМ; Aspout=0; Gly=0,1 М. б) Gluout=20 мМ; KGout=0; Malout=10 мМ; Aspout=0; Gly=0,1 М Рис. 8. а) Влияние изменения концентраций внемитохондриальных метаболитов на стационарный поток в цикле Кребса при воздействии салицилата; б) Зависимость величины потоков через аспартат-глутаматный переносчик (AGC), кетоглутаратдегидрогеназу (-KGDH) и -кетоглутаратмалатный переносчик (-KMC) от концентрации салицилата при реактивации их изменением концентраций метаболитов. Глава 5. Описание с помощью модели цикла Кребса E.coli in vivo данных и предсказания модели Значения тех свободных параметров модели цикла Кребса E. coli, функционирующего на ацетате в аэробных условиях, которые не могли быть определены из имеющихся экспериментальных данных, а именно константы равновесия для фермента, киназы/фосфатазы фосфоенолпируваткарбоксилазы, изоцитратдегидрогеназы и малик- константы скорости для суммарной биосинтетической нагрузки, были оценены из in 23 vivo данных по распределению стационарных потоков в цикле Кребса [Л21] как такие значения перечисленных параметров, которые позволили воспроизвести на модели экспериментальные данные. Для этого также пришлось увеличить концентрации изоцитратлиазы, фосфоенолпируваткарбоксилазы и глутаматдегидрогеназы, которые ранее были вычислены по специфической активности клеточного экстракта E. coli. На рис. 9 представлено сравнение распределения потоков, измеренного in vivo (Рис. 9,а) и рассчитанного с помощью модели (Рис. 9,б). Видно, что модель хорошо воспроизводит экспериментально измеренное распределение потоков в цикле Кребса. Также с помощью модели было изучено поведение системы в зависимости от ATPазной и биосинтетической нагрузок и от интенсивности работы дыхательной цепи. а) 24 б) Рис. 9. Распределение потоков в цикле Кребса E. coli, растущей на ацетате аэробно, измеренное экспериментально [Л21] (a) и полученное на модели цикла Кребса E. coli (б). Обозначения как на Рис. 2. ВЫВОДЫ Построены 1. кинетические модели сегмента цикла Кребса митохондрии, потребляющей глутамат и малат, и цикла Кребса Escherichia coli, растущей на ацетате аэробно, основанные на детальном описании функционирования отдельных ферментов. 2. С помощью модели сегмента цикла Кребса митохондрии получены ответы системы на изменение внешних условий - энергетической и биосинтетической нагрузки клетки. На модели цикла Кребса E.coli получены зависимости стационарных концентраций метаболитов и стационарных потоков от концентрации субстрата ацетата, уровня биосинтетической и ATPазной нагрузки и от уровня потребления NADH комплексом I дыхательной цепи. 25 3. Показано, что при воздействии салицилата подавление функции митохондрий обусловлено сукцинатдегидрогеназы и главным образом ингибированием кетоглутаратдегидрогеназы. Модель позволяет предсказать, что одновременное увеличение концентрации малата и глутамата и снижение концентрации кетоглутарата в цитозоле, а также снижение концентрации внутримитохондриального глицина приводит к значительному восстановлению стационарной скорости потребления глутамата в результате перераспределения потоков в сегменте цикла Кребса. 4. Модель цикла Кребса E.coli, построенная на основе in vitro экспериментальных данных из различных источников, позволяет описать экспериментальное увеличения распределение концентраций потоков в цикле ферментов при условии изоцитратлиазы, фосфоенолпируваткарбоксилазы и глутаматдегидрогеназы относительно концентраций этих ферментов, рассчитанных из соответствующих специфических активностей экстракта клеток E.coli, выращенных на ацетате. СПИСОК РАБОТ, ОПУБЛИКОВАННЫХ ПО ТЕМЕ ДИССЕРТАЦИИ 1. Могилевская (Зобова) Е.А. Исследование гепатотоксичности салицилатов и поиск способов ее предотвращения с помощью кинетической модели цикла трикарбоновых кислот // Открытый всероссийский конкурс на лучшую научную студенческую работу в разделе медицинские науки. – М. – 2002. – С. 59. 2. Могилевская (Зобова) Е.А. Исследование механизмов гепатотоксичности салицилатов с помощью кинетической модели цикла трикарбоновых кислот / Могилевская (Зобова) Е.А., Демин О.В. // Тезисы докладов 2-го съезда токсикологов России. – М. – 2003. – С. 453. 26 3. Могилевская Е.А. Кинетическая кетоглутаратдегидрогеназы модель Escherichia coli функционирования / Могилевская 2- Е.А., Лебедева Г.В., Демин О.В. // Российский Биомедицинский журнал 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 27 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. ЛИТЕРАТУРА [Л1] Edwards J.S., Ibarra R.U., Palsson B.O. In silico predictions of Escherichia coli metabolic capabilities are consistent with experimental data // Nat. Biotechnol. – 2001. – V.19. – P.125-130. [Л2] Demin O.V., Lebedeva G.V., Kolupaev A.G., Zobova E.A., Plyusnina T.Yu., Lavrova A.I., Dubinsky A.Yu., Goryacheva E.A., Tobin F., Goryanin I.I. Kinetic Modelling as a Modern Technology to Explore and Modify Living Cells // Modelling in Molecular Biology. Natural Computing Series. - Springer, 2004. - P.59-103. [Л3] Krebs H.A., Johnson W.A. The role of citric acid in intermediate metabolism in animal tissues // Enzymologia. - 1937. - N.4. - P.148-156. 28 [Л4] Кондрашова М.Н. Структурно-кинетическая организация цикла трикарбоновых кислот при активном функционировании митохондрий // Биофизика. – 1989. - Т. 34, вып. 3. - С.450-457. [Л5] Chance B.C., Hollunger G. The Interaction of Energy and Electron Transfer Reactions in Mitochondria. I. GENERAL PROPERTIES AND NATURE OF THE PRODUCTS OF SUCCINATE-LINKED REDUCTION OF PYRIDINE NUCLEOTIDE // J. Biol. Chem. - 1961. - V. 236. - P.1534-1543. [Л6] Parlo R.A., Coleman P.S. Enhanced Rate of Citrate Export from Cholesterol-rich Hepatoma Mitochondria // J. Biol. Chem. – 1984. – V. 259, N. 16. – P.9997-10003. [Л7] Reitzer L.J., Wice B.M., Kennell D. Evidence that glutamine, not sugar, is the major energy source for cultured HeLa cells // J. Biol. Chem. – 1979. - V. 254. – P.2669-2676; Yudkoff M., Nelson D., Daikhin Y., Erecinska M. Tricarboxylic acid cycle in rat brain synaptosomes. Fluxes and interactions with aspartate aminotransferase and malate/aspartate shuttle // J. Biol. Chem. - 1994. – V. 269, N.44. – P. 27414-27420. [Л8] Cronan J.E., LaPorte D. Tricarboxylic acid cycle and glyoxylate bypass // E.coli and Salm.typhimurium: Cellular and Molecular Biology. - ASM Press, 1996. - P.206216. [Л9] Kornberg H.L. and Krebs H.A. Synthesis of cell constituents from C2-units by a modified tricarboxylic acid cycle // Nature. – 1957. - V. 179. - P.988–991. [Л10] Temple A.R. Acute and chronic effects of aspirin toxicity and their treatment // Archives of Internal Medicine. – 1981. – V. 141. – P. 364-369. [Л11] Prescott L.F. Effects of non-narcotic analgesics on the liver // Drugs. – 1986. – V. 32. - P. 129-147. [Л12] Fromenty B., Pessayre D. Inhibition of mitochondrial beta-oxidation as a mechanism of hepatotoxicity // Pharmacological Therapeutics. – 1995.- V. 67. – P. 101-154. [Л13] Vessey D.A., Hu J., Kelly M. Interaction of salicylate and ibuprofen with the carboxylic acid: CoA ligases from bovine liver mitochondria // Journal of Biochemical Toxicology. – 1996. – V. 11. – P. 73-78. 29 [Л14] Kaplan E.H., Kennedy J., Davis J. Effects of salicylate and other benzoates on oxidative enzymes of the tricarboxylic acid cycle in rat tissue homogenates // Archives of Biochemistry. - 1954. – V. 51. – P. 47-61. [Л15] Haas R., Parker W.D., Stumpf Jr.D., Erugen L.A. Salicylate-induced loose coupling: protonmotive force measurements // Biochemical Pharmacology. – 1985. – V. 34. - P. 900-902. [Л16] Bohnensack R., Sel’kov E.E. Stoichiometric regulation in the citric acid cycle. I. Linear interactions of intermediates // Studia biophysica. – 1977. - V.65. – P. 161-173; Bohnensack R., Sel’kov E.E. Stoichiometric regulation in the citric acid cycle. II. Nonlinear interactions // Studia biophysica. – 1977. - V.66. – P. 47-63. [Л17] Ramakrishna R., Edwards J.S., McCulloch A., Palsson B.O. Flux-balance analysis of mitochondrial energy metabolism: consequences of systemic stoichiometric constraints // Am. J. Physiol. Regul. Integr. Comp. Physiol. – 2001. - V.280, N.3. – P. R695-R704. [Л18] Дынник В.В., Темнов А.В. Математическая модель окисления пирувата в митохондриях печени. Регуляция цикла Кребса адениновыми и пиридиновыми нуклеотидами // Биохимия. – 1977. - Т.42, вып.6. - С.1030-1044; Дынник В.В., Хайнрих Р., Сельков Е.Е. Математическая модель углеводного энергетического обмена. Взаимодействие гликолиза, цикла Кребса и Н-транспортных челноков при изменении нагрузки ATPазы // Биохимия. – 1980. – Т. 45, вып. 5. – С. 771782; Дынник В.В. Механизмы регуляции мышечного энергетического обмена при окислении глюкозы и жирных кислот. Математическая модель // Биохимия. – 1982. – Т. 47, вып. 8. – С. 1278-1288; Дынник В.В., Маевский Е.И., Григоренко Е.В., Ким Ю.В. Субстратное ингибирование в цикле трикарбоновых кислот // Биофизика. – 1984. – Т. 29, вып. 6. – С. 954-958. [Л19] Kohn M.C., Achs M.J., Garfinkel D. Computer simulation of metabolism in pyruvate-perfused rat heart. II. Krebs cycle // Am.J.Physiol. – 1979. – V. 273, N. 3. – P. R159-R166; Джафаров Р.Х. Теоретическое исследование механизмов ингибирования цикла трикарбоновых кислот избытком субстратов. Диссертация на соискание ученой степени кандидата физико-математических наук. - Пущино, 1988; Cortassa S., Aon M.A., Marban E., Winslow R.L., O’Rourke B. An integrated Model of cardiac mitochondrial energy metabolism and calcium dynamics // 30 Biophysical Journal. – 2003. – V. 84. – P. 2734-2755; K.Yugi, M. Tomita. A general computational model of mitochondrial metabolism in a whole organelle scale // Bioinformatics. – 2004. - V. 20. - P. 1795-1796. [Л20] E.M.T.El-Mansi, G.C.Dawson and C.F.A.Bryce. Steady-state modelling of metabolic flux between the tricarboxylic acid cycle and the glyoxylate bypass in Escherichia coli // Comput.Applic.Biosci. – 1994. – V. 10, N. 3. – P. 295-299; Singh V. K., Ghosh I. Kinetic modeling of tricarboxylic acid cycle and glyoxylate bypass in Mycobacterium tuberculosis, and its application to assessment of drug targets // Theoretical Biology and Medical Modelling. - 2006. – V. 3. – P. 27. [Л21] Walsh K., Koshland D.E. Determination of flux through the branch point of two metabolic cycles // J. Biol. Chem. – 1984. - V. 259, N. 15. – P. 9646-9654. [Л22] Демин О.В., Горянин И.И., Холоденко Б.Н., Вестерхофф Х.В. Кинетическое моделирование энергетического метаболизма и генерации активных форм кислорода в митохондриях гепатоцита // Молекулярная биология. – 2001. – Т. 35, вып. 6. – С. 1095-1104. [Л23] Miller S.P. et al. Locations of regulatory sites for Isocitrate Dehydrogenase Kinase/Phosphatase // J. Biol. Chem. – 2000. – V. 275, N. 2. – P. 833-839. [Л24] Корниш-Боуден Э. Основы ферментативной кинетики. М., 1979. [Л25] Goryanin I., Hodgman T.C., and Selkov E. Mathematical simulation and analysis of cellular metabolism and regulation // Bioinformatics. – 1999. - V. 15. - P. 749-758. [Л26] Hooke R. and T.A. Jeeves. “Direct search” solution of numerical and statistical problems // J. of the Association for Computing Machinery. - 1961. – V. 8. – P. 212229. [Л27] Waskiewicz D.E., Hammes G.G. Elementary steps in the reaction mechanism of the alpha-ketoglutarate dehydrogenase multienzyme complex from Escherichia coli: kinetics of succinylation and desuccinylation // Biochemistry. – 1984. – V. 23, N. 14. P. 3136-3143. [Л28] Amarasingham C.R., Davis B.D. Regulation of alpha-ketoglutarate dehydrogenase formation in Escherichia coli // J. Biol. Chem. – 1965. – V. 240, N. 9. – P. 3664-3668.