Scientific journal

International Journal of Applied and fundamental research

ISSN 1996-3955

ИФ РИНЦ = 0,593

Исследование в области виртуального скрининга занимает длительное время, что негативно сказывается на проведении поиска новых потенциальных лекарственных средств. Одним из путей реализации цели виртуального скрининга является построение количественной зависимости «структура – активность».

Точное предсказание биологической активности больших наборов разнообразных соединений является чрезвычайно сложной задачей. Исследования в моделировании биологической активности помогают сделать вычислительный прогноз, они необходимы для открытия новых лекарств. Каждая функциональная зависимость биологической активности от структуры принимает заранее определенную функциональную форму, основанную на теории, для взаимосвязи между переменными, характеризующими комплекс, которые также включают параметры, соответствующие экспериментальным или имитационным данным, и их прогнозируемое сродство к привязке. Неотъемлемая проблема этого функционального подхода заключается в том, что он приводит к плохой предсказуемости для тех соединений, которые не соответствуют допущениям моделирования. Более того, стратегии повторной выборки, такие как перекрестная проверка, все еще систематически не используются для защиты от чрезмерного соответствия калибровочных данных при оценке параметров модели «структура – активность» [1].

Проведение компьютерного моделирования в области исследования зависимости «структура – активность» является актуальным, позволяет создавать прогнозные модели молекулярных систем, как, например, прогнозирования биологической активности. При проведении исследований количественной зависимости биологической активности от структуры необходимо учитывать параметры, оказывающие влияние на проникновение веществ в организм. К таким параметрам относятся липофильность, численно обозначается константой липофильности (log P) и ионизация с помощью констант ионизации (кислотности (рКа) и основности (рКв)).

Исследование соединений ряда производных антраниловой кислоты (а. к.) позволяет обнаружить вещества с противомикробной [2, 3], противовоспалительной и анальгетической активностями. Выделены а. к. из морских грибов Penicillium paneum SD-44 с антимикробной активностью [4].

Исследование зависимости «структура – активность» по прогнозированию противомикробной активности по отношению к E. coli позволит проводить целенаправленный синтез веществ с выраженной противомикробной активностью. Использование при прогнозировании фармакокинетических параметров: logP, pKa и рКв – позволяет проводить построение моделей биологической активности с учётом особенностей всасывания и распределения соединений в организме [5].

Цель данной работы связана с количественным моделированием антимикробной активности по отношению к E. coli (АМА К.П.) от структуры «структура – активность» с использованием констант липофильности (log P рассч) и констант ионизации (pKa рассч и pKв рассч).

Объектом проводимого исследования являются N-замещенные производные антраниловой кислоты (18 соединений).

X = Br, Y = H, R = C6H5 (1); X = Br, Y = H, R = 2–OCH3C6H4 (2); X = Br, Y = H, R = 3–NO2C6H4 (3); X = Br, Y = H, R = 4–BrC6H4 (4); X = Br, Y = H, R = CH2C6H5 (5); X = I, Y = H, R = CH = CHCOOH (6); X = I, Y = H, R = C2H4Cl (7); X = I, Y = H, R = C3H6Cl (8); X = I, Y = H, R = 3–NO2C6H4 (9); X = I, Y = H, R = 4–BrC6H4 (10); X = I, Y = H, R = 4–CH3C6H4 (11); X = I, Y = H, R = CH2C6H5 (12); X = Br, Y = Br, R = CH2C6H5 (13); X = Br, Y = Br, R = CH2CH2CH2Cl (14); X = Br, Y = H, R = CH2CH2CH2Cl (15); X = Br, Y = H, R = адамантил (16); X = Br, Y = Br, R = адамантил (17); X = I, Y = H, R = адамантил (18)

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

Математическое регрессионное моделирование проводили программой Statistica 6.

Теоретические значения log Pрассч, pKa рассч и pKв рассч 24 соединений получены с использованием программы AK_QSAR (logP, pKa, pKb) [6] (табл. 1, 3).

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

