Научная статья на тему 'Квантово-химическое моделирование механизма действия лигандов nr 3c 4 рецептора'

Квантово-химическое моделирование механизма действия лигандов nr 3c 4 рецептора Текст научной статьи по специальности «Химические науки»

CC BY
180
76
i Надоели баннеры? Вы всегда можете отключить рекламу.
Журнал
Auditorium
Область наук
Ключевые слова
КВАНТОВО-ХИМИЧЕСКАЯ МОДЕЛЬ / ЛИГАНДЫ NR 3C 4-РЕЦЕПТОРА / ПРОТОНИРОВАНИЕ / АГОНИСТЫ И АНТАГОНИСТЫ РЕЦЕПТОРА / ПЕРЕНОС ЭЛЕКТРОНА / quantum-chemical model / ligands of NR3C4 receptor / protonisation / receptor agonists and antagonists / transportation of electron

Аннотация научной статьи по химическим наукам, автор научной работы — Кузнецова Н. Б., Кузнецов П. Е.

На основе рентгеноструктурных данных из базы Protein Data Base и квантово-химических расчетов предложена модель взаимодействия лигандов NR 3C 4 рецептора. Ключевую роль в проявлении агонистической активности, вероятно, играет протонирование лигандов. Предлагаемый механизм действия заключается в том, что лиганды протонируются, а затем происходит перенос электрона с рецептора на агонист. Согласно предложенной модели, антагонисты не протонируются и не участвуют в переносе электрона. Данный механизм подтверждается расчетами взаимодействия с рецептором для 22 агонистов и антагонистов. При этом объясняются экспериментальные данные о типе и степени их активности.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

QUANTUM-CHEMICAL MODELLING OF MECHANISM OF ACTION OF NR

A new model of interaction of ligands of NR3C4 receptor is offered on the basis of X-ray structural data from Protein Data Base and quantum-chemical computation. It seems that protonisation of ligands plays the key role in development of agonist activities. The mechanism of action we propose is as follows: ligands are protonized and this is followed by transportation of electron from receptor to agonist. In accordance with the suggested model, antagonists are not protonized and they do no take part in transportation of electron. This mechanism is supported by computation of interaction with the receptor for 22 agonists and antagonists. Experimental data on the type and degree of activities of agonists and antagonists are being explained.

Текст научной работы на тему «Квантово-химическое моделирование механизма действия лигандов nr 3c 4 рецептора»

УДК 544.18/577.29:577.17

КВАНТОВО-ХИМИЧЕСКОЕ МОДЕЛИРОВАНИЕ МЕХАНИЗМА ДЕЙСТВИЯ

ЛИГАНДОВ NR3C4 РЕЦЕПТОРА

© 2014 Н. Б. Кузнецова1, П. Е. Кузнецов2

1канд. хим. наук, доцент каф. химии e-mail: moscow 1978@mail.ru

2

докт. хим. наук, проф. каф. химии e-mail: kuznetsovpe@mail. ru

Курский государственный университет

На основе рентгеноструктурных данных из базы Protein Data Base и квантовохимических расчетов предложена модель взаимодействия лигандов NR3C4 - рецептора. Ключевую роль в проявлении агонистической активности, вероятно, играет протонирование лигандов. Предлагаемый механизм действия заключается в том, что лиганды протонируются, а затем происходит перенос электрона с рецептора на агонист. Согласно предложенной модели, антагонисты не протонируются и не участвуют в переносе электрона. Данный механизм подтверждается расчетами взаимодействия с рецептором для 22 агонистов и антагонистов. При этом объясняются экспериментальные данные о типе и степени их активности.

Ключевые слова: квантово-химическая модель, лиганды NR3C4-рецептора,

протонирование, агонисты и антагонисты рецептора, перенос электрона.

NR3C4-рецептор относится к стероидным ядерным рецепторам. Известно, что его лиганд индуцирует трансформацию NR3C4-рецептора, после чего комплекс лиганд -рецептор транслоцируется в ядро и связывается с ДНК, регулируя экспрессию генов. Цитозольный NR3C4-рецептор является переносчиком активной формы гормона в ядро, являясь там активным ядерным рецептором. При этом ядерный комплекс связывается с ядерной фракцией и полианионами, а цитозольный - не связывается [Щитовидная железа 1998: 150].

