ENSEMBLE DISCRIMINANT ANALYSIS OF RELATIONSHIPS BETWEEN NHE-1 INHIBITORYACTIVITY OF CYCLIC GUANIDINS AND THE ELECTRONIC PARAMETERS OF GUANIDINE FRAGMENT

Full Text

При анализе химической структуры известных ингибиторов натрий-водородного обменника 1 изоформы (NHE-1) выявлено, что большинство из них содержат гуанидиновую группу. Она представлена в структурах амилорида и его производных (DMA, EIPA, MIBA и HMA); производных бензоилгуанидина (НОЕ-642 (ка-рипорид), НОЕ-694, энипорид (EMD-85131), сабипорид (BIIB-722). Известны моно- (S1611, S3226) и бицикли-ческие производные гуанидина, содержащие хинолин (зонипорид (CP-597396), MS 31038), индол (SM 20220, SM 20550, SMP-300), бензоксазинон (№-R9032), дигид-робензофуран (BMS 284640), тетрагидронафтален (T-162559), циклогептапиридин (TY-12533) и другие [2, 5, 7]. Международная база Thomson Reuters Integrity (2011) содержит 262 гуанидинсодержащих производных, обладающих NHE-ингибирующей активностью. Кроме того, опубликованы данные о наличии у циклических гуанидинов - производных бензимидазолов - способности ингибировать NHE-1 [3, 7], которые содержат не только полностью, но и частично встроенную гуанидиновую группу. В ранее проведенных исследованиях на кафедре фармакологии ВолгГМУ было показано, что высокая NHE-1-ингибирующая активность наиболее характерна для производных аминобензимидазола и определяется наличием аминогруппы в положении С2, диэтилами-ноэтильного, пиперидиноэтильного или морфолиноэ-тильного заместителя в N1 [3]. Кроме того, было установлено, что введение в структуру производных имидазо[1,2-а]бензимидазола «активных» для производных аминобензимидазола радикалов в положение N9 или N1 не приводило к высокому уровню NHE-1-ингибирую-щей активности. В последующем, при проведении фармакофорно-го анализа соединений с высоким NHE-1 -ингибирующим уровнем активности, были выявлены потенциальные 2D-фармакофоры, в том числе гуанидиновая группировка -частично или полностью встроенная в гетероцикл. ЦЕЛЬ РАБОТЫ Было выдвинуто предположение, что на реализацию данного вида активности влияют электронные Выпуск 2 (58). 2016 87 ЩШгорСз [ЩсмеТКЩ параметры гуанидинового фрагмента и атомов его ближайшего окружения. Поэтому для исследования были выбраны соединения двух химических классов: производных 2-аминобензимидазола (частично встроенная в цикл гуанидиновая группа) и производные №-имида-зо[1,2-а]бензимидазола (полностью встроенная в цикл гуанидиновая группа), содержащие в одном и том же положении по отношению к азоту в гуанидиновой группе или диэтиламиноэтильный, или пиперидиноэтильный, или морфолиноэтильный заместители. МЕТОДИКА ИССЛЕДОВАНИЯ Расчет центрированных на гуанидиновый фрагмент электронных параметров и анализ их влияния на NHE-1-ингибирующую активность изучаемых веществ проводили в несколько этапов. 1. Изучение NHE-1-ингибирующей активности 25 соединений производных циклических гуанидинов под лабораторными шифрами «РУ», из которых 21 производных №-имидазо[1,2-а]бензимидазола и 4-2-ами-нобензимидазола (табл. 1, 2) в сравнении с селективным ингибитором зонипоридом in vitro на тромбоцитах кролика по методу D. Rosskopf [6] и K. Kusumoto [4], модифицированному [1] в концентрациях (1 х 10-6М; 1 х 10-7М; 1 х 10-8М). 2. Формирование базы данных по структуре соединений с NHE-1 ингибирующей активностью (программа ChemFinfer 9.0, CambridgeSoft Corporation, USA, 2004). 3. Построение методами молекулярной механики первичных 3D-моделей (программа MarvinSketch 15.6.15.0, ChemAxon, USA, 2014). 4. Оптимизация исследуемых конформаций по-луэмпирическим квантово-химическим методом PM7 (пакет MOPAC-2012, Stewart Computational Chemistry, USA, 2012). 5. Выбор лучшей 3D-конформации по минимальной величине общей энергии молекулы (BestConf, Вол-гГМУ, Россия, 2015) и конвертация out-файлов в hin-файлы (программа OpenBabel-2.3.2, Free Software Foundation, USA, 2011). 6. Расчет центрированных на гуанидиновый фрагмент электронных параметров (программа HiperChem 8.0.8, Hypercube, USA, 2009). 7. Формирование обучающих выборок (по результатам кластерного анализа). 8. Пошаговый дискриминантный анализ, выявление значимых систем дискриминатных функций (программа Statistica 8.0, StatSoft, USA, 2007). 9. Расчет граничных значений гуанидинцентриро-ванных электронных параметров, соответствующих высокоактивным соединениям. 10. Анализ соответствия электронных параметров изученных соединений граничным условиям NHE-1 активности ансамбля дискриминантных функций. 11. Расчет средних значений электронных параметров и статистическое сравнение (t-критерий Стьюдента, р < 0,05) групп соединений с полностью встроенной и частично встроенной в цикл гуанидиновой группировкой. РЕЗУЛЬТАТЫ ИССЛЕДОВАНИЯ И ИХ ОБСУЖДЕНИЕ В условиях in vitro все изучаемые соединения проявили NHE-1 -ингибирующую активность различной степени выраженности. Было показано дозозависимое действие. Все изученные соединения в соответствии с граничными значениями активности, полученными в результате кластерного анализа, были разделены на 2 класса: «активные» (а) в концентрации 1 х 10-8М Д(10-8) > 42,0 %; б) в концентрации 1 х 10-7М Д(10-7) > 30,0%; в) в концентрации 1 х 10-6М Д(10-6) > 41,0 %) и «неактивные». В концентрации 1 х 10-8М было выявлено «активных» 32 % производных циклических гуанидинов. В концентрациях 1 х 10-7М и 1 х 10-6М количество таких соединений увеличилось и составило 58 и 68 % соответственно. В исследованиях in silico методами молекулярной механики и квантовой химии были рассчитаны центрированные на гуандиновый фрагмент электронные параметры различного порядка - суммы зарядов на неводородных атомах (в электронных единицах, эл. ед.), начиная с центрального атома углерода гуанидинового фрагмента (Q0) и далее для атомов от первого до седьмого окружения (Q1_Q7). Методом пошагового дискриминантного анализа были рассчитаны зависимости уровня активности в концентрациях веществ 1 х 10-6, 1 х 10-7, 1 х 10-8 М от электронных параметров изученных соединений. В результате было получено 9 систем дискриминантных уравнений (табл. 1). Для всех систем дискриминантных функций были рассчитаны граничные условия гуани-динцентрированных электронных параметров, соответствующих активным соединениям. На активность в концентрации 1 х 10-8М статистически значимо влияют характеристики от нулевого до четвертого уровня (Q0-Q4) (табл. 1). На активность в концентрации веществ 1 х 10-7М статистически значимо оказывают влияние от нулевого до второго уровня (Q0-Q2). При оценке циклических гуанидинов в концентрации 1 х 10-6М было выявлено, что статистически значимо влияет только один электронный параметр Q2 сумма зарядов для атомов второго окружения. При сопоставлении дискриминантных зависимостей активности соединений в диапазоне концентраций 1 х 10-8-1 х 10-6М от центрированных на гуанидиновый фрагмент электронных параметров был найден общий показатель, предположительно обеспечивающий устойчивое наличие NHE-1-ингибирующей активности - это сумма зарядов на атомах второго окружения Q2 от центрального атома углерода (первое окружение атомов азота в гуанидиновой группировке). Обращает на себя тот факт, что с уменьшением концентрации число статистически значимых электронных 88 Выпуск 2 (58). 2016 параметров возрастает от одного до трех и до пяти соответственно. Это может свидетельствовать о том, что именно в низких концентрациях в большей степени проявляется специфичность взаимодействия изучаемых соединений с мишенью, тогда как в высоких концентрациях наблюдается менее селективное взаимодействие. С учетом граничных условий был определен прогнозный спектр уровней активности по данным дискриминантного анализа для каждого из испытанных циклических гуанидинов в концентрации 1 х 10-6, 1 х 10-7, 1 х 10-8 М. Индекс активности соединения принимается равным 1, если значение его электронного параметра соответствовало пороговому условию; в противном случае индекс активности полагали равный 0. Прогнозный спектр по данным ансамблевого дискриминантного анализа циклических гуанидинов, испытанных на NHE-1 ингибирующую активность, позволил выявить соединения с высокой активностью (табл. 2). К ним относятся производные 2-аминобезимидазола - соединения РУ-1355, РУ-1390, РУ-1199, РУ-1200 с частично встроенной в цикл гуанидиновой группой. Таблица1 Статистически значимые системы дискриминатных функций для NHE-1-ингибирующей активности некоторых производных циклических гуанидинов Элект ронный параметр Системы дискриминатных функций для 3 исследуемых концентраций X О со 1 х 10'7 М 1 х 10'6 М Qû ginact-17,54+1 06,41*Qû (1) ■ gact=-29,79+137,53*Qü N=25; F(1,23) = 15,96; p = 5,7*10'4; Fo=0,840; Finact = 1,000; Fact = 0,500; Г раничные условия для активных соединений Q0 > 0,394 ginact=-16,44+84,27*Q0 (6) ■ gact=-12,36+71,63*Q0 N=25; F(1,23)=4,17; p=5,3*10'2; F0=0,792; Finact=600; Fact=0,929; Граничные условия для активных соединений Q0>0,322 Q1 ginact=-30,18-58,52*Q1 (2) • gact=-45,08-71,07*Q1 N=25; F(1,23)=14,92; p=7,9*10'4; F0=0,840; Finact=1,000; Fact=0,500; Г раничные условия для активных соединений Q1 < -1,187 ginact=-26,23-44,92*Q1 (7) • gact=-21,19-39,95*Q1 N=25; F(1,23)=3,67; p=6,8*10'2; F0=0,625; Finact=0,500; Fact=0,714; Граничные условия для активных соединений Q1 < -1,0130 Q2 'ginact=-0,47+6,26*Q2 (3) ' gact=-1,35-9,91*Q2 Vl=25; F(1,23)=6,01; p = 2,2*10'2; F0=0,800; Finact=0,880; Fact=0,630; Граничные условия для активных соединений Q2 < 0,054 'ginact-0,55-3,28*Q2 (8) ■ gact=-1,11 + 11,29*Q2 ■ N=25; F(1,23) = 4,69; p = 4,1*10-2; F0=0,625; Finact=0,700; Fact=0,571; Граничные условия для активных соединений Q2>0,038 finact = -1,37+10,16*Q2 (9) gact = -0,41-3,27*Q2 N = 25; F(1,23) = 4,39; p = 4,7*10'2; F0 = 0,640; Finact = 0,000; Fact=0,941; Граничные условия для активных соединений Q2 > 0,0715 Q3 ginact=-4,1 4-14,81*Q3 (4) ' gact=-2,11-7,52*Q3 N=25; F(1,23)=9,90; p = 4,5*10'3; F0=0,800; Finact=0,941; Fact = 0,500; Граничные условия для активных соединений Q3>-0,279 Q4 ginact=-5,53-12,85*Q4 (5) ■ gact=-11,60-18,31*Q4 n=25; F(1,23)=10,15; p=4,1*10-3; F0=0,840; Finact=1,00; Fact=0,500; Граничные условия для активных соединений Q4<-1,110 Примечание. ginact - дискриминатная функция для класса неактивных соединений; gact - дискриминатная функция для класса активных соединений; N - количество исследуемых соединений; F (1,23) - значение критерия Фишера для числа степеней свободы; Fo - общая точность прогноза; Finact - точность прогноза для неактивных соединений; Fact - точность прогноза для активных соединений; p - статистическая значимость системы дискриминантных функций. Выпуск 2 (58). 2016 89 ÜSempfä tgOMrf ШШ Таблица 2 Общий прогнозный спектр циклических гуанидинов, испытанных на NHE-1-ингибирующую активность (данные ансамблевого дискриминантного анализа) № Соеди нение Индекс активности соединения (в исследуемой концентрации, М) X О à X О X О о: м ^-имидазо[1,2-а]бензимидазола 1 РУ-0185 2 2 0 4 2 РУ-0376 1 3 0 4 3 РУ-87 0 3 0 3 4 РУ-35 1 2 0 3 5 РУ-67 0 3 0 3 6 РУ-284 1 2 0 3 7 РУ-285 1 2 0 3 8 РУ-292 1 2 0 3 9 РУ-293 0 2 1 3 10 РУ-294 1 2 0 3 11 РУ-374 1 2 0 3 12 РУ-474 1 2 0 3 13 РУ-1205 1 2 0 3 14 РУ-64 1 1 0 2 15 РУ-17 0 2 0 2 16 РУ-66 1 1 0 2 17 РУ-273 0 2 0 2 18 РУ-1193 0 2 0 2 19 РУ-1204 0 2 0 2 20 РУ-1189 0 1 0 1 21 РУ-1191 1 0 1 2-аминобензимидазола 22 РУ-1199 5 2 0 7 23 РУ-1200 5 2 0 7 24 РУ-1355 5 2 0 7 25 РУ-1390 5 2 0 7 При сравнении средних значений электронных параметров двух групп изученных соединений с полностью и частично встроенной в цикл гуанидиновой группировкой выявлены статистически значимые различия для всех уровней электронного окружения (рис.). я электронных зарядов, э о о г L/Ï о '-л >- У -1-1-1-1-1 * n- 1 -г_ 1 1 и g -1 - ^ -1,5 -1 , - Соединения, содержащие в своей структуре гунидиновую группу □ частично встроенную в цикл □ полностью встроенную в цикл ° QO Q1 Q2 Q3 Q4 Q5 Q6 Q7 Уровень электронного окружения *Данные достоверны по отношению к группе, содержащей полностью встроенную в цикл гуанидиновую группу, t-критерий Стьюдента (p < 0,05). Рис. Сравнительный анализ центрированных на гуанидиновый фрагмент электронных параметров соединений, содержащих частично и полностью встроенную в цикл гуанидиновую группу, испытанных на NHE-1-ингибирующую активность Можно предположить, что центрированные на гуанидиновый фрагмент электронные параметры соединений с частично встроенной гуанидиновой группировкой наиболее оптимальны для взаимодействия с NHE-1. ЗАКЛЮЧЕНИЕ В исследованиях in silico методом молекулярной механики и квантовой химии рассчитаны гуанидинцентри-рованные электронные параметры для соединений с частично (производные 2-аминобензимидазола) и полностью (производные №-имидазо[1,2-^бензимидазола) встроенной в цикл гуанидиновой группой. Методом ансамблевого дискриминантного анализа проведена оценка зависимости NHE-ингибирующий активности от данных показателей. Было выявлено, что электронные параметры второго порядка Q2, характеризующие ближайшее окружение атома азота гуанидинового фрагмента, оказывают статистически значимое влияние на уровень активности во всех изученных концентрациях и предположительно определяют наличие устойчивой NHE-1 -ингибирующий активности. Наиболее статистически значимые результаты были получены для концентрации 1 х 10-8М, что, скорее всего, связано с усилением роли селективных взаимодействий в формировании NHE-1-ингибирующей активности. В результате прогноза методом ансамблевого дискриминантного анализа было выявлено, что соединения, имеющие частично встроенную гуанидиновую группу, обладают более высокой NHE-1-ингибирующий активностью. Центрированные на гуанидиновый фрагмент электронные параметры соединений с полностью встроенной в цикл гуанидиновой группой статистически значимо отличаюся от веществ с частично встроенной.

About the authors

P. M. Vassiliev

Volgograd State Medical University

Email: pvassiliev@mail.ru
department of pharmacology

A. A. Spasov

Volgograd State Medical University

department of pharmacology

V. U. Muravyova

Volgograd State Medical University

department of pharmacology

N. A. Gurova

Volgograd State Medical University

department of pharmacology

L. R. Yanalieva

Volgograd State Medical University

department of pharmacology

V. V. Vorfolomeeva

Volgograd State Medical University

department of pharmacology

V. A. Anisimova

Research institute of physical and organic chemistry of the Southern Federal University

M. V. Bogoslavtseva

Research institute of physical and organic chemistry of the Southern Federal University

References

Supplementary files