Научный журнал

Международный журнал прикладных и фундаментальных исследований

ISSN 1996-3955

ИФ РИНЦ = 0,593

• Авторы
• Резюме
• Файлы
• Ключевые слова
• Литература

Андрюков К.В. 1

---

1 ФГБОУ ВО «Пермская государственная фармацевтическая академия» Министерства здравоохранения Российской Федерации

Методом молекулярного моделирования с использованием молекулярного докинга построены трехмерные модели исследования взаимодействия амидов и гидразидов N-ароилзамещенных антраниловых кислот с мутантной формой 5-липоксигеназы. Выполнен расчет новых конформаций для двадцати гомологов производных антраниловой кислоты ингибированием активного участка фермента. Расчет и построение новых конформаций проводили с помощью программного комплекса MGL Tools. Молекулярное конструирование производных амидов и гидразидов N-ароилзамещенных антраниловых кислот выполняли в активную ячейку фермента 5-липоксигеназы (PDB ID: 3v99), содержащей арахидоновую кислоту. Анализ результатов расчета новых конформаций проводили по RMSD относительно арахидоновой кислоты. Энергия связывания с 3V99 по величине Be5-ЛОГ получена в интервале от -3,24 до -5,92 ккал/моль и характеризует низкое сродство исследуемых веществ к данной мутантной форме фермента. Оценку взаимодействия с использованием конформаций проводили так же по межмолекулярной энергии (Ime5-ЛОГ) и ингибированию константе (Ki5-ЛОГ). Максимальное значение энергии связывания с 5-липоксигеназой, равное -5,92 ккал/моль, найдено по конформации № 9 с RMSD = 7,17 для соединения IV, межмолекулярная водородная связь при этом образуется по аминокислоте аргинин (ARG596). Осуществлена проверка результатов моделирования на экспериментальных данных торможения каррагенинового отека с помощью множественного линейного регрессионного анализа.

Статья в формате PDF

645 KB

антраниловая

молекулярный докинг

5-липоксигеназа

противовоспалительная активность

1. Hošek J., Leláková V., Bobal P., Pizova H., Gazdová M., Malaník M. Prenylated stilbenoids affect inflammation by inhibiting the NF-κB/AP-1 signaling pathway and cyclooxygenases and lipoxygenase // Journal of natural products. 2019. Vol. 82, Is. 7. Р. 1839–1848.

2. Abdelgawad M.A., Labib M.B., Ali W.A.M., Kamel G., Azouz A.A., El-Nahass E.S. Design, synthesis, analgesic, anti-inflammatory activity of novel pyrazolones possessing aminosulfonyl pharmacophore as inhibitors of COX-2/5-LOX enzymes: Histopathological and docking studies // Bioorganic chemistry. 2018. Vol. 78. P. 103–114.

3. Walker C., Biasucci L.M. Cardiovascular safety of non-steroidal anti-inflammatory drugs revisited // Postgraduate Medicine. 2018. Vol. 130, Is. 1. P. 55–71.

4. Koundouros N., Poulogiannis G. Reprogramming of fatty acid metabolism in cancer // British journal of cancer. 2020. Vol. 122, Is. 1. P. 4–22.

5. Han S.H., Suh H.S., Jo H., Oh Y., Mishra N.K., Han S., Kim I.S. Synthesis and anti-inflammatory evaluation of N-sulfonyl anthranilic acids via Ir (III)-catalyzed C–H amidation of benzoic acids // Bioorganic & Medicinal Chemistry Letters. 2017. Vol. 27, Is. 10. P. 2129–2134.

6. Du L., Du S., Li J., Wang H. Design, synthesis, and biological evaluation of dual-target COX-2/5-LOX inhibitors for the treatment of inflammation // Medicinal Chemistry Research. 2023. Vol. 32, Is. 2. P. 218–238.

7. Kausar N., Muratza S., Raza M.A., Rafique H., Arshad M.N., Altaf A.A., Shafqat S.R. Sulfonamide hybrid schiff bases of anthranilic acid: synthesis, characterization and their biological potential // Journal of Molecular Structure. 2019. Vol. 118. P. 8–20.

8. Gilbert N.C., Rui Z., Neau D.B., Waight M.T., Bartlett S.G., Boeglin W.E., Newcomer M.E. Conversion of human 5-lipoxygenase to a 15-lipoxygenase by a point mutation to mimic phosphorylation at Serine-663 // The FASEB Journal. 2012. Vol. 26, Is. 8. P. 3222–3229.