Определены экспериментальные значения противомикробной активности (E. coli) 18 соединений из ряда N-замещенных антраниловых кислот [7]. Для изучаемых производных определены различные значения МПК, от 2 до 1500 мкг/мл (табл. 1). Полученный результат по противомикробной активности: МПК (мкг/мл) переводили в логарифм: log(1/МПК) E. coli [8] для выравнивания полученных экспериментальных результатов и получения дескриптора биологической активности – антимикробная активность по кишечной палочке (log АМА К.П.).

Известно много программ прогнозирования констант липофильности, например программа Ligand Scout, которую использовали для изучения структуры производных N-замещённых антраниловых кислот (рис. 1).

В результате изучения структуры N-замещённого производного антраниловой кислоты определены гидрофобные фрагменты: атом галогена (бром) и два ароматических кольца. Далее проводили квантово-химические расчёты и полученные результаты использовали при проведении прогноза свойств программой AK_QSAR (logP, pKa, pKв).

Таблица 1

МПКэксп. (E. coli), log АМА К.П. эксп., log P рассч, pKa рассч., pKв рассч. N-замещенных антраниловых кислот (18 соединений)

С помощью программного комплекса AK_QSAR (logP, pKa, pKв) рассчитаны значения log P рассч, pKa рассч и pKв рассч 18 соединений (табл. 1).

Рассчитанные значения log P рассч находятся в интервале от 0,25 у соединения 6 до 3,06 (соединение 2). Значения pKaрассч (констант кислотности) находятся в пределах от 3,26 (соединение 17) до 7,33, для соединения 9. Константы основности (pKврассч) в интервале от 12,43 у соединения 6 до 14,63, для соединения 17.

Выполнены исследования зависимости биологической активности от структуры «структура – активность», в качестве изучаемой активности использовали log АМА К.П. эксп. Взаимосвязь структурных дескрипторов (log P рассч, pKa рассч, pKв рассч) с log АМА К.П. эксп. изучали множественным линейным регрессионным анализом программой Statistica 6.

По результатам проведённого анализа было получено свыше 25 уравнений, из которых выбрали два наиболее значимых для прогнозирования log АМА К.П. эксп. (табл. 2).

Значимость полученных уравнений регрессии проверяли значениями коэффициента множественной регрессии (R), среднеквадратической ошибки (S), критерия Фишера (F) и Q2LOO.

Коэффициент множественной регрессии зависимости структурных параметров от log АМА К.П. эксп. лежит в диапазоне 0,766–0,769. Критерий Фишера составляет от 6,76 до 10,70. Ошибка прогнозирования находится в диапазоне 0,47–0,48. Оценка с использованием программы Statographics, по значению коэффициента детерминации предсказаний Q2LOO по уравнению 2, выше и равна 0,33.

Проведена валидация уравнения 2, для расчёта АМА К.П. на примерах шести соединений 19–24 из ряда N-замещенных антраниловых кислот (табл. 3).

X = Br, Y = H, R = 2–фурил (19); X = I, Y = H, R = C3H7 (20); X = I, Y = H, R = C6H5 (21); X = I, Y = H, R = 4–ClC6H4 (22); X = Br, Y = Br, R = CH2Cl (23); X = Br, Y = H, R = CH2Cl (24)

Результаты определения АМА К.П. эксп. соединений 19–24 приведены в табл. 3.

Оценку точности прогнозирования МПКE. coli по полученному уравнению 2, проводили с использованием коэффициента корреляции (Rпрогн.) прогнозируемых значений (log АМА К.П. рассч.) с экспериментальными (log АМА К.П. эксп.) (табл. 3), при этом log АМА К.П. переводили в МПК (10log АМА К.П.).

На рис. 2 приведена линейная диаграмма с Rпрогн. = 0,758.

Заключение

Проведены математические исследования в области количественного моделирования антимикробной активности с использованием зависимости предложенного дескриптора биологической активности log АМА К.П. от log P рассч, pKa рассч, pKв рассч N-замещенных производных антраниловой кислоты.

Из найденных двух уравнений регрессионного моделирования log АМА К.П. выделено уравнение 2 с наибольшими значениями статистических параметров, характеризующих значимость модели.

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

Высокие результаты проверки позволяют предложить использование найденной количественной модели «структура – активность» для математического моделирования противомикробной активности при поиске соединений лидеров с антимикробной активностью в отношении E. coli производных антраниловой кислоты.