В Protein Data Base [[Сайт] Protein Data Base] содержатся определенные рентгеноструктурным анализом координаты атомов лиганд связанных NR3C4-рецепторов. Среди них комплексы с дигидротестостероном, а также с селективными модуляторами NRзC4-рецептора BMS-564929, S-24, C-23.

Рассмотрим механизм действия лигандов NR3C4-рецептора исходя из

экспериментальных данных. Анализ структуры четырех комплексов NR3C4-рецептора с дигидротестостероном, BMS-564929, S-24 и C-23 показывает, что все лиганды содержат группы, которые могут протонироваться и у которых есть неподеленные электронные пары, а именно: у дигидротестостерона - это группа С=О, у BMS-564929 - группа C=N, у S-24 -группа С=^ а у C-23 -группа N02. Геометрически эти группы в комплексах с рецептором находятся рядом с Arg82 на расстоянии 2-3А. Можно предположить, что протонирование, возможно, протекает в среде рецептора с переносом протона от Arg82 на протонируемую группу лиганда.

ХИМИЧЕСКИЕ НАУКИ

В пользу протонирования указывает и тот факт, что в электронной структуре агонистов на протонируемых группах фиксируется один электрон вблизи ВЗМО-уровня (верхняя занятая молекулярная орбиталь). Так, на ВЗМО-уровне атома кислорода С=О группы дигидротестостерона находится электронная плотность 1,04 ат. ед. Вблизи C=N группы BMS-564929 на ВЗМО ароматического кольца группируется > 1 ат. ед. Кроме того, протонированная производная BMS-564929 по C=N группе на 13 ккал/моль устойчивее других протонированных производных этого соединения. То же можно сказать и о селективных модуляторах NRзC4-рецептора S-24 и С-23.

Ранее мы уже рассматривали роль протонирования в биохимических процессах [Кузнецова и соавт. 2006: 18-24; Kuznetsova 2006: 288-294]. При этом были выявлены условия протонирования, в частности наличие границы раздела фаз. Существует и экспериментальное подтверждение протонирования лигандов NR3C4-рецептора. Так, в экспериментах ВЭЖХ-МС формируются условия, которые способствуют протонированию агонистов NR3C4-рецептора, - химическая ионизация на границе раздела фаз в хроматографической колонке между твердой фазой (например, С 18) и элюентом (0,1 % муравьиной кислоты в воде или ацетонитриле). Отметим, что для дигидротестостерона, который выходит из колонки в виде массы (М+1)+, не существует других центров протонирования, кроме группы С=О с неподеленными электронными парами на атоме кислорода.

Таким образом, вполне возможно протонирование лигандов NR3C4-рецептора в среде рецептора или на его поверхности.

На рисунке 1 представлен комплекс дигидротестостерона с NR3C4 рецептором.

Рис. 1. Рентгеноструктурные данные комплекса дигидротестостерона с NR^Q-рецептором

Анализ экспериментальных данных показывает, что ОН-группа дигидротестостерона образует две водородные связи молекулы с рецептором, а именно с NH-группой бокового остатка аминокислоты Asn34 и ОН-группой аминокислоты Thr206. С=О-группа дигидротестостерона находится на расстоянии 2,4 А от атома водорода Arg82.

Активный центр NR3C4 рецептора с зарядом минус 1 представляет собой две аминокислотные цепи, расположенные на разных спиралях рецептора и сближенные в

Auditorium: электронный научный журнал Курского государственного университета. 2014. № 1

Кузнецова Н. Б., Кузнецов П. Е. Квантово-химическое моделирование механизма действия лигандов МК3С4-рецептора

пространстве. Одна - нейтральная аминокислотная последовательность с активной ОН-группой аминокислоты Thr206, она образует водородную связь с атомом водорода ОН-группы дигидротестостерона. Другая аминокислотная последовательность содержит активные аминокислоты Asn34 и Glu35 с зарядом минус 1. Asn34 образует водородную связь с атомом кислорода ОН-группы дигидротестостерона. Можно предположить, что активный центр является донором электронной плотности.

Расчеты проводились с использованием программ Hyper Chem Release 7 for Windows, Discovery Studio Client v2.5.0.9164. Рассматривалась упрощенная модель, которая описывает взаимодействие дигидротестостерона с NR3C4-рецептором в виде комплекса дигидротестостерона с активным центром рецептора, состоящим, как отмечалось ранее, из двух цепей аминокислот (см. рис. 2).