9. Шакирова А.Б., Подчезерцева А.В., Коркодинова Л.М., Васильева М.Ю., Колла В.Э., Марданова Л.Г., Вахрина Н.И., Данилов Ю.Л. Синтез, противовоспалительная и анальгетическая активности некоторых ариламидов N-замещенных антраниловых кислот // Химико-фармацевтический журнал. 2001. Т. 35, № 4. С. 193–195.

Воспаление является нормальной реакцией на любой вредный раздражитель, угрожающий организму, и может варьироваться от локализованной до генерализованной реакции. Нестероидные противовоспалительные препараты (НПВП) достигают своего противовоспалительного действия за счет ингибирующего действия на фермент циклооксигеназу (ЦОГ), белок необходим для биосинтеза простагландинов из арахидоновой кислоты. ЦОГ существует в двух изоформах, а именно ЦОГ-1 и ЦОГ-2. В общих чертах, ЦОГ-1 отвечает за защиту слизистой оболочки желудка и поддержание гомеостаза, тогда как ЦОГ-2 индуцируется провоспалительными стимулами в очагах воспаления. Арахидоновая кислота является одним из основных компонентов фосфолипидов клеточных мембран и действует как нормальный субстрат для различных ферментов, таких как 5-липоксигеназа (5-ЛОГ) [1] и ЦОГ [2].

5-ЛОГ метаболизирует арахидоновую кислоту с образованием лейкотриенов, которые являются медиаторами при астме и воспалении, и инициирует биосинтез провоспалительных лейкотриеновых липидных медиаторов и вместе с 15-ЛОГ также необходима для синтеза противовоспалительных липоксинов. 5-липоксигеназа (5-ЛОГ) является одним из трех основных изоферментов, участвующих в метаболическом процессе арахидоновой кислоты, и существует в виде компенсаторного механизма с ЦОГ [3]. Лейкотриены (ЛТ), в конечном счете продуцируемые ЛОГ-путем, рассматриваются как важные медиаторы воспалительных и аллергических заболеваний и, как полагают, тесно связаны с раком и сердечно-сосудистыми заболеваниями [4]. Каталитическая активность 5-ЛОГ регулируется с помощью множества механизмов, включая Ca(2+)-направленное связывание с мембраной и фосфорилирование по специфическим сериновым остаткам. Следовательно, можно снизить потенциальный риск сердечно-сосудистых заболеваний, связанных с НПВП, путем ингибирования 5-ЛОГ и снижения синтеза ЛТ.

Биологическая активность производных антраниловой кислоты характеризуются широким спектром: анальгетической, противомикробной, в том числе и противовоспалительной [5, 6] и антиоксидантной [7].

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

Материалы и методы исследования

Молекулярное моделирование исследуемых соединений ряда антраниловой кислоты осуществляли с помощью: ChemDraw и Gaussian 03. Расчет анализируемых конформаций с 5-липоксигеназой с PDB ID Code: 3v99 [8] проводили программным комплексом MGL Tools и AutoDock 4.2.. Подготовку исследуемых производных для молекулярного докинга AutoDock 4.2 осуществляли переведением веществ в трехмерный формат (.pdb) программой Marvin Sketch. Проверку результатов моделирования осуществляли на экспериментальных данных торможения каррагенинового отека программой Statistica 6.

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

Молекулярное моделирование двадцати гомологов (I–XX) амидов и гидразидов N-ароил антраниловых кислот проводили программой ChemDraw.

R1= H, R2= 4–BrC6H4, R3= NHCH2CH=CH2 (I); R1= H, R2= 2–фурил, R3= NHCH2CH=CH2 (II); R1= Br, R2= 2–COOH C6H4, R3= NH2 (III); R1= I, R2=2–фурил, R3=NHCH2C6H5 (IV); R1=H, R2=4–NO2C6H4, R3=NHCH2CH=CH2 (V); R1=Br, R2=2–фурил, R3= NHNHCOCH2Cl (VI); R1=Br, R2=2–фурил, R3= NHCH2CH2CH(CH3)2 (VII); R1= Br, R2= 2–фурил, R3= NHCH3 (VIII); R1= I, R2= –фурил, R3=NHCH3 (IX); R1= H, R2= 3, 4, 5–(OCH3)3C6H2, R3= NHCH2CH=CH2 (X); R1=Br, R2=2–фурил, R3=NHNHCOCH2CH2CH3 (XI); R1=Br, R2=2–фурил, R3=NHNHCO(2-фурил) (XII); R1= Br, R2= 2–фурил, R3= NHNHCOCH2C6H5 (XIII); R1= Br, R2= 2–фурил, R3=NHNHCO 2–COOHC6H4 (XIV); R1=I, R2=2–фурил, R3=N(CH3)2 (XV); R1=Br, R2=2–OCH3C6H4, R3=NH2 (XVI); R1=Br, R2=3–NO2C6H4, R3=NH2 (XVII); R1=Br, R2=2–фурил, R3=NHNH2 (XVIII); R1=H, R2= 4–OCH3C6H4, R3=NHCH2CH=CH2 (XIX); R1=I, R2=C6H5, R3=NHCH2CH2OH (XX).

