УДК 577.338:577.123.36
ТЕОРЕТИЧЕСКОЕ ИССЛЕДОВАНИЕ АДЕНОЗИНА И ЕГО ИЗОСТЕРНЫХ АНАЛОГОВ. ВОЗМОЖНЫЙ МЕХАНИЗМ ИХ АКЦЕПТИРОВАНИЯ В АКТИВНОМ ЦЕНТРЕ АДЕНОЗИНДЕЗАМИНАЗЫ МЛЕКОПИТАЮЩИХ
© 2003 Ю.П.Зарубин, Д.В. Миргородский, П.П. Пурыгин,1 И.А.Ильичева, В.Л. Флорентьев2
Методами компьютерного моделирования были изучены аденозин и его изостерные аналоги, включая их Nl-протонированные формы (кроме 1-деазааналогов), с целью нахождения корреляций между их молекулярной структурой и субстратными свойствами для аденозиндезами-назы млекопитающих. Для этих соединений методом ab initio STO-3G были рассчитаны распределение зарядов на атомах и изоповерхности электростатического потенциала вблизи их ван-дер-Ваальсовых радиусов. Конформационные исследования были выполнены методом молекулярной механики ММ+. Обсужден механизм акцептирования, определяющий субстратную селективность аденозиндезаминазы млекопитающих, и предсказаны возможные субстратные свойства для ранее неизученных аналогов аденозина.
Введение
Поиск новых лекарственных средств, устойчивых к действию различных ферментов, в частности ферментов катаболизма пуринов и пирими-динов, является одной из актуальных проблем современной медицинской химии и фармакологии [1—3]. Особое место среди таких ферментов занимает аденозиндезаминаза (АДА, КФ 3.5.4.4), осуществляющая превращение (2'-дезокси)аденозина в (2'-дезокси)инозин путем гидролитического замещения NH2-группы на OH-группу в положении 6 остатка пурина. В результате этого различные фармакологически активные аналоги аденозина (Ado) превращаются в неактивные или малоактивные аналоги инозина [4]. АДА
1 Юpий Павлович Зарубин, Денис Виктоpович Миpгоpодский, Петp Петpович Пурыгин, кафедра органической химии Самаpского госудаpственного унивеpситета, 443011, г. Самара, ул. Акад. Павлова, 1.
2Ирина Алексеевна Ильичева, Владимир Леонидович Флорентьев, лаборатория химии белкового синтеза Института молекулярной биологии имени В.А. Энгельгардта РАН, 119991, г. Москва, ул. Вавилова, 32.
играет также исключительно важную роль в регуляции иммунитета у млекопитающих [5, 6].
К настоящему времени синтезировано значительное число аналогов Ado, часть из которых нашла применение в медицине и в биохимических исследованиях процессов метаболизма на уровне клеточных культур, тканей и органов [2]. Во многих случаях изучалась устойчивость этих соединений к ферментативному дезаминированию.
Однако существует немного работ, в которых бы исследовалась взаимосвязь между структурой, конформациями и субстратными свойствами для Ado и его аналогов в отношении АДА методами компьютерной химии [7-18]. Несмотря на наличие данных по рентгеноструктурному анализу (РСА) комплексов АДА-ингибитор [19-25], возможный механизм акцептирования субстрата или его аналогов в активном центре этого фермента в литературе до недавнего времени не был описан [14, 17, 18]. Поэтому представляется перспективным теоретическое изучение взаимосвязи "структура— субстратные и ингибиторные свойства" у различных аналогов Ado для АДА не только с целью описания возможного механизма акцептирования субстрата или его аналогов в активном центре данного фермента, но и предсказания возможных субстратных и ингибиторных свойств у плохо изученных или неизученных аналогов Ado.
Особое место среди них занимают изостерные аналоги, в которых имеются немодифицированный остаток ß-D-рибофуранозы и замещения атомов в различных положениях остатка исходного пуринового гетероцикла на атомы C и N, на CH- и NH-группы при сохранении локализации NH2-группы [26-54]. Для большинства этих соединений в литературе описаны субстратные и ингибиторные свойства в отношении АДА млекопитающих [26-42]. Так как стерические изменения структуры гетероцикла в молекулах этих соединений минимальны, а характер их взаимодействий с функциональными группами остатков аминокислот в активном центре АДА сильно различается, данные аналоги Ado являются удобными "химическими инструментами" для изучения механизма акцептирования субстрата и его аналогов в активном центре этого фермента.
Структурные формулы всех исследованных молекул Ado, его изостер-ных аналогов и их 1h+-катионов (кроме 1-деазапроизводных) приведены ниже на рис. 1 с пояснениями в табл. 1.
NH-
Z iL ü, HO- O4' V'^V
4'C У1'
31 T2' OH OH
Рис. 1. Обобщенная структурная формула молекул Ado, его изостерных аналогов и их 1^-катионов
Таблица 1
Расположение заместителей в остатке гетероцикла для Ado, его изостерных аналогов и их 1й+-катионов на рис. 1
Молекула Заместители в соответствующих положениях
нуклеозида остатка гетероцикла по пуриновой нумерации
Z1 Z2 Z3 Z4 Z5 Z7 Z8 Z9
Ado N CH N C C N: CH N
z2Ado N N: N C C N: CH N
z8Ado N CH N C C N: N: N
c1 Ado CH CH N C C N: CH N
c3Ado N CH CH C C N: CH N
c7Ado N CH N: C C N: CH N
c9Ado N CH N: C C NH CH C
z2c3Ado N N: CH C C N: CH N
zVAdo^ N CH N: C N N: CH C
z8c1 Ado CH CH N: C C N: N N
z^Ado^ N CH CH C C N: N N
z8c7Ado N CH N C C CH N N
7H-z8c9Ado N CH N C C NH N C
8H-z8c9Ado N CH N C C N: NH C
z^Ado^1 N N: N C C N: N: N
c1-3 Ado CH CH CH C C N: CH N
c1-' Ado CH CH N: C C CH CH N
c^'Ado^ N CH CH C C CH CH N
z2,sc3Ado N N: CH C C N: N: N
z2,sc'Ado N N: N: C C CH N: N
z^-^Ado^ N N: CH C C CH CH N
zV-'Ado CH CH N: C C CH N: N
c^-'Ado^ CH CH CH C C CH CH N
z4c3,yAdo N CH CH N C N: CH C
zV^Ado^ N CH N N C CH CH C
z^^Ado^ N CH N N C N: N C
z^c^Ado^1 CH CH N N C N: N C
z4,sc3,yAdo N: CH CH N C N: N C
z^c^Ado^1 CH CH CH N C N: N C
Обозначения:
zA — A-аза-, cB — B-деаза-, zA,C — А,С-диаза-, cB,D — Б,Б-дидеаза-, cB,D,E — Б,Б,Е-тридеаза-.
а) для данных соединений субстратные свойства в отношении АДА не исследовались;
б) Z1 = NH для соответствующих 1й+-катионов Ado (1A+-Ado) и его аналогов.
1. Расчетная часть
Исходные структуры молекул Ado, его изостерных аналогов (29 структур), а также их 1h+-катионов (21 структура) были построены из шаблонов структур остатка аденозин-5'-монофосфата в программном пакете HyperChem 7.0 [55] со стартовыми конформациями рибофуранозного цикла С2'-эндо (S-область) и С3'-эндо (N-область).
Заряды на атомах структур молекул нуклеозидов и их ^+-катионов рассчитаны неэмпирическим методом STO-3G в режиме single point.
Для определения наиболее устойчивых конформаций у изучаемых структур методом молекулярной механики MM+ использовались модифицированные скрипты [56], позволяющие проводить расчеты потенциальной энергии по трем торсионным углам — гликозидному O4'C1'N(C)9C(N)4 (х) в диапазоне от -180° до +180° с шагом 5° и 1° (для определения точных значений параметров энергетических минимумов) и двум экзоциклическим — O5'C5'C4'C3' (у) и HO5'O5'C5'C4' (ß) — со следующими стартовыми значениями: -60° (—синклинальная, —ск), +60° (+синклинальная, + ск) и 180° (антиперипланарная, ап). Эти скрипты отличаются от ранее использованных [17, 18] тем, что позволяют анализировать значительно большее число конформационных областей молекул нуклеозидов и их ^+-катионов, при этом ниже вероятность перехода рибофуранозного цикла из N - в S-область псевдовращения.
Расчет сродства к протону у оптимальных для акцептирования в активном центре АДА конформаций молекул нуклеозидов проведен полуэмпирическим методом AM1 [57] в режиме single point.
Расчет изоповерхностей электростатических потенциалов в устойчивых конформациях молекул нуклеозидов и их ^+-катионов, оптимальных соответственно для акцептирования и удерживания в активном центре АДА, проведен неэмпирическим методом STO-3G в режиме single point.
Определение конформационных характеристик всех изученных структур осуществлялось с помощью специально созданной нами программы NuclConf. Все условные обозначения конформаций нуклеозидов и расчетные формулы для определения их параметров взяты из [58, 59].
2. Результаты и их обсуждение
Расчет парциальных зарядов (по Малликену) на атомах молекул Ado и его различных изостерных аналогов и их ^+-катионов показывает, что протонирование атома азота N1 приводит к заметному изменению электронной плотности в остатке гетероцикла и мало влияет на ее распределение в остатке ß-D-рибофуранозы. Для всех ^+-катионов нуклеозидов наблюдается возрастание положительного заряда на атоме углерода C6 (с 0.25. . . 0.29 до 0.33... 0.38 дол. эл.). Это подтверждает экспериментальные данные о необходимости протонирования атома N1 гетероцикла для протекания ферментативной реакции дезаминирования [60, 61], когда в качестве нуклеофи-
ла выступает цинк-активированная молекула воды в виде гидроксид-иона [21]. Аномально высокие значения зарядов на атоме C6 в молекуле z5c9Ado (0.41 дол. эл.) и катионе 1h+-z5c9Ado (0.51 дол. эл.) объясняются наличием эндоциклических электроноакцепторных атомов азота в положениях 1, 3 и 5 триазинового фрагмента гетероцикла.
Изучение конформационных характеристик наиболее низкоэнергетических структур молекул нуклеозидов в N - и S-областях псевдовращения рибофуранозного цикла показало, что анти-конформация вокруг глико-зидной связи C1'-N(C)9 характерна для наиболее устойчивых минимумов большинства структур как в N-, так и в S-областях. Для ряда 8-азапроиз-водных Ado характерна устойчивая син-конформация вокруг гликозидной связи C1'-N(C)9 в случае нахождения рибофуранозного цикла в S-области. Объясняется это в основном отталкиванием отрицательных зарядов на атомах азота N8 гетероцикла и кислорода O4' рибофуранозы. При экранировании атома азота N8 водородом у 8H-z8c9Ado наблюдается классическая анти-конформация. S-конформеры молекул нуклеозидов более устойчивы, чем их соответствующие N-конформеры, они являются глобальными энергетическими минимумами соответствующих структур. Тип конформа-
ции вокруг экзоциклической связи C4—C5'--+ск, он не зависит от области
псевдовращения рибофуранозного цикла, при этом значения торсионного угла у лежат в очень узком диапазоне (~60°).
Анализ водородных связей в активном центре АДА для различных структур фермент-ингибиторных комплексов по данным РСА [20-22, 25] показал, что наиболее оптимальная ориентация 5'-ОН-группы в системе водородных связей между водородом у атома азота N1 имидазольной группы остатка His17 и одним из атомов кислорода COO--группы остатка Asp19 находится в области -ск-конформации торсионного угла ß. Но наиболее устойчивой областью конформации торсионного угла ß для молекул исследуемых нуклеозидов является + ск, что связано, по всей вероятности, с образованием водородной связи между атомами водорода HO5' экзоциклической 5'-OH-группы и кислорода O4' рибофуранозного цикла в случае + ск-конформации торсионного угла у.
В N-области для большинства молекул характерны близкие кон-формации рибофуранозного цикла группы C3'^Hdo (в основном типа C3'^Hdo-C4'^K3o). В S-области конформации рибофуранозного цикла относятся к C2'-эндо-гpуппе (в основном типа C2'^Hdo-C1'^K30).
Таким образом, среди исследуемых молекул нуклеозидов нет ни одной, которая имела бы в наиболее устойчивой конформации параметры, характерные для аналогов субстрата в активном центре АДА, известные из данных РСА [20-22, 25], а именно: х = -118... -101° (анти-конформации), у = = 40... 63° (+ск-конформация), ß = -80 ...-90° (-ск-конформация), конформация рибофуранозного цикла типа C3'-эндо [20-22] или C3'-эндо-C4'-экзо [25]. Поэтому конформационная область (х = -120... - 100°, у = 40... 65°) была дополнительно детально просканирована с целью нахождения для мо-
лекул нуклеозидов всех вариантов устойчивых структур, причем как в N-, так и в S-областях. Было найдено, что, независимо от области псевдовращения рибофуранозного цикла, значения торсионных углов х и у лежат в очень узких диапазонах и попадают в область указанных выше параметров. Это позволяет АДА акцептировать в активном центре АДА молекулу субстрата или его аналога как в N-, так и в S-областях. Рибонуклеозиды акцептируются в N-области с конформацией C3'-эндо [62]. Значения торсионного угла ß для данных устойчивых конформаций лежат в основном в +ск- и в ±ап-областях, что может дополнительно замедлять акцептирование аналога субстрата в активном центре АДА. Переход в C3'-эндо-[20-22] или C3'-эндо-C4'-экзо-конформацию [25] может происходить также в активном центре фермента. Рассчитанные конформационные напряжения в молекулах не превышают 5 ккал/моль относительно энергий их глобальных минимумов, принятых за 0.
Протонирование атома N1 гетероцикла в молекулах нуклеозидов приводит к существенным изменениям конформационных характеристик образующихся 1h+-катионов нуклеозидов в N - и S-областях. Однако в S-области типы их конформаций остаются в основном прежними. Для 8-азапроизвод-ных 1h+-Ado в S-области характерна устойчивая син-конформация вокруг гликозидной связи, подобно 8-азапроизводным Ado. Значения торсионного угла ß для данных устойчивых конформаций ^+-катионов нуклеозидов лежат в ±ап-областях, что объясняется ослаблением водородной связи между атомами водорода HO5' и кислорода O4' вследствие протонирования.
Конформационные параметры наиболее низкоэнергетических структур 1h+-катионов, оптимальных для удерживания в активном центре АДА, практически не отличаются от таковых для непротонированных молекул нуклеозидов. Для этих конформаций ^+-катионов нуклеозидов значения торсионного угла ß также лежат в ±ап-областях.
Значения амплитуды складчатости рибофуранозного цикла для всех молекул и ^+-катионов нуклеозидов лежат в узком диапазоне и близки к ранее полученным [11, 17, 18].
Рассчитанные величины сродства к протону (см. [63]) по положению N1 гетероцикла для молекул нуклеозидов (кроме 1-деазапроизводных) показывают (табл. 2), что их значения находятся в пределах от 199.9 ккал/моль (z2'8Ado) до 229.1 ккал/моль (c3'7Ado). С увеличением числа электроноак-цепторных атомов азота в гетероцикле сродство к протону снижается, а с увеличением числа электронодонорных CH-групп — возрастает. Для изомерных аналогов Ado величина сродства к протону по сравнению с Ado может как снизиться, так и возрасти.
Поэтому только конформационными особенностями и неодинаковым сродством к протону изученных аналогов Ado нельзя объяснить значительные различия в их субстратных свойствах в отношении АДА, известные из экспериментальных данных (табл. 3).
В связи с этим нами было выдвинуто предположение не только о важ-
Таблица 2
Значения сродства к протону для молекул нуклеозидов в областях конформаций, активных по отношению к АДА
Молекула нуклеозида Сродство к протону, ккал/моль
Ado 219.1
z2Ado 209.5
z8Ado 209.6
c1 Ado -
c3Ado 224.3
c'Ado 224.0
c9Ado 219.7
z2 c3Ado 219.2
z5 c9Ado 204.2
z8 c1 Ado -
z8 c3Ado 208.4
z8 c'Ado 215.3
7Я-z8cyAdo 213.9
8Я-z8c9Ado 219.4
z2,8Ado 199.9
c1-3Ado -
c1-7Ado -
c3-'Ado 229.1
z2,8c3Ado 204.2
z2-8 c'Ado 205.5
z2c3,7Ado 223.9
z8c1,7Ado -
c1-3-7Ado -
z4c3,yAdo 215.7
z4c'-9Ado 208.2
z4,8cyAdo 205.5
z4-8 c1-9Ado -
z4-8 c3-9Ado 202.4
z4-V-3-yAdo -
ной роли взаимодействий тех или иных структурных фрагментов в молекулах исследуемых изостерных аналогов Ado с различными функциональными группами остатков аминокислот в активном центре АДА, но и об определенном порядке их акцептирования. В качестве определяющего критерия была выбрана изоповерхность электростатического потенциала (ИЭП), так как изменение характера или отсутствие тех или иных водородных связей в системе фермент—лиганд должно приводить к различиям в субстратных и ингибиторных свойствах соответствующих нуклеозидов для АДА.
Анализ данных РСА (рис. 2), кинетических данных (табл. 3) и результатов расчетов ИЭП для молекул Ado и его изостерных аналогов показывает, что для оптимального акцептирования молекулы нуклеозида функциональными группами остатков аминокислот в активном центре фермента необ-
Таблица 3
Кинетические параметры дезаминирования (Km, Vmax) и ингибирования (Kj) для Ado и его изостерных аналогов в отношении АДА
Соединение Km, мкМ Kj, мкМ Vnax (отн.), % Литература
Ado 67 - - [29]
29 - 100 [26]
25 - 100 [30]
58 - 100 [31]
33 - 100 [32]
20 - 100 [33]
z2Ado 670 н/и 102 [26]
z8Ado 130 н/и 310 [30]
96 145 217 [34]
c1 Ado - 2.0 - [29]
- 0.662 - [35]
- 0.66 - [36]
- 0.18 - [37]
c3Ado - н/д - [29]
- 359 - [35]
- 360 - [36]
н/д 257 н/д [37]
c'Ado - - - [29, 31, 35, 36]
c9Ado - - - [38]
z2c3Ado н/д н/д н/д [28]
z8c1 Ado - 0.903 - [27]
z8c7Ado 250 н/и 6.4 [32]
125 н/и 6.4 [33]
z8c9Ado 1000 н/и 750-850 [30]
c1-3 Ado - 110 - [35, 36]
c1-7 Ado - - - [35, 36, 39]
z2-8c3Ado н/д н/д н/д [28]
z2-8c7Ado 440 н/и 0.2 [26]
z8c1,7Ado - - - [40]
z4c3-9Ado - ? - [41]
z4,8c3,9Ado н/д - н/д [42]
Примечания:
н/д — нет данных;
н/и — не исследовалось;
? — не установлено достоверно наличие или отсутствие ингибиторной активности для данного соединения
ходимо наличие областей отрицательного электростатического потенциала в положениях 1, 3 и 7 остатка гетероцикла (атомы азота N1, N3 и N7) и положениях 3' и 5' остатка -рибофуранозы (атомы кислорода 03' и 05'), а также областей положительного электростатического потенциала у атомов водорода 5'-0Н-группы (Н05') и 3'-0Н-группы (Н03').
Снижение величин парциальных отрицательных зарядов на атомах азо-
His238
"В
Asp295
ге
Asp296
C
HO/4O
His17
Or
Glu217—CJ3 O
HN 6
VN
O
e,'| -o^= C4
H2O
Asp19
His238
Asp295.
\
C0 O
O
V в/
. Hp NH-
Ж
OH
Glu217—C
у
Asp296
HO O
His17
H2O
Asp19
Рис. 2. Схема взаимодействий по данным РСА 6(£ )-гидрокси-1,6-дигидро-пуринрибозида [20] (слева) и с'А^ [21] (справа) с функциональными группами остатков аминокислот в активном центре АДА. Связи иона цинка ^п2+) показаны сплошными линиями, водородные связи — пунктирными
H
та N1, N3 или N7 вследствие электроноакцепторного действия атомов азота в соседних положениях гетероцикла (табл. 4) приводит к ухудшению субстратных свойств, а нахождение атомов азота или CH-групп в неканонических положениях гетероцикла (1-5, 7 и 8) приводит также к появлению ингибиторных свойств у изостерных аналогов Ado (табл. 3).
Молекулы нуклеозидов со сходным распределением областей электростатических потенциалов в гетероцикле проявляют или должны проявлять сходные субстратные и ингибиторные свойства для АДА.
Ключевую роль в акцептировании молекулы субстрата играет образование водородной связи между атомами азота N7 гетероцикла в субстрате и водорода COOH-группы остатка Asp296 в активном центре АДА. Поэтому c7Ado [29, 31, 35, 36], c9Ado [37] и c1-7Ado [35, 36, 39] не являются ни субстратами, ни ингибиторами для данного фермента (табл. 1). Можно предположить, что c3-7Ado [46], z2c3'7Ado [47], c1-3-7Ado [49] и z4c7-9Ado [50, 51], подобно c7Ado, c9Ado и c1,7Ado, также не должны быть ни субстратами, ни ингибиторами для АДА.
Не менее важную роль в процессе акцептирования субстрата играют взаимодействия остатка ß-D-рибофуранозы в нуклеозидах с функциональными группами остатков аминокислот в активном центре АДА (рис. 2). По литературным данным, ни 5'-дезоксиаденозин [64, 65], ни 5'-O-метиладенозин [31] не являются субстратами или ингибиторами для АДА, а 3'-дезоксиаденозин дезаминируется немногим медленнее, чем Ado или 2'-дезоксиаденозин [66]. Из этого следует, что 5'-ОН-группа в остатке ß-D-рибофуранозы играет такую же важную роль в акцептировании субстрата АДА, как и атом азота N7 в остатке гетероцикла.
Таким образом, процесс начального акцептирования субстрата в активном центре фермента носит кооперативный характер с участием функциональных групп аминокислотных остатков His17, Asp19 и Asp296 (рис. 2).
Таблица 4
Значения парциальных зарядов на атомах азота N1, N3 и N7 гетероцикла для молекул Ado и его некоторых изостерных аналогов
Молекула Заряды на Заряды на Заряды на
нуклеозида атоме N1, атоме N3, атоме N7,
дол. эл. дол. эл. дол. эл.
Ado -0.34 -0.32 -0.28
z2Ado -0.22 -0.20 -0.27
z8Ado -0.34 -0.32 -0.18
c1 Ado - -0.33 -0.28
c3Ado -0.34 - -0.28
z2c3Ado -0.24 - -0.27
z5c9Ado -0.34 -0.30 -0.22
z8c1 Ado - -0.32 -0.19
z8c3Ado -0.34 - -0.18
8H-z8c9Ado -0.35 -0.32 -0.18
z2'8Ado -0.22 -0.19 -0.17
z2'8c3Ado -0.23 - -0.18
z2'VAdo -0.22 -0.20 -
z4'8c9Ado -0.34 -0.19 -0.19
z4'8c1'9Ado - -0.19 -0.20
z4'8c1'3'9Ado - - -0.21
3'-0Н-Группа играет второстепенную роль, хотя она и важна для акцептирования субстрата в активном центре АДА.
Далее должно происходить последовательное образование двух водородных связей: между атомами азота N1 гетероцикла и водорода COOH-груп-пы остатка Glu217 и между атомами азота N3 гетероцикла и водорода CONH-группы остатка Gly184. В случае c1Ado и в образовании водородной связи участвует атом N3 гетероцикла, а в случае c3Ado — атом N1 гетероцикла. Молекула c1Ado изоэлектронна и изоструктурна катиону 1h+-Ado, а молекула c3 Ado неспособна протонироваться по атому азота N1 из-за увеличения расстояния между ним и атомом водорода COOH-группы остатка Glu217 (рис. 3). Поэтому c1 Ado — более сильный конкурентный ингибитор, чем c3Ado3. Оба соединения не являются субстратами для АДА (табл. 3).
Подобно c1Ado и c3Ado, c1'3Ado также не субстрат для АДА, но в то же время он ингибитор, причем более сильный, чем c3Ado, и более слабый, чем c1Ado (табл. 3). Неожиданное усиление ингибиторных свойств c1,3Ado по сравнению с c3Ado объясняется, вероятно, тем, что молекула c1'3Ado структурно сходна с катионом 1h+-Ado, несмотря на замену в положении 1 гетероцикла атома азота на CH-группу и отсутствие положительного заряда. Это подтверждает ранее выдвинутое предположение о том, что атом азота N7 гетероцикла является ключевым для акцептирования
3В [37], однако, сообщается, что c3Ado является субстратом для АДА, но очень
плохим, нет данных по Km и Vmax.
Г
-Asp296
Glu217—C
4
Л
Gy184
-OH
Рис. 3. Предполагаемая система водородных связей молекул c3Ado (X = C,
Y = CH, Z = N), z4c3-9Ado (X = N, Y = CH, Z = C) и z4-8c3-9Ado (X = Y =
= N, Z = C) с функциональными группами остатков аминокислот в активном центре АДА
субстрата или его аналога активным центром данного фермента. В молекуле z8c1 Ado атом азота N8 влияет на распределение парциальных зарядов в 1,2,3-триазольной части гетероцикла, что приводит к ослаблению электростатического взаимодействия между атомами азота N7 гетероцикла и водорода COOH-группы остатка Asp296. В результате акцептирование молекулы z8c1Ado протекает медленнее, что приводит к снижению ингибитор-ной способности z8c1Ado по сравнению с c1Ado (табл. 1). z4'8c3'9Ado является очень плохим субстратом для АДА и не ингибирует дезаминирование Ado. Только взаимодействие z4-8c3-9Ado со 100-кратным количеством фермента по сравнению со стандартными условиями эксперимента привело к 70% превращению z4'8c3'9Ado в 4,8-диаза-3,9-дидеазаинозин через 13 ч [42]. Поэтому данные по дезаминированию c3Ado [37], скорее всего,—реально протекающая реакция, а не артефакт. Для z4c3-9Ado наблюдалась картина, аналогичная для z4'8c3'9Ado при стандартной концентрации АДА [41]. Экспериментально наблюдаемое отсутствие ингибиторных свойств у z4c3-9Ado и z4'8c3'9Ado объясняется, по всей вероятности, влиянием электроноакцеп-торного атома азота в положении 4 (мета-положение относительно атома C6) остатка гетероцикла.
У z2Ado и z8Ado в соответствующих положениях имеются электроноак-цепторные атомы азота вместо CH-групп, но их субстратные свойства резко различаются. В молекуле z2Ado атом азота N2 в мета-положении к атому углерода C6 сильно снижает способность гетероцикла к дезаминированию по механизму SnAi\ Поэтому z2Ado дезаминируется намного медленнее, чем Ado (табл. 3).
В молекуле z8Ado атом азота N8 влияет на распределение парциальных зарядов в 1,2,3-триазольной части гетероцикла. В результате акцептирование молекулы z8Ado протекает медленнее, что приводит к уменьшению скорости ферментативного дезаминирования z8Ado по сравнению с Ado в одинаковых условиях V (отн.) при увеличении Km и Vmax (отн.) (табл. 3).
В отличие от c7Ado, z8c7Ado или z2'8c7Ado дезаминируются АДА, хо-
тя и с меньшей скоростью, чем Ado. В данном случае роль атома азота N7 в Ado выполняет атом азота N8 в z8c7Ado или z2'8c7Ado. Начальное акцептирование молекул z8c7Ado или z2'8c7Ado активным центром фермента носит, по-видимому, более сложный характер (рис. 4). Это приводит к значительному увеличению Km и резкому снижению Vmax для z8c7Ado по сравнению с Ado (табл. 1).
,.Asp296 С 180° '
NH.
O^O I
H H
/
Г T/N
nh2
OH
N ■ Ii x*.
O
\
H
/
/ /
,N—--H- O
N
Asp296
С
OH
ОН ОН он он
Рис. 4. Возможный механизм образования водородной связи между атомами азота N8 в молекулах z8c7Ado (X = СН), z2'8c7Ado (X = М) и водорода в-СООН-группы остатка Авр296 в активном центре АДА на начальной стадии акцептирования
Для таутомеров 7Н^8е9Аёо и 8Н^8с9Лёо характер начальной стадии акцептирования будет различным. В случае 8И-Т8е9Аёо она должна протекать аналогично Аёо, но медленнее из-за наличия в положении 8 гете-роцикла электроноакцепторного атома азота N8. Так как в водном растворе имеются оба таутомера z8c9Ado [67], то возможно акцептирование также и 7H-z8c9Ado. Начальная стадия этого процесса должна сопровождаться меньшими конформационными изменениями бокового радикала остатка Авр296 в активном центре АДА (рис. 5) по сравнению с z8c7Ado или z2'8с^о (рис. 4).
Asp296 180°
Г
Asp296
O
N"
<f\
H H
C
ОН он он он
Рис. 5. Возможный механизм образования водородной связи между атомом водорода Н7 в молекуле 7H-z8c9Ado и карбонильным атомом кислорода в-СООН-группы остатка Азр296 в активном центре АДА на начальной стадии акцептирования
Это приводит к значительному увеличению Кт и Vmax для z8c9Ado по сравнению с Ado (табл. 3). При этом скорость дезаминирования V для
N
Ado и его изостерных аналогов при одинаковых молярных концентрациях будет уменьшаться в такой последовательности: V(Ado) > V(z8c9Ado) > > V(z8c7Ado).
В отличие от c3Ado, z2c3Ado и z2'8c3Ado являются очень медленно реагирующими субстратами и слабыми ингибиторами АДА [28]4. Сродство к протону у атома азота N1 в z2c3Ado и z2'8c3Ado ниже, чем у c3Ado. Наличие электроноакцепторного атома азота в положении 2 гетероцикла в молекулах z2c3Ado и z2'8c3Ado вместо CH-группы должно было бы привести к ослаблению их ингибиторных свойств по сравнению с c3Ado и не способствовать дезаминированию. Однако z2c3Ado и z2'8c3Ado очень медленно дезаминируются, что можно объяснить только образованием водородной связи между атомами азота N2 в z2c3Ado и z2'8c3Ado и водорода CONH-группы остатка Gly184 в активном центре АДА (рис. 6).
^Asp296
оЧ
I
nh2 H
Glu217—C
O_H'" I .OH
» I >
■H^
N oh oh
Gly184
Рис. 6. Предполагаемая система водородных связей молекул z2c3Ado (X = CH) и z2-8c3Ado (X = N) с функциональными группами остатков аминокислот в активном центре АДА
Она значительно слабее, чем в случае наличия атома азота N3, но тем не менее способствует, хотя и в очень малой степени, переносу протона на атом азота N1 гетероцикла z2c3Ado или z2'8c3Ado от COOH-группы остатка Glu217, что необходимо для протекания ферментативной реакции дезами-нирования.
Подобно z8c7Ado, z2-8c7Ado также дезаминируется АДА, но значительно медленнее. При этом взаимодействие с функциональными группами остатков аминокислот в активном центре фермента для атомов азота N1 и N3 в 1,2,3-триазиновой части гетероцикла подобно z2Ado, а для атома азота N8 в пиразольной части гетероцикла — z8c7Ado (рис. 4). Это значительно увеличивает Km и резко снижает Vmax для z2-8c7Ado по сравнению с z8c7Ado (табл. 3).
Как и c1,7Ado, z8c1,7Ado не является не только субстратом, но и ингибитором для АДА (табл. 3). Это представляется парадоксальным, так как структурно родственные соединения z8c1Ado и z8c7Ado — инги-
4Для z2c3Ado Кт и К имеют очень большие значения, а Утах очень мала, поэтому
данные отсутствуют.
битор и субстрат соответственно. Но если предположить, что в случае 8-аза-7-деазапроизводных Ado образование двух водородных связей — между атомами азота N8 субстрата и водорода COOH-группы остатка Asp296 и между атомами азота N1 субстрата и водорода COOH-группы остатка Glu217 — кооперативный процесс, то замена атома азота в положении 1 ге-тероцикла на CH-группу приводит к полной потере и субстратной, и ин-гибиторной активности z8c1'7Ado для данного фермента. По этой причине изоструктурность молекулы z8c17Ado катиону 1h+-z8c7Ado уже не играет той роли, как и в случае других 1-деазапроизводных Ado.
Анализ распределения ИЭП [14, 17, 18] и изменения величин парциальных отрицательных зарядов на атомах азота N1, N3 или N7 в молекулах изостерных аналогов Ado (табл. 4) с изученными и неизученными субстратными (Km, Vmax и V (отн.))и ингибиторными (Ki) свойствами для АДА позволил предложить следующие ряды возможных изменений этих свойств у различных нуклеозидов.
Km (Ado) « Km (z5 c9 Ado) < Km(z8c9Ado);
Vmax(Ado) « Vmax(z5c9 Ado) < Vmax(z8c9Ado);
V(Ado) > V (z5c9 Ado) > V(z8c9Ado).
Km (z8Ado) « Km (z2Ado) < Km(z2'8Ado);
V(z8Ado) >> V(z2Ado) > V (z2'8Ado).
Vmax(z8Ado) >> Vmax(z2'8Ado) > Vmax(z2Ado).
Km(z2c3Ado) < Km(z2'8c3Ado);
Ki(z2c3Ado) < Ki(z2'8c3Ado);
V (z2c3Ado) > V(z2'8c3Ado).
Km(z8Ado) < Km(z4'8c9Ado);
V(z8Ado) > V (z4'8c9Ado).
Ki(c1 Ado) < Ki(z8c1Ado) < Ki(c1'3Ado) < Ki(z4'8c9Ado) <
< Km(z4'8c1'3'9Ado).
Km(c3Ado) < Km(z8c3Ado) < Km(z4c3'9Ado) < Km(z4'8c3'9Ado);
Ki(c3Ado) < Ki(z8c3Ado) < Ki(z4c3'9Ado) < Ki(z4'8c3'9Ado);
V(c3Ado) > V(z8c3Ado) > V(z4c3'9Ado) > V(z4'8c3'9Ado).
Для катионов 1h+-Ado и его изостерных аналогов области положительного электростатического потенциала в положениях 1 и 6 гетероцикла сходны с таковыми у молекул c1Ado [14, 17], z8c1Ado [18], c1'3Ado, z4'8c1'3'9Ado и z4'8c1'9Ado. Можно предположить, что молекулы c1Ado, z8c1Ado и c1'3Ado проявляют ингибиторные свойства в отношении АДА не только в силу структурного сходства с катионом 1h+-Ado, но и сходства ИЭП, не имея при этом положительного заряда.
На основании литературных данных по изучению субстратных и инги-биторных свойств для АДА различных изостерных аналогов Ado и результатов квантово-химических и молекулярно-механических расчетов, можно
предположить, что процесс акцептирования молекулы Ado или его изостер-ного аналога функциональными группами остатков аминокислот в активном центре АДА включает следующие стадии.
Стадия 1 — переход молекулы субстрата в активное конформационное состояние со следующими параметрами:
X = -120...-100° (анти-конформация);
у = 40... 65° (+ск-конформация);
ß = -80... -90° (—ск-конформация);
равновесная конформация рибофуранозного цикла типа C3'-эндо.
Стадия 2 — кооперативное образование трех водородных связей: между атомами кислорода O5' субстрата и водорода при атоме N1 ими-дазольной группы остатка His17;
между атомом водорода HO5' субстрата и атомами кислорода COO--группы остатка Asp19;
между атомами азота N7 субстрата и водорода COOH-группы остатка Asp296.
В случае наличия у молекулы субстрата в положении 7 остатка гете-роцикла CH-группы, а в положении 8 — атома азота, образование водородной связи происходит между атомами азота N8 субстрата и водорода COOH-группы после конформационных изменений бокового радикала остатка Asp296 в активном центре АДА.
Стадия 3 — образование водородной связи между атомами азота N1 субстрата и водорода COOH-группы остатка Glu217.
В случае 8-аза-7-деазапроизводных Ado образование двух водородных связей на стадиях 2 и 3 в активном центре АДА (между атомами азота N8 субстрата и водорода COOH-группы остатка Asp296 и между атомами азота N1 субстрата и водорода COOH-группы остатка Glu217) является кооперативным процессом.
Стадия 4 — образование водородной связи между атомами азота N3 субстрата и водорода CONH-группы остатка Gly184.
Если у молекулы в положении 3 остатка гетероцикла находится CH-группа, а в положении 2 — атом азота, то происходит образование очень слабой водородной связи между атомами азота N2 субстрата и водорода CONH-группы остатка Gly184 в активном центре АДА вследствие их взаимного смещения (рис. 6).
Стадия 5 — протонирование атома азота N1 гетероцикла в молекуле путем переноса протона от COOH-группы остатка Glu217.
Если в положениях 2 и 3 гетероцикла находятся CH-группы, то в этом случае перенос протона крайне затруднен по причине неоптимальной геометрии соответствующей водородной связи, образующейся на стадии 3.
В случае аденозина или 2'-дезоксиаденозина такая система образующихся водородных связей субстрата в активном центре АДА максимально способствует нуклеофильной атаке цинк-активированной молекулы воды (в ви-
де гидроксид-иона) на атом углерода C6 гетероцикла и, следовательно, протеканию ферментативной реакции дезаминирования по механизму SNAr.
Работа выполнена при поддержке РФФИ (проект 03-04-48938).
Литература
[1] Glazer R.I. Adenosine deaminase inhibitors: their role in chemotherapy and immunosuppression // Cancer Chemother. and Pharmacol. 1980. V. 4. No. 4. P. 227-235.
[2] Agarwal R.P. Inhibitors of adenosine deaminase // Pharmacol. and Ther. 1982. V. 17. No. 3. P. 399-429.
[3] Cristalli G. Costanzi S., Lambertucci C., Lupidi G., Vittori S., Volpini R., Camaioni E. Adenosine deaminase: functional implications and different classes of inhibitors // Med. Res. Rev. 2001. V. 21. No. 2. P. 105-128.
[4] Montgomery J.A. Studies on the biologic activity of purine and pyrimidine analogs // Med. Res. Rev. 1982. V. 2. No. 3. P. 271-308.
[5] Resta R., Thompson L.F. SCID: the role of adenosine deaminase deficiency // Immunol. Today. 1997. V. 18. No. 8. P. 371-374.
[6] Franco R., Valenzuela A., Lluis C., Blanco J. Enzymatic and extraenzymatic role of ecto-adenosine deaminase in lymphocytes // Immunol. Rev. 1998. V. 161. P. 27-42.
[7] Orozco M., Canela E.I., Franco R. Theoretical study of the protonation and tautomerization of adenosine, formycin, and their 2-NH2 and 2-F derivatives: functional implications in the mechanism of reaction of adenosine deaminase // Mol. Pharmacol. 1989. V. 35. No. 2. P. 257-264.
[8] Orozco M., Lluis C., Mallol J., Canela E.I., Franco R. Quantum chemical study of the electronic and conformational characteristics of adenosine and 8-substituted derivatives: functional implications in the mechanism of reaction of adenosine deaminase //J. Pharm. Sci. 1990. V. 79. No. 2. P. 133-137.
[9] Orozco M., Canela E.I., Franco R. A quantum chemical study of the enzymatic deamination of benzoadenine derivatives. A theoretical model of the interactions occurring between nucleosides and the active site of adenosine deaminase // Eur. J. Biochem. 1990. V. 188. No. 1. P. 155-163.
[10] Orozco M., Canela E.I., Franco R. Theoretical study of the hydroxyl nucleophilic attack on the 6-aminopyrimidine molecule: functional implications in the reaction mechanism of nucleoside deaminative enzymes //J. Org. Chem. 1990. V. 55. No. 9. P. 2630-2637.
[11] Orozco M., Velasco D., Canela E.I., Franco R. Determination of the conformational preferences of adenosine at the active site of adenosine deaminase // J. Amer. Chem. Soc. 1990. V. 112. No. 23. P. 8221-8229.
[12] Hansen L.M., Kollman P.A. Free energy perturbation calculations on models of active sites: Applications to adenosine deaminase inhibitors // J. Comput. Chem. 1990. V. 11. No. 8. P. 994-1002.
[13] Marrone T.J., Straatsma T.P., Briggs J.M., Wilson D.K., Quiocho F.A., McCammon J.A. Theoretical study of inhibition of adenosine deaminase by (8^)-coformycin and (8^)-deoxycoformycin //J. Med. Chem. 1996. V. 39. No. 1. P. 277-284.
[14] Пурыгин П.П., Зарубин Ю.П., Ильичева И.А., Флорентьев В.Л. Кван-тово-химическое исследование аденозина и его аза- и деазаанало-гов как антагонистов и ингибиторов аденозиндезаминазы млекопитающих // Вестник Самарского гос. ун-та. 1999. №4(12). С. 111-135.
[15] Ford H., Jr., Dai F., Mu L., Siddiqui M.A., Nicklaus M.C., Anderson L., Marquez V.E., Barchi J.J., Jr. Adenosine deaminase prefers a distinct sugar ring conformation for binding and catalysis: kinetic and structural studies // Biochemistry. 2000. V. 39. No. 10. P. 2581-2592.
[16] Bojack G., Earnshaw C.G., Klein R., Lindell S.D., Lowinski C., Preuss R. Design and synthesis of inhibitors of adenosine and AMP deaminases // Org. Lett. 2001. V. 3. No. 6. P. 839-842.
[17] Зарубин Ю.П., Ильичева И.А., Пурыгин П.П., Флорентьев В.Л. Теоретическое исследование антагонистов и ингибиторов аденозиндезамина-зы млекопитающих. I. Аденозин и его аза- и деазааналоги // Биоорган. химия. 2002. Т. 28. №4. С. 315-323.
[18] Зарубин Ю.П., Ильичева И.А., Пурыгин П.П., Флорентьев В.Л. Теоретическое исследование антагонистов и ингибиторов аденозиндезамина-зы млекопитающих. II. Изомерные аза-деазааналоги аденозина // Биоорган. химия. 2002. Т. 28. №5. С. 447-454.
[19] Wilson D.K., Rudolph F.B., Quiocho F.A. Atomic structure of adenosine deaminase complexed with a transition-state analog: understanding catalysis and immunodeficiency mutations // Science. 1991. V. 252. No. 5010. P. 1278-1284.
[20] Sharff A.J., Wilson D.K., Chang Z., Quiocho F.A. Refined 2.5 Е structure of murine adenosine deaminase at pH 6.0 //J. Mol. Biol. 1992. V. 226. No. 4. P. 917-921.
[21] Wilson D.K., Quiocho F.A. A pre-transition-state mimic of an enzyme: X-ray structure of adenosine deaminase with bound 1-deazaadenosine and zinc-activated water // Biochemistry. 1993. V. 32. No. 7. P. 1689-1694.
[22] Wilson D.K., Quiocho F.A. Crystallographic observation of a trapped tetrahedral intermediate in a metalloenzyme // Nat. Struct. Biol. 1994. V. 1. No. 10. P. 691-694.
[23] Sideraki V., Mohamedali K.A., Wilson D.K., Chang Z., Kellems R.E., Quiocho F.A., Rudolph F.B. Probing the functional role of two conserved active site aspartates in mouse adenosine deaminase // Biochemistry. 1996. V. 35. No. 24. P. 7862-7872.
[24] Sideraki V., Wilson D.K., Kurz L.C., Quiocho F.A., Rudolph F.B. Site-directed mutagenesis of histidine 238 in mouse adenosine deaminase: substitution of histidine 238 does not impede hydroxylate formation // Biochemistry. 1996. V. 35. No. 47. P. 15019-15028.
[25] Wang Z., Quiocho F.A. Complexes of adenosine deaminase with two potent inhibitors: X-ray structures in four independent molecules at pH of maximum activity // Biochemistry. 1998. V. 37. No. 23. P. 8314-8324.
[26] Bennett L.L., Jr, Allan P.W., Carpenter J.W., Hill D.L. Nucleosides of 2-aza-purines — cytotoxicities and activities as substrates for enzymes metabolizing purine nucleosides // Biochem. Pharmacol. 1976. V. 25. No. 5. P. 517-521.
[27] Franchetti P., Cappellacci L., Grifantini M., Lupidi G., Nocenti-ni G., Barzi A. 8-Aza analogues of deaza purine nucleosides. Synthesis and biological evaluation of 8-aza-1-deazaadenosine and 2'-deoxy-8-aza-1-deazaadenosine // Nucleosides and Nucleotides. 1992. V. 11. No. 5. P. 1059-1076.
[28] Bussolari J.C., Ramesh K., Stoeckler J.D., Chen S.F., Panzica R.P. Synthesis and biological evaluation of N4-substituted imidazo- and v-triazolo-[4,5-d]pyridazine nucleosides // J. Med. Chem. 1993. V. 36. No. 25. P. 4113-4120.
[29] Ikehara M., Fukui T. Studies of nucleosides and nucleotides. LXIII. Deamination of adenosine analogs with calf intestine adenosine deaminase // Biochim. Biophys. Acta. 1974. V. 338. No. 2. P. 512-519.
[30] Agarwal R.P., Sagar S.M., Parks R.E., Jr. Adenosine deaminase from human erythrocytes: purification and effects of adenosine analogs // Biochem. Pharmacol. 1975. V. 24. No. 6. P. 693-701.
[31] Krajewska E., De Clercq E., Shugar D. Nucleoside-catabolizing enzyme activities in primary rabbit kidney cells and human skin fibroblasts // Biochem. Pharmacol. 1978. V. 27. No. 10. P. 1421-1426.
[32] Bennett L.L., Jr, Allan P.W., Smithers D., Vail M.H. Resistance to 4-aminopyrazolo (3,4-d) pyrimidine // Biochem. Pharmacol. 1969. V. 18. No. 4. P. 725-740.
[33] Hecht S.M., Frye R.B., Werner D., Fukui T., Hawrelak S.D. Synthesis and biological activity of pyrazolo[3,4-d]pyrimidine nucleosides and nucleotides related to tubercidin, toyocamycin, and sangivamycin // Biochemistry. 1976. V. 15. No. 5. P. 1005-1015.
[34] Simon L.N., Bauer R.J., Tolman R.L., Robins R.K. Calf intestine adenosine deaminase. Substrate specificity // Biochemistry. 1970. V. 9. No. 3. P. 573-577.
[35] Lupidi G., Riva F., Cristalli G., Grifantini M. Inhibition of adenosine deaminase by deaza derivatives of adenosine and purine riboside // Ital. J. Biochem. 1982. V. 31. No. 6. P. 396-403.
[36] Lupidi G., Cristalli G., Marmocchi F., Riva F., Grifantini M. Inhibition of adenosine deaminase from several sources by deaza derivatives of adenosine and EHNA // J. Enzyme. Inhib. 1985. V. 1. No. 1. P. 67-75.
170
fö.n. 3apy6uH, ff.B. MnpropogCKHH, n.n.nypbiruH h gp.
[37] Kurz L.C., Moix L., Riley M.C., Frieden C. The rate of formation of transition-state analogues in the active site of adenosine deaminase is encounter-controlled: implications for the mechanism // Biochemistry. 1992. V. 31. No. 1. P. 39-48.
[38] Zimmerman T.P., Deeprose R.D., Wolberg G., Stopford C.R., Duncan G.S., Miller W.H., Miller R.L., Lim M.-I., Ren W.-Y., Klein R.S. Inhibition of lymphocyte function by 9-deazaadenosine // Biochem. Pharmacol. 1983. V. 32. No. 7. P. 1211-1217.
[39] Cristalli G., Vittori S., Eleuteri A., Volpini R., Cola D., Camaioni E., Gariboldi P.V., Lupidi G. Synthesis of 1,7-dideazapurine ribonucleosides and deoxyribonucleosides // Nucleos. Nucleot. 1993. V. 12. No. 1. P. 39-53.
[40] Sanghvi Y.S., Larson S.B., Willis R.C., Robins R.K., Revankar G.R. Synthesis and biological evaluation of certain C-4 substituted pyrazolo[3,4-b]pyridine nucleosides // J. Med. Chem. 1989. V. 32. No. 5. P. 945-951.
[41] MacCoss M., Meurer L.C., Hoogsteen K., Springer J.P., Koo G., Peterson L.B., Tolman L., Emini E. Synthesis and biological evaluation of nucleosides containing 8-aminoimidazo[1,2-a]pyrazine as an isosteric replacement for adenine // J. Heterocycl. Chem. 1993. V. 30. No. 5. P. 1213-1220.
[42] Schneller S.W., Thompson R.D., Cory J.G., Olsson R.A., De Clercq E., Kim I.-K., Chiang P.K. Biological activity and a modified synthesis of 8-amino-3-ß-D-ribofuranosyl-1,2,4-triazolo[4,3-a]pyrazine, an isomer of formycin //J. Med. Chem. 1984. V. 27. No. 7. P. 924-928.
[43] May J.A., Jr., Townsend L.B. Synthesis of v-triazolo(4,5-c)pyridine nucleosides and 4-(ß-D-ribofuranosyl)amino-1,2,3-thiadiazolo(5,4-b)pyridine via a rearrangement //J. Org. Chem. 1976. V. 41. No. 8. P. 1449-1456.
[44] Franchetti P., Messini L., Cappellacci L., Grifantini M., Nocentini G., Guarracino P., Marongiu M.E., La Colla P. 8-Aza derivatives of
3-deazapurine nucleosides. Synthesis and in vitro evaluation of antiviral and antitumor activity // Antiviral. Chem. Chemother. 1993. V. 4. No. 6. P. 341-352.
[45] Tam S. Y.-K., Hwang J.-S., De Las Heras F.G., Klein R.S., Fox J.J. Nucleosides. CV. Synthesis of the 8-(ß-D-ribofuranosyl)-pyrazolo[1,5-a]-1,3,5-triazine isosteres of adenosine and inosine (1) //J. Heterocycl. Chem. 1976. V. 13. No. 6. P. 1305-1308.
[46] Cristalli G., Franchetti P., Grifantini M., Nocentini G., Vittori S. 3,7-Dideazapurine nucleosides. Synthesis and antitumor activity of 1-deazatubercidin and 2-chloro-2'-deoxy-3,7-dideazaadenosine // J. Med. Chem. 1989. V. 32. No. 7. P. 1463-1466.
[47] Meade E.A., Wotring L.L., Drach J.C., Townsend L.B. Synthesis, antiproliferative, and antiviral activity of certain
4-aminopyrrolo[2,3-d]pyridazine nucleosides: an entry into a novel
series of adenosine analogues // J. Med. Chem. 1992. V. 35. No. 3. P. 526-533.
[48] Montgomery J.A., Thomas H.J. Nucleosides of 2-azapurines and certain ring analogs // J. Med. Chem. 1972. V. 15. No. 2. P. 182-187.
[49] Walton E., Holly F.W., Jenkins S.R. Indole and 4-aminoindole nucleosides //J. Org. Chem. 1968. V. 33. No. 1. P. 192-197.
[50] Pathil S.A., Otter B.A., Klein R.S. 4-Aza-7,9-dideazaadenosine, a new cytotoxic synthetic C-nucleoside analogue of adenosine // Tetrahedron Lett. 1994. V. 35. No. 30. P. 5339-5342.
[51] Nishimura N., Kato A., Maeba I. Synthesis of pyrrolo[2,1-f][1,2,4]triazine C-nucleosides. Isosteres of sangivamycin, tubercidin, and toyocamycin // Carbohydr. Res. 2001. V. 331. No. 1. P. 77-82.
[52] Ramasamy K., Ugarkar B.G., McKernan P.A., Robins R.K., Revankar G.R. Synthesis and antitumor activity of certain 3-|3-D-ri-bofuranosyl-1,2,4-triazolo[3,4-f ]-1,2,4-triazines related to formycin prepared via ring closure of a 1,2,4-triazine precursor // J. Med. Chem. 1986. V. 29. No. 11. P. 2231-2235.
[53] Kang Y., Larson S.B., Robins R.K., Revankar G.R. Synthesis and biological evaluation of certain 3-|3-D-ribofuranosyl-1,2,4-tri-azolo[4,3-b]pyridazines related to formycin prepared via ring closure of pyridazine precursors //J. Med. Chem. 1989. V. 32. No. 7. P. 1547-1551.
[54] Kobe J., Brdar B., Soric J. Formycin analogs. II. Antiviral and cytotoxic s-triazolo[4,3-a]- and [1,5-a] pyridine derivatives // Nucleos. Nucleot. 1986. V. 5. No. 2. P. 135-151.
[55] HyperChem™ 7.0 for Windows 95/98/ME/NT/2000/XP. Hypercube, Inc. 2002.
[56] http://www.hyper.com/support/software/Scripts/script_index.html (раздел "Conformational Analysis")
или ftp://ftp.hyper.com/pub/scripts/tor_mm.scr .
[57] Dewar M.I.S., Zoebisch E.G., Healy E.F., Stewart J.J.P. AM1: a new general purpose quantum mechanical molecular model //J. Amer. Chem. Soc. 1985. V. 107. No. 13. P. 3902-3909.
[58] IUPAC-IUB Joint Commission of Biochemical Nomenclature: Abbreviations and symbols for the description of conformations of polynucleotide chains // Eur. J. Biochem. 1983. V. 131. P. 9-15.
[59] Зенгер В. Принципы структурной организации нуклеиновых кислот // М.: Мир, 1987. С. 28-29.
[60] Ronca G., Zucchelli G. Competitive inhibition of adenosine deaminase by purine and pyrimidine bases // Biochim. et biophys. acta. 1968. V. 159. No. 1. P. 203-205.
[61] Maguire M.H., Sim M.K. Studies on adenosine deaminase. 2. Specificity and mechanism of action of bovine placental adenosine deaminase // Eur. J. Biochem. 1971. V. 23. No. 1. P. 22-29.
[62] Калиниченко Е.Н., Бейгельман Л.Н., Михайлов С.Н., Михайлопу-ло И.А. Субстратная специфичность аденозиндезаминазы. Роль ме-тильных ^упп при 2', 3' и 5'-aтомах углерода аденозина // Биоорган. химия. 1988. Т. 14. №9. С. 1157-1161.
[63] Laboratory Exercises Using HyperChem®. Chapter 2. Molecular Geometry and Properties. 12. Proton Affinity. P. 50-51.
[64] Bloch A., Robins M.J., McCarthy J.R., Jr. The role of the 5'-hydroxyl group of adenosine in determining substrate specificity for adenosine deaminase //J. Med. Chem. 1967. V. 10. No. 5. P. 908-912.
[65] Hampton A., Harper P.J., Sasaki T. Substrate properties of cycloadenosines with adenosine aminohydrolase as evidence for the conformation of enzyme-bound adenosine // Biochemistry. 1972. V. 11. No. 25. P. 4736-4739.
[66] Mikhailopulo I.A., Wiedner H., Cramer F. Substrate specificity of adenosine deaminase: the role of the substituents at the 2'- and 3'-carbons of adenine nucleosides, of their configuration and of the conformation of the furanose ring // Biochem. Pharmacol. 1981. V. 30. No. 9. P. 1001-1004.
[67] Dodin G., Bensaude O., Dubois J.-E. Tautomerism of formycin. Mechanism of interconversion // J. Amer. Chem. Soc. 1980. V. 102. No. 11. P. 3897-3899.
Поступила в редакцию 19/X1/2003; в окончательном варианте — 19/X1/2003.
THEORETICAL STUDY OF ADENOSINE AND ITS ISOSTERIC ANALOGUES. A POSSIBLE MECHANISM OF THEIR BINDING IN AN ACTIVE SITE OF MAMMALIAN ADENOSINE DEAMINASE
© 2003 Y.P. Zarubin, D.V. Mirgorodsky, P.P. Puryginf I.A. Il'icheva,
V.L. Florent'ev6
Adenosine and its isosteric analogues including their Nl-protonated forms (except for that of 1-deazaanalogues) are studied by computer modelling to find a relationships between their molecular structures and substrate properties for the mammalian adenosine deaminase. The atomic charge distribution and maps of the electrostatic potential around their van der Waals molecular surface are calculated for these compounds by using the ab initio STO-3G method. The conformational studies are carried out by the MM+ method of molecular mechanics. A mechanism of binding that determines the substrate selectivity of mammalian adenosine deaminase is discussed. The potential substrate properties are predicted for earlier unstudied adenosine analogues.
Paper received 19/X1/2003; Paper accepted 19/X1/2003.
5Zarubin Yury Pavlovich ([email protected]), Mirgorodsky Denis Viktorovich, Purygin Pyotr Petrovich, Dept. of Organic Chemistry, Samara State University, Samara, 443011, Russia.
6Il'icheva Irina Alekseevna, Florent'ev Vladimir Leonidovich Laboratory of protein synthesis chemistry, V.A. Engelhardt Institute of Molecular Biology of Russian Academy of Sciences, Moscow, 119991, Russia.