с активным центром NR3C4-рецептора

Модель комплекса сформирована в соответствии с рентгеноструктурными данными: активный центр рецептора находится в виде цвиттер-иона, причем расположение концевых групп активного центра рецептора соответствует экспериментальным данным; ОН-группа дигидротестостерона образует две водородные связи с Asn34 и Thr206.

Как предполагалось ранее, дигидротестостерон протонирован по атому кислорода С=О-группы. Также на расстоянии 2,4 А от него находится подвижный атом водорода Arg82 (см. рис. 1), который может диффундировать к атому кислорода СО-группы. Более того, расчеты показали, что этот процесс термодинамически разрешен. Рентгеноструктурные данные свидетельствуют о том, что комплекс содержит довольно много молекул воды, что тоже способствует протонированию.

Основное состояние комплекса на рисунке 2 триплетное, так как энтальпия его синглетного состояния на 40,4 ккал/моль выше.

ХИМИЧЕСКИЕ НАУКИ

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

Таким образом, протонированный гормон как акцептор электрона, катализирует перенос электрона с рецептора и этим вызывает его трансформацию. Расчеты показали, что перенос электрона термодинамически выгоден, разность энтальпий образования конечного и исходного состояний равна -52,5 ккал/моль. Кроме того, энтальпия рассматриваемого комплекса с переносом электрона на 15,4 ккал/моль устойчивее комплекса в триплетном состоянии без переноса электрона.

Дигидротестостерон является эндогенным лигандом КК3С4-рецептора. Рассмотрим взаимодействие экзогенных лигандов с рецептором на примере BMS-564929.

На рисунке 3 представлен комплекс BMS-564929 с NR3C4 рецептором из Protein Data Base.

\

\

\

Рис. 3. Рентгеноструктурные данные комплекса BMS-564929 с NR^-рецептором

Анализ экспериментальных данных показывает, что ОН-группа BMS-564929 образует единственную водородную связь молекулы с рецептором, а именно с NH-группой аминокислоты Asn35. Мы видим, что в данном случае, как и в случае других экзогенных агонистов, образуется одна водородная связь, а не две, как в случае дигидротестостерона, но с тем же активным центром.

Рассмотрим электронную структуру BMS-564929. На фенильном кольце и цианогруппе находятся п-электроны. Поэтому очень высока вероятность протонирования атома азота C=N группы. Также на расстоянии 3,1 А от него находится подвижный атом водорода Arg82 (см. рис. 3).

Auditorium: электронный научный журнал Курского государственного университета. 2014. № 1

Кузнецова Н. Б., Кузнецов П. Е. Квантово-химическое моделирование механизма действия лигандов МК3С4-рецептора

Активный центр NR3C4 рецептора с зарядом минус 1 представляет собой аминокислотную последовательность, содержащую активные аминокислоты Asn35 и Glu36 с зарядом минус 1. Она образует водородную связь с атомом кислорода ОН-группы BMS-564929. Можно предположить, что активный центр является донором электронной плотности.

Рассматривалась упрощенная модель, которая описывает взаимодействие BMS-564929 с NR3C4 рецептором в виде комплекса BMS-564929 с активным центром рецептора. Модель комплекса сформирована в соответствии с рентгеноструктурными данными: активный центр рецептора находится в виде цвиттер-иона, ОН-группа BMS-564929 образует водородную связь с №13С4-рецептором.

Основное состояние комплекса на рисунке 4 триплетное, так как энтальпия его синглетного состояния на 14,4 ккал/моль выше.

Распределение спиновой и зарядовой плотности в комплексе после взаимодействия показывает, что в каждой из областей, отмеченных на рисунке 4 пунктиром, находится на валентной оболочке по одному электрону с параллельными спинами. До взаимодействия электроны с параллельными спинами находились на рецепторе. Следовательно, произошло перераспределение электронов. Так, заряд на протонированном BMS-564929 до взаимодействия был равен плюс 1, а после него - нулю. Следовательно, электрон перераспределился с рецептора на BMS-564929. Соответственно, на активном центре рецептора до взаимодействия заряд был равен минус 1, а после - нулю. Таким образом, взаимодействие BMS-564929 с активным центром рецептора инициирует перенос электрона с рецептора на BMS-564929.