Спектральные характеристики и результаты экспериментального определения противовоспалительной активности ПВА (ПВАэксп.) на модели торможения каррагенинового отека объекта исследования (I–XX) опубликованы в работе [9].

Для проведения расчета новых конформаций исследуемых амидов и гидразидов N-ароил антраниловых кислот выполняли молекулярный докинг с мутантной формой 5-ЛОГ с 10 мутациями [8].

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

Рис. 1. Диаграмма взаимодействия арахидоновой кислоты с 5-ЛОГ

Таблица 1

Результаты молекулярного докинга по 5-ЛОГ: номер конформации, скоринговые функции и остатки аминокислот с водородной связью (Н-связь) и ПВАэксп.

Кроме того, активный центр модели 3V99 находится в открытой конфигурации, доступной для второй небольшой молекулы – исследуемого лиганда.

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

Локализацию сайта связывания проводили на основе комплекса внутри исследуемой мутантной формы 5-ЛОГ: арахидоновой кислоты с активным участком (рис. 1).

На рис. 1 изображен комплекс арахидоновой кислоты с активным участком 5-ЛОГ.

Диаграмма, представленная на рис. 1, показывает наличие ионизационных и гидрофобных взаимодействий. Ионизационная связь образуется по карбоксильной группе арахидоновой кислоты с остатком фенилаланина (PHE177). Структурные фрагменты арахидоновой кислоты, содержащие двойные связи и алкильный радикал, гидрофобно ориентируются рядом с остатками фенилаланина (PHE177, PHE555, PHE610), изолейцина (ILE406), лейцина (LEU607), аланина (ALA672). Локальный участок связывания внутри 5-ЛОГ мы определили в координатах: x= 16,91, y= – 72,12, z= – 47,24.

Далее проводили построение 20 Grid-карт (x= 16,91, y= – 72,12, z= – 47,24), с координатами точек (60 × 60 × 60) внутри исследуемой ячейки 5-ЛОГ. Выполнен расчет 10 конформаций по каждому исследуемому амиду и гидразиду N-ароил антраниловой кислоты. Анализ результатов расчета новых конформаций проводили по RMSD относительно арахидоновой кислоты. Полученные конформации по каждому исследуемому соединению представлены в табл. 1.

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

Наибольшее значение энергии связывания с 5-ЛОГ, равное -5,92 ккал/моль, найдено по конформации № 9 для соединения IV (рис. 2), межмолекулярная связь показывает взаимодействие карбонильной группы в NH-ацильном фрагменте по аминокислоте аргинин (ARG 596).

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

Проверку осуществляли исследованием зависимости ПВАэксп. от Be5-ЛОГ, Ime5-ЛОГ, Ki5-ЛОГ. Полученные результаты представлены уравнениями с наибольшими значениями коэффициента множественной регрессии (R), критерия Фишера (F) и минимальным значением средней квадратичной ошибки (S) (табл. 2).

Рис. 2. Диаграмма взаимодействия соединения IV с 5-ЛОГ в активном участке фермента

Таблица 2

Уравнения регрессии зависимости ПВАэксп. от Be5-ЛОГ, Ime5-ЛОГ, Ki5-ЛОГ

Выводы

1. Методом молекулярного моделирования построены трехмерные модели исследования взаимодействия амидов и гидразидов N-ароил антраниловых кислот с мутантной формой 5-ЛОГ.

2. Выполнен расчет новых двадцати конформаций взаимодействия исследуемых производных с мишенью 3v99 мутантной формы 5-липоксигеназы.

3. Осуществлена проверка результатов поиска соединений с ПВА на экспериментальных данных торможения каррагенинового отека, в результате получены две двухфакторные модели зависимости с высоким значением R в интервале от 0,761 до 0,762.

Библиографическая ссылка