WWW.NEW.Z-PDF.RU
БИБЛИОТЕКА  БЕСПЛАТНЫХ  МАТЕРИАЛОВ - Онлайн ресурсы
 

Pages:   || 2 |

«Коркош Вячеслав Сергеевич Молекулярное моделирование потенциал-управляемых натриевых каналов эукариот и их взаимодействия с лигандами ...»

-- [ Страница 1 ] --

РОССИЙСКАЯ АКАДЕМИЯ НАУК

Федеральное государственное бюджетное учреждение наук

и

Институт эволюционной физиологии и биохимии им. И.М. Сеченова РАН

На правах рукописи

Коркош Вячеслав Сергеевич

Молекулярное моделирование

потенциал-управляемых натриевых

каналов эукариот и их взаимодействия с

лигандами

Специальность 03.03.01 – физиология

Диссертация

на соискание ученой степени

кандидата биологических наук

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

доктор биологических наук Тихонов Д.Б .

Санкт-Петербург Оглавление Список сокращений

Общая характеристика работы

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

1.1. Ионные каналы

1.1.1. Общие сведения

1.1.2. Потенциал-управляемые ионные каналы

1.1.3. Лиганд-управляемые ионные каналы

1.1.4. Каналы, активируемые растяжением и теплом

1.2. Молекулярная структура ионных каналов как основа классификации

1.2.1. P-loop каналы

1.2.2. Потенциал-управляемые натриевые каналы

1.3. Фармакология ионных каналов

1.3.1. Фармакология потенциал-управляемых Na-каналов.

1.3.2. Конотоксины

1.3.3. µ-конотоксины и Na-каналы

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

2.1. Методы исследования ионных каналов

2.1.1. Экспериментальные методы

2.1.2. Рентгеноструктурный анализ

2.1.3. Теоретические методы

2.1.4. Молекулярная динамика

2.1.5. Метод Монте-Карло

2.2. Параметры и методы расчётов, использованные в работе

2.2.1. Расчёт энергетически оптимальных конформеров методами молекулярной механики

2.2.2. Гомологическое моделирование и системы ограничителей

2.2.3. Анализ лиганд-рецепторного взаимодействия

2.3. Программное обеспечение

2.3.1. ZMM

2.3.2. Прочие

Глава 3 . Результаты и их обсуждение .

3.1. Проблемы гомологического моделирования

3.1.1. Выбор выравнивания порообразующих спиралей Р-loop каналов

3.1.2. Анализ стабилизаторов в структурах каналов

3.1.3. Выбор выравнивания S5 и S6 спиралей в гомологических моделях

3.1.4. Заключение

3.2. Стерические и электростатические факторы блокады ионных каналов

3.2.1. Изучение возможных способов связывания лиганда

3.2.2. Распределение ионов в области селективного фильтра

3.2.3. Влияние Pl1 на энергию связи ионов

3.2.4. Заключение

3.3. Блокада натриевых каналов µ-конотоксинами

3.3.1. Выбор модели

3.3.2. Связывание GIIIA

3.3.3. Связывание PIIIA

3.3.4. Связывание KIIIA

3.3.5. Механизм неполного блока

Общее заключение

Выводы

Список литературы

Список сокращений TRP – транзиторный рецепторный потенциал nAChR – никотиновый ацетилхолиновый рецептор TTX – тетродотоксин STX – сакситоксин µCTX – мю-конотоксин ЛА – локальные анестетики PyR – пиретроидный рецептор BTX – батрахотоксин МСМ – Монте-Карло минимизация МД – молекулярная динамика NMDA – N-метил-D-аспартат AMPA - -амино-3-гидрокси-5-метил-4-изоксазолпропионовая кислота СО – среднеквадратическое отклонение Pl1 - 2-(4-бромфенил)-1-(5-(4-хлорфенил)-1Н-имидазол-2-ил)этанамин) Общая характеристика работы Актуальность проблемы .

Белки ионных каналов, обеспечивающие пассивный транспорт ионов, являются важными компонентами клеточных мембран. Активность потенциал-чувствительных натриевых и калиевых каналов лежит в основе генерации и распространения потенциала действия в нейронах. Синаптические ионные каналы, активируемые нейромедиаторами, необходимы для передачи возбуждающих и тормозных сигналов между нейронами. Многие модуляторы клеточного метаболизма также оказывают свое действие через ионные каналы. Поскольку каналообразующие белки являются трансмембранными, их внеклеточные участки и сама пора канала доступны для действия внеклеточных химических агентов, как естественных, так и искусственных .

Изучение молекулярных механизмов ионных каналов необходимо для медицины, фармакологии. Несмотря на наличие огромного арсенала лекарственных препаратов, смертность от заболеваний сердечно-сосудистой системы, например, артериальной гипертонии, постоянно растет. Одной из основных проблем современной фармакологии является недостаток информации о молекулярных основах действия препаратов, то есть о клеточных и субклеточных структурах, участвующих в патологических процессах .

Большая проблема заключается еще и в том, что «фармакологическая мишень» так же неисчерпаема, как «атом». Как только исследователь установит, что мишенью является, например, потенциал-управляемый натриевый канал, встает вопрос о механизмах действия, зависимости активности препарата от уловий работы канала и т.д .

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

Современные научные подходы, такие как рентгеноструктурный анализ, методы молекулярной биологии (например, точечные мутации), а также теоретические методы, включая молекулярное моделирование, дают возможность проводить такой анализ Натриевые каналы играют важную роль во многих процессах в организме, например, в генерации потенциала действия. Неудивительно, что фармакология ионных каналов является одним из наиболее важных предметов изучения на данный момент. Хорошо известны такие токсины как тетродотоксин, локальные анестетики, батрахотоксин, однако механизмы действия лигандов изучены не до конца. Отсутствие рентгеновских структур эукариотических натриевых каналов можно восполнить при помощи теоретических методов. Молекулярное моделирование позволяет объединить разрозненные экспериментальные данные в одной модели, которая может не только объяснить имеющиеся экспериментальные данные, но и имеет предсказательную силу. Такие модели позволяют исследователям точнее направлять свои усилия и создавать подходящие фармакологические агенты .

Одним из наиболее интересных классов лигандов натриевых каналов являются µконотоксины. Они особенно интересны тем, что некоторые нативные конотоксины, а также мутанты, созданные на основе конотоксинов, вызывают эффект неполного блока в натриевых потенциал-управляемых каналах. Данный эффект интересен с точки зрения разработки новых лигандов, поскольку не приводит к необратимым последстви ям для организма. Однако на данный момент не существует структурного объяснения для этого эффекта, что затрудняет работу по созданию новых агентов .

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

В рамках этой цели были поставлены следующие задачи:

Изучение роли межсегментных конткатов в стабилизации пространственной укладки ионных каналов Изучение роли стерических и электростатических факторов при блокаде натриевых каналов Построение гомологической модели поры эукариотического натриевого канала Моделирование связывания µ-конотоксинов с эукариотическим натриевым каналом Моделирование эффекта неполного блока эукариотических натриевых каналов µ-конотоксинами Научная новизна. В результате проведённых исследований получена модель потенциал-управляемого натриевого ионного канала Nav1.4, позволяющая объяснить экспериментальные данные по взаимодействию с µ-конотоксинами. Впервые предложено объяснение возникновению остаточных токов при неполной блокаде поры канала токсинами. Также в рамках работы впервые показаны эволюционные изменения в структурных мотивах ионных каналов, приводящие к пространственным различиям .

Научно-практическая значимость работы. Полученные в работе данные о строении, механизмах проводимости и блокады натриевых каналов объясняют широкий спектр экспериментальных данных. Модели, полученные в результате работы, могут использоваться для дизайна новых препаратов, имеющих медицинское значение. Модель ионной проводимисти позволяет целенаправленно подходить к созданию качественно новых лигандов, обеспечивающих неполный блок канала, что открывает новые перспективы в области фармаколонии .

Положения, выносимые на защиту:

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

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

Апробация работы. Результаты работы докладывались на следующих конференциях: IV съезд физиологов СНГ (Дагомыс, Россия, 2014), Biophysical Society 58 Annual Meeting (San Francisco), 11-й конгресс «Нейронаука для Медицины и Психологии »

(Судак, Россия, 2015) .

Публикации. По теме диссертации опубликовано 4 статьи в рецензируемых журналах и 3 тезисов докладов .

Структура и объём диссертации. Диссертация изложена на 135 страницах машинописного текста и состоит общей характеристики работы, обзора литературы по исследуемой теме – глава 1, описания методики – глава 2, описания результатов исследования и их обсуждения – глава 3, заключения – глава 4, выводов и списка литературы, который включает 204 источника (из них 202 иностранных).

Работа иллюстрирована 53 рисунками и 7 таблицами .

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

1.1. Ионные каналы 1.1.1. Общие сведения Согласно теории Ходжкина-Хаксли модель возбудимой мембраны предполагает регулируемый перенос ионов через мембрану (Hille, 1992, Sperelakis, 1995, Soria and Cena, 1998, Nicholls, 2003, Camerino, Tricarico et al., 2007, Siegel, 2011). Однако непосредственный переход иона через липидный бислой весьма затруднен. Очевидно, что в мембране должны существовать некоторые специальные структуры, проводящие ионы. Такие структуры были найдены и названы ионными каналами (рисунок 1). Каналы могут быть выделены из различных объектов: плазматической мембраны клеток, постсинаптической мембраны мышечных клеток и т.д. Известны также ионные каналы, образованные антибиотиками (Nicholls, 2003) .

К основным свойствам ионных каналов относятся:

селективность;

независимость работы отдельных каналов;

дискретный характер проводимости;

зависимость параметров каналов от мембранного потенциала .

Селективностью называют способность ионных каналов избирательно пропускать ионы какого-либо одного типа. Еще в первых опытах на аксоне кальмара было обнаружено что ионы Na+ и К+ по-разному влияют на мембранный потенциал. Ионы К + регулируют потенциал покоя, а ионы Na+ - ответственны за формирование потенциала действия (Hodgkin, Huxley et al., 1952). В модели Ходжкина-Хаксли это описывается путем введения независимых калиевых и натриевых ионных каналов. Предполагалось, что первые пропускают только ионы К +, а вторые - только ионы Na+. Измерения показали, что ионные каналы обладают очень высокой селективностью по отношению к катионам (катион селективные каналы) либо к анионам (анион-селективные каналы) (Hodgkin and Huxley, 1952). В то же время через катион-селективные каналы способны проходить различные катионы. Например, для Na+-кaнала калиевый ток в 20 раз меньше натриевого (Hille, 1992) .

Независимость работы отдельных каналов. Прохождение тока через отдельный ионный канал не зависит от того, идет ли ток через другие каналы. Например, К +-каналы могут быть включены или выключены, но ток через Nа +-каналы не меняется (Aldrich, 1985) .

Влияние каналов друг на друга происходит опосредованно: изменение проницаемостей каких-либо каналов (например, натриевых) меняет мембранный потенциал, а уже он влияет на токи прочих ионных каналов .

Дискретный характер проводимости ионных каналов. Ионные каналы представляют собой субъединичный комплекс белков, пронизывающий мембрану (K ew and Davies 2010). В центре его существует пора, сквозь которую могут проходить ионы .

Результаты экспериментов, выполненных на различных ионных каналах, показали, что проводимость ионного канала дискретна, и он может находиться в двух состояниях:

открытом или закрытом. Переходы между состояниями происходят в случайные моменты времени и подчиняются статистическим закономерностям (Kew and Davies, 2010) .

Зависимость параметров канала от мембранного потенциала. Ионные каналы нервных волокон чувствительны к мембранному потенциалу, например, натриевый и калиевый каналы аксона кальмара. Это проявляется в том, что после начала деполяризации мембраны соответствующие токи начинают изменяться. Ион-селективный канал имеет сенсор - некоторый элемент своей конструкции, чувствительный к действию электрического поля. При изменении мембранного потенциала меняется величина действующей на него силы, в результате эта часть ионного канала перемещается и меняет вероятность открывания или закрывания ворот - своеобразных заслонок, действующих по закону "всё или ничего" (Hille, 1978) .

Изучение ионных каналов в мембранах яляется одной из важнейших задач современной биофизики. Начиная с 90-х, исследования в рамках молекулярной генетики, в сочетании с методом локальной фиксации потенциала (пэтч-кламп) и другими методиками, значительно расширили наши знания об ионных каналах. Были выделены и клонированы гены, кодирующие натриевые, калиевые и многие другие типы каналов. Было обнаружено более чем 400 различных генов, кодирующих ионные каналы, чего никак нельзя было ожидать исходя из ранних исследований функций ионных каналов .

Чтобы понять значимость такого количества различных генов ионных каналов, проводились многочисленные исследования с помощью пэтч-клампа и других физиологических методов. Подобные исследования обнаружили большое количество потенциал-управляемых ионных каналов, которые отвечают на мембранный потенциал во многом так же, как натриевые и калиевые каналы, которые лежат в основе потенциала действия. Но существуют и другие типы каналов. Существуют каналы, активируемые химическими агентами (лигандами) с внеклеточной или с внутриклеточной стороны мембраны (Nicholls, 2003). Другие каналы могут реагировать на механические воздействия или изменения температуры .

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

1.1.2. Потенциал-управляемые ионные каналы На данный момент обнаружены потенциал-управляемые ионные каналы, которые выборочно проницаемы для каждого из основных физиологических ионов – натрия, калия, кальция и хлора (Purves, 2012). Кроме того, было обнаружено множество различных генов для каждого типа потенциал-управляемых ионных каналов. Примером может служить обнаружение 10 человеческих генов натриевых каналов. Эти гены натриевых каналов (SCN гены) кодируют белки, различающиеся по структуре, свойствам и распространению в определённых тканях .

Кроме быстро инактивируемых натриевых каналов, обнаруженных Ходжкином и Хаксли в аксоне кальмара, были обнаружены потенциал-чувствительные натриевые каналы, которые не инактивируются (Bendahhou, Cummins et al., 1999). Как и следовало ожидать, этот канал производит потенциалы действия большой продолжительности и является мишенью для местных обезболивающих, таких как бензокаин и лидокаин (Khodorova, Meissner et al., 2001) .

В некоторых нейронах, потенциал-управляемые кальциевые каналы порождают потенциалы действия, как и потенциал-чувствительные натриевые каналы (Catterall, PerezReyes et al., 2005). В других нейронах кальциевые каналы контролируют форму потенциалов действия, сгенерированных преимущественно натриевыми каналами (Van Goor, Krsmanovic et al., 1999). В целом, воздействуя на внутриклеточную концентрацию кальция, активность кальциевых каналов регулирует огромный диапазон биохимических процессов внутри клеток (Catterall, Perez-Reyes et al., 2005). Одним из важных процессов, регулируемых потенциал-чувствительными кальциевыми каналами является секреция нейромедиаторов в синапсах .

Но всё же наиболее крупным и разнообразным классом потенциал-управляемых ионных каналов являются калиевые каналы. На данный момент известно около 100 генов калиевых каналов, которые разбиты на несколько групп, значительно различающихся по активационным, пропускным и инактивационным свойствам. Некоторые инактивируются в течение минут, как в случае калиевых каналов аксона кальмара, изученного Ходжкином и Хаксли (Hodgkin and Huxley, 1952, Hodgkin, Huxley et al., 1952). Другие инактивируются в течение миллисекунд, что типично для большинства потенциал-управляемых натриевых каналов. Пожалуй, наиболее важную роль калиевые каналы играют в регуляции мембранного потенциала покоя. Как минимум две семьи калиевых каналов, которые открыты при значительном отрицательном потенциале на мембране, вносят вклад в установление мембранного потенциала покоя. Наконец, было обнаружено несколько типов потенциал-управляемых хлорных каналов. Эти каналы находятся в нейронах, контролируя возбудимость и внося вклад в мембранный потенциал покоя (Jentsch, Gnther et al., 1995) .

1.1.3. Лиганд-управляемые ионные каналы Другой класс ионных каналов реагирует на химические сигналы (лиганды), но не реагирует на изменения в мембранном потенциале (Purves, 2012). Наиболее важным из этих лиганд-управляемых ионных каналов в нервной системе является класс каналов, активируемых связыванием нейромедиаторов (Sobolevsky, Rosconi et al., 2009, Lee, Lu et al., 2014). В то время как потенциал-управляемые ионные каналы обычно пропускают только один вид ионов, каналы, активируемые внеклеточными лигандами, обычно менее избирательны, позволяя двум и более типам ионов проходить через пору канала (Green, Lambert et al., 1996) .

Другие лиганд-управляемые каналы чувствительны к химическим сигналам, возникающим в цитоплазме нейронов и могут быть избирательны к конкретным ионам, таким как калий или хлор, или же могут быть проницаемы для всех физиологических катионов. Примерами каналов, отвечающих на внутриклеточные воздействия являются кальций-активируемые калиевые каналы, катионные каналы, управляемые циклическими нуклеотидами и кислотно-чувствительные ионные каналы (Bassler, Ngo-Anh et al., 2001, Craven and Zagotta, 2006). Основная функция таких каналов - это преобразование внутриклеточных химических сигналов в электрические сигналы. Классическим примером являются каналы, управляемые циклическими нуклеотидами в светочувствительных клетках сетчатки (Becirovic, Nguyen et al., 2014) .

Хотя многие из лиганд-управляемых ионных каналов находятся во внешней клеточной мембране, некоторые находятся в мембранах внутриклеточных органелл, таких как митохондрии или эндоплазматический ретикулум. Некоторые из этих внутриклеточных каналов выборочно проницаемы для кальция и регулируют выход кальция из цистерн эндоплазматического ретикулума в цитоплазму, где этот вторичный передатчик может в свою очередь запускать набор клеточных ответов (Pinto, Kihara et al., 2015) .

1.1.4. Каналы, активируемые растяжением и теплом Существуют также классы ионных каналов, которые отвечают на изменение температуры или на деформацию мембраны. Температурно-чувствительные ионные каналы, такие как некоторые члены семейства каналов TRP (Hellmich and Gaudet, 2014, Yue, Xie et al., 2015), вносят вклад в чувствительность к боли и температурам, а также участвуют в воспалительных процессах. Эти каналы зачастую специализируются на определённых температурах, а некоторые активируются при охлаждении. Механо-чувствительные ионные каналы реагируют на механические изменения в плазматической мембране и лежат в основе рецепторов растяжения и нейромышечных рефлексов растяжения (Delmas and Coste, 2013). Одна из специфических форм таких каналов играет ключевую роль в механизме восприятия звука (Delmas and Coste, 2013) .

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