Таким образом, протонированный BMS-564929 как акцептор электрона катализирует перенос электрона с рецептора и этим вызывает его трансформацию. Расчеты показали, что перенос электрона термодинамически выгоден, разность энтальпий образования конечного и исходного состояний равна -60,2 ккал/моль.

ХИМИЧЕСКИЕ НАУКИ

Для селективных модуляторов S-24, C-23 были получены аналогичные исследования и расчеты, с той лишь разницей, что ОН-группа этих соединений образовывала водородную связь с не с ЫН-, а с С=О-группой Asn35.

Таким образом, можно выделить следующие особенности механизма действия агонистов КК3С4-рецептора.

1. ОН-группа дигидротестостерона образует водородные связи с NH-группой бокового остатка аминокислоты Asn34 и ОН-группой аминокислоты Т^20б-рецептора. При этом происходит перенос электрона с рецептора на СОН+ группу протонированного дигидротестостерона.

2. Селективные агонисты протонируются по группам С=О, NO2, NH2, CN, ароматический N. Трехвалентный N или -ОН-группа образуют водородные связи с NH-или С=О-группами аминокислоты Asn35. При этом также происходит перенос электрона с рецептора на агонист.

3. Агонисты образуют водородные связи с разными атомами Asn35.

4. Все агонисты содержат активные группы, обладающие неподеленными электронными парами. Поэтому они могут протонироваться.

5. Агонисты обладают определенной жесткой геометрией. Так, в комплексе с рецептором группа ОН агониста должна быть связана водородной связью с Asn35, а протонируемая группа агониста должна находиться на расстоянии примерно 2,5А от Arg82.

6. Активный центр КК3С4-рецептора, по-видимому, является донором электронов. Протонированный агонист является акцептором электронов. Таким образом, протонирование является необходимым условием переноса электронов с рецептора на агонист.

Были исследованы электронные характеристики и комплексы с рецептором для 22 соединений агонистов и антагонистов КК3С4-рецептора. Анализировались следующие параметры:

- наличие групп, которые могут протонироваться;

- наличие электронной плотности на ВЗМО на этих группах;

- величина АЕ в эВ (разность энергий уровней ВЗМО и предВЗМО), отражающая скорость протонирования [Кузнецова 2006: 18-24];

- lgP;

- НСМО (нижняя свободная молекулярная орбиталь) в протонированной форме, отражающая возможность взаимодействия с рецептором.

Все соединения были разбиты по структурной формуле на 4 группы. В таблице 1 представлены структурные и электронные параметры этих соединений.

Таблица 1

Электронные параметры лигандов КК3С4-рецептора

№ Активные группы Кь нМ АЕ, эВ НСМО Протонир ованной формы Протонирование lgP Примечание

Соединения первой группы

1 ОН, С=О 0,2 0,26 -6,53 По ат. О СО группы, 1 ё на ВЗМО ат.О 3,93 Агонист

2 OH, 0,25 0,28 Гидролизуется - Агонист

Auditorium: электронный научный журнал Курского государственного университета. 2014. № 1

Кузнецова Н. Б., Кузнецов П. Е. Квантово-химическое моделирование механизма действия лигандов МК3С4-рецептора

СО

3 ОН 1542 0,94 -6,69 По бензольному кольцу 5,3 Агонист

Продолжение таблицы 1

4 NH 1485 0,9 -5,53 По ат^ 5,7 Агонист

5 ОН, CO 1300 0,45 -6,36 Нет 5,1 Антагонист

Соединения второй группы

6 OH, NO2 3,98 0,69 -6,44 По ат. О NО2 группы 2,04 Агонист

7 OH, NO2 6,11 0,26 -6,53 По ат. О NО2 группы 3,33 Агонист

8 OH, NO2 11,6 0,39 -6,52 По ат. О NО2 группы 3,98 Агонист

Соединения третьей группы

9 OH, CN 0,2 0,25 -5,3 По атому N CN группы 1,76 Агонист

10 OH, CN 20 0,18 -5,18 По атому N CN группы 1,76 Агонист

11 OH, CN 2,1 0,28 -5,37 По атому N CN группы 1,54 Агонист

12 CN 24 0,75 -5,32 По ат. N C=N-группы и по N 2,77 Агонист

13 Нет ОН, CN 0,01 0,32 - - 3,42 Антагонист

14 Нет ОН, есть NO2 0,5 0.59 -6,31 По ат. О NО2 группы 2,72 Антагонист

15 Нет ОН, есть NO2, N 1 0,57 -6,38 По ат. О NО2 группы 2,69 Антагонист

16 Нет ОН, есть NO2 3 0,55 По ат. О NО2 группы 2,88 Антагонист

Соединения четвертой группы

17 N,N 2 0,51 -5,68 По ат. N,N 3,7 Агонист

18 N,N 6 0,57 -5,67 По ат. N,N 3,4 Агонист

19 N,N 21 0,38 -5,74 По ат. N,N 3,9 Агонист

20 N,N 115 0,44 -5,63 По ат. N 3,2 Антагонист

21 N,N 358 0,56 -5,54 По ат. N 5,1 Антагонист

22 N,N 456 0,57 -5,56 По ат. N 4,6 Антагонист

Здесь Ki - константа ингибирования;

ДЕ - разность энергий уровней ВЗМО и предВЗМО в эВ; НСМО - нижняя свободная молекулярная орбиталь; lgP - коэффициент липофильности.

Таблица 2

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Названия лигандов КР.3С4-рецептора в таблице 1

№ Наименование

1 Дигидротестостерон

2 Propionic acid (8R,9S,10R,13S,14S,17S)-10,13-dimethyl-3-oxo-2,3,6,7,8,9,10,11, 12,13,14,15,16,17-tetradecahydro-1H-cydopenta[a]phenanthren-17-yl ester

3 5-Bromo-2-(2,2-dimethyl-6-methylene-cyclohexylmethyl)-4-methoxy-phenol

4 2,2,4-Trimethyl-5-(3-methyl-butyl)-2,5-dihydro-1H-6-oxa-1-aza-chrysene

Продолжение таблицы 2

5 (4aS,10aR)-4a-benzyl-7-hydroxy-3,4,4a,9,10,10a-hexahydrophenanthren-2(1H)-one::4a-Benzyl-7- hydroxy-3,4,4a,9,10,10a-hexahydro-1H-phenanthren-2-one

6 (2S)-3-(4-acetamidophenoxy)-2-hydroxy-2-methyl-N-[4-nitro-3-

ХИМИЧЕСКИЕ НАУКИ

(trifluoromethyl)phenyl]propanamide

7 (2S)-3-(4-fluorophenoxy)-2-hydroxy-2-methyl-N-[4-nitro-3-(trifluoromethyl)phenyl]propanamide

8 (2S)-3-(4-bromophenoxy)-2-hydroxy-2-methyl-N-[4-nitro-3-(trifluoromethyl)phenyl]propanamide

9 2-chloro-4-((3aR,4S)-4-hydroxy-3-methyl-1,1-dioxo-tetrahydro-1lambda6-thia-2,6a-diaza-pentalen- 2-yl)-3 -methyl-benzonitrile

10 2-chloro-4-((3aR,4S)-4-hydroxy-3-methyl-1,1-dioxo-tetrahydro-1lambda6-thia-2,6a-diaza-pentalen- 2-yl)-3 -methyl-benzonitrile

11 4-[(7R,7aS)-7-hydroxy-1,3-dioxo-hexahydro-1H-pyrrolo[1,2-a]imidazolidin-2-yl]-2-chloro-3- methylbenzonitrile

12 4-[(1R,6R,7S)-3,5-dioxo-2,4-diazatricydo[5.2.1.0A{2,6}]decan-4-yl]naphthalene-1-carbonitrile

13 (2R,6S)-4-(4-Bromo-3-methyl-phenyl)-4-aza-tricyclo[5.2.2.0*2,6*]undec-8-ene-3,5-dione

14 (2S,6R)-4-(4-Nitro-naphthalen-1-yl)-4-aza-tricyclo[5.2.1.0*2,6*]dec-8-ene-3,5-dione

15 (1R,6R,7S)-4-(4-nitromphthalen-1-yl)-2,4-diazatricydo[5.2.1.0A{2,6}]decane-3,5-dione

16 (3a alpha,4beta,8beta,8a alpha)-4,5,6,7,8,8a-hexahydro-4,8-etheno-1H-cydohepta[c]furan-1,3(3aH)- di-one

17 Cis-(4aR,12bS)-11-Trifluoromethyl-2,3,4,4a,5,6,8,12b-octahydro-1H-6,8-diaza-benzo[a]anthracen-9- one

18 (2R,3R)-2,3-Dimethyl-5-trifluoromethyl-1,2,3,8-tetrahydro-pyrrolo[3,2-g]quinolin-7-one

19 (5bR,9aR)-5b-Methyl-4-trifluoromethyl-1,5b,6,7,8,9,9a,10-octahydro-pyrido[2,3-b]carbazol-2-one

20 6,8,8-Trimethyl-4-trifluoromethyl-8,9-dihydro-1H-pyrido[3,2-g]quinolin-2-one

21 6,7,7,8,8-Pentamethyl-4-trifluoromethyl-6,7,8,9-tetrahydro-1H-pyrido[3,2-g]quinolin-2-one

22 (6R,7S)-6,7,8,8-Tetramethyl-4-trifluoromethyl-6,7,8,9-tetrahydro-1H-pyrido[3,2-g]quinolin-2-one

Соединения первой группы - гетероциклические, близкие по структуре к дигидротестостерону (К;=0,2 нМ) и содержащие группу С=О. Соединение № 2 (К;=0,25 нМ) гидролизуется и превращается в дигидротестостерон. Соединения № 3, 4 (Ki=1542 нМ, Ki=1485 нМ) протонируются по бензольному кольцу и поэтому являются слабыми агонистами.

Соединение № 5 (Ki=1300 нМ) является антагонистом. В его электронной структуре группа С=О не содержит на ВЗМО электронной плотности, поэтому протонирование и последующий перенос электрона невозможны.

Исходя из расчетных данных, прослеживается связь структура-активность по следующим параметрам: ДБ, lgP.

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

Вторая группа соединений представлена в экспериментальных данных Protein Data Base селективным модулятором C-23 и его аналогами № 6, 7, 8 (Ki=3,98 нМ; 6,11 нМ; 11,6 нМ) и содержит NO2-группу, которая может протонироваться. Во всех соединениях возможно протонирование по NO2-группе и последующее взаимодействие с NR3C4-рецептором с переносом электрона на агонист, что также не противоречит предлагаемому механизму действия. Исходя из расчетных данных, прослеживается связь структура -активность по параметру lgP.

Третья группа соединений представлена в экспериментальных данных Protein Data Base селективным модулятором S-24 и его аналогами № 9, 10, 11, 12 (Ki=0,2 нМ; 20 нМ; 2,1 нМ; 24 нМ (агонисты)); № 13, 14, 15, 16 (Ki=0,01 нМ; 0,5 нМ; 1 нМ; 3 нМ (антагонисты)) и содержит активную С^группу.

Анализ электронных и структурных параметров показал, что все агонисты содержат ОН-группу, а антагонисты не содержат, поэтому невозможны образование водородных связей антагонистов с рецептором и последующий перенос электрона. Агонист № 12 (Ki=24 нМ) может протонироваться по атому трехвалентного азота и затем образовывать

Auditorium: электронный научный журнал Курского государственного университета. 2014. № 1

Кузнецова Н. Б., Кузнецов П. Е. Квантово-химическое моделирование механизма действия лигандов МЯзС^рецептора

водородную связь с Asn^-рецептора. Но его связывание с рецептором более слабое, так как он не содержит группу ОН, которая может подстраиваться под рецептор, а содержит трехвалентный протонированный атом азота, встроенный в гетероциклическую структуру. Таким образом, более активные соединения содержат ОН-группу, а менее активные -протонированный атом азота. Кроме того, важна полярность соединений. Высокая полярность протонированной производной способствует правильному геометрическому расположению агонистов вблизи Arg82 и последующему протонированию. Кроме того, важна геометрия лиганда. Так, в соединении № 10 (К;=20 нМ) метильная группа будет стерически препятствовать взаимодействию с рецептором, так как она расположена рядом с ОН-группой. С другой стороны аналоги соединения bms и № 9 (К;=0,2 нМ) не содержат вблизи от ОН группы метильную группу, поэтому являются более активными.

Прослеживается связь структура-активность по следующим параметрам: lgP для агонистов, ДЕ - для агонистов и антагонистов.

Четвертая группа № 17-22 (К;=2 нМ, 6 нМ, 21 нМ, 115 нМ, 358 нМ,456 нМ) соединений не представлена в экспериментальных данных Protein Data Base и содержит трициклическую структуру и два активных атома азота, которые могут протонироваться.

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

Квантово-химическое исследование комплекса соединения № 17 (^=2 нМ) с NR3C4-рецептором показало, что с протонированием происходит перенос электрона с рецептора на трициклическую структуру соединения № 17. Данный процесс термодинамически разрешен.

Из табличных данных видно, что все активные агонисты содержат в своей структуре -ОН-группу; активные группы, по которым возможно протонирование -С=О, -CN, -NO2 и на которых содержится практически один электрон.

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

Из анализа механизма действия и структурных данных агонистов были сформулированы условия, предъявляемые к структурной формуле агонистов:

1) значения lgP находятся в интервале от 1,4 до 4,4;

2) наличие ОН группы для наиболее активных агонистов или NH группы - для менее активных;

3) наличие групп, которые могут протонироваться: -С=О, -CN, -NO2;

4) ДЕ < 0,5 эВ;

5) НСМО протонированной формы составляет [-5,1; ... -6,5] эВ;

ХИМИЧЕСКИЕ НАУКИ

6) расстояние между атомом, который протонируется, и атомом кислорода ОН-группы в протонированной форме составляет [7, ... 11] А;

7) наличие вакансий для электрона на НСМО уровне на атомах - центрах протонирования в протонированной форме;

8) размер лигандов, сопоставимый с активным центром рецептора;

9) жесткая геометрия и определенное взаиморасположение протонируемой и -ОН-групп, причем неподеленная электронная пара протонируемой группы должна находиться вблизи Arg82 и взаимодействовать с ним;

10) термодинамическая разрешенность процесса переноса электрона.

Механизм действия лигандов ККзС4-рецептора заключается в следующем.

1. Неподеленная электронная пара протонируемой группы лиганда взаимодействует с Arg82, заряженным положительно.

2. Одновременно происходит перенос протона с Arg82 на лиганд и ориентация протонированного лиганда в пространстве с образованием водородных связей с Asn35.

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

4. В результате переноса электрона происходит трансформация рецептора и его переход из цвиттер-иона в катионную форму.

5. Активные агонисты характеризуются наличием в структуре:

- ОН-группы,

- активных групп, по которым возможно протонирование и на ВЗМО которых содержится практически один электрон.

Таким образом, активная (трансформированная) форма NR^-рецептора представляет собой катион-радикал в комплексе с протонированным лигандом. Основным во взаимодействии лиганда с рецептором является протонирование. Если оно возможно, то лиганд - агонист, если невозможно, то антагонист. В случае агониста после протонирования происходит дальнейший перенос электрона с рецептора на лиганд.

Проведенные расчеты показывают, что протонирование играет ключевую роль в проявлении агонистической активности лигандов КК3С4-рецептора.

Библиографический список

Кузнецова Н.Б., Кузнецов П.Е., Согуренко И.А. Возможный механизм токсического действия экотоксикантов класса дибензо-п-диоксина // Вестник Саратовского государственного аграрного университета. 2006. № 2. С. 18-24

Щитовидная железа. Фундаментальные аспекты / под ред. проф. А.И. Кубарко и проф. S. Yamashita. Минск - Нагасаки, 1998. 366 с.

Discovery Studio Client v2.5.0.9164 [Сайт]. URL: http//www.accelrys.com/ (дата обращения: 20.06.2011).

Hyper Chem Release 7 for Windows. 2002. HyperCube, Inc. 860 p.

Kuznetsova N.B., Vlasov I.A., Kuznetsov P.E., Mihirev D.A. Quantum-chemical and Experimental Modelling of Hydrogen Peroxide Formation in Heterogeneous Water Surroundings/Milieus Which Contain Dibenzo-p-dioxines or Their Analogues // SPIE Proceedings, 2006. V. 5067. P. 288-294.

Auditorium: электронный научный журнал Курского государственного университета. 2014. № 1

Кузнецова Н. Б., Кузнецов П. Е. Квантово-химическое моделирование механизма действия лигандов МК3С4-рецептора

Protein Data Ба8е[Сайт]. URL: www.rcsb.org/pdb/home/home.do (дата обращения:

14.07.2013).

i Надоели баннеры? Вы всегда можете отключить рекламу.