Все реагенты были аналитического качества и использовались без дальнейшей очистки. Ход реакции отслеживался с помощью тонкослойной хроматографии (TLC) на силикагеле с использованием n-гексана/этилацетата (7:3–6:4, v/v), при этом пятна визуализировались под ультрафиолетовым светом. Если не указано иное, все реакции проводились в окружающей атмосфере. Синтезированные соединения были охарактеризованы инфракрасной спектроскопией преобразования Фурье (FTIR), 1H ЯМР и 13C ЯМР-спектроскопией. ЯМР-спектры 1H и 13C были зафиксированы в CDCl₃ на частотах 600 МГц и 150 МГц соответственно, а химические сдвиги фиксировались в ppm относительно остаточного сигнала растворителя. Синтетический путь показан на рисунке 1. Полные данные физико-химической характеристики всех синтезированных соединений приведены в дополнительном файле 1.
Синтез этилбензоата (1)
Бензойная кислота (50 г, 0,41 моль) была помещена в колбу с круглым дном объемом 500 мл, а также добавлен абсолютный этанол (150 мл). Концентрированная серная кислота (6 мл) добавлялась по каплям в качестве катализатора при непрерывном магнитном перемешивании. Реакционная смесь нагревалась под рефлюксом при 78 °C в течение 4–6 часов с использованием водяного конденсатора. Ход реакции отслеживался TLC. После завершения смесь охлаждалась до комнатной температуры и нейтрализовалась насыщенным водным карбонатом натрия до достижения pH 7–8. Реакционная смесь была перенесена в сепаративную воронку, органический слой собирался, промывался дистиллированной водой и сушился на безводном сульфате натрия. После фильтрации этилбензоат (1) использовался непосредственно на следующем этапе без дальнейшей очистки (рисунок 1).
Синтез бензогидразида (2)
Этилбензоат (1; 0,41 моль) растворялся в абсолютном этаноле (100 мл), а гидразингидрат (80%, 25 мл, избыток) добавлялся при непрерывном перемешивании. Реакционная смесь нагревалась под рефлюксом при 78 °C в течение 4–5 часов под воздухом с использованием водяного конденсатора. Ход реакции отслеживался TLC. После завершения работы белый кристаллический осадок бензогидразида (2) был собран фильтрацией, промытый холодным этанолом, высушен до постоянного веса при комнатной температуре и использован на последующих синтетических этапах (см. рисунок 1).
Синтез 5-фенил-1,3,4-оксадиазол-2-тиол (3)
Бензогидразид (2; 10 г, 0,073 моль) и гидроксид калия (4,12 г, 0,073 моль) были помещены в круглодонную колбу объёмом 250 мл с этанолом (30 мл). Смесь нагревалась под рефлюксом при 78 °C в течение 10 минут с использованием водяного конденсатора. Затем был добавлен дисульфид углерода (CS₂; 5 мл, 0,083 моль), и рефлюкс продолжался ещё 3–4 часа. Ход реакции отслеживался TLC. По завершении реакционная смесь перемещалась в дистиллированную воду и нагревалась с непрерывным перемешиванием. Разбавлённую соляную кислоту добавляли по каплям до pH до 5, что приводило к осадке соединения (3) (рисунок 1). Осад собирали фильтрацией, тщательно промывали холодной дистиллированной водой, сушили до постоянного веса и взвешивали для определения выхода. Соединение характеризулось спектроскопией FTIR, 1H ЯМР и 13C ЯМР.
Синтез 4,5-дифенил-4H-1,2,4-триазол-3-тиол (4)
Бензогидразид (2; 5 г, 0,037 моль) и фенилизотиоцианат (4 мл, 0,037 моль) были помещены в колбу с круглым дном объемом 100 мл, оснащённую конденсатором с водяным охлаждением. Реакционная смесь нагревалась под рефлюксом при 78 °C до полного образования промежуточного продукта тиосемикарбазида, что было подтверждено TLC с использованием n-гексана/этилацетата (7:3, v/v). Затем добавлялся гидроксид калия (2 г, 0,036 моль), растворённый в этаноле (20 мл), и рефлюкс продолжался до образования соединения (4), что подтвердил TLC. Реакционная смесь переливалась в дистиллированную воду, нагревалась с непрерывным перемешиванием и подкислялась разбавленной соляной кислотой до pH 5 для осаждения соединения (4) (рисунок 1). Осадок собирался фильтрацией, тщательно промывался холодной дистиллированной водой, высушен до постоянного веса и характеризулся спектроскопией FTIR, 1H ЯМР и 13C ЯМР.
Синтез промежуточных продуктов бромоацетамида 7(a–f)
Каждая заменённая анилиновая производная 5(a–f) (1 эквивалент) растворялась в хлороформе и охлаждалась в ледяной ванне для поддержания температуры реакции 0–5 °C. Бромоацетилбромид (6, ≥98%, 1 эквивалент) добавлялся по капле с непрерывным перемешиванием. Во время реакции по необходимости добавляли 18% водный раствор карбоната натрия для поддержания pH реакции между 9 и 10 путём нейтрализации образующейся броминой кислоты. После завершения добавления бромоацетилбромида реакционную смесь давали нагреться до комнатной температуры и перемешывали до тех пор, пока TLC не подтвердил полное потребление исходного материала. Полученные промежуточные продукты бромоацетамида 7(a–f) собирались фильтрацией, промывались холодной дистиллированной водой, сушили до постоянного веса и использовались непосредственно на последующем синтетическом этапе без дальнейшей очистки (см. рисунок 1).
Синтез производных ацетамида на основе триазола 8(a–c)
Производное триазола (4; 0,5 г, 0,002 моль) растворялось в N-диметилформамиде (10 мл) в колбе с круглым дном объемом 50 мл. Соответствующий промежуточный продукт бромоацетамида (7a, 7b или 7c; 0,002 моль, 1 эквивалент) добавлялся с непрерывным перемешиванием, затем гидрид натрия (0,05 г, 0,002 моль) для стимулирования S-алкилирования. Смесь для реакции перемешивалась при комнатной температуре в течение 24–48 часов, при этом ход реакции отслеживался с помощью TLC. После завершения реакционной смеси заливали в ледяную воду для выпадения осадков. Осадок был собран фильтрацией, тщательно промыт холодной водой, высушен и характеризуется FTIR, 1H NMR и 13CЯМР-спектроскопией (рисунок 1).
Синтез производных ацетамида на основе оксадиазола 9(d–f)
Производное оксадиазола (3; 0,5 г, 0,002 моль) растворялось в N-диметилформамиде (10 мл) в колбе с круглым дном объемом 50 мл. Был добавлен соответствующий промежуточный продукт бромоацетамида — 7d, 7e или 7f (0,002 моль, эквивалент), а затем гидрид натрия (0,05 г, 0,002 моль) для облегчения S-алкилирования. Реакционная смесь перемешивалась при комнатной температуре в течение 24–48 часов, при этом прогресс отслеживался TLC. После завершения смесь заливали в ледяную воду для осадка. Осадок был собран фильтрацией, тщательно промыт холодной водой, высушен и характеризуется FTIR, 1H NMR и 13CЯМР-спектроскопией (рисунок 1).
Анализ ингибирования α-глюкозидазы
Ингибирующая активность α-глюкозидазы оценивалась путём смешивания 50 мкл раствора образца с 100 мкл раствора α-глюкозидазы (0,35 U/мл), приготовленного в фосфатном буфере (pH 6,8). Реакционная смесь инкубировала при 37 °C в течение 10 минут, затем в качестве субстрата добавляли 100 мкл 1,5 мМ p-нитрофенил-α-D-глюкопиранозида. После дальнейшей инкубации при 37 °C в течение 20 минут реакция завершалась добавлением 1 мл карбоната натрия (Na₂CO₃). Поглощение измерялось на уровне 400 нм14. Акарбоза использовалась в качестве положительного контроля. Реакционная смесь без тестового соединения служила отрицательным контролем, тогда как реакционная смесь без фермента использовалась в качестве заготовки.
Все эксперименты проводились в тройных экземплярах, а результаты выражаются в виде среднего ± стандартного отклонения (SD). Процентное ингибирование активности α-глюкозидазы рассчитывалось с использованием следующего уравнения(1)14:
(1)
где As обозначает поглощение раствора пробного образца, Ab — поглощение реагентного бланка (без α-глюкозидазы), а A 0 — поглощение отрицательного контроля (без образца). Значения IC50 определялись путём тестирования концентраций соединений в диапазоне от 3,125 до 100 мкм. Нелинейный регрессионный анализ проводился с использованием стандартного программного обеспечения для построения графиков.
Анализ ингибирования уреазы
Ингибирующая активность уреазы была определена с помощью модифицированного метода Бертело. Реакционная смесь состояла из 300 мкл раствора пробного образца, 200 мкл раствора уреазы джека и 200 мкл мочевины 25 мМ. Смесь инкубировала в водяной бане при 37 °C в течение 30 минут. После инкубации добавляли 500 мкл раствора с цветообразующим форму, содержащего 10 г/л фенола и 50 г/л нитропруссида натрия, а затем 500 мкл второго раствора, содержащего 5% (v/v) гипохлорита натрия и 5 г/л гидроксида натрия. Реакционная смесь затем инкубировала дополнительно 15–30 минут для развития цвета.
Поглощение раствора измерялось на уровне 625 нм по сравнению с реагентом15. Тиомочевина использовалась в качестве положительного контроля, реакционная смесь без тестового соединения — в отрицательном контроле, а смесь без фермента уреазы — в качестве бланка. Все эксперименты проводились в тройном экземпляре (n = 3), а результаты выражаются в средних значениях ± стандартном отклонении (S.D.). Процент ингибирования рассчитывался по следующей математической формуле (2)15:
(2)
где A0 — это поглощение контроля (без ингибитора), а A s — поглощение тестового образца (с ингитором). Построение графиков и нелинейный регрессионный анализ для определения значений IC50 проводились с использованием стандартного программного обеспечения для построения графиков данных.
BSA-тест связывания
Связывающие взаимодействия между сывороточным альбумином крупного рогатого скота (BSA) и синтезированными соединениями изучались с помощью флуоресцентной спектроскопии тушения. В буфере фосфатов 20 мМ (pH 7,4) был приготовлен раствор BSA с массой 10 мкм. Для каждого эксперимента с титрованием 1,0 мл раствора BSA смешивалось с 2,0 мл фосфатного буфера в кварцевом кюветте. Раствор титрировался последовательно добавлением 5 мкл аликвотов тестового соединения (1,5 мМ метанола). Спектры флуоресцентного излучения были зафиксированы в диапазоне длин волн от 300 до 400 нм после возбуждения на 295 нм, при этом максимальное излучение наблюдалось примерно на 336 нм. Измерения флуоресценции проводились с помощью флуоресцентного спектрофотометра, оснащённого стандартным программным обеспечением для получения флуоресценции, с шириной возбуждения и эмиссионных щельей соответственно 10 нм и 2,5 нм. Для определения термодинамических параметров связывания, включая свободную энергию Гиббса (ΔG), энтальпию (ΔH) и энтропию (ΔS), эксперименты по титрированию были повторены при 298, 308 и 313 K.
Для изучения механизма флуоресцентного затухания экспериментальные данные анализировались с использованием уравнения Штерна–Вольмера (уравнение 3)16:
(3)
где F0 и F обозначают флуоресцентные интенсивности BSA при отсутствии и присутствии тушителя соответственно; Kq — это константа скорости бимолекулярного закаления; τ0 — средний срок службы BSA при отсутствии тушителя; KSV — это постоянная гасения Штерна-Вольмера; и [Q] — концентрация тушителя. Константа связывания (Kb) и количество сайтов связывания (n) впоследствии определялись путём построения двойного логарифма согласно уравнению (4)16:
(4)
Кроме того, стандартное изменение свободной энергии Гиббса (ΔG°), стандартное изменение энтальпии (ΔH°) и стандартное изменение энтропии (ΔS°), регулирующее взаимодействие связывания, были рассчитаны с помощью взаимосвязей Ван'т Хоффа и Гиббса–Гельмгольца, выраженных в уравнениях (5) и (6)17:
ΔG∘ = -RT ln Kb (5)
(6)
где R — универсальная газовая константа, а T — абсолютная температура в Кельвине.
Молекулярная стыковка
Все молекулярное моделирование и симуляции стыковки проводились с использованием комплекса Шрёдингера (интерфейс Maestro), а стыковка осуществлялась с помощью модуля Glide. Кристаллические структуры α-глюкозидазы (PDB ID: 5NN8) и уреазы (PDB ID: 4H9M) были извлечены из белкового банка RCSB. Белковые структуры были подготовлены с помощью Белкового мастера подготовки путём удаления сокристаллизованных лигандов, молекул воды и несущественных гетероатомов, добавления недостающих остатков и петлей, присвоения правильного порядка связей и добавления полярных атомов водорода. Протонирование и таутомерные состояния определялись при физиологическом pH (7,0 ± 2,0) с помощью модуля Epik. Подготовленные белковые структуры впоследствии подвергались сдержанной минимизации энергии с помощью силового поля OPLS для облегчения стерических столкновений.
Двумерные структуры соединений 3, 4, 8(a–c) и 9(d–f) вместе с эталонными ингибиторами акарбоза и тиочевница были подготовлены с помощью модуля LigPrep. Были оптимизированы порядок связей, стереохимия и ионизационные состояния при физиологическом pH, а полученные трёхмерные конформации лиганда минимизированы по энергии с помощью силового поля OPLS.
Сетки рецепторов генерировались с помощью инструмента генерации сетки рецепторов, при этом коробки сетки были центрированы вокруг кокристаллизованных нативных лигандов для определения активных участков. Размеры сетки были установлены на 25 × 25 × 25 Å для обоих белков. Уреазная сетка была центрирована по координатам x = −38,76, y = −44,14 и z = 66,88, тогда как сетка α-глюкозидазы была центрирована в x = 21,14, y = −7,56 и z = 24,37.
Расчёты стыковки проводились с использованием режима повышенной точности (XP) Glide с гибким лигандным выборкой. Для каждого лиганда было сгенерировано несколько поз связывания, и для дальнейшего анализа выбиралась поза с наименьшим баллом стыковки Glide. Протокол стыковки был проверен путём повторного стыковки сокристаллизованных лигандов с соответствующими рецепторными сетками. Значения отклонения по среднеквадратичному корню (RMSD) ≤2.0 Å подтверждали надёжность и воспроизводимость протокола стыковки.
Конформации связывания стыкованных комплексов визуализировались с помощью PyMOL (версия 3.1)18. Взаимодействия белков и лиганда, включая водородные связи, гидрофобные контакты и электростатические взаимодействия, были проанализированы с помощью Discovery Studio Visualizer19.
Анализ теории функционала плотности (ДПФ)
Все расчёты по теории функционала плотности (DFT) выполнялись с использованием гауссианского 09, а оптимизированные молекулярные геометрии и граничные молекулярные орбитали визуализированы с помощью GaussView 5. Синтезированные производные ацетамида на основе триазола и оксадиазола, 8(a–c) и 9(d–f), подверглись геометрической оптимизации в водной фазе с использованием гибридного функционала Бекка с тремя параметрами Ли–Янга–Парра (B3LYP) в сочетании с базисным набором 6-311+G(d,p) 20. Эффекты растворителя были реализованы с помощью модели поляризируемого континуума (CPCM)21. Впоследствии на том же уровне теории были проведены расчёты колебательных частот, чтобы подтвердить, что все оптимизированные структуры соответствуют истинным локальным минимумам на поверхности потенциальной энергии, что подтверждается отсутствием мнимых частот.
Для понимания электронных переходов и кинетической устойчивости синтезированных соединений были рассчитаны энергии самой высокой занятой молекулярной орбитали (EHOMO) и самой низкой свободной молекулярной орбитали (ELUMO). Для каждой производной был определен энергетический разрыв между молекулярными орбитами (ΔE = ELUMO - EHOMO) для оценки её химической реактивности и кинетической стабильности, при этом меньший энергетический разрыв указывал на более высокую реактивность.
Кроме того, глобальные дескрипторы химической реактивности — включая химическую твёрдость (η), мягкость (S), химический потенциал (μ), электроотрицательность (X) и индекс электрофильности (w) — вычислялись на основе орбитальных энергий границы с использованием следующих стандартных операционныхуравнений 22:
(7)
(8)
(9)
(10)
Вычисленные электронные свойства и глобальные дескрипторы систематически сравнивались по всей серии для установления связей структуры и активности (SAR), отражающих наблюдаемые профили α-глюкозидазы и уреазы.
Кроме того, были созданы и визуализированы карты молекулярного электростатического потенциала (MESP) с помощью VESTA для картирования распределения заряда по оптимизированным молекулярным каркастам. Поверхности MESP позволяли выявлять электрофильные (электронодефицитные, предпочитаемые нуклеофильными атаками) и нуклеофильные (электронно-богатые, благоприятные для электрофильных атак) областей. Это распределение электростатического потенциала в конечном итоге коррелировало с специфическими нековалентными связывающими взаимодействиями, такими как водородные связи и электростатические контакты, наблюдаемые в молекулярных стыковочных симуляциях.
Анализ ADME
Для оценки похождости на препарат и биофармацевтических свойств синтезированных производных на основе триазола и оксадиазола на основе ацетамида 8(a–c) и 9(d–f) был проведён анализ абсорбции, распределения, метаболизма и экскреции in silico (ADME). Упрощённые молекулярно-входные системы линейного ввода (SMILES) представления всех синтезированных соединений, вместе с эталонными ингибиторами акарбозой и тиоуреи, были сгенерированы с помощью программного обеспечения для молекулярного моделирования и отправлены на веб-сервер SwissADME23.
Рассчитанные физико-химические параметры включали молекулярную массу (MW), липофильность (logP), топологическую полярную поверхность (TPSA), количество доноров водородных связей (HBD), количество акцепторов водородных связей (HBA) и количество вращающихся связей (RB).
Фармакокинетические свойства оценивались путём прогнозирования всасывания желудочно-кишечного тракта (ГИ) и проницаемости гематоэнцефалического барьера (BBB). Кроме того, соединения оценивались на соответствие правилу пяти Липинского и другим установленным критериям подобия препарата для оценки их пероральной биодоступности. Наконец, для каждого соединения рассчитывался балл синтетической доступности (SA) по шкале от 1 (легко синтезировать) до 10 (очень сложно синтезировать). Прогнозируемые профили ADME и подобия препаратов синтезированных соединений сравнивались с профилями референсных ингибиторов для оценки их потенциала как ведущих кандидатов для разработки лекарств.