1.2. Молекулярная структура ионных каналов как основа классификации До последнего времени большая часть информации о структуре каналов была получена не напрямую, а посредством изучения свойств этих белков. Существующие классификации каналов основаны именно на таких признаках. Один из подходов к классификации основан на том, чтобы объединять ионные каналы в зависимости от способа регуляции их действия. Именно так систематизированы ионные каналы в предыдущих разделах. Другой вариант заключается в разделении ионных каналов в соответствии с основным проводимым ионом, однако оба варианта имеют очевидные недостатки .

Например, nAChRs (Itier and Bertrand, 2001) проводит ионы натрия, как и потенциалуправляемые натриевые каналы. Однако он в корне отличается от них по механизму активации, а также по аминокислотной последовательности. Другой пример, который иллюстрирует необходимость создания новой системы классификации, состоит в том, что некоторые каналы, управляемые циклическими нуклеотидами (Kaupp and Seifert, 2002), настолько гомологичны по последовательности потенциал-управляемым каналам, что могут быть эволюционно ближе к ним, чем к лиганд-управляемым каналам. Ещё одним примером служит то, что глутаматные рецеторы (Traynelis, Wollmuth et al., 2010), которые часто относят к семейству лиганд-управляемых ионных каналов, не имеют никакого сходства последовательностей и, следовательно, вероятно, не связаны эволюционно с другими каналами этой группы .

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

1.2.1. P-loop каналы Суперсемейство так называемых “P-loop” каналов включает потенциал-управляемые калиевые, натриевые, кальциевые каналы, каналы, активируемые циклическими нуклеотидами, каналы ионотропных рецепторов глутамата, и TRPV каналы (Zhorov and Tikhonov, 2004, Yue, Xie et al., 2015). Несмотря на огромное разнообразие свойств, функций и механизмов активации, все эти каналы имеют главную отличительную особенность - домен, образующий ионную пору в мембране, который формируется четырьмя субъединицами (или псевдосубъединицами в случае эукариотических натриевых и кальциевых каналов). От каждой субъединицы в этом домене присутствует две трансмембранных -спирали (S5 и S6 у потенциал-управляемых каналов), соединенных так называемой Р-петлей (P-loop), которая включает аминокислоты, образующие селективный фильтр (рисунок 1). Наличие этой петли и послужило основанием для классификации – все разнородное семейство называют p-loop каналами (Zhorov and Tikhonov, 2004). Именно в петле, которая формирует самую узкую часть канала расположены остатки, образующие так называемый селективный фильтр – структуру, обеспечивающую выборочное пропускание конкретных ионов. Как и следовало ожидать, аминокислотный состав селективных фильтров различен для каналов, которые пропускают разные ионы .

Рисунок 1. Трёхмерная структура p-loop канала Kv1 .

2. Цилиндры отображают спиральные участки канала. Оранжевым показаны участки p-спирали (p-loop), голубым – внутренние спирали (inner helices, s6), жёлтым – внешние спирали (outer helices, s5) .

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

Рисунок 2. Схематическое представление канала из суперсемейства p-loop .

Четыре домена с S1 до S4 отвечает за потенциал-чувствительность, в то время как S5 и S6, вместе с соединяющей их p-петлёй образуют непосредственно пору канала .

В последние десятилетия более точная информация о структурных особенностях ионных каналов получается посредством рентгеновской кристаллографии. Фактически, для ряда представителей этого суперсемейства сама идея гомологии пришла не из сопоставления последовательностей, а из структурных данных. К примеру, каналы TRPV1 и Kv1.2 относятся к суперсемейству Р-loop каналов, однако, как видно из выравнивания, они обладают весьма низкой гомологией в области селективного фильтра (рисунок 3). В то же время, наложение рентгеновских структур показывает, что каналы обладают схожей топологией в области поры (рисунок 3) .

Для удобства работы с белками, имеющими различную нумерацию аминокислотных остатков, использовалась универсальная индексация, предложенная Жоровым и Тихоновым. В данной индексации для обозначения аминокислотного остатка используется буквенное обозначение сегмента (“o” – outer helix, соответствующая сегменту S5, “i" – inner helix, соответствующая сегменту S6, “p” – Р-loop, внутримембранная петля между S5 и S6) и номер остатка внутри сегмента .

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

Однако, несмотря на общую структурную организацию, внутри суперсемейства Рloop каналов существуют и свои важные отличия. Особенно интересны в этом плане структуры прокариотических натриевых каналов, полученные недавно (Payandeh, Scheuer et al., 2011, Zhang, Ren et al., 2012). Их отличительная особенность заключается в том, что область Р-петли состоит не из одной спирали, как в калиевых каналах (Doyle, Morais Cabral et al., 1998), а из двух спиралей, разделённых селективным фильтром (рисунок 4) .

Различия в структуре порообразующих сегментов также объясняют наличие уникальных сайтов связывания лекарств и различных нейротоксинов, которые действуют на конкретные подклассы Р-loop каналов .

Таким образом, суперсемейство p-loop каналов объединяет каналы с похожей пространственной организацией поры. Это позволяет получать дополнительные данные о структуре слабо изученных каналов .

Рисунок 3. Пример наложения двух p-loop каналов: TRPV1 и Kv1 .

2. Сверху: выравнивание последовательностей между TRPV1 и Kv1.2. Снизу: трёхмерное наложение структур Kv1.2 (зелёным) и TRPV1 (лиловым); несмотря на низкую гомологию последовательностей, трёхмерное выравнивание показывает близость структур .

Рисунок 4. Трёхмерная структура p-loop канала NavMs .

Цилиндры отображают спиральные участки канала. Оранжевым показаны участки p-спирали (p-loop), голубым – внутренние спирали (inner helices, s6), жёлтым – внешние спирали (outer helices, s5). Примечательно, что в области pпетли находятся две спирали: p1 и p2 .

1.2.2. Потенциал-управляемые натриевые каналы Одной из важных проблем является изучение структуры потенциал-управляемых натриевых каналов. Натриевые каналы обнаружены практически во всех возбудимых и многих невозбудимых клетках. В возбудимых клетках потенциал-активируемые Na-каналы участвуют в формировании потенциала действия и обеспечивают достаточно высокую скорость его распространения по мембране (Зефиров, 1997). Самая высокая плотность Naканалов определяется в перехватах Ранвье миелинизированных нервных волокон, более низкая - в мембранах сомы нервной клетки, нервных окончаний, мышечных и эндокринных клеток. В пресинаптических структурах Na-каналы, участвуя в формировании потенциала действия, регулируют количество входящих ионов Са и, соответственно, количество освобождаемых квантов медиатора (Katz, 1971), а также синхронизируют секрецию медиатора во времени .

Детальный анализ функционирования потенциал-активируемых Na-каналов методами фиксации потенциала на гигантских аксонах беспозвоночных и миелинизированных нервных волокнах позвоночных животных прив ёл к созданию модели работы Na-канала (Hille, 1992). Воротный механизм Na-канала характеризуется 4 процессами: активация при деполяризации, инактивация при длительной деполяризации, деактивация после реполяризации и реактивация канала при его выходе из инактивированного состояния. Na-каналы закрыты при нормальных значениях мембранного потенциала покоя и открываются на очень короткое время (1 мс и менее) при деполяризации .

Первым шагом в выявлении структуры Na-каналов было выделение молекул каналов с использованием нейротоксинов, специфически связывающихся с каналом. В результате были выделены три субъединицы, которые присутствуют в равных пропорциях в Na-канале: одна представляет собой большой гликопротеин с молекулярной массой 280 кДа (1) и два малых полипептида с молекулярными массами 36 кДа (1) и 33 кДа (2) (Catterall 2000). 1-субъединица является основной и универсальной, тогда как малые субъединицы различаются в разных типах тканей .

Биохимические исследования обнаружили, что -субъединица состоит из четырёх повторов, нумеруемых от I до IV, каждый из которых содержит 6 трансмембранных сегментов, обозначаемых от S1 до S6 (рисунок 5). Сегмент S4 отвечает за чувствительность канала к изменению мембранного потенциала, и потому является одним из наиболее высококонсервативных участков. Кроме того, показано, что уменьшение общего положительного заряда в сегменте S4 путем сайт-направленного мутагенеза снижает потенциал-зависимость активационных ворот, а антитела, связывающиеся с тем же участком канала, замедляют инактивацию Na-каналов .

Рисунок 5. Строение и субъединиц Na-канала .

Трансмембранные спиральные сегменты изображены в виде цилиндров. S4-сегмент содержит положительно заряженные аминокислоты и является сенсором потенциала; h - остатки, необходимые для быстрой инактивации; P в кружочках

– сайты фосфорилирования цАМФ-зависимой протеинкиназой; P в ромбике - сайты фосфорилирования протеинкиназой С. Маленькие белые кружочки – остатки аминокислот с указанием их заряда, необходимые для высокоаффинного связывания с тетродотоксином (Siegel, 2011) .

Положительно заряженные остатки аминокислот на одной стороне S4-спирали позволяют ей определять изменения электрического потенциала на мембране .

Деполяризация мембраны воздействует на заряженные аминокислоты таким образом, что происходит перемещение спирали, которое, в свою очередь, приводит к открытию поры канала. Предполагается, что спираль вращается для того, чтобы пора могла открыться (Yang and Horn, 1995) .

Однако, в отличие от многих типов каналов, потенциал-управляемые натриевые каналы эукариот до сих пор не закристаллизованы и их точная пространственная структура остаётся неизвестной .

1.3. Фармакология ионных каналов 1.3.1. Фармакология потенциал-управляемых Na-каналов .

Потенциал-управляемые натривые каналы являются мишенями для многих токсинов и медицински важных препаратов. Несмотря на десятилетия интенсивных исследований, механизмы действия препаратов все ещё недостаточно хорошо изучены на структурном уровне. Основной причиной является именно отсутствие структур высокого разрешения эукариотических каналов и их комплексов с лигандами. Для решения этих проблем в частности используются методы гомологического моделирования, которое позволяет изучать действие лигандов на молекулярном уровне, используя имеющиеся рентгеновские структуры в качестве шаблонов. Из-за ограничений, присущих этому подходу, модели должны быть проверены на соответствующих экспериментальных данных. С другой стороны, такие модели дают уникальную возможность интегрировать различные экспериментальные данные в рамках одной модели .

Использование гомологического моделирования позволило лучше изучить действие сразу нескольких классов лигандов. Это в первую очередь тетродотоксин, который блокирует внешнюю пору, местные анестетики, блокирующие внутреннюю пору, батрахотоксин, связывающийся во внутренней поре, но, как это ни парадоксально, активирующий канал, пиретроидные инсектициды, которые активируют канал, связываясь со стороны липидного слоя в интерфейсе между субъединицами и альфа и бета-токсины скорпионов, которые связываются в промежутке между порой и сенсором мембранного напряжения, влияя на активацию (рисунок 6). Рассмотрим подробнее достижения в области изучения механизмов действия этих лигандов .

Рисунок 6. Схематическое расположение сайтов связывания различных токсинов в потенциал-управляемом натриевом ионном канале .

Сверху – схематическое изображение структуры, снизу – трёхмерное схематическое изображение. Оранжевым показаны участки S1-S3 спиралей, синим – сенсор напряжения S4, голубым – линкер между S4 и S5, жёлтым – внешняя спираль S5, зелёным – петлевая часть поры p-loop, лиловым – внутренняя спираль S6 (Tikhonov, Zhorov, unpublished) .

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

К лигандам, действующим на внешнюю пору, в первую очередь относятся тетродотоксин (TTX), сакситоксин (STX) и их аналоги. Интенсивные исследования методами сайт-направленного мутагенеза показали, что связывание этих токсинов происходит в основном за счёт остатков селективного фильтра и нескольких остатков рядом с ним. Среди них ключевая роль принадлежит кольцу внешних карбоксилатов, т.е .

кислотных остатков в трёх-четырёх позициях вправо по последовательности от остатков селективного фиьтра. Детальный анализ позволил обнаружить контакты между некоторыми остатками каналов и отдельными химическими группами TTX и STX (Chang, French et al., 1998, Xue, Ennis et al., 2003). Для некоторых контактов была определена свободная энергии взаимодействия (Li, Ennis et al., 2001, Choudhary, Aliste et al., 2007). Хотя сайты связывания различных токсинов не идентичны, по существу они перекрывают друг друга. В частности, важную роль играют внешние карбоксилаты, которые также служат сайтами связывания для проникающих ионов натрия (Khan, Romantseva et al., 2002) .

Существует интригующая разница между действием TTX / STX, которые полностью блокируют ток через канал, и µ-конотоксинами, которые могут приводить к неполному блоку .

Первая модель поры натриевого канала в комплексе с TTX и STX была предложена Липкиндом и Фоззардом до получения кристаллической структуры P-loop каналов (Lipkind and Fozzard, 1994). Важную роль в разработке модели сыграли экспериментальные данные, согласно которым действие токсинов резко ослабевало при мутациях аспартата и глутамата в селективном фильтре (I и II субъединицы), а также внешних карбоксилатов в повторах I, II, и IV (Terlau, Heinemann et al., 1991). Согласно предложенной модели, гуанидиновая группа TTX связывалась с селективным фильтром, а карбоксилат повтора II формировал водородную связь с гидроксильной группой TTX. Также модель показала, что при рассмотрении канала с внеклеточной стороны, субъединицы I, II III и IV расположены по часовой стрелке, и это важное предсказание было подтверждено позднее методами циклического мутационного анализа взаимодействия канала с µ-конотоксином GIIIA (Dudley, Chang et al., 2000) .

С появлением рентгеновской структуры KcsA, Липкинд и Фоззард разработали гомологическую модель Nav1.4 на основе KcsA и использовали экспериментальные данные о токсин-канальных взаимодействиях для посадки TTX и STX в модель (Lipkind and Fozzard, 2000). Поскольку ширина внешней поры KcsA недостаточна, чтобы вместить крупные полужёсткие токсины, авторы поместили остатки селективного фильтра на Cконцах P-спиралей и предположили, что P-спирали сильнее смещены от оси поры, по сравнению с калиевыми каналами. В этой модели кольцо DEKA находится у границы между центральной полостью и внешней порой. Чтобы обеспечить структурное обоснование многочисленных экспериментально определённых контактов между большими пептидными µ-конотоксинами и натриевыми каналами, авторы вынуждены были дополнительно модифицировать модель натриевого канала путём сдвига P-спиралей дальше от оси поры и увеличения их наклона (Choudhary, Aliste et al., 2007) .

Альтернативные гомологические модели Nav1.4 с TTX и STX на основе KcsA предполагали сильное сходство между трёхмерными структурами натриевых и калиевых каналов (Tikhonov and Zhorov, 2005, Tikhonov and Zhorov, 2011). Эти модели были построены, чтобы объяснить те же экспериментальные данные по взаимодействию токсинов с каналом, но остатки селективного фильтра находились в середине восходящих участков, а не на изгибе Р-петли. Это позволило посадить TTX и STX во внешнюю пору без смещения P-спирали относительно калиевых каналов (рисунок 7) .

Появление рентгеновских структур бактериальных натриевых каналов (Payandeh, Scheuer et al., 2011, Zhang, Ren et al., 2012) свидетельствует в пользу моделей (Tikhonov and Zhorov, 2005, Tikhonov and Zhorov, 2011). Положение глутаматов селективного фильтра согласуется с предложенным в моделях, а P-спирали располагаются аналогично калиевым каналам. Тем не менее, обе модели (Lipkind and Fozzard, 2000, Tikhonov and Zhorov, 2005) не смогли предсказать пространственную укладку С-концевых половин P-петли. Структура NavAb показывает, что, в отличие от калиевых каналов, С-концевая часть P-петли содержит P2 спираль, направленную во внеклеточную сторону, в результате чего внешняя пора оказывается достаточно широкой, чтобы вместить большие токсины .

Остается неясным, в какой степени рентгеновские структуры гомотетрамерических прокариотических натриевых каналов близки к структурам эукариотических каналов. В частности, симуляция методами молекулярной динамики взаимодействия TTX с моделью Nav1.4 на основе прямого выравнивания с NavAb предсказывает, что длинная ось TTX перпендикулярна относительно оси поры и TTX образует водородные связи с внешними карбоксилатами в повторах I, II и III (Chen and Chung, 2014). Эти результаты плохо согласуются с экспериментальными данными по взаимодействию TTX с натриевым каналом. В частности, модель не объясняет значительный вклад D4p54 и слабый вклад D3p53 в действие токсина .

В модели Nav1.4 на основе NavAb, которая была построена с помощью прямого выравнивания последовательностей Р-петель (Tikhonov and Zhorov, 2012), внешние карбоксилаты также не в состоянии взаимодействовать с TTX так эффективно, как следует из мутационных данных. Возможное решение проблемы было найдено с помощью скорректированного выравнивания последовательности между бактериальными и эукариотическими натриевыми каналами. Предлагаемые делеции в последовательности вблизи кольца DEKA (Tikhonov and Zhorov, 2012) обеспечивают ориентацию внешних карбоксилатов в согласии с экспериментальными данными (рисунок 7) .

Рисунок 7. Модель натриевого канала Nav1 .

4 в комплексе с TTX. Слева: модель на основе KcsA, данные из работы (Tikhonov and Zhorov 2005); справа: модель на основе NavAb, данные из работы (Tikhonov and Zhorov, 2012) .

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

Особенно интересен механизм блокады натриевых каналов местными анестетиками (Catterall, 2000, Ruetsch, Boni et al., 2001, Nau and Wang, 2004, Fozzard, Lee et al., 2005, Clare, 2010). Несмотря на простую структуру этих лигандов и общий принцип механизма блока поры, феноменология их действия весьма сложна. Достаточно слабый "тонический блок" наблюдается при редкой ативации канала. Активация посредством длинной деполяризации усиливает блокирующее действие, этот эффект называется "частотно-зависимый блок". Для объяснения этих эффектов, Бертиль Хилль предложил "гипотезу модулированного рецептора" (Hille, 1977). Гипотеза предполагает слабое сродство ЛА с каналом в закрытом состоянии и значительное увеличиние сродства, когда канал переходит в открытое или инактивированное состояние (Hondeghem and Katzung, 1977, Hondeghem and Katzung, 1984). Эта гипотеза также предлагает существование двух путей доступа молекул ЛА к сайту связывания: гидрофильный путь через открытые ворота активации и гидрофобный путь через мембрану. Исследования методами мутагенеза позволили обнаружить остатки в S6 сегменте субъединиц I, III и IV, которые контролируют действие ЛА (Wang, Quan et al., 1998, Yarov-Yarovoy, Brown et al., 2001, Yarov-Yarovoy, McPhee et al., 2002). Это наводит на мысль, что связывание ЛА происходит во внутренней области поры. Интересно, что влияние отдельных замен на тонический блок, частотно-зависимый блок, и гидрофобный доступ неравны (Li, Galue et al., 1999) .

Антагонизм между проникающим ионом и заряженной частицей, блокирующей пору, не удивителен. Проникающие ионы могут влиять на блокаду канала либо путём конкуренции с блокаторами за сайт связывания внутри поры, или за счёт электростатического отталкивания от близко расположенных сайтов (Hille and Schwarz, 1978, Zhorov and Tikhonov, 2004, Zhorov and Tikhonov, 2013). Различные типы лигандионных взаимодействий серьезно осложняют картину зависимости действия лиганда на ионные каналы от состояния и напряжения (Hille and Schwarz, 1978, Triggle, 2007, Bruhova, Tikhonov et al., 2008, Wulff, Castle et al., 2009) .

Основная проблема с моделированием ЛА и родственных препаратов заключается в получении структурного обоснования сложного механизма действия, который описан в "гипотезе модулированного рецептора" (Hille, 1977). Гидрофильный путь доступа, очевидно, включает в себя открытые активационных ворот, которые находятся с внутриклеточной стороны относительно сайта связывания ЛА. Тем не менее до сих пор не обнаружен точно гидрофобный путь, который обеспечивает доступ ЛА в закрытый канал .

Остатки, которые контролируют этот путь были найдены в середине S6 и в P-петле (Qu, Rogers et al., 1995, Sunami, Dudley et al., 1997, Wang, Quan et al., 1998, Sunami, Glaaser et al., 2000) .

Так как ЛА имеют большее сродство к открытому/инактивированному каналу, чем к каналу в закрытом состоянии, Липкинд и Фоззард использовали первую опубликованную структуру открытого калиевого канала, MthK, в качестве шаблона для моделирования открытого состояния натриевого канала (Lipkind and Fozzard, 2005). Так как степень сходства структуры между MthK и натриевыми каналами неизвестна, авторы не пытались найти энергетически оптимальную конформацию, поскольку такое исследование критически чувствительно к структурным деталям, а исследовали набор различных конформаций, которые соответствуют экспериментальным данным по мутациям и взаимодействиям. В частности, Липкинд и Фоззард предположили, что алкиламинные и ароматические фрагменты ЛА находятся, соответственно, около фенилаланина и тирозина в D4S6, чувствительных к молекулам ЛА, ранее обнаруженым Уильямом Катераллом и соавторами (Ragsdale, McPhee et al., 1994). Были получены конформации, которые позволяли объяснить экспериментальные данные. Это позволило сделать важный вывод, что ЛА не полностью занимают пространство вунтри поры (рисунок 8). А электростатические расчёты позволили предположить, что блок должен иметь в первую очередь электростатическую природу (Lipkind and Fozzard, 2005). Авторы пришли к выводу, что лидокаин и ему подобные молекулы слишком велики, чтобы вписаться в модель закрытого канала на основе KcsA аналогично тому, как они располагаются в модели открытого канала на основе MthK .

Рисунок 8. Модель открытого натриевого канала в комплексе с мепивакаином на основе MthK .

Данные из работы (Lipkind and Fozzard, 2005) .

В последующем была предложена модель открытого Nav1.4 на основе KvAP (Tikhonov and Zhorov, 2007), в которой оптимальное положение лидокаина схоже с моделью на основе MthK (Lipkind and Fozzard, 2005). В модели канала в закрытом состоянии, ЛА не могут поместиться в центральную полость и их ароматическое кольцо проникает в интерфейс между P-спиралью и соседними сегментами S6 (Bruhova, Tikhonov et al., 2008) .

Было отмечено, что некоторые остатки, контролирующие гидрофобный путь, смотрят в тот же интерфейс. Результаты моделирования позволили предположить, что интерфейс между субъединицами III и IV служит в качестве гидрофобного пути доступа для ЛА в закрытом натриевом канале (Tikhonov and Zhorov, 2005, Tikhonov and Zhorov, 2007, Bruhova, Tikhonov et al., 2008) .

Модели натриевых каналов, которые основаны на шаблонах калиевых каналов, имеют два слабых места. Во-первых, интерфейсы между субъединицами калиевых каналов достаточно узкие и возможность доступа ЛА через соответствующие интерфейсы в моделях натриевых каналов находится под вопросом. Во-вторых, некоторые остатки, которые влияют на доступ препаратов в закрытом канале, смотрят в сторону селективного фильтра, а не в интерфейс (Sunami, Dudley et al., 1997, Sunami, Glaaser et al., 2000) .

В работе 2008 года Брухова, Тихонов и Жоров, используя гомологическое моделирование, предположили, что гидрофобный путь пролегает через боковое ответвление между S6 спиралями III и IV субъединиц канала, в области p-спирали III субъединицы (Bruhova, Tikhonov et al., 2008). Это предположение позволило создать модель, в которой молекулы локальных анестетиков связываются с закрытым натриевом каналом в горизонтальной ориентации, при этом ароматическая часть лиганда оказывается в интерфейсе между III и IV субъединицами, а аминогруппа попадает в пору. При таком связывании молекула анестетика должна вытеснить ион натрия из центральной полости и электростатически воздействовать на заселённость селективного фильтра .

Недавно полученные рентгеновские структуры натриевых каналов служат подтверждением предложенного в моделях гидрофобного пути доступа. Субъединичные интерфейсы этих каналов значительно шире, чем в калиевых каналах, что может позволить большим ЛА проникать в поры. Особый интерес представляет структура бактериального натриевого канала, NavMs, с бромфенил содержащим блокатором (Bagneris, DeCaen et al., 2014). Из-за ограниченного разрешения в структуре видно только атом брома лиганда, но его положение в интерфейсе субъединицы совпадает с атомом азота, в модели 2008 года (Bruhova, Tikhonov et al., 2008). Кроме того, в присутствии блокатора ионная заселённость селективного фильтра уменьшается, что говорит о вытеснении ионов лигандом. Это согласуется с гипотезой о "взаимном влиянии" органических и неорганических катионов, в соответствии с которой молекула ЛА, которая входит в закрытый канал через межсубъединичный интерфейс, должна вытеснить ион из центральной полости, который может покинуть канал только через селективный фильтр (Bruhova, Tikhonov et al., 2008) .

Эта модель объясняет данные, согласно которым мутации в наружной поре или блок TTX приводят к предотвращению блока закрытого канала ЛА (Bruhova, Tikhonov et al., 2008) .

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

Природные вещества батрахотоксин, аконитин, вератридин и грайанотоксин известны как стероидные активаторы натриевых каналов (агонисты). Сайты и механизмы действия этих токсинов, возможно, одна из самых интригующих загадок, связанных с лигандами натриевых каналов. В присутствии таких соединений порог активации смещается в сторону отрицательных значений так, что каналы открываются вблизи мембранного потенциала покоя. Кроме того, это делает невозможной инактивацию и каналы остаются постоянно открытыми. Столь сложное действие этих липидорастворимых молекул на активацию и проницаемость изначально вдохновило предложение об аллостерическом характере действия лигандов, которые связываются на границе липидного слоя (Catterall, 1979, Hille, 2001). Тем не менее, последующие данные мутационного анализа не подтвердили это предположение. Остатки, чувствительные к BTX были найдены в S6 сегментах всех четырёх субъединиц канала. Более того, схема остатков, чувствительных к BTX сильно перекрывается со схемой остатков, чувствительных к ЛА, тем самым предполагая, что сайт связывания находится в поре. Некоторые прочие характеристики активаторов и блокаторов также схожи (Wang and Wang, 2003). Например, крупные блокаторы предотвращают закрытие канала (Yeh and Narahashi, 1977, Yamamoto, 1986) .

В разработанных моделях BTX, вератридин и аконитин предположительно связываются во внутренней поре, не блокируя проникновение ионов (Tikhonov and Zhorov, 2005). В этих моделях активаторы не перекрывают путь ионов полностью. Напротив, они обращены катионофильной стороной (так называемые триады кислорода) в пору, что позволяет ионам проникать между агонистом и внутренней стенкой поры. Кроме того, связывание этих относительно больших токсинов предотвращает закрытие канала. Эта модель и её последующее развитие (Du, Garden et al., 2011) предполагают, что BTX действует в ионном канале аналогично тому, как хируригический стент действует в кровеносном сосуде .

На основании предсказаний модели были проведены мутационные исследования, в которых были определены несколько новых остатков, чувствительных к BTX (Wang, Mitchell et al., 2006, Wang, Tikhonov et al., 2007, Du, Garden et al., 2011). Особый интерес представляет собой замена аспарагина в D2S6 на лизин, который преобразует действие BTX от активации канала к необратимой блокаде (Wang, Tikhonov et al., 2007). Это показывает общее сходство в сайтах и характере связывания стероидных молекул, которые, в зависимости от их структуры и аминокислотного состава внутренних спиралей канала, может либо блокировать, либо активировать натриевые каналы .

Другой тип активаторов натриевых каналов включает в себя пиретроиды, большой класс синтетических аналогов природных пиретринов из цветков pyrethrum daisy (Dong, Du et al., 2014). Пиретроиды, которые имеют относительно низкую токсичность для млекопитающих и относительно благоприятные экологические свойства, широко используются в качестве одной из наиболее эффективных мер контроля в глобальной борьбе против вредных в сельском хозяйстве вредителей и заболеваний, включая малярию и лихорадку денге. Другим хорошо известным агонистом натриевого канала является ДДТ, инсектицид, который, как полагают, имеет сайт связывания, перекрывающийся с пиретроидным рецептором. Интенсивное использование ДДТ, а затем и пиретроидов в течение нескольких десятилетий привело к развитию резистентности во многих популяциях насекомых. Основной механизм резистентности, известный как KDR, был зарегистрирован почти во всех медицински и агрономически важных членистоногих вредителях и в настоящее время более 50 мутаций натриевых каналов связаны с сопротивлением к перитроидам в различных членистоногих вредителях (Silver, Du et al., 2014). Поэтому до сих пор активно ведётся поиск пиретроидных сайтов связывания с целью разработки новых инсектицидов .

Новаторская модель (O'Reilly, Khambay et al .

, 2006) интегрировала мутационные данные, в том числе резистентные мутации, и предсказала связывание пиретроидов в интерфейсе между повторами II и III. Анализ дополнительных экспериментальных данных привел к моделированию рецептора пиретроидов между повторами I и II, и на основании этой модели были проведены исследования методами мутагенеза, что позволило обнаружить новые остатки, отвечающие за связывание (Du, Nomura et al., 2013). К последним относятся остатки, которые отличаются между натриевыми каналами насекомых и млекопитающих, давая тем самым молекулярное основание для избирательной токсичности пиретроидов. Недавно разработана модель натриевого канала комаров с двумя пиретроидными рецепторами, PyR1 и PyR2, показавшая наличие квазисимметрии вокруг оси поры (Du, Nomura et al., 2015). Мутагенез на основе данных модели привёл к открытию семи новых остатков, связывающихся с пиретроидами. Это позволяет предположить, что каждый рецептор расположен между четырьмя спиралями (L45, S5 и два S6), принадлежащими соседним повторам, причём четыре остатка D2S6 относятся к PyR1 и ещё четыре остатка относятся к PyR2 (Du, Nomura et al., 2015). Предполагается, что связывание лигандов как с PyR1, так и с PyR2 необходимо для активации канала .

Из рентгеновских структур потенциал-управляемых гомотетрамерных каналов известно, что внеклеточные участки S1-S2 и S3-S4 петель данного повтора близки к S5-S6 петле смежного повтора. Пептидный токсин, связавшийся с этими петлями между областью поры и областью сенсора напряжения может зафиксировать последний в определённом состоянии и, таким образом, модифицировать гейтинг канала. К таким токсинам относятся токсины скорпионов (Wang, Yarov-Yarovoy et al., 2011, Zhang, Yarov-Yarovoy et al., 2011, Zhang, Yarov-Yarovoy et al., 2012), пауков (Bosmans and Swartz, 2010, Minassian, Gibbs et al.,

2014) и актиний (Xiao, Blumenthal et al., 2014). Альфа- и бета-токсины скорпионов связываются с различными участками и вызывают различные эффекты. Альфа-токсины замедляют инактивацию, а бета-токсины вызывают отрицательный сдвиг активации путём захвата сенсора напряжения II повтора в активированном состоянии. Внеклеточные S1-S2 и S3-S4 участвуют в образовании рецепторов альфа- и бета-токсинов, но связывающие области расположены на противоположной стороне канала и включают IV повтор для альфа токсинов и II повтор для бета токсинов (Thomsen and Catterall, 1989, Rogers, Qu et al., 1996, Cestele, Yarov-Yarovoy et al., 2006). Удивительно, но остатки, чувствительные к бета-токсину найдены не только в IIS1-S2 и IIS3-S4, но и в IIIP2-S6 (Zhang, Yarov-Yarovoy et al., 2012) .

Сочетание мутационных и молекулярных исследований с моделированием позволило определить картину действия токсина скорпионов на молекулярном уровне (Catterall, Cestele et al., 2007). Главная проблема заключалась в моделировании гибких участков каналов, которые не образуют регулярную вторичную структуру. Для решения этой проблемы авторы сгенерировали набор энергетически возможных моделей и выбрали модель, в которой ключевые остатки канала, связанные с токсином, образовывали кластеры (Zhang, Yarov-Yarovoy et al., 2012). Клиновидная форма бета-токсина, CssIV, предположительно связывается глубоко в области между петлями IIS1-S2 и IIS3-S4 .

Посадка токсина в этой модели воспроизводит большинство экспериментально известных взаимодействий. Однако один из важных экспериментально обнаруженных остатков, D1445 (Zhang, Yarov-Yarovoy et al., 2012), находится за пределами предлагаемого сайта связывания. Это может свидетельствовать либо об ограниченной точности модели, либо о косвенном влиянии замещения остатков на действие токсина .

Аналогичный анализ был проведен для альфа-токсинов, которые по-видимому связываются квази-симметрично относительно бета-токсинов (рисунок 9) (Wang, YarovYarovoy et al., 2011). Моделирование показывает, что они проникают в коническую щель между петлями IVS1-S2 и IVS3-S4. Разница между основными эффектами альфа-токсинов (ингибирующими быструю инактивацию) и бета-токсинов (усиливающими активацию) может быть связана с функциональной и структурной асимметрией канала. В частности, известно, что сенсоры напряжения в повторах I и II быстро реагируют в ответ на деполяризацию, в то время как сенсор напряжения в IV повторе реагирует медленно, и его время активации аналогично времени быстрой инактивации (Chanda and Bezanilla, 2002) .

Рисунок 9. Схематическое изображение сайтов связывания для - и -токсинов скорпиона (Tikhonov, Zhorov, unpublished) .

Связывание лиганда в одной и той же области может привести к различным, а иногда и противоположным эффектам. Например, как ЛА, так и BTX действуют на внутреннюю область поры, но ЛА блокируют канал, тогда как BTX активирует его. Другие примеры могут включать токсины, которые связываются с внешними петлями. Так, поначалу предполагалось, что токсины пауков, такие как Huwentoxin-IV из Ornithoctonus huwena и некоторые Hainantoxins связываются во внешней поре, аналогично TTX, так как они ингибируют Na+ проводимость (Li, Xiao et al., 2003). Однако, остатки, чувствительные к Huwentoxin и остатки, влияющие на бета-токсины скорпионов, находятся в одной области (Leipold, DeBie et al., 2007). Моделирование преполагает, что Huwentoxin-IV связывается в области между IIS1-S2 и IIS3-S4 (Minassian, Gibbs et al., 2014), что сходно с моделью связывания альфа-токсина скорпиона. По-видимому, главное отличие в том, что бетатоксины скорпиона фиксируют сенсор напряжения в активированном состоянии, в то время как Huwentoxin-IV и связанные с ним токсины фиксируют сенсор напряжения в состоянии покоя. Подробная информация об этих специфических взаимодействиях требует дальнейшего анализа .

Таким образом, существует большое количество токсинов, воздействующих на натриевые каналы. Их действие весьма различно, как и механизмы воздействия. Какие-то являются хорошо изученными, а какие-то изучены хуже. Одним из наиболее интересных семейств токсинов с точки зрения наличия экспериментальных данных является семейство конотоксинов .

1.3.2. Конотоксины Конотоксины представляют собой природные нейротоксины, выделенные из яда хищных морских моллюсков рода Conus. Яд моллюска обычно содержит смесь различных конотоксинов, и применяется моллюском для защиты и нападения на жертву. По характеру действия, конотоксины подобны постсинаптическим токсинам змей, но могут значительно превосходить их по токсичности. Большинство конотоксинов способны избирательно и с высокой специфичностью связываться и модулировать работу различных ионных каналов воздействуя на передачу импульса по нервному волокну (Terlau and Olivera, 2004). При этом детальный механизм действия для большинства соединений остается не до конца ясным .

Состав яда может различаться в зависимости от рациона улитки, таким образом, детальный состав яда может сильно различаться даже между особями одного вида .

По причине высокой специфичности блокирования нервной проводимости, конотоксины широко используются в качестве избирательно действующих агентов при электрофизиологических и клинических исследованиях механизмов передачи возбуждения в нервной системе, в качестве основы для производства обезболивающих препаратов (Wang and Chi, 2004, Olivera and Teichert, 2007). В связи с высоким терапевтическим и диагностическим потенциалом конотоксинов, большой интерес представляет более детальное изучение механизмов действия подобных биологических активных веществ и производство на их основе новых соединений с помощью современных методов био и нанотехнологий .

Структурно конотоксины представляют собой пептиды, состоящие из 10-30 аминокислотных остатков, и имеющие обычно одну, две или более дисульфидных связей (рисунок 10). Структурное разнообразие конотоксинов определяется их высокой специфичностью, и разнообразием возможных механизмов действия. Мишенями для различных групп конотоксинов могут быть различные семейства ионных каналов, помимо этого конотоксины могут иметь разные сайты связывания с одним и тем же типом ионного канала .

Рисунок 10. Структура конотоксинов. Строение -конотоксина MVIIA (слева) и -конотоксина PnIB (справа). Структура окрашена от N-конца (синий) к C-концу (красный), жёлтым показаны дисульфидные связи Яд улиток содержит смесь конотоксинов (до 50), при этом, состав яда зависит от рациона и таким образом может различаться даже у представителей одного вида. Таким образом, 500-700 существующих видов Conus, могут синтезировать до 50000 различных конотоксинов (Olivera, 2006). Общепринятая классификация подобного разнообразия представлена на рисунке 11 .

Рисунок 11. Классификация конотоксинов .

Конотоксины делятся на две основные группы: с одной дисульфидной связью, и богатые дисульфидными связями пептиды (2 и более). Большинство известных конотоксинов принадлежат второй группе, являющейся наиболее многочисленной .

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

Изучение конотоксинов показало значительное разнообразие их воздействия:

некоторые блокируют передачу сигнала от нервов к мышцам, другие способны вызывать судороги, сон или снимать болевой синдром. Детальный механизм их действия не всегда известен, хотя известно, что большинство способно модулировать работу ионных каналов клеточных мембран (Terlau and Olivera, 2004). В некоторых случаях, удается не только установить мишень токсина, но и определить сайт связывания лиганда и построить структуру комплекса токсина с лигандом. Так, например, для -конотоксина PnIA была определена детальная структура (Celie, Kasheverov et al., 2005) (рисунок 12) .

Рисунок 12. Структура комплекса -конотоксина PnIA с сайтом связывания в надмембранной части никотинового-ацетилхолинового рецептора. Токсин встраивается между 2-мя субъединицами рецептора, вызывая конформационные изменения в надмембранной части рецептора, которые влияют на проводимость ионного канала

Конотоксины делятся согласно соответствующему белку-мишени на следующиегруппы:

-конотоксины способные с высокой специфичностью блокировать различные типы никотиновых-ацетилхолиновых рецепторов;

-конотоксины ингибирующие инактивацию потенциал-зависимых натриевых каналов;

-конотоксины ингибирующие калиевые каналы;

-конотоксины способные блокировать потенциал-зависимые натриевые каналы;

-конотоксины способные блокировать потенциал-зависимые кальциевые каналы .

Высокое разнообразие по механизму действия делает конотоксины весьма перспективными объектами с медицинской точки зрения. Уже сейчас на их основе разрабатываются лекарства для лечения некоторых неврологических заболеваний, таких, например, как эпилепсия (Alonso, Khalil et al., 2003). Способность -конотоксинов ингибировать потенциал-зависимые кальциевые каналы N-типа, отвечающие за чувствительность к боли, позволила использовать их в качестве основы для производства обезболивающих средств, не вызывающих привыкания, для лечения хронических болей (Miljanich, 2004) .

Сравнительно небольшой размер, а также чрезвычайно высокая специфичность воздействия конотоксинов на различные компоненты нервной системы, делает их ценным инструментом при исследовании функций ионных каналов, механизмов передачи возбуждения в нервной системе на молекулярном уровне (Nicke, Wonnacott et al., 2004) .

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

Как видно из вышесказанного, фармакологические агенты, такие как конотоксины, играют ключевую роль в создании и разработке лекарств. Однако недостаток экспериментальных данных приводит к неполному пониманию механизмов действия подобных лигандов, что, в свою очередь, затрудняет исследования. Именно поэтому необходимо использовать теоретические методы, такие как молекулярное моделирование, для углублённого изучения механизмов связывания и поиска молекулярных детерминантов связывания лигандов с мишенями, в частности, с ионными каналами .

1.3.3. µ-конотоксины и Na-каналы Натриевые каналы эукариот, которые играют ключевую роль в физиологии возбудимых клеток, были предметом экспериментальных исследований в течение последних десятилетий (Trimmer, Cooperman et al., 1989, Priestley, 2004, Wood, Boorman et al., 2004, Rogers, Tang et al., 2006). Действие природных токсинов и синтетических препаратов на Nav1 каналы интенсивно изучается уже много лет. Эти токсины, как было рассмотренно ранее, включают в себя широко известные антагонисты натриевых каналов тетродотоксин, сакситоксин и µ-конотоксины, агонисты натриевых каналов батрахотоксин, вератридин, аконитин и грайанотоксин, местные анестетики, и многие другие лиганды (French, Worley et al., 1984, Moczydlowski, Hall et al., 1984, Linford, Cantrell et al., 1998, French and Terlau, 2004, Eijkelkamp, Linley et al., 2012, Waszkielewicz, Gunia et al., 2013) .

Выше была рассмотрена работа Тихонова и Жорова, в который было смоделировано действие TTX на модель натриевого канала Nav1.4, построенной на основании структуры NavAb (Tikhonov and Zhorov, 2012). Тем не менее, TTX взаимодействует только с областью селективного фильтра (DEKA-кольца) и кольцом внешних карбоксилатов, расположенным ближе к внеклеточному пространству относительно DEKA-кольца. Пока неясно, может ли структура прокариотических натриевых каналов быть использована для объяснения действия гораздо более крупных токсинов, которые нацелены на внешнюю пору .

Особенно интересным является действие -конотоксинов (Olivera and Cruz, 2001, Lewis, Dutertre et al., 2012). -конотоксины представляют собой довольно небольшие пептиды (16-26 аминокислот) с несколькими дисульфидными поперечными связями, которые стабилизируют типичную пространственную укладку, так называемую III (Lewis, Dutertre et al., 2012). С точки зрения структуры, конотоксины за счёт своего размера нацелены на другой участок канала, нежели TTX и STX (Zhang, Gruszczynski et al., 2010) .

Они не так сильно взаимодействуют с селективным фильтром, но зато взаимодействуют с внешними карбоксилатами и несколькими соседними с ними остатками .

Семья -конотоксинов включает в себя 3 токсина, которые особенно хорошо изучены: GIIIA, PIIIA и KIIIA (Cruz, Kupryszewski et al., 1989, Shon, Olivera et al., 1998, Zhang, Green et al., 2007). Наиболее изученным из них является GIIIA, токсин с высоким сродством к натриевому каналу Nav1.4 (Sato, Ishida et al., 1991, Becker, Prusak-Sochaczewski et al., 1992, Dudley, Todt et al., 1995, French, Prusak-Sochaczewski et al., 1996, Chang, French et al., 1998, Dudley, Chang et al., 2000, Li, Ennis et al., 2001, Cummins, Aglieco et al., 2002, Hui, Lipkind et al., 2002, Hui, McIntyre et al., 2003, Xue, Ennis et al., 2003, Choudhary, Aliste et al., 2007). PIIIA представляет собой близкий гомолог GIIIA, но с некоторыми отличительными особенностями, которые делают его интересным объектом для изучения (Safo, Rosenbaum et al., 2000, McArthur, Ostroumov et al., 2011, McArthur, Singh et al., 2011, Chen and Chung, 2012, Tietze, Tietze et al., 2012). И, наконец, KIIIA, самый короткий из них (всего 16 остатков), представляет особый интерес из-за своего небольшого размера и способности связываться с каналом одновременно с ТТХ (Zhang, McArthur et al., 2009, Van Der Haegen, Peigneur et al., 2011, Wilson, Yoshikami et al., 2011, Khoo, Gupta et al., 2012, Stevens, Peigneur et al., 2012) .

Предыдущие модели связывания -конотоксинов в поре канала (Choudhary, Aliste et al., 2007, McArthur, Singh et al., 2011, McArthur, Singh et al., 2011) были разработаны для объяснения экспериментальных данных, используя гомологические модели, основанные на структуре калиевого канала KcsA. Но так как NavAb имеет отличную от KcsA укладку в области селективного фильтра и наружной поры, то появляется возможность изучить альтернативную модель на предмет объяснения экспериментальных данных. Модель натриевого канала Nav1.4, построенная на основе NavAb, служит хорошей отправной точкой (Tikhonov and Zhorov, 2012), чтобы изучить связывание -конотоксинов GIIIA, PIIIA и KIIIA в поре канала .

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

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

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

биохимические методы;

фармакологические методы;

генетические методы;

рентгеноструктурный метод .

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

2.1.2. Рентгеноструктурный анализ Наиболее важным подходом, с точки зрения структурного анализа ионных каналов, несомненно, является метод рентгеноструктурного анализа. Использование метода рентгенокристаллографии совместно с биохимическими, генетическими, электрофизиологическими методами позволило выявить принципы структуры и функционирования канальных белков. Рентгенокристаллографический метод позволяет непосредственно оценить трехмерную структуру ионного канала. Этот метод стал относительно недавно использоваться для структурного анализа каналов, так как они содержат трансмембранные гидрофобные домены, что затрудняет их кристаллизацию .

Мак-Киннон и его коллеги впервые проанализировали молекулярную структуру ионного канала с помощью рентгенокристаллографического метода с высоким разрешением на примере бактериального К-канала, (Doyle, Morais Cabral et al., 1998) (рисунок 13) .

Рисунок 13. Структура калиевого канала, полученная с помощью рентгенокристаллического метода. А) Вид на канал с внешней стороны мембраны. Б) Вид на канал сбоку. В) Вид на канал сбоку, показаны две из четырех субъединиц, селективный фильтр образован карбоксильными группами трех остатков аминокислот – глицин (G), тирозин (Y) и глицин (G). Г) Схема канала в мембране, иллюстрирующая три калий-связывающих сайта внутри канала. Спирали поры образуют отрицательный диполь, который помогает стабилизировать ион К внутри поры (Doyle, Morais Cabral et al., 1998) .

Особенностью этих каналов является относительно небольшой размер и простая трансмембранная топология. Проведённый анализ трехмерной структуры показал, что канал образован четырьмя идентичными субъединицами, окружающими центральную пору (рисунок 13, А). Каждая субъединица состоит из двух трансмембранных участков, соединяющихся петлей – областью Р, которая формирует селективный фильтр канала (рисунок 13, Б). Четыре внутренних сегмента каждой субъединицы образуют область поры с цитоплазматической стороны. Форма и структура поры определяют ион-проводящие свойства. Внутреннее и внешнее устья поры выстланы кислыми аминокислотами, чьи отрицательные заряды притягивают катионы из внеклеточной среды .

Наиболее узкий участок поры формирует селективный фильтр, который ограничивает скорость прохождения ионов. Селективный фильтр образован тремя карбонильными атомами белкового скелета каждой из четырех субъединиц, отрицательный заряд которых обеспечивает полярное окружение для ионов К, проходящих через канал (Jiang, Lee et al., 2003) (рисунок 13В). Рентгенокристаллографический анализ показал, что c порой канала могут связываться три иона К, один из которых находится в широкой внутренней полости, и до двух ионов К могут одновременно занимать селективный фильтр (рисунок 13, Г) (Jiang, Lee et al., 2003) .

На данный момент получены рентгеновские структуры многих ионных каналов:

например калиевые каналы в открытом (MthK, Kv1.2, KvaP) (Jiang, Lee et al., 2002, Jiang, Lee et al., 2003, Long, Campbell et al., 2005) и закрытом (KcsA, KirBac1.1, KirBac3.1) (Doyle, Morais Cabral et al., 1998, Clarke, Caputo et al., 2010, Zubcevic, Bavro et al., 2014) состояниях, каналы глутаматного рецептора (NMDA и AmpA) (Sobolevsky, Rosconi et al., 2009, Lee, Lu et al., 2014), а недавно были получены структуры прокариотических натриевых каналов (NavMs, NavAb, NavRh) (Payandeh, Scheuer et al., 2011, Zhang, Ren et al., 2012, Bagneris, Decaen et al., 2013). Обилие доступных структур позволяет использовать теоретические методы компьютерного моделирования для изучения молекулярных детерминант и функциональных особенностей строения каналов .

2.1.3. Теоретические методы Экспериментальные данные, полученные из разных источников и разными методами, могут быть объединены в рамках модели. Это позволяет более корректно интерпретировать экспериментальные данные и дает возможность более тонкого планирования новых исследований, например, выбора участков для направленного мутагенеза. При этом отклонение получаемых экспериментальных результатов от предсказанных свидетельствует о расхождении между использованным «шаблоном» или моделью и реальным объектом исследования. Анализ этих отклонений позволяет корректировать модель и, используя новые экспериментальные подходы, получать более точные сведения об изучаемом ионном канале .

Основным подходом к изучению неизвестных канальных белков методами компьютерного моделирования является, безусловно, гомологическое моделирование (Rigden, 2009). Эмпирически установлено, что если последовательности двух белков идентичны друг другу более чем на 30%, то белки почти наверняка являются «родственниками» и степень эволюционной дивергенции ещё не столь велика, чтобы их структуры утратили общность. Эти наблюдения и являются основой методики предсказания пространственной структуры, называемой моделированием на основании гомологии .

На настоящий момент моделирование по гомологии позволяет установить структуру более половины белков, чьё строение ещё неизвестно. Если же выбирать мишени для экспериментального определения структуры таким образом, чтобы в результате для каждого белка был получен хотя бы один структурный гомолог (с идентичностью последовательностей свыше 30%), то окажется, что достаточно получить всего 16 000 структур (Vitkup, Melamud et al., 2001), а «степень покрытия» при этом составит более 90%, включая и мембранные белки. Моделирование по гомологии в этом случае поможет предсказать структуры большей части оставшихся белков .

Процесс моделирования по гомологии (Hillisch, Pineda et al., 2004, Ginalski, 2006) включает несколько шагов:

Идентификация структурного шаблона — белка с известной пространственной структурой, гомологичного моделируемому (идентичность последовательностей свыше 30%). Поиск производится с помощью серверов FASTA или PSI-BLAST (или их аналогов) в базе структур белков PDB (едином депозитарии структурных данных для биомакромолекул);

Выравнивание аминокислотных последовательностей шаблона и модели .

Парное выравнивание служит «инструкцией» программам, осуществляющим моделирование. Множественное выравнивание может быть полезно для выявления консервативных остатков во всём семействе или отдельных подсемействах белков .

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

Построение модели заключается, главным образом, в «натягивании»

последовательности моделируемого белка на «остов» шаблона согласно выравниванию. «Петлевые» участки (не имеющие гомологии с шаблоном) достраиваются независимо, положение боковых цепей оптимизируется с помощью методов эмпирических силовых полей .

Оценка качества, оптимизация и использование модели. Самый сложный этап моделирования по гомологии — оптимизировать модель с учётом всей доступной биологической информации по моделируемому белку. Вообще, моделирование структуры по гомологии с белком, выполняющим иную функцию, не способно автоматически дать модель, пригодную для практически важных задач. Обязательно требуется аккуратная оптимизация, превращающая «заготовку» (которой, по сути, является модель «нулевого приближения») в рабочий инструмент, — задача, зависящая скорее от интуиции и опыта исследователя, чем от конкретных компьютерных методик .

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

Рисунок 14. Схема моделирования по гомологии на примере рецептора мелатонина MT 1 человека (Chugunov, Farce et al., 2006) .

Одним из наиболее важных шагов в гомологическом моделировании является построение выравнивания. На данный момент разработано и реализовано множество алгоритмов, предназначенных для выравнивания первичных структур. Разделяют два основных подхода: парное выравнивание, при котором ищутся сходные участки для двух последовательностей, и множественное выравнивание, которое применяется для нахождения консервативных участков в наборе последовательностей .

Различают глобальное и локальное выравнивание. Глобальное выравнивание

предполагает, что последовательности подобны по всей длине .

В глобальное выравнивание включаются входные последовательности целиком. Локальное выравнивание применяется, если последовательности содержат как родственные (гомологичные), так и негомологичные участки. Результатом локального выравнивания является выбор участка в каждой из последовательностей и выравнивание между этими участками. Локальное выравнивание, в частности, широко применяется для многодоменных белков в том случае, если родственными предполагаются не все домены сравниваемых белков. Глобальное выравнивание применяется в случае исследования двух родственных белков, зачастую из одного семейства. В этом случае можно ожидать близкого сходства, и большинство изменений носят характер точечных замен в аминокислотной последовательности .

Наиболее известным и часто используемым инструментом для выравнивания последовательностей является программа Clustal (Goujon, McWilliam et al., 2010, Sievers, Wilm et al., 2011, McWilliam, Li et al., 2013) .

Для построения и оптимизации моделей применяются два основных метода:

молекулярная динамика и метод Монте-Карло .

2.1.4. Молекулярная динамика Метод классической молекулярной динамики (метод МД) — это метод, в котором временная эволюция системы взаимодействующих атомов или частиц отслеживается интегрированием их уравнений движения (Rigden, 2009). Метод молекулярной динамики является одним из наиболее важных инструментов теоретического изучения структуры и динамики биологических макромолекул, жидкостей, газов и других молекулярных систем .

Изначально разработанный в теоретической физике, метод молекулярной динамики получил большое распространение в науке о веществе и, начиная с 1970-х годов, в биохимии и биофизике. Метод молекулярной динамики играет важную роль в определении структуры белка и уточнении его свойств. Взаимодействие между объектами может быть описано силовым полем (классическая молекулярная динамика), квантово-химической моделью или смешанной теорией, содержащей элементы двух предыдущих .

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

Молекулярно-динамические расчёты позволяют получить информацию о флуктуациях и конформационных изменениях белков, нуклеиновых кислот и других молекул. В настоящее время метод молекулярной динамики повсеместно используется для исследования структурных, динамических и термодинамических характеристик биологических молекул и их комплексов. К задачам, решаемым методом молекулярной динамики, относится изучение стабильности белков, исследование конформационных превращений систем, сворачивания белков, докинга лигандов, изучение процессов ионного транспорта. Важную роль метод молекулярной динамики играет и при установке конечной структуры молекул и комплексов, определённой методами ЯМР и рентгеноструктурного анализа, а также в фармакологии при разработке новых лекарственных препаратов .

Основные положения метода:

Для описания движения атомов или частиц применяется классическая механика .

Уравнение движения частиц находят при помощи аналитической механики;

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

Точное знание траекторий движения частиц системы на больших промежутках времени не является необходимым для получения результатов макроскопического, (термодинамического) характера;

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

Методы молекулярной динамики играют одну из ключевых ролей в развитии теоретической биохимии. Исследования методами МД позволили получить огромное количество моделей белков, как ионных каналов, так и прочих. Однако и этот метод связан с определёнными ограничениями, что не позволяет использовать его для любых задач. Во многом это связано с тем, что расчёт требует слишком большого времени и компьютерных мощностей (Rigden, 2009) .

2.1.5. Метод Монте-Карло Наряду с методом молекулярной динамики используется ещё одна группа квазидинамических методов, объединённых под общим названием «метод Монте-Карло» .

Метод Монте-Карло отличается от метода молекулярной динамики тем, что каждая следующая конформация молекулы определяется не путём решения уравнений Ньютона, а как результат случайного процесса. Суть этого подхода состоит в сочетании случайных бросков конформации с последующим градиентным спуском в ближайший энергетический минимум.

При этом полная конформационная энергия системы (Eк) представляется в виде суммы компонент:

Eк = Енев + Еэл + Етор + Евал + Евс + Егидр

Енев – ван дер Ваальсово взаимодействие пар валентно не связанных атомов .

Евал – энергия деформации валентных углов .

Евс – энергия водородных связей .

Eэл – элетростатическое (Кулоновское) взаимодействие пар валентно не связанных атомов .

Етор – торсионная энергия в четверке валентно связанных атомов .

Егидр – энергия гидратации .

На первом этапе проводится минимизация энергии стартовой конформации (Е0) .

Затем одна из обобщённых координат системы изменяется случайным образом, и проводится новая минимизация (Е1). Новая конформация принимается в траекторию, если её энергия меньше предыдущей Е1Е0. В противном случае, вероятность принятия новой конформации равна Больцмановскому фактору от разности энергий сравниваемых конформаций E0 и E1. Если новая конформация принята, следующий шаг делается уже из неё. Температура обычно принимается 600 или 1000К. Таким образом, за счёт случайного изменения конформации система может перейти в новый энергетический минимум. Кроме того, возможность принятия конформера с повышением энергии обеспечивает более широкое исследование всего конформационного пространства. При МСМ поиске все низкоэнергетические конформации накапливаются для последующего анализа. Расчет прекращается после того, как определенное число минимизаций не приводит к понижению энергии системы (Тихонов, 2004) .

Таким образом, в настоящее время только комплексное исследование способно значительно углубить и развить существующие представления о строении и функционировании ионных каналов. Сочетание экспериментальных и теоретических подходов представляется особенно продуктивным. Простое накопление опытных данных без теоретического анализа в рамках молекулярных или кинетических моделей не сулит существенного прогресса. В то же время моделирование, не являющееся обобщением всего набора доступных данных, остаётся ограниченным и не способствует дальнейшему прогрессу в экспериментальных исследованиях. Если же экспериментальные и теоретические подходы используются параллельно, в рамках единого проекта, их взаимодействие происходит непрерывно, и модельнотеоретическая работа становится непосредственным продолжением классических экспериментов, а подчас и исходным пунктом для планирования новых экспериментов и подходов .

2.2. Параметры и методы расчётов, использованные в работе 2.2.1. Расчёт энергетически оптимальных конформеров методами молекулярной механики Во всех расчётах использовалось приближение атом-атомных потенциалов, при котором квантовые эффекты игнорировались. При этом подходе валентная структура системы принимается неизменной, то есть такие процессы, как протонирование/депротонирование, переходы между энантиомерами, разрыв и образование ковалентных связей не могут быть смоделированы. Однако для обратимого связывания лигандов, что и служило главной расчётной задачей, приближение атоматомных потенциалов является приемлемым компромиссом между точностью и скоростью .

Ван дер Ваальсово взаимодействие рассчитывалось по потенциалу Морзе (6/12) с использованием силовых полей AMBER (Weiner, Kollman et al., 1984) .

Поскольку расчёт энергии гидратации является времяёмким, гидратация игнорировалась. Для учета влияния растворителя, состоящего в экранировании электростатических взаимодействий, учитывалось введение зависимости диэлектрической постоянной от расстояния (Guenot and Kollman, 1992, Guenot and Kollman, 1993). Для того, чтобы учесть влияние среды, диэлектрическая проницаемость рассчитывалась по формуле = d · [0 + (1 - s) · (1 - 0 )], где d – расстояние между взаимодействующими атомами, а s – величина, принимающая значения от 0 до 1 в зависимости от того, если взаимодействующие атомы находятся в водной среде (0) или в белковой поре (1), 1 = 8, 0 =2 (Garden and Zhorov, 2010) .

Для ускорения расчёта, взаимодействия между атомами, удаленными друг от друга более чем на 8, не учитывались. Такая отсечка вносит заметные искажения при расчете дальнодействующих электростатических взаимодействий (Brooks, Brunger et al., 1985) .

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

Для пептидных цепей бралась стандартная валентная геометрия и заряды на атомах .

В процессе расчётов валентные углы пептидов удерживались жесткими (Zimmerman, Pottle et al., 1977) .

Поиск оптимальной геометрии проводился методом Монте-Карло с минимизацией энергии (МСМ) (Li and Scheraga, 1987). Для оценочных расчётов этот параметр обычно брался равным 200 – 500, а для сложных задач определения возможных конформаций моделей - от 5000 до 10000 .

При построении трёхмерного выравнивания между различными каналами, проводилась минимизация расстояния между углеродами основной цепи двух витков pспирали, образованными остатками начиная с p39 и заканчивая p47. При этом минимизация проводилась одновременно для всех субъединиц .

2.2.2. Гомологическое моделирование и системы ограничителей Использование рентгеновских структур белков для создания моделей белков родственных в настоящее время широко распространено. Самый простой способ – зафиксировать главные цепи модели в соответствии с рентгеновской структурой «шаблона». В этом случае свойства моделируемого белка будут отличаться от свойств шаблона только за счет локального влияния замещенных аминокислот. Очевидно, что такой подход продуктивен только для белков, демонстрирующих очень высокую гомологию с шаблоном. В противном случае, необходимо учитывать деформации главной цепи, вызванные замещением аминокислот .

Альтернативным подходом является задание начальной геометрии в соответствии с рентгеноструктурными данными и её дальнейшая оптимизация. При этом, как правило, возникают две главные проблемы. В ряде случаев моделируется только часть белка, конформация которой без внешних стабилизирующих факторов не является стабильной .

Очевидно, что свободная оптимизация такой незамкнутой системы приведет к её распаду .

Вторая проблема состоит в том, что в начальной геометрии модели замещённые аминокислоты часто формируют запрещённые (близкие) контакты. Естественно, такие контакты обладают очень высокой энергией, и при оптимизации система буквально «взрывается», то есть теряет гомологию с шаблоном. Из сказанного очевидно, что создание корректной гомологической модели требует применения специфических приемов .

В рамках данной работы при создании моделей использовались системы «ограничителей» (constraints). Ограничители представляют собой дополнительные потенциальные функции для удержания конформации системы в заданных пределах .

Потенциальные функции ограничителей имели вид парабол с плоским дном (flat-bottom parabolic functions). Энергия ограничителя равна 0, если конформация системы находится в заданных пределах (плоское дно). Если же конформация отклоняется от предельных значений, энергия ограничителя растет пропорционально квадрату отклонения. Наиболее часто использовались два типа ограничителей: ограничитель отклонения декатровых координат атома от заданных значений и ограничитель расстояния между атомами .

Ограничитель отклонения декартовых координат накладывался на С-альфа атомы аминокислот. В качестве центра разрешённой области брались рентгеновские координаты .

Радиус плоского дна ограничителя обычно составлял от 1 до 2. Ограничитель расстояния между атомами применялся для фиксации специфических взаимодействий, например, водородных связей или солевых мостиков .

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

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

Простое ограничение декартовых координат С-альфа атомов аминокислот вблизи рентгеновских значений не всегда продуктивно. Простейший пример – вставки и делеции в цепи, вблизи которых применять ограничители декартовых координат С-альфа атомов просто невозможно .

В некоторых случаях, как например при моделировании натриевых рецепторов по шаблону калиевого канала, изначально было ясно, что N-конец мембранной петли имеет иную конформацию, чем в шаблоне, поскольку свойства селективных фильтров каналов резко различаются. В этих случаях использовалась фиксация специфических взаимодействий. Например, ограничители расстояний, фиксирующие -спиральные водородные связи, не препятствуют изменению ориентации спиралей и их деформации, однако удерживают вторичную структуру. Ценным методом выявления специфических связей является анализ первичных структур гомологичных белков. Например, в структуре прокариотического натриевого канала NavAb существует контакт между Q p45 и Rp58, которому можно найти аналоги в модели Nav1.4, построенной на его основе .

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

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

Конформационная энергия комплекса блокатор-канал может быть подразделена на три основные составляющие: внутренняя энергия блокатора, внутренняя энергия канала и энергия взаимодействия. Раздельный анализ этих компонент часто оказывается полезен. В свою очередь, энергия взаимодействия может быть разложена по отдельным аминокислотным остаткам модели, то есть может быть определён вклад каждом аминокислоты в энергию взаимодействия. Такой анализ необходим в тех случаях, когда блокатор взаимодействует одновременно с большим числом аминокислот, что очень актуально в случае связывания конотоксинов .

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

Так как основным объектом работы являлись µ-конотоксины, которые полностью или частично блокируют ток ионов через канал, то необходимо было разработать подход, который бы позволил оценить величину ионного тока через канал, связанный с конотоксином .

Для этого был предложен следующий подход: создавалось 10 000 стартовых позиций, в которых ион случайным образом размещался в области селективного фильтра и наружной поры, после чего проводилась короткая минимизация энергии (10 минимизаций для каждой из 10 000 структур). Для анализа выбирались только те структуры, где энергия взаимодействия иона с остальной системой была отрицательной (то есть превалировали силы притяжения). Такой подход позволял найти области, которые могут быть заполнены ионами. В качестве количественной характеристики тока использовалось число ионов на определённой точки оси координат .

2.3. Программное обеспечение Сложные протоколы гомологического моделирования и анализа связывания блокаторов предъявляют высокие требования к программному обеспечению расчетов .

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

2.3.1. ZMM Такой пакет расчётных программ (ZMM) был создан Б.С.Жоровым. Он поддерживается и постоянно усовершенствуется (Zhorov, 1981, Zhorov, 1982, Zhorov, 1993, Zhorov and Ananthanarayanan, 1994, Zhorov and Ananthanarayanan, 1995, Zhorov and Ananthanarayanan, 1996, Zhorov and Ananthanarayanan, 1998, Zhorov and Ananthanaraya na n, 2000, Zhorov and Lin, 2000). В ZMM реализованы алгоритмы Монте-Карло с минимизацией энергии. Пакет поддерживает функции, необходимые для гомологического моделирования, в том числе системы ограничителей. Версия 2015 г. включает 11 отдельных программ, позволяющих задать практически любые условия расчёта (свыше 100 численных параметров расчёта!) .

Одно из основных преимуществ этого пакета – описание молекулярных систем и минимизация энергии во внутренних координатах.

Внутренними координатами являются:

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

2.3.2. Прочие При анализе сложных белковых структур, эффективные и удобные средства визуализации становятся важным компонентом расчётных систем. В рамках данной работы использовалась программа молекулярной графики (MVM) для персональных компьютеров .

Программа работает в операционной системе Windows. Для создания изображений используется библиотека OpenGl. Программа рассчитана на визуализацию молекулярных систем размером до 256 тысяч атомов. Кроме набора разнообразных средств управления и визуализации, MVM обеспечивает пользователю интерактивную среду для управления расчетным пакетом ZMM. При этом возможно управление расчетными программами и вывод информации в ходе работы. В то же время графическая (MVM) и расчетная (ZMM) части являются независимыми и могут использоваться раздельно .

Для проведения статистического анализа и построения корреляционных зависимостей использовалась свободно распространяемая программная среда R, и её графическая оболочка RStudio .

Также, в рамках изучение специфических контактов и молекулярных детерминантов фолдинга каналов, строились графы контактов внутри канала. При этом использовалась свободно распространяемая программа yEd graph editor .

Глава 3. Результаты и их обсуждение .

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

Поэтому дальнейшую работу можно условно разбить на четыре взаимосвязанные части:

изучение проблем гомологического моделирования, изучение механизмов блокады, создание модели канала и, непосредственно, изучение действия µ-конотоксинов при помощи модели .

3.1. Проблемы гомологического моделирования Как было сказано ранее, основным методом, позволяющим исследовать ионные каналы, структура которых ещё не известна, является гомологическое моделирование .

Одним из ключевых моментов в гомологическом моделировании является именно выравнивание аминокислотных последовательностей, позволяющее найти совпадающие по положению остатки в структуре-шаблоне для остатков моделируемого белка. Очевидно, что неверное выравнивание делает всю гомологическую модель неадекв атной. Например, аминокислотный остаток, который смотрит в пору канала, в результате неверного выравнивания окажется сориентирован в другом направлении. Хотя в гомологическом моделировании используются методы оптимизации пространственных структур, результаты неверного выравнивания как правило не могут быть устранены, особенно в области регулярных вторичных структур .

У современных методов выравнивания есть существенные недостатки. Например, макромолекулы, даже родственные, но находящиеся на достаточном удалении друг от друга в эволюционном дереве, зачастую очень сильно различаются по своим последовательностям. Алгоритмы же, которые используются в программах для выравнивания последовательностей, ищут определённые сходства, ориентируясь на консервативность аминокислот. Таким образом, построенное выравнивание носит вероятностный характер и не застраховано от ошибок .

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

3.1.1. Выбор выравнивания порообразующих спиралей Р-loop каналов Для начала построим выравнивание между последовательностями канала TRPV1 и трёх калиевых каналов: Kv1.2, KcsA и MlotiK при помощи средств линейной (одномерной) биоиноформатики (в данном случае это сервер Clustal Omega), и сравним это линейное выравнивание (1D выравнивание) с результатами 3D выравнивания известных пространственных структур этих белков (рисунок 15) (здесь и далее используется универсальная нумерация, предложенная (Zhorov and Tikhonov 2004), см. п. 1.2.1) .

Рисунок 15. Выравнивание структур p-loop каналов. Сверху – линейное выравнивание. Снизу – расстояние между атомами углеродов основной цепи разных каналов относительно MlotiK .

Рисунок 16. Трёхмерное наложение p-спиралей p-loop каналов. Сверху вниз: наложение структур Kv1.2 и MlotiK; наложение структур TRPV1 и MlotiK, исходя из выравнивания Clustal Omega (A); наложение структур TRPV1 и MlotiK, с модифицированным выравниванием (B) .

Подобие последовательностей между TRPV1 и калиевыми каналами весьма низкое (порядка 21%). В таких случаях белки часто имеют совершенно разные пространственные структуры .

3D выравнивание структур каналов Kv1.2, KcsA, MlotiK и TRPV1 (PDB коды 2a79, 1bl8, 3beh и 3j5p соответственно) было произведено по процедуре минимизации RMSD для p-спиралей, описанной в разделе «Методы» (стр 52). Результаты представлены на рис. 16Видно, что в случае, к примеру, TRPV1 и MlotiK пространственная укладка сохраняется .

Что же касается калиевых каналов, то они достаточно консервативны как по последовательности, так и по пространственной укладке .

Рисунок 17. Трёхмерное наложение структур каналов TRPV1 и MlotiK. Структура MlotiK показана оранжевым, структура TRPV1 – зелёным .

Обратим более пристальное внимание на построенное нами одномерное выравнивание. Clustal-Omega правильно предсказывает выравнивание в области S5спиралей всех исследуемых каналов, ориентируясь на положение заряженного остатка аргинина в позиции о5 и наличие характерных ароматических остатков в позициях о14 и о19 (рисунок 15) .

В Р-спирали, используя мотив TХGХGD, свойственный калиевым каналам, ClustalOmega выравнивает области селективного фильтра. Как видно из рисунка 16, в случае калиевых каналов трёхмерное наложение полностью подтверждает этот вывод: к примеру, триптофаны в позиции p42 явно находятся в одной и той же позиции как в MlotiK, так и в Kv1.2 (то же самое верно и для KcsA). Но в то же время, в связи с низкой консервативностью TRPV1 относительно калиевых каналов, выравнивание Р-спирали неднозначно. Для количественной оценки различий в трёхмерном выравнивании была рассчитана величина среднеквадратического отклонения (СО) для -углеродов двух центральных витков Рспирали (от остатка p39 до p47). Если использовать выравнивание, основанное на совпадении сигнатуры селективного фильтра (GXG), то мы получим большое значение СО, равное 2.27 (TRPV/A). В то же время между калиевыми каналами мы получаем величину СО равную 0.58 для пары MlotiK/Kv1.2 и 0.91 для пары MlotiK/KcsA. Из рассмотрения трёхмерных структур видно, что остатки MlotiK_Lp45 и TRPV_K p45, занимающие одно и то же положение в выравнивании по Clustal-Omega, имеют совершенно разную пространственную ориентацию (рисунок 16). Визуальное сравнение трехмерных структур позволило предложить альтернативное выравнивание последовательностей, представленное на рис 15 (TRPV/B). В этом выравнивании мотив GMG сдвинут относительно GYG на две позиции. Расчет 3D выравнивания, проведенный на основе этого предположения показал, что СО падает более чем в полтора раза, и становится равным 1.32. Хотя это по-прежнему выше, чем соответствующие значения для калиевых каналов, но это неудивительно, так как консервативность последовательностей между каналами достаточно низкая. Более того, 3D выравнивание ставит в одно положение остатки MlotiK_Lp42 и TRPV_Ep42, ориентации которых совпадают (рисунок 16) .

Таким образом, сопоставление структур атомного разрешения, позволило однозначно установить роль и пространственное расположение каждого конкретного аминокислотного остатка в цепи .

3.1.2. Анализ стабилизаторов в структурах каналов Построение трёхмерного выравнивания позволило нам увидеть сходства и различия в структурах изучаемых ионных каналов, но само по себе оно не даёт ответа на вопрос, почему эти отличия возникают. Для этого необходимо было дополнительно изучить внутренние взаимодействия в канале, а именно: внутрисегментные и межсегментные контакты между остатками аминокислот .

Анализ контактов внутри каналов проводился по следующей методике: контактом считались все взаимодействия, полученные в модели, при расстоянии между взаимодействующими группами не более 5 и энергией взаимодействия менее -0.3 ккал/моль (что соответствует более сильному притяжению). Такой порог был выбран эмпирически, путём рассмотрения контактов в трёхмерной структуре и сравнения их с общим списком взаимодействий. Следует дополнительно отметить, что энергия, рассчитываемая в модели, не является свободной энергией, а служит оценкой конформационной энергии системы .

Были рассмотрены следующие четыре канала: TRPV1, MlotiK, Kv1.2 и KcsA. Эти каналы имеют схожую структуру, свойственную калиевым каналам, что позволяет изучить их на предмет сходств и различий в контактах между аминокислотными остатками .

Структуры этих ионных каналов были предварительно оптимизированы согласно методике минимизации, описанной в разделе «Методы» (стр 51). На рисунке 18 показаны рассчитанные энергии внутрисегментных и межсегментных контактов и количество значимых контактов для каждого остатка в p-спирали. Верхняя половина графика соответствует суммарному количеству контактов, образованных соответствующим остатком, а нижняя часть показывает суммарную энергию его контактов, взятую по модулю. В калиевых каналах консервативные триптофаны в позициях p41 и p42 образуют больше контактов, чем любой другой остаток, и суммарная энергия их контактов значительно ниже (более предпочтительна), чем для остальных остатков .

Рисунок 18. Количество и энергия контактов, образуемых остатками в p-спирали. Верхняя часть графика показывает суммарное количество контактов, образовываемых соответствующим остатком, а нижняя – суммарную энергию этих контактов .

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

Боковая цепь триптофана Wp41 вытянута в сторону соседней субъединицы и взаимодействует с N-концом S6-спирали, особенно с остатками в позициях i4 и i8, а также образует плотный контакт с консервативным аспартатом в позиции p54. Эти контакты стабилизируют пространственную укладку селективного фильтра канала. Боков ая цепь второго триптофана (в позиции p42) вытянута в сторону полости между S5 и S6 спиралями субъединицы, в которой находится триптофан, и образует несколько контактов. Хотя в разных каналах Wp42 образует контакты с разными остатками в S6-спирали (это могут быть i6, i7 или i10), тем не менее, общая ориентация триптофана и его контакты в целом консервативны среди каналов. Такие же контакты можно найти и в других калиевых каналах, таких как KirBac1.1, KvAP и MthK .

В канале KcsA есть ещё один остаток, играющий важную роль в стабилизации структуры: глутамат в позиции p45, боковая цепь которого согласно рентгеноструктурным данным (Doyle, Morais Cabral et al., 1998) близка к аспартату в позиции p54. Если оба остатка ионизированы, то электростатическое отталкивание не позволило бы им образовать контакт. Для решения этого парадокса (Berneche and Roux, 2002) было предположено, что такая пара имеет один общий протон, что позволяет им образовывать плотный контакт. В нашей модели Ep45 находился в протонированном состоянии, что позволило получить отрицательную энергию взаимодействия, что соответствует притяжению .

В канале TRPV1 отсутствуют триптофаны в позициях p41 и p42. Вместо них в последовательности находятся лейцин и глутамат, которые, согласно расчётам, не формируют сильных контактов. Тем не менее, остатки F p44, K p45 и Fp46 образуют много плотных контактов. Фенилаланин в позиции p44 взаимодействует с S6 спиралью той же субъединицы и дотягивается до боковой цепи тирозина o14 в центре S5 спирали соседней субъединицы. Фенилаланин в позиции p46 взаимодействует с S6 спиралью той же субъединицы, в частности, с тирозином o10. Стоит отметить, что o14 и o10 формируют контакт друг с другом. Также интересен контакт, возникающий между аспартатом в позиции p54 и лизином в позиции p45 (который соответствует Ep45 в KcsA). По-видимому, между остатками образуется солевой мостик, что даёт сильный контакт .

Таким образом, калиевые каналы и TRPV1 образуют цепь контактов, формирующих циклическую структуру, состоящую из контактов между остатками Р-петли и S6 спиралей различных субъединиц (рисунок 20). Если смотреть с внеклеточной стороны, то общий вид стабилизаторов похож в калиевых каналах и TRPV1. И там, и там пара ароматических остатков образует кольцо с участием остатков S6 спирали .

Рисунок 19. Схематическое и трёхмерное изображение наиболее важных контактов в pспирали. Цвет вершины графа соответсвует цвету субъединицы (жёлтый – I, зелёный – II, лиловый – III, голубой – IV), цифры рядом со стрелочками показывают энергию взаимодействия между вершинами .

Рисунок 20. Схематическое и трёхмерное изображение колец, образуемых контактами остатков p-спирали. Цвет вершины графа соответсвует цвету субъединицы (жёлтый – I, зелёный

– II, лиловый – III, голубой – IV) .

3.1.3. Выбор выравнивания S5 и S6 спиралей в гомологических моделях Все приведённые выше данные относятся преимущественно к Р-петлям и Рспиралям. Однако, не менее важную роль в формировании поры играют внешняя и внутренняя спирали канала: и здесь тоже существуют свои сложности в построении выравнивания. Ранее разными авторами было предложено несколько различных вариантов .

Липкинд и Фоззард (Lipkind and Fozzard, 2000) в своей модели натриевого канала Na(V)1.4, построенной на шаблоне бактериального калиевого канала KcsA предложили выравнивание порообразующих спиралей (S6 в Na(V)1.4 и M2 в KcsA), исходя из следующих соображений. В KcsA ключевую роль в образовании спирали M2 играет триптофан в позиции i2, поэтому были выбраны наиболее подходящие крупные гидрофобные остатки в каждой из субъединиц Na(V)1.4, чтобы выравнивать последовательности по ним. Это позволило боковым цепям фенилаланина и тирозина в сегменте S6 четвертого домена (F4i15 and Y4i22 ), участвующим в связывании Na(V)1.4 с местными анестетиками (Chahine, Chen et al., 1992, Fozzard, Sheets et al., 2011), занять позиции, соответствующие изолейцину-i15 и треонину-i22, которые обращены в пору KcsA, согласно ретгеновской структуре. Аспарагин-1i19 и лейцин-1i22 в первом домене также оказались обращенными в пору, что согласуется с результатами мутационного анализа воздействия местных анестетиков (рисунок 21). Выравнивание между собой S5 различных субъединиц основывалось на наличии консервативных фенилаланинов на Cконце спирали (и тирозина в четвёртой субъединице) .

Хубер и соавторы (Huber, Wappl et al., 2000) для выравнивания последовательностей кальциевого канала L-типа и калиевого канала KcsA использовали процедуру «поиска посредника». Авторы использовали парные выравнивания между кальциевым каналом Cav1.2 и натриевым каналом Na(V)1.4, между KcsA и Shaker калиевым каналом (Schrempf, Schmidt et al., 1995), и между Na(V)1.4 и Shaker калиевым каналом (Guy and Conti, 1990) .

Однако ни одно из предложенных выравниваний не было достаточно убедительным, чтобы принять его как единственно верное. Поэтому в своей работе Жоров и соавторы (Zhorov, Folkman et al., 2001) предложили проверить эти варианты, используя их для построения модели кальциевого канала L-типа, и анализируя данные по взаимодействию каналов с дигидропиридинами. Построив различные гомологичные модели кальциевого канала, авторы показали, что лучше всего объяснить экспериментальные данные можно только в определённой модели (рисунок 21), где остатки, отвечающие за связывание лиганда будут находиться в поре канала, а S5 спирали – согласно выравниванию из работы (Huber, Wappl et al., 2000). Последовательности S6 спиралей выравниваются аналогично с выравниванием, предложенным Липкиндом и Фоззардом, однако с изменениями в сегменте (IS6), которое в 3D моделях приводит к противоположным ориентациям ключевого остатка, аспарагина N1i20, который не обращён в пору канала и участвует в межсегментных взаимодействиях, чем и объясняется его влияние на свойства каналов и связывание лигандов (Tikhonov, Bruhova et al., 2015). При таком выравнивании все высококонсервативные остатки аспарагина в S6 сегментах натриевых и кальциевых каналов занимают одно и то же положение (рисунок 21) .

–  –  –

Появление структуры бактериального натриевого канала NavAb (Payandeh, Scheuer et al., 2011) дает ценную информацию для решения проблемы выравнивания. Отметим два важных аспекта структуры NavAb. Во-первых, поскольку тип укладки NavAb такой же, как и в калиевых каналах, можно построить однозначное выравнивание последовательностей, основанное не на методах линейной биоинформатики, а на расчете выравнивания трёхмерных структур, по процедуре, описанной в методах (стр. 52). Результаты расчета, представленные на рис 22. позволяют получить однозначное выравнивание с калиевыми каналами, представленное на рис 21 .

Рисунок 22. Трёхмерное наложение структур каналов NavAb и MlotiK. Структура NavAb показана оранжевым, структура MlotiK – зелёным .

С другой стороны, NavAb, как близкий родственник эукариотических натриевыхканалов, однозначно выравнивается по последовательностям с натриевыми и кальциевыми каналами (для выравнивания использовались алгоритмы Clustal Omega) .

NavAb и Nav1.4 гораздо легче выравнивать относительно друг друга, за счёт наличия высококонсревативных остатков в разных сегментах. Так, на N-конце S5 спирали находится мотив из аспарагина, лейцина и фенилаланина, полностью повторяющийся в двух доменах Nav1.4 и в NavAb, но отсутствующий в калиевых каналах, что позволяет однозначно построить выравнивание между S5 спиралями. Кроме того, такое выравнивание подтверждается консервативными ароматическими остатками в середине спирали S5. В S6 спиралях также есть несколько ключевых позиций: консервативные аспарагины в середине спирали, консервативный набор фенилаланинов, которые, вероятно, являются стабилизаторами структуры натриевых каналов и поэтому у них должен быть аналог в родственных каналах. Разумеется, тот факт, что Nav1.4 является гетеротетрамером, затрудняет изучение структуры, но всё же сходство между последовательностями достаточно высоко, чтобы с высокой степенью вероятности говорить о возможности построения гомологической модели по такому выравниванию .

Таким образом проведенный анализ одномерного и трехмерного выравнивания NavAB с известными последовательностями и трехмерными структурами показал, что этот канал может служить своеобразным посредником для построения выравнивания между калиевыми каналами и натриевыми каналами эукариот. Выравнивание S5 и S6 спиралей, полученное при использовании структуры NavAB совпало с выравниванием, предложенным в работе (Tikhonov and Zhorov, 2012). Именно это выравнивание было использовано в дальнейшем при построении модели .

3.1.4. Заключение Выравнивание последовательностей является ключевым шагом в построении гомологических моделей. Поскольку традиционные методы выравнивания, опирающиеся исключительно на анализ аминокислотных последовательностей, не несут информации о функциональных свойствах белков, то включение в рассмотрение трёхмерных структур и других данных является крайне важным .

Проведенный анализ и расчеты позволили построить однозначное выравнивание междупоследовательностями и пространственными структурами p-loop каналов, включая натриевые каналы эукариот, моделирование который является основной задачей работы .

3.2. Стерические и электростатические факторы блокады ионных каналов Гомологическое моделирование представляет собой лишь первый этап работы по изучению механизмов лиганд-рецепторных взаимодействий. Не менее сложной частью является моделирование посадки лигандов в канале. И здесь, в свою очередь, существует большое количество факторов, которые необходимо принимать во внимание при проведении комплексного исследования .

Взаимодействие между ионами, блокирующими канал и ионами, обеспечивающими ток через канал, является одной из фундаментальных проблем в области фармакологии ионных каналов. Блокирующая частица создает пространственный/гидрофобный/электростатический барьер для нормального проникновения ионов (Catterall, 2002). Однако вклад этих компонент в блокирующий эффект остается предметом научной дискуссии (см. стр. 26) .

Недавно была опубликована структура NavMs канала, связанного с препаратом, родственным локальным анестетикам: 2-(4-бромфенил)-1-(5-(4-хлорфенил)-1Н-имидазолил)этанамин) (Bagneris, DeCaen et al., 2014). Несмотря на то, что в структуре связанного лиганда виден лишь один тяжелый атом лиганды (бром), были обнаружены две очень важные особенности. Во-первых, атом брома, являющийся терминальным атомом ароматической группы использованного лиганда, находится в интерфейсе субъединицы .

Во-вторых, ионная заселённость участка вблизи центральной полости существенно снижается в присутствии блокирующего препарата. На основании этих результатов авторы предложили модель, в которой ароматический фрагмент, содержащий атом брома, находится в интерфейсе субъединицы, в то время как аминогруппа выступает в пору и противоположный конец препарата располагается в селективном фильтре, тем самым стерически смещая находящийся там ион .

Некоторые ключевые особенности этой модели совпадают с предыдущими предложениями, высказанными в работе Бруховой и соавторов, которые основаны на модельной интерпретации косвенных данных (Bruhova, Tikhonov et al., 2008). На рисунке 23 показано, что единственное существенное отличие состоит в том, что используемые в моделях молекулы местных анестетиков (такие как тетракаин и лидокаин) не имеют длинных групп на амино-конце, поэтому они не в состоянии войти в селективный фильтр и стерически вытеснить из него ионы. Поэтому, чтобы объяснить экспериментально наблюдаемые взаимодействия между местными анестетиками и ионами натрия, в этой модели предполагается, что природа этих взаимодействий в первую очередь электростатична. Таким образом, модель, предложенная в работе (Bagneris, DeCaen et al.,

2014) может быть верна для блокады канала NavMS лигандом Pl1, но, по-видимому, не может объяснить, взаимодействие между локальными анестетиками и ионами в эукариотических натриевых каналах .

Рисунок 23. Посадка локальных анестетиков в канале. Сверху: наложение модели из работы (Bruhova, Tikhonov et al. 2008) со структурой, полученной в работе (Bagneris, DeCaen et al. 2014);

видно, что хвост молекулы блокатора находится в интерфейсе между субъединицами канала;

снизу: химические структуры локальных анестетиков .

В связи с этим возникает вопрос: возможно ли предложить альтернативную интерпретацию данных, полученных исследователями в этой работе? Как будет показано далее, экспериментальные данные зачастую можно объяснить несколькими моделями, каждая из которых имеет свои плюсы и минусы .

3.2.1. Изучение возможных способов связывания лиганда Первым шагом, необходимым для построения модели связывания лиганда с каналом всегда является создание модели самого канала. Для этого использовалась структура высокого разрешения, которая была дополнительно оптимизирована для получения лучшей оценки энергетики системы. Оптимизация проводилась в два этапа: на первом этапе проводилась минимизации с ограничениями на альфа-углероды структуры, чтобы они не отклонялись от оригинальной структуры более чем на 1. На втором этапе все ограничения были сняты и модель канала была оптимизирована до тех пор, пока 5000 шагов минимизации не уменьшили энергию кажущегося глобального минимума .

С точки зрения моделирования, наиболее интересным из представленных в работе (Bagneris, DeCaen et al., 2014) данных является полученная структура с известной локализацией атома брома, входящего в состав лиганда – блокатора канала. Действительно, размеры интерфейса и размеры молекулы лиганда практически однозначно определяют положение ароматического фрагмента последнего. В то же время противоположный конец молекулы, который находится в достаточно широкой полости канала, обладает значительной конформационной свободой. Для того, чтобы изучить возможные ориентации лиганда, использовался следующий подход: лиганд помещался по центру вестибюля канала, после чего проводились случайные изменения его ориентации. Всего было создано 100 000 различно ориентированных стартовых позиций, каждая из которых была дополнительно оптимизирована. На рисунке 24 показан ансамбль низкоэнергетических структур, полученных в результате минимизаций .

Рисунок 24. Ансамбль низкоэнергетических структур, полученных для посадки Pl1. Жёлтым показан атом хлора в составе Pl1, оранжевым – атом брома, а синим – аминогруппа. Структура канала окрашена в цвета различных участков канала: p-спираль (лиловый), s5 (зелёный) и s6 (голубой) .

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

Аминогруппа лиганда в большинстве случаев подходит к свободному кислороду Tp48, что согласуется с предположением, высказанным в работе (Bagneris, DeCaen et al., 2014). Если же рассматривать положение дальнего от брома конца лиганда, то в целом можно выделить два основных варианта его расположения (рисунок 25). В первом случае лиганд проникает в область селективного фильтра и смотрит в сторону внеклеточного пространства, что, опять же, согласуется с предложенным авторами вариантом посадки .

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

Таблица 1. Энергия взаимодействия лиганда с остатками канала в двух конформациях .

–  –  –

Таким образом, в одной возможной конформации лиганд стерически вытесняет ионы из селективного фильтра, что объясняет полученные экспериментальные данные, а в другой стерическое замещение ионов в фильтре невозможно. Но можно ли предложить другое объяснение для полученных данных, исходя из модели горизонтальной посадки лиганда? Чтобы ответить на этот вопрос изучим поведение ионов в канале, в особенности в области селективного фильтра .

Рисунок 25. Два варианта расположения лиганда Pl1 в канале. Слева: лиганд проникает в область селективного фильтра и стерически препятствует току ионов через канал; справа: лиганд располагается между различными субъединицами таким образом, что заряженная группа оказывается на оси канала, оказывая электростатическое воздействие на ионы .

3.2.2. Распределение ионов в области селективного фильтра Точное число ионов натрия, которые могут одновременно находиться в поре канала, неизвестно. Поэтому для поиска оптимальных вариантов распределения ионов в канале было решено провести поэтапное изучение распределения в случае с разным количеством ионов в системе. В первом приближении в систему был добавлен только один ион. Для оценки его вероятного местоположения было создано 40 стартовых позиций иона равномерно по всей оси канала в области селективного фильтра и проведена МК минимизации его положения. Если посмотреть на полученное распределение возможных положений иона вдоль оси канала, то видно, что наблюдается одна предпочтительная позиция для иона – в области глутаматов селективного фильтра (рисунок 26). Эта позиция хорошо соответствует данным о карте электронной плотности в селективном фильтре (Bagneris, DeCaen et al., 2014) .

Зная расположение одного иона в канале, можно усложнить задачу. Добавим второй ион в систему, распределив его стартовые положения так же, как и для одного иона, вдоль оси канала в области селективности фильтра. При этом первый ион находится в системе в оптимальном положении. Но при этом мы минимизировали положение не только «нового»

иона, но и первого, чтобы учесть ион-ионные взаимодействия, которые несомненно играют важную роль в распределении ионов в канале .

Расчёты показали, что первый ион по-прежнему занимает свою позицию возле EP50 (будем называть такую позицию иона сайтом I), в то время как второй ион находится глубже в канале (рисунок 26), рядом с входом в центральную полость (назовём это сайт II). Опять же, существование такого сайта хорошо согласуется с данными по электронной плотности ионов в рентгеновской структуре. На следующем этапе был добавлен третий ион и расчеты были повторены, как описано выше (рисунок 26). Увеличение ион-ионного отталкивания привело к сглаживанию пиков оптимальных позиций I и II и к появлению третьего населенного сайта с внеклеточной стороны селективного фильтра (III сайта). При добавлении еще одного, четвёртого иона, ситуация существенно не изменилась, за исключением появления пика распределения в глубине центральной полости и увеличение популяции на I и III сайтах. Однако ионы, находящиеся вдалеке от селективного фильтра, нас не интересуют, поскольку лиганд воздействует на них слишком слабо .

Таким образом, расчёты показали, что распределения ионов в селективном фильтре существенно не различаются для систем, содержащих три и четыре иона, достаточно остановиться на трёх ионах для анализа действия лиганда .

Рисунок 26. Распределение ионов в канале NavMs. Сверху: распределение ионов вдоль оси канала при наличии одного, двух и трёх ионов в канале. Увеличение количества ионов в модели приводит к незначительному расширению пиков и появлению дополнительных пиков вне канала (на дальних участках оси). Снизу: усреднённая энергия ионов в системе в отсутствии (красным) и при наличии (чёрным) лиганда. Зелёным показана разность энергий .

3.2.3. Влияние Pl1 на энергию связи ионов Для того, чтобы оценить вклад электростатического воздействия, оказываемого лигандом на ионы, необходимо рассмотреть две различные системы: в присутствии лиганда и без него. Энергетика для каждой конформации была пересчитана при наличии в системе лиганда PL1. В качестве модели посадки лиганда использовалась полученная ранее модель с горизонтальной посадкой, оптимизированная для получения наилучшей моды связывания с каналом. На рисунке 26 показаны графики распределения средней энергии ионов для различных положений вдоль оси канала в отсутствии лиганда, в присутствии лиганда и разность между энергиями, что позволяет оценить вклад Pl1 в энергетику ионов в разных позициях. Для ионов, расположенных вблизи кольца Ep5, электростатический эффект присутствия лиганда был незначителен. Но для ионов в сайте II (около полости) расстояние до центра аминогруппы Pl1 составляла порядка 3 – 4 .

Соответственно, электростатическое отталкивание увеличивается до 1,3 – 1,5 ккал / моль .

Абсолютное значение этого отталкивания зависит от диэлектрической функции и ее параметров. Поэтому необходимо сравнить это отталкивание с другими электростатическими взаимодействиями. В таблице 2 показана энергия взаимодействия ионов с различными остатками канала в модели. Основной вклад в удержание ионов в канале вносят Тp48, Q p45 и Ep50. Из представленных данных видно, что увеличение электростатического отталкивания, путём включения PL1 в модель, может значительно изменить общую энергию иона в положении II. В целом суммарная энергия электростатического взаимодействия изменяется от -4,8 до -3,4 (изменение энергии на 31%) .

–  –  –

Следует отметить, что результаты находятся в соответствии с многочисленными наблюдениями. Например, было показано, что наличие ионов К+ внутри и вне клетки приводит к снижению блокады тетраэтиламмонием (Spassova and Lu, 1999, Thompson and Begenisich, 2003). Стоит также отметить, что эффект ионов четвертичного аммония отличается для разных калиевых каналов (например, Kir, Shaker и KcsA), что может быть объяснено различным сродством сайтов в каналах к ионам калия, что приводит к различиям в заселённости сайтов в течение цикла проводимости (Kutluay, Roux et al., 2005). Позднее было показано, что связывание тетрабутиламмония в канале KcsA зависит от заселённости селективного фильтра (1/3 или 2/4 сайтов) ионами калия (Faraldo-Gomez, Kutluay et al., 2007). Полученные же данные по NavMs и Pl1 предлагают аналогичное объяснение антагонизма между ЛА и ионами натрия .

3.2.4. Заключение Подытоживая, проведённое исследование методами молекулярного моделирования позволило продемонстрировать, что связывание лиганда PL1 в канале NavMs может повлиять на заселённость селективного фильтра ионами не только стерически (как предлагается авторами работы (Bagneris, DeCaen et al., 2014)), но также электростатически (когда ни одна из частей лиганда не входит непосредственно в селективный фильтр) .

Это не означает, что модель (Bagneris, DeCaen et al., 2014) является ошибочной. Как видно из приведённых рассуждений, для конкретного препарата и конкретного канала стерическое смещение действительно возможно. Тем не менее, даже в данном случае мы видим, что возможно предложить альтернативное объяснение экспериментальным данным .

Поэтому вряд ли стоит предполагать, что стерический механизм воздействия лиганда является универсальным для молекул локальных анестетиков. В частности, подобный механизм не может быть реализован для молекул местных анестетиков в эукариотических натриевых каналах, потому что в этих молекулах обычно отсутствуют длинные участки, идущие от аминогруппы. Поэтому подобные лиганды просто слишком малы для проникновения в селективный фильтр и для них возможен только электростатический механизм воздействия. Как видно из примера Pl1 и NavMs, электростатический вклад лиганда может быть достаточно значительным, по сравнению с другими электростатическими взаимодействиями в канале .

Полученные результаты хорошо согласуются с данными, ранее полученными для гомологических моделей эукариотических натриевых каналов. Однако точность гомологических моделей всегда ограничена и выводы, сделанные из подобных исследований, всегда носят более или менее гипотетический характер. В данном случае использовалась реальная структура атомарного разрешения, в которой существуют чёткие ограничения, накладываемые на положение лиганда, и, что особенно важно, на положение его аминогруппы, которая играет ключевую роль в электростатическом отталкивании .

Именно использование подобных экспериментальных данных позволяет получать результаты, которые более надежны, нежели чисто гомологическая модель .

3.3. Блокада натриевых каналов µ-конотоксинами Все предыдущие расчёты и полученные результаты подводят нас к наиболее важной части работы: изучению взаимодействия µ-конотоксинов с эукариотическим натриевым каналом Nav1.4. Для этого необходимо построить гомологическую модель канала, провести моделирование связывания токсинов и сравнить предсказания модели с экспериментальными данными .

3.3.1. Выбор модели В первую очередь, необходимо обратить внимание на некоторые особенности выбора структуры для моделирования. Доступные рентгеновские структуры прокариотических натриевых каналов (NavAb, NavMs, NavRh) имеют одинаковое строение в области поры и селективного фильтра (рисунок 27). В частности, как видно из рисунка, при наложении рентгеновских структур, P1 и P2 спирали разных каналов оказываются очень близки .

Рисунок 27. Наложение структур прокариотических натриевых каналов NavAb, NavMs и NavRh. NavAb показан зелёным, NavMs показан лиловым и NavRh показан красным .

Дополнительно выделены остатки, формирующие селективный фильтр канала .

Так как при построении модели проводится минимизация энергии, то выбор конкретного шаблона не является существенным для моделирования: расхождения в структурах меньше, чем допустимые отклонения модели от шаблона. Следующий важный этап работы – построение выравнивания между каналами. Выравнивание для спиральных сегментов рассматривалось ранее (см. стр. 68). Однако для моделирования связывания лигандов в области селективного фильтра канала критически важно правильно построить выравнивание полседовательностей именно в этой области .

Последовательности прокариотических и эукариотических каналов достаточно сильно различаются. В области P-петли идентичность между последовательностью NavAb и первой субъединицей Nav1.4 составляет ~ 21%. В работе Тихонова и Жорова было показано, что прямое выравнивание последовательностей (без вставок-удалений, основанное на положении остатков селективного фильтра) не позволяет построить гомологическую модель Nav1.4, способную объяснить экспериментальные данные по связыванию TTX с каналом (Tikhonov and Zhorov, 2012). Для решения этой проблемы был предложен скорректированный вариант выравнивания с модификациями в области селективного фильтра (рисунок 28) .

Рисунок 28. Скорректированный вариант выравнивания последовательностей субъединиц Nav1.4 с NavAb. Так как выравнивание последовательностей по остаткам, образующим селективный фильтр, не дало нужных результатов, в работе (Tikhonov and Zhorov 2012) было предложено внести дополнительные изменения .

Такое выравнивание позволило построить 3D-модель, в которой остатки, взаимодействующие с TTX ориентированы внутрь поры, в то время как длинные полярные остатки, находящиеся на участках до и после селективного фильтра, образуют межспиральные водородные связи, которые стабилизируют структуру. При этом особенно стоит отметить, что в модели с прямым и модифицированным выравниваниями заметно отличается расположение кольца внешних карбоксилатов, что крайне важно для связывания µ-конотоксинов, как будет видно далее (рисунок 29) .

Рисунок 29. Трёхмерное наложение моделей Nav1.4, построенных на основании разных выравниваний. Зелёным показана модель с модифицированным выравниванием, а лиловым – модель с прямым выравниванием. Дополнительно показаны внешние карбоксилаты .

Однако, тот факт, что модель хорошо позволяет воспроизвести экспериментальн ые данные по взаимодействию с TTX, ещё не значит, что её можно использовать для изучения действия других лигандов, в данном случае – µ-конотоксинов. Для проверки было проведено сравнение моделей с прямым и модифицированным выравниваниями с точки зрения их способности описать экспериментальные данные по наиболее изученному из µконотоксинов: GIIIA. Наиболее важным для нас результаты циклического мутационного анализа, который позволили определить пары взаимодействующих остатков между токсином и каналом (таблица 3). Эти данные были получены в работе (Choudhary, Aliste et al., 2007), и ими можно воспользоваться, чтобы определить набор дистанционных ограничений, которые накладываются на водородные связи или солевые мостики между боковыми цепями соответствующих остатков в канале и токсине .

–  –  –

Известно достаточно большое количество взаимодействующих пар (12 пар), но в этом заключается и определённая сложность. В некоторых случаях невозможно однозначно описать взаимодействие между остатком токсина и остатком канала. Например, в соответствии с экспериментальными данными, внешние карбоксилаты III и IV субъединиц канала участвуют в контактах сразу с 3 остатками токсина (Choudhary, Aliste et al., 2007) .

При этом экспериментально полученная энергия взаимодействия для этих контактов различается. Поэтому при расчётах используются только те взаимодействия, которые отмечены в таблицы как «сильные» .

Для того, чтобы посадить токсин в пору, не внося никаких дополнительных искажений в структуру канала, GIIIA поместили во внеклеточное пространство над каналом. При этом токсин обращён к каналу стороной, которая содержит основные остатки (K9, K11, R13, и R19 K16), центр масс находится на оси поры и расстояние между плоскостью, образованной альфа углеродами внешних карбоксилат и ближайшим к ним альфа углеродом токсина составляет 11. На этом расстоянии токсин не взаимодействует с каналом. Но если использовать только одну стартовую позицию для токсина, то может оказаться, что его изначальная ориентация относительно оси канала скажется на его итоговой посадке. Чтобы исключить такую возможность, были созданы 36 стартовых конформаций токсина, последовательно вращая его вокруг оси поры. Используя описанные выше взаимодействия в качестве ограничений, проводилась МК-минимизация для каждого варианта. На рисунке 30 представлено изменение энергии системы в процессе минимизации для одного из вариантов .

Рисунок 30. Изменение энергии системы в процессе минимизации при посадке GIIIA в канал. Показаны общая конформационная энергия системы, энергия на ограничителях и энергия взаимодействия токсина с каналом .

После этого был проведён второй этап минимизаций, на котором все ограничения на взаимодействия были сняты. Это необходимо, чтобы убедиться, что полученные структуры достаточно стабильны. В целом, для различных стартовых точек была получена одна и та же ориентация токсина как в случае прямого, так и в случае модифицированного выравнивания. Рисунок 31 показывает наложение всех моделей, полученных при различных ориентациях GIIIA в модели с модифицированным выравниванием. Этот результат показывает, что система ограничений, предложенная выше, однозначно определяет ориентацию токсина в канале .

Рисунок 31. Наложение всех конформаций, полученных при различных ориентациях GIIIA в модели с модифицированным выравниванием .

И в одной, и во второй моделях крупный токсин не может поместиться в узкий селективный фильтр. Скорее, он зависает над сужающейся порой таким образом, что боковые цепи основных остатков связывались с внешними карбоксилатами, которы е расположены выше селективного фильтра, ближе к внеклеточному пространству. Важно отметить, что центр масс токсина смещается от оси канала в сторону повтора III (рисунок 32). При этом С-конец спирали токсина сильно взаимодействует с Р2-спиралью повтора III .

А N-конец токсина направлен во внеклеточное пространство. Центр же токсина, где находится аргинин R13, оказывается в поре около повтора II .

Рисунок 32. Ассиметричная посадка GIIIA в канал. Спираль токсина (оранжевый) смещена относительно оси канала в сторону повтора III .

Несмотря на общее сходство в ориентации токсина, результаты показали значительное преимущество модели с модифицированным выравниванием. Для этой модели, СО альфа углеродов в области, где связывается токсин, от структуры NavAb составляло 1.39. После снятия ограничений и повторной минимизации токсин слегка передвинулся (СО - углеродов между моделью с ограничениями и без них составляет всего лишь 0,5 ). Для модели с прямым выравниванием, в свою очередь, наложение ограничений привело к искажению структуры канала, и СО между -карбонами модели и NavAb составляло 2.04. При этом повторная минимизация без ограничений приводила к значительному смещению токсина ~ 4.5 (рисунок 33). Энергия взаимодействия между токсином и каналом в модели с прямым выравниванием была на 33% выше («хуже»), чем в модели с модифицированным выравниванием .

Рисунок 33. Модели посадки GIIIA полученные для варианта с прямым и модифицированным выравниваниями. Жёлтым показана ориентация GIIIA в модели с прямым выравниванием, а оранжевым – ориентация в модели с модифицированным выравниванием .

Для количественной оценки модели дополнительно была рассчитана энергия взаимодействий между отдельными остатками токсина и канала в двух расслабленных моделях (рисунок 34). Используя экспериментально подсчитанные значения энергии парного взаимодействия, приведённые в таблице, можно построить график, отображающий связь между энергией, полученной в модели для данной пары, и экспериментальными данными. Следует отметить, что не все экспериментальные контакты были использованы в качестве ограничений во время процедуры стыковки. Для оценки связи между экспериментом и моделью был использован коэффициент корреляции Пирсона. Для модели с прямым выравниванием корреляция слаба (коэффициент корреляции Пирсона (Р) -0.44, достоверность р = 0,17). А для модели с модифицированным выравниванием коэффициент корреляции -0.71 и корреляция оказалась статистически значимой (р = 0,02). Это во многом неожиданный результат, поскольку модель хоть и ставит перед собой задачу объяснить происходящие в системе процессы, но далеко не всегда можно говорить о численной интерпретации данных .

Рисунок 34. Графики корреляций между экспериментально полученными энергиями взаимодействия остатков и расчётными значениями в моделях. Сверху: модель с прямым выравниванием; Снизу: модель с модифицированным выравниванием. Синим показаны контакты, использованные при построении модели .

Подводя итоги, можно констатировать, что, модель Nav1.4 с модифицированным выравниванием лучше подходит для изучения взаимодействия µ-конотоксинов с каналом. Поэтому для дальнейшего исследования используется только эта модель .

3.3.2. Связывание GIIIA Рассмотрим подробнее результаты, полученные для взаимодействия GIIIA с Nav1.4 .

Хотя модель легко воспроизводит многие экспериментально установленные специфические взаимодействия между токсином и каналом, есть и такие, которые не наблюдаются в модели. Наприме, контакт между гидроксипролином h P17 и M1240 (Dudley, Chang et al., 2000). h P17 в модели находится недалеко от Р2-спирали третьего повтора, где находится M1240, но не было найдено никаких специфических взаимодействий между двумя остатками. Гидроксильная группа h P17 при этом находится в непосредственной близости с D1248 на С-конце P2-спирали III повтора. Чтобы изучить возможный контакт между этими остатками, были введены дополнительные ограничения на расстояние между гидроксильной группой h P17 и кислородом D1248, что привело к модели, показанной на рисунке 35. При этом контакты, полученные в предыдущей модели, сохраняются. Такая модель обеспечивает структурное объяснение роли h P17 в действии токсина. Важно отметить, что новый контакт, который предсказывают расчеты, пока не был описан в экспериментальных исследованиях. D1248 не имеет аналогов в гомологичных позициях других повторов, и, таким образом, новый контакт может быть важным фактором, определяющим ориентацию токсина в канале .

Рисунок 35. Модель посадки GIIIA с учётом ограничителей на hP17 .

Таблица 4. Энергия взаимодействия между остатками токсина и остатками канала, рассчитанная в модели .

Ключевые остатки канала выделены жёлтым, ключевые остатки токсина – голубым. Красным обозначены экспериментально обнаруженные контакты .

–  –  –

Рисунок 36. Наиболее важные контакты, образуемые остатками токсина GIIIA с остатками канала. Показаны различные пары взаимодействующих остатков .

В работе (Xue, Ennis et al., 2003), авторы обнаружили, что замена D762 и E765 во втором повторе канала влияет на связывание GIIIA. В полученной модели эти остатки расположены далеко от GIIIA и не могут непосредственно взаимодействовать с ним. Тем не менее, в соответствии с кристаллической структурой NavAb, D762 является важным стабилизатором пространственной структуры поры, образуя солевой мостик с R399. E765 расположен на той же стороне P2-спирали, что и D762, будучи, таким образом, обращён к P1-спирали I повтора. Согласно модели, как D762, так и E765 образуют контакты с R399, что способствует стабилизации пространственной укладки поры. Таким образом, можно предположить, что описанные в работе замещения D762 и E765 аланином (Xue, Ennis et al.,

2003) не влияют на связывание токсина непосредственно, а вызывают аллостерические эффекты .

Рисунок 37. Межсубъединичные контакты в Nav1.4 играют важную роль в стабилизации структуры. Показаны остатки, образующие устойчивую связь, стабилизирующую структуру поры канала .

В полученном комплексе, тело GIIIA, содержащее -спираль, в основном взаимодействует с P2-спиралью III повтора. Примечательно, что в позиции 1248 находится аспартат, который является привлекательной мишенью для электростатического взаимодействия с токсином. Другие повторы канала имеют нейтральные или положительно заряженные остатки в соответствующих позициях выравнивания последовательностей (рисунок 28). Другой возможной причиной асимметричного связывания токсина является асимметрия внешнего карбоксилатного кольца. В то время как внешние карбоксилаты в повторах I, II, IV находятся в одинаковых позициях, в повторе III D1241 сдвигается на одну позицию по отношению к C-концу .

Если взглянуть ещё раз на график соответствия энергии взаимодействия в модели экспериментально рассчитанным энергиям, то видно, что, для некоторых «средних»

контактов между внешними карбоксилатами повторов III и IV с токсином точки графика сильно отклоняются от прямой (рисунок 34). Возможной причиной является то, что в расчётах учитывается только одна структура, таким образом, нельзя учесть энтропийный вклад в энергию. Из-за большой конформационной свободы длинных боковых цепей основных остатков токсина и кислотных остатков канала, экспериментально определенные энергии парного взаимодействия складываются из многих микросостояний комплекса канал-белок. Чтобы исследовать эту возможность, были проведены дополнительные расчеты. Используя модель, описанную выше, были сгенерированы 20 000 стартовых точек со случайными конформациями боковых цепей основных остатков GIIIA. Для каждой сгенерированной конформации была проведена минимизация энергии с учётом ограничений, накладываемых на модель. Результаты показаны на рисунке 38 .

Рисунок 38. Ансамбль низкоэнергетических комплексов, полученных для боковых цепей токсина. Шариками показаны центры заряженных групп остатков токсина .

В ансамбле самых низкоэнергетических комплексов, заряженные концы гибких боковых цепей остатков токсина колебались в пределах 5. Это позволяет предположить, что остатки токсина могут формировать контакты с различными остатками токсина, «переключаясь» между ними, не нарушая при этом общей пространственной ориентации токсина (рисунок 39) .

Рисунок 39. Пример переключения контактов между токсином и каналом. Сверху: R13 образует контакты с III и IV субъединицами, а K11 – с I и II; Снизу: в другой конформации остатки меняются партнёрами .

К основным остаткам токсина относятся K11, R13, K16 и R19. Согласно полученным результатам, К11 взаимодействует преимущественно с повторами I, III и IV, но не взаимодействует с повтором II в связи с наличием R13 в непосредственной близости от неё .

R13 взаимодействует в основном с I и II, но может взаимодействовать с любым повтором из-за его длинной боковой цепи и пространственной ориентации токсина. K16 в середине спиральной части токсина может взаимодействовать с любым повтором канала, в зависимости от контактов, образуемых другими остатками токсина. R19 может достичь только повтора III, потому что он расположен близко к C-концу токсина, сразу после спиральной части, и никак не может дотянуться до повторов I и II. Таким образом, полученная модель предсказывает, что в рамках одной модели посадки тела токсина в пору канала его остатки могут образовывать различные контакты с остатками канала .

3.3.3. Связывание PIIIA Теперь рассмотрим проблему взаимодействия с каналом другого конотоксина – PIIIA. Согласно экспериментальным исследованиям, главную роль в связывании PIIIA с каналом играют остатки R12, R14, K17, и R20 (McArthur, Singh et al., 2011). Тем не менее, в отличие от GIIIA, парные контакты между PIIIA и Nav1.4 неизвестны. Таким образом, методология посадки, которая применялась для GIIIA, не подходит для PIIIA. Так как пространственная структура PIIIA подобна структуре GIIIA, можно предположить, что в комплексе с каналом -спираль и ключевые основные остатки PIIIA ориентированы так же, как у GIIIA (рисунок 40) .

Исходя из этого предположения была сгенерирована стартовая позиция для расчётов. Модель, полученная в результате минимизации показана на рисунке 41. В полученной модели ключевые остатки PIIIA участвуют в следующих взаимодействиях: R12 / E403, R14 / E758, К17 / D1241, К17 / D1532, R20 / D1532 (рисунок 41 и таблица 5) .

Рисунок 40. Трёхмерные структуры и выравнивание последовательностей между конотоксинами GIIIA, PIIIA и KIIIA. Слева направо: GIIIA (оранжевый), PIIIA (синий) и KIIIA (зелёный) .

Рисунок 41. Модель посадки и контакты, образованные токсином PIIIA с каналом. Показаны различные взаимодействующие пары остатков .

Таблица 5. Энергия взаимодействия между остатками токсина PIIIA и остатками канала, рассчитанная в модели .

Ключевые остатки канала выделены жёлтым, ключевые остатки токсина – голубым .

–  –  –

Рисунок 42. Трёхмерное наложение моделей посадки GIIIA и PIIIA. Модель GIIIA показана жёлтым, модель PIIIA – оранжевым .

Хотя полученные конформации GIIIA и PIIIA весьма схожи (СО для -углерода ключевых основных остатков составляет 1,8 ), есть определенные различия в парных взаимодействиях между токсином и каналом. Например, для GIIIA важную роль играет специфическое взаимодействие между D12 и T759 (Choudhary, Aliste et al. 2007), что отражается и в модели. В PIIIA отсутствует подходящий остаток в гомологичной позиции .

В результате, в модели R12 в PIIIA образует контакты отличные от гомологичного K11 в GIIIA. R12 в PIIIA сильно взаимодействует с I повтором, в то время как K11 в GIIIA взаимодействует с III и IV потворами, потому что D12 стерически препятствует возникновению контактов между K11 и I повтором. R14 формирует тесные контакты с I и II повторами, но слабо взаимодействует с III повтором, и ещё слабее с IV повтором. К17 взаимодействует с повторами III и IV, как и его аналог в GIIIA, К16. Наконец, R20 (аналог R19 в GIIIA) взаимодействует только с IV повтором. Как и в случае с GIIIA, показанное макросостояние может быть представлено множеством микросостояний, которые существуют благодаря переключению контактов между длинными гибкими боковыми цепями .

В работе (McArthur, Singh et al., 2011) авторы исследовали потенциал-зависимое действие PIIIA и его мутантов и определяли глубину () остатков в электрическом поле мембраны. Значения располагались по убыванию следующим образом: R14 K17 R20 R12 S13 R2 G6.

В полученной модели, относительная глубина заряженных остатков (определяется расположением их заряженных групп) располагается следующим образом:

R14 K17 R20 R12 R2; S13 находится примерно на той же глубине, что и R12, в то время как G6 и R2 расположены на N-конце токсина во внеклеточном пространстве .

Таким образом, ориентация токсина PIIIA в модели согласуются с экспериментальными данными. Следует отметить, что эти данные не использовались при получении модели, что говорит о хорошем качестве модели .

3.3.4. Связывание KIIIA KIIIA представляет собой самый маленький µ-конотоксин в семействе, и не содержит аналога N-концевого сегмента, который присутствует в GIIIA и PIIIA .

Поэтому его аминокислотная последовательность содержит только три основных остатка в -спирали .

Данные по парным взаимодействиям (McArthur, Singh et al., 2011) показывают, что все эти остатки взаимодействуют с внешним карбоксилатом III повтора (D1241). Используя эти данные, можно наложить ограничения по расстоянию на водородные связи между боковыми цепями токсина и карбоксилатом. Последующие минимизации привели к модели, где положение и ориентация -спирали KIIIA аналогична ориентации спиралей PIIIA и GIIIA. Весь набор контактов, образуемых основными остатками KIIIA показан на рисунке 43 и в таблице 6. Важно отметить, что из-за своего небольшого размера и отсутствия N-концевого сегмента, KIIIA не может перекрыть пору канала своим телом .

Таблица 6. Энергия взаимодействия между остатками токсина KIIIA и остатками канала, рассчитанная в модели .

Ключевые остатки канала выделены жёлтым, ключевые остатки токсина – голубым .

–  –  –

Рисунок 43. Модель посадки и контакты, образованные токсином KIIIA с каналом. Показаны различные взаимодействующие пары остатков .

Наиболее интересно то, что в присутствии KIIIA TTX может достигать своего сайта связывания в селективном фильтре, хотя кинетика связывания замедляется относительно свободного канала (Zhang, McArthur et al., 2009). Чтобы найти объяснение этих экспериментальных данных, можно построить модель канала, в комплексе с KIIIA и TTX .

Модель взаимодействия TTX с каналом была предложена ранее в работе (Tikhonov and Zhorov, 2012). Поэтому для построения общей модели использовалась модель TTX из этой работы и модель посадки KIIIA, полученная нами ранее .

Рисунок 44. Модель одновременной посадки TTX и KIIIA в канале. Слева: стрелочкой показан возможный путь выхода TTX из канала при наличии KIIIA в системе; справа: тело токсина смещено в сторону III субъединицы, что даёт больше места для посадки TTX .

Как видно из рисунка 44, полученная суммарная модель не приводит к стерическим конфликтам в положении токсинов. Компактная молекула TTX глубоко связана в узкой части поры между кольцом DEKA и кольцом внешних карбоксилатов, в то время как более крупный конотоксин не проникает глубоко в пору, а приближается к внешним карбоксилатам с внеклеточной стороны, образуя солевые мостики с ними. Каждый карбоксилат может служить акцептором для нескольких водородных связей, что позволяет образовывать одновременные контакты с TTX и KIIIA. Внешний карбоксилат в III повторе, который, в отличие от других внешних карбоксилатов, не смотрит в сторону оси канала, не взаимодействует с TTX, но формирует солевой мостик с R10 в KIIIA .



Pages:   || 2 |
Похожие работы:

«ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ГОРОДА МОСКВЫ ГИМНАЗИЯ № 1562 ИМЕНИ АРТЕМА БОРОВИКА 109341., г.Москва, ул. Братиславская, д. 4. Е-mail:gimn1562@уаndex.ru тел., факс: (495) 349-00-11 Рабочая программа Биология (внеурочная деятельность) "Р...»

«КАЗАНСКИЙ ФЕДЕРАЛЬНЫЙ УНИВЕРСИТЕТ ИНСТИТУТ ЭКОЛОГИИ И ПРИРОДОПОЛЬЗОВАНИЯ Кафедра прикладной экологии О.В. НИКИТИН КОНТРОЛЬ ИСТОЧНИКОВ ЗАГРЯЗНЕНИЯ АТМОСФЕРНОГО ВОЗДУХА Конспект лекций Казань – 2015 УДК 504.064:504.3.054 Принято на заседании кафедры прикладной экологии Протокол №...»

«ЭКОЛОГИЯ МЕТОД ОПРЕДЕЛЕНИЯ ИНДЕКСА ОРНИТОЛОГИЧЕСКОЙ ПРИВЛЕКАТЕЛЬНОСТИ ХОЗЯЙСТВЕННЫХ ОБЪЕКТОВ И.Р. Еналеев Кафедра биоэкологии Татарский государственный гуманитарно-педагогический университет...»

«НЕКОММЕРЧЕСКОЕ ПАРТНЕРСТВО "Центр экологии топливно-энергетического комплекса"ПОЛЕВОЕ РУКОВОДСТВО ПО ЛИКВИДАЦИИ РАЗЛИВОВ НЕФТИ: НАСОСНОЕ ОБОРУДОВНИЕ ДЛЯ ЛРН МОСКВА Перед Вами демонстрационная версия (проект) одной из частей "ПОЛЕВОГО РУКОВОДСТВА ПО ЛИКВИДАЦИИ РА...»

«Евроазиатская региональная ассоциация зоопарков и аквариумов Правительство Москвы Московский государственный зоологический парк ДРОФИНЫЕ ПТИЦЫ ПАЛЕАРКТИКИ: РАЗВЕДЕНИЕ И ОХРАНА МОСКВА – 2008 ЕВРОАЗИАТСКАЯ РЕГИОНАЛЬНАЯ АССОЦИАЦИЯ ЗООПАРКОВ И АКВАРИУМОВ EURASIA...»

«ВЕСТНИК УДМУРТСКОГО УНИВЕРСИТЕТА 55 БИОЛОГИЯ. НАУКИ О ЗЕМЛЕ 2011. Вып. 4 УДК 582.4 Т.А. Полянская ОНТОГЕНЕТИЧЕСКАЯ СТРУКТУРА ЦЕНОПОПУЛЯЦИЙ БОРЕАЛЬНЫХ ДЛИННОКОРНЕВИЩНЫХ РАСТЕНИЙ Экологиче...»

«УДК 619:616:34.008.314.4:636.2-053.2 РЕЗУЛЬТАТЫ ИЗУЧЕНИЯ НАРУШЕИЯ ОБМЕНА ВЕЩЕСТВ У КОРОВ И ЗАБОЛЕВАЕМОСТИ ТЕЛЯТ ДИСПЕПСИЕЙ Пахомов Г.А. – к.в.н., доцент Казанская государственная академия ветеринарной медицины имени Н.Э. Баумана тел.: (843...»

«Сценарий внеклассного мероприятия Питание и здоровье Автор: учитель биологии высшей категории Новикова Н.Н. ЦЕЛИ: Пропагандировать здоровый образ жизни; активизировать познавательную деятельность учащихся; заинт...»

«МУНИЦИПАЛЬНОЕ БЮДЖЕТНОЕ УЧРЕЖДЕНИЕ КУЛЬТУРЫ "МЕЖПОСЕЛЕНЧЕСКАЯ БИБЛИОТЕКА СОВЕТСКОГО РАЙОНА" КРАЕВЕДЧЕСКИЙ КАЛЕНДАРЬ: ЗНАМЕНАТЕЛЬНЫЕ И ПАМЯТНЫЕ ДАТЫ СОВЕТСКОГО РАЙОНА НА 2013 ГОД Советский 91.9:63 К 78 Составители: Яблочкова В. В. – зав. эколого-краеведческим отделом МБУК МБСР Верушкина...»

«RU 2 464 801 C2 (19) (11) (13) РОССИЙСКАЯ ФЕДЕРАЦИЯ (51) МПК A23K 3/02 (2006.01) ФЕДЕРАЛЬНАЯ СЛУЖБА ПО ИНТЕЛЛЕКТУАЛЬНОЙ СОБСТВЕННОСТИ (12) ОПИСАНИЕ ИЗОБРЕТЕНИЯ К ПАТЕНТУ (21)(22) Заявка: 2010136773/13, 03.09.2010 (72) Автор(ы): Отрошко Сергей Алексеевич (RU), (24) Дата начала отсчет...»

«Научный журнал НИУ ИТМО. Серия "Экономика и экологический менеджмент" № 3, 2016 УДК 338.984 Управление финансовыми потоками в логистических комплексах Канд. экон. наук Негреева В.В. v.negreeva@mail.ru Цимбалист-Колесникова И.А. tsimbala@mail.ru Шевченко Я.В. folein@mail.ru Университет ИТМО 191002, Россия, Санкт-Петербург, ул. Ломоносо...»

«Яфарова Гузель Гульусовна ФУНКЦИОНАЛЬНОЕ СОСТОЯНИЕ ДВИГАТЕЛЬНЫХ ЦЕНТРОВ СПИННОГО МОЗГА В УСЛОВИЯХ ИЗМЕНЕНИЯ СУПРАСПИНАЛЬНЫХ ВЛИЯНИЙ 03.00.13физиология Диссертация на соискание ученой степени кандида...»

«1.Результаты освоения курса внеурочной деятельности "Комнатное цветоводство" 5 класс Познавательные УУД Ученик будет знать: Биологические и морфологические особенности декоративных растений, их роль в жизни челове...»

«УДК 577.357, 591.31 Володяев Илья Владимирович СВЕРХСЛАБОЕ ИЗЛУЧЕНИЕ И ОПТИЧЕСКОЕ ВЗАИМОДЕЙСТВИЕ ЯЙЦЕКЛЕТОК И ЗАРОДЫШЕЙ ШПОРЦЕВОЙ ЛЯГУШКИ 03.00.30-03 — биология развития, эмбриология Автореферат диссертации на соискание ученой степени канди...»

«ИННОВАЦИОННЫЙ ЦЕНТР РАЗВИТИЯ ОБРАЗОВАНИЯ И НАУКИ INNOVATIVE DEVELOPMENT CENTER OF EDUCATION AND SCIENCE Сельскохозяйственные науки в современном мире Выпуск II Сборник научных трудов по итогам международной научно-практической конференц...»

«ИТОГИ ОТБОРОЧНОГО ТУРА ПО БИОЛОГИИ. В ФИНАЛ ПРОХОДЯТ ПРИМЕРНО 25% УЧАСТНИКОВ (ВЫДЕЛЕНЫ ЦВЕТОМ). Чтобы найти свой результат, пользуйтесь поиском (Ctrl-F)! 7 класс (проходной балл-16) Максимальный балл 100 8 класс (проходной балл-16) Максимальный балл 100 9 класс (проходной балл-16) Мак...»

«Глобальный Экологический Фонд (ГЭФ), объединяющий 182 страны, в партнерстве с международными организациями, гражданским обществом и частным бизнесом предоставляет гранты развивающимся странам и странам с переходной экономикой для совместного решения местных, национальных и глобальных экологических проблем с целью достижения устойчивого ра...»

«ВЕСТНИК БУРЯТСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА Спецвыпуск С/2012 докринологии. 1956. №6. С. 75-78.7. Скакун Н.П., Олейник А.Н . Сравнительное действие атропина и метацина на внешнесекреторную функцию печени // Фармакология и ток...»

«УРОК ПО БИОЛОГИИ В 8 КЛАССЕ ГИГИЕНА ЗРЕНИЯ "ДОРОЖЕ АЛМАЗА СВОИ ДВА ГЛАЗА" ЦЕЛЬ УРОКА: изучить основные правила гигиены зрения.ЗАДАЧИ УРОКА: Раскрыть причины дальнозоркости и близорукости, показать примеры их предупреждения; Разъяснить требования, предъявляемые гигиеной зрения;...»

«РОССИЙСКАЯ ФЕДЕРАЦИЯ (19) (11) (13) RU 2 543 278 C2 (51) МПК A01G 1/00 (2006.01) A01B 79/02 (2006.01) B82Y 40/00 (2011.01) ФЕДЕРАЛЬНАЯ СЛУЖБА ПО ИНТЕЛЛЕКТУАЛЬНОЙ СОБСТВЕННОСТИ (12) ОПИСАНИЕ ИЗОБРЕТЕНИЯ К ПАТЕНТУ 2012151771/13, 03.12.2012 (21)(22) Заявка: (72) Автор(ы): Яппаров Ахтам Хусаинович (RU), (24) Дата начала...»

«Annotation Грегори Бейтсон выдающийся мыслитель XX века, философ, эколог, кибернетик и системный теоретик, внесший значительный вклад в антропологию, психиатрию и теорию коммуник...»

«РОССИЙСКАЯ ФЕДЕРАЦИЯ (19) (11) (13) RU 2 543 812 C1 (51) МПК A01N 61/00 (2006.01) ФЕДЕРАЛЬНАЯ СЛУЖБА ПО ИНТЕЛЛЕКТУАЛЬНОЙ СОБСТВЕННОСТИ (12) ОПИСАНИЕ ИЗОБРЕТЕНИЯ К ПАТЕНТУ 2013153460/13, 03.12.2013 (21)(22) Заявка...»

«ISSN 0869-4362 Русский орнитологический журнал 2015, Том 24, Экспресс-выпуск 1181: 3053-3065 Данные по биологии восточного подвида дрофы Otis tarda dybowskii в Даурии О.А.Горошко Второе издание. Первая публикация в 2008* Дрофа Otis tarda занесена в Красный спис...»

















 
2018 www.new.z-pdf.ru - «Библиотека бесплатных материалов - онлайн ресурсы»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 2-3 рабочих дней удалим его.