Научная статья на тему 'QSAR-МОДЕЛИРОВАНИЕ НЕКОТОРЫХ ПРОИЗВОДНЫХ УРАЦИЛА С ВЫРАЖЕННОЙ ИНГИБИТОРНОЙ АКТИВНОСТЬЮ В ОТНОШЕНИИ ДЕЗОКСИУРИДИНФОСФОРИЛАЗЫ'

QSAR-МОДЕЛИРОВАНИЕ НЕКОТОРЫХ ПРОИЗВОДНЫХ УРАЦИЛА С ВЫРАЖЕННОЙ ИНГИБИТОРНОЙ АКТИВНОСТЬЮ В ОТНОШЕНИИ ДЕЗОКСИУРИДИНФОСФОРИЛАЗЫ Текст научной статьи по специальности «Фундаментальная медицина»

CC BY
67
23
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ИНГИБИТОРЫ ДЕЗОКСИУРИДИНФОСФОРИЛАЗЫ / МОДЕЛИРОВАНИЕ QSAR / ПРОГРАММА GUSAR 2013 / MNA- И QNA-ДЕСКРИПТОРЫ / DEOXYURIDINETRIPHOSPHATASE INHIBITORS / QSAR MODELLING / GUSAR 2013 SOFTWARE / QNA-DESCRIPTORS / MNA-DESCRIPTORS

Аннотация научной статьи по фундаментальной медицине, автор научной работы — Хайруллина В.Р., Акбашева Ю.З., Янгирова К.Р., Гимадиева А.Р., Мустафин А.Г.

С использованием программы GUSAR 2013 выполнен количественный анализ взаимосвязи между структурой и эффективностью 135 ингибиторов дезоксиуридинтрифосфатазы на основе производных урацила в интервале значений IC50= 30÷185000 нмоль/л. Построено шесть статистически значимых валидных консенсус-моделей, характеризующихся высокой точностью прогноза IC50 для структур обучающих и тестовых выборок.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по фундаментальной медицине , автор научной работы — Хайруллина В.Р., Акбашева Ю.З., Янгирова К.Р., Гимадиева А.Р., Мустафин А.Г.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

QSAR-SIMULATION OF SOME URACIL DERIVATIVES WITH EXPRESSED INHIBITOR ACTIVITY IN DESOXYURIDINPHOSPHYLASE

Due to the widespread prevalence, deoxyuridinetriphosphatase (UTPase) is considered by modern biochemists and physicians as a promising substance for development of drugs with a wide range of activities. The therapeutic effect of such drugs is conditioned by suppression of DNA biosynthesis in various viruses, bacteria, and protozoa. In order to rationalize the search for new dUTPase inhibitors, Russian and foreign researchers are actively using the QSAR methodology at the stage of selection of compounds. However, the practical application of this methodology is impossible without the existence of valid QSAR models. With the use of the GUSAR 2013 software, a quantitative analysis of the relationship between the structure and efficacy of 135 deoxyuridinetriphosphatase inhibitors based on uracil derivatives was performed in the IC50 range of 30 ÷ 185 000 nmol/l. Six statistically significant valid consensus models characterized by high descriptive ability and moderate prognostic ability on the structures of training and test samples were constructed. To design valid QSAR models for dUTPase inhibitors, one can use QNA or MNA descriptors and their combinations in a consensus approach.

Текст научной работы на тему «QSAR-МОДЕЛИРОВАНИЕ НЕКОТОРЫХ ПРОИЗВОДНЫХ УРАЦИЛА С ВЫРАЖЕННОЙ ИНГИБИТОРНОЙ АКТИВНОСТЬЮ В ОТНОШЕНИИ ДЕЗОКСИУРИДИНФОСФОРИЛАЗЫ»

ISSN 1998-4812

987

раздел ХИМИЯ

УДК: 544.165

QSAR-МОДЕЛИРОВАНИЕ НЕКОТОРЫХ ПРОИЗВОДНЫХ УРАЦИЛА С ВЫРАЖЕННОЙ ИНГИБИТОРНОЙ АКТИВНОСТЬЮ В ОТНОШЕНИИ ДЕЗОКСИУРИДИНФОСФОРИЛАЗЫ

© В. Р. Хайруллина*1, Ю. З. Акбашева1, К. Р. Янгирова1, А. Р. Гимадиева2, А. Г. Мустафин1'2

1Башкирский государственный университет Россия, Республика Башкортостан, 450076 г. Уфа, ул. Заки Валиди, 32.

2Уфимский институт химии УФИЦРАН, Россия, Республика Башкортостан, 450054 г. Уфа, просп. Октября, 71.

* Veronikal979@yandex. ru

С использованием программы GUSAR 2013 выполнен количественный анализ взаимосвязи между структурой и эффективностью 135 ингибиторов дезоксиуридинтрифосфатазы на основе производныхурацила в интервале значений IC50= 30^185000 нмоль/л. Построено шесть статистически значимых валидных консенсус-моделей, характеризующихся высокой точностью прогноза IC50 для структур обучающих и тестовых выборок.

Ключевые слова: ингибиторы дезоксиуридинфосфорилазы, моделирование QSAR, программа GUSAR 2013, MNA- и QNA-дескрипторы.

Введение

Фермент дезоксиуридинтрифосфатаза (dUT-Pase, ЕС3.6.1.23) относится к семейству гидролаз и катализирует реакцию гидролиза дезоксиуридин-трифосфата с образованием дезоксиуридинмоно-фосфата и пирофосфат-аниона. Он присутствует в различных живых системах, включая организм животных и человека, а также простейшие, бактерии и вирусы. Благодаря вышеописанной биохимической реакции dUTPase выполняет две важные функции в живых системах: 1) регулирует уровень дезоксиу-ридинтрифосфата с целью снижения уровня мутаций в молекулах ДНК; 2) способствует нормальному функционированию другого жизненно важного фермента - тимидилатсинтазы, непосредственно участвующей в синтезе дезокситимидинмонофос-фата из дезоксиуридинмонофосфата [1-11]. Затем дезокситимидинмонофосфат в результате фосфо-рилирования превращается в дезоксиуридинтри-фосфат, представляющий собой одно из мономерных звеньев молекул ДНК [3-11].

В связи с повсеместной распространенностью dUTPase рассматривается современными биохимиками и медиками как перспективная мишень для разработки противоопухолевых, противовирусных (в т.ч. и против вируса СПИДа), антибактериальных, противомалярийных лекарственных препаратов, обладающих широким спектром активности. Терапевтическое действие этих препаратов будет осуществляться благодаря подавлению ими биосинтеза ДНК в различных вирусах, бактериях и простейших. Кроме того, актуальность поиска эффективных и вместе с тем специфических ингибиторов dUTPase обусловлена и тем фактом, что этот фермент способен гидролизовать метаболит противоопухолевого лекарственного препарата 5-фтор-урацила - 5-фтордезоксиуридинтрифосфата с образованием 5-фтордезоксиуридинмонофосфата. Эта нежелательная реакция со временем вызывает резистентность к 5-фторурацилу у онкологических

больных, значительно снижая эффективность химиотерапии вплоть до летального исхода. Таким образом, поиск эффективных ингибиторов dUTPase является не только фундаментальной, но и практически важной задачей медицинской химии и биохимии [1; 12-14].

Однако рациональный поиск новых потенциальных лекарств невозможен без привлечения различных методов виртуального скрининга. Одними из перспективных подходов в этом направлении рассматриваются методы QSAR, которые позволяют осуществлять экспрессный отбор структур с заданным комплексом свойств для доклинических испытаний. Однако практическая реализация такого подхода не представляется возможной без наличия статистически значимых моделей QSAR, ориентированных на прогноз ферментативной и рецеп-торной специфичности потенциальных прототипов лекарств [15-21]. В связи с этим, целью настоящей работы было построение и валидация моделей QSAR, ориентированных на прогноз численных значений параметра ГС50 для ингибиторов дезоксиуридинтрифосфатазы (dUTPase) в ряду некоторых производных урацила с общей структурной

I

Рис. 1. Общая структурная формула ингибиторов UPhase, вошедших в обучающие массивы ОВ1 и ОВ2, где п = 13; Х = O, Б, N Y - линейная или разветвленная алкиле-новая группа с С[-С8; 7 = -БО^'Я2 или -Ж3802-Я4, где Я1 и Я2 - алкил с С^С6 или незамещенный аралкил; Я3 - алкил с СГС6; Я4 - арил или ненасыщенная гетероциклическая группа.

Анализ влияния структурных фрагментов на эффективность ингибирования каталитической активности уридинфосфорилазы и дезоксиуридинтрифосфатазы

Рис. 2. Общая идеология построения моделей QSAR М1-М6.

Рис. 3. Формирование обучающих и тестовых

Методическая часть

Изучение количественной связи «структура-активность ингибиторов dUTPase» проводили с использованием компьютерной программы GUSAR 2013 (General Unrestricted Structure Activity

выборок для построения моделей QSAR М1-М6.

Relationships) [15-21]. Модели QSAR строили в несколько этапов, краткое перечисление которых представлено на рис. 2. Идеология формирования и отбора дескрипторов для построения моделей QSAR описана в работах [15-21].

1. Формирование обучающих и тестовых выборок. Обучающие ОВ1-ОВ2 и тестовые выборки ТВ1-ТВ2 для построения моделей М1-М6 и проверки их валидности формировали в соответствии со схемой, представленной на рис. 3, на основе массива структур 135 соединений MS1, предварительно сформированного на основе структур производных урацила и соответствующих им значений IC50 в отношении UPhase, взятых из работ [12-14]. Как видно из рис. 3, структуры моделируемых соединений обладали разной степенью ингибиторно-го действия в отношении dUPhase человека.

Предварительно перед построением моделей, значения IC50 из нмоль/л пересчитывали в моль/л и переводили в р1С50 по формуле:

рГС50 = -lg IC50 (моль/л),

Проверку валидности итоговых моделей М1-М6 проводили на структурах тестовых выборок ТВ1-ТВ2.

2. Построение моделей QSAR. Модели QSAR М1-М6 для количественного прогноза ингибитор-ной активности производных урацила в отношении dUThase человека строили на основе двух типов подструктурных дескрипторов атомных окрестностей: QNA (Quantitative Neighbourhoods of Atoms) и MNA (Multilevel Neighbourhoods of Atoms) [15-21]. Расчет этих типов дескрипторов в программе GUSAR 2013 осуществлялся автоматически из структурных формул химических соединений с учетом валентности и парциальных зарядов атомов. Специфика типов связи при расчетах не учитывалась. Идеология расчета QNA- и MNA-дескрипторов детально описана в статьях [15-21]. В дополнение к этому для повышения описательной и прогностической способности моделей использовали три дескриптора, описывающих молекулу целиком (топологическую длину, топологический объем и липофильность), расчет которых также осуществлялся в выбранной программе автоматически.

Для сокращения дескрипторного пространства и отбора наиболее значимых дескрипторов использовали метод RBF-SCR, который сочетает в себе преимущества как интерполяции RBF (Radial Basis Function), так и методов самосогласованной регрессии. Детальное описание этого метода представлено в работе [25]. Дополнительно проверку адекватности построенных моделей осуществляли с использованием процедуры скользящего контроля с 20-тикратным рандомизированным выбросом 20% соединений из обучающих выборок ОВ1 и ОВ2. Обе процедуры в программе GUSAR 2013 реализованы автоматически [15-21].

Каждая из четырех итоговых консенсус моделей QSAR М1-М2, М4-М5 включала в себя по 20 частных регрессионных зависимостей, объединенных автоматически исходя из принципа их общего сходства. Несмотря на то, что консенсус модели М3 и М6 строились по тому же принципу, они включали в себя по 320 парциальных регрессионных зависимостей. При этом каждая из этих 320 частных моделей, входящих в консенсус-

модели М3 и М6, строилась независимо друг от друга на основании либо QNA-, либо MNA-дескрипторов. Однако регрессионные уравнения, построенные на основе QNA- и MNA-дескрпторов, не поддаются однозначной физической интерпретации, и по этой причине регрессионные уравнения, построенные на их основе, явном виде в программе GUSAR 2013 не отображаются. Итоговый прогноз численного значения pIC50 для конкретного соединения формировался на основе результатов усреднения прогнозных значений pIC50 частных регрессионных моделей QSAR, входящих в используемую консенсус-модель.

3. Оценка описательной и прогностической способности моделей QSAR. Описательную способность моделей QSAR М1-М6, а также наличие систематических ошибок при их расчетах оценивали по результатам прогноза параметра IC50 для структур обучающих выборок ОВ1-ОВ2 с использованием метрики на основе коэффициентов детерминации R2 (R2test, R20, R^js;, ССС). Прогностическую способность этих же моделей оценивали по результатам прогноза параметра IC50 для структур внешней и внутренней тестовых выборок TB1 и TB2 с использованием двух категорий метрик:

1) по метрикам на основе коэффициентов детерминации R2 (R2test, R20, Q2(f:), Q2(f2), Ri^s , ССС); и

2) по метрикам, оценивающим ошибки прогноза значений pIC50 (среднеквадратичная ошибка (RMSEP), средняя абсолютная ошибка (MAE), стандартное отклонение (S.D.)) [22-24]. Расчет этих статистических параметров проводили с использованием программы Xternal Validation Plus 1.2 [26].

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

На основе консенсус-подхода, реализованного в программе GUSAR 2013, проведено моделирование количественной взаимосвязи между строением и активностью ингибиторов dUTPase с общей структурной формулой I (рис. 1), содержащимися в обучающих выборках OB1-OB2. В результате, в зависимости от типа использованных при расчетах дескрипторов, для каждой из этих обучающих выборок построено по три консенсус-модели QSAR. Описательная и прогностическая способность этих моделей оценена на структурах обучающих выборок ОВ1-ОВ2 при скользящем контроле с 20-тикратным рандомизированным исключением 20% соединений, а также на структурах тестовых выборок ТВ1-ТВ2. Критерии, на основе которых оценивали описательную и прогностическую способность моделей М1-М6, представлены в табл. 1. Следует отметить, что различные коэффициенты детерминации, величины стандартного отклонения и критерия Фишера, представленные в табл. 1, представляют собой усредненные значения, полученные с учетом всех частных регрессионных моделей, входящих в консенсус модель M,-.

Таблица 1.

Статистические характеристики и оценка показателей точности предсказаний значений р1С50 для ингибиторов UPhase по

консенсус-моделям М1-М6 при использовании метода БСЯ-ЯВЕ

Обучающая выборка

Модель

N

Nmod

RT

QO

RT

S.D.,

V

ОВ1 ОВ2

ОВ1 ОВ2

ОВ1 ОВ2

Модели QSAR на основе QNA-дескрипторов М1 113 T0 0.975 0.661

М4 94 T0 0.985 0.6T9

Модели QSAR на основе MNA-дескрипторов М2 113 T0 0.965 0.7T4

М5 94 T0 0.985 0.671

Модели QSAR на основе QNA- и MNA-дескрипторов

М3 М6

113

94

3T0 3T0

0.973 0.985

0.718 0.681

0.7T18 0.7791

0.6494 0.7T11

0.68TT 0.781

16.053 13.684

11.975 9.988

1T.665 11.633

0.477 0.494

0.4T0 0.459

0.4T4 0.45T

15 1T

T1

16

19 14

F

Рис. 4 Статистические параметры, оцененные по результатам прогноза параметра р1С50 для ингибиторов dUTPase, содержащихся в тестовых выборках ТВ1 (а) и ТВ2 (б), по консенсус-моделям QSAR М1-М6, построенных с использованием метода отбора дескрипторов БСЯ-ЯВЕ.

Обозначения: N - число структур в обучающей выборке; ^mod - число частных регрессионных моделей, входящих в консенсус модель М,-; Я2ов - коэффициент детерминации для соединений из обучающей выборки; Я2тв - коэффициент детерминации для соединений из тестовой выборки; Q2^ - коэффициент корреляции, рассчитанный на обучающей выборке при скользящем контроле с исключением по одному; F - критерий Фишера; S.D. - стандартное отклонение; V -число переменных в конечном регрессионном уравнении.

В соответствии с общепринятыми рекомендациями ведущих исследователей в области QSAR-моделирования [22-24] полагали, что модель М,-обладала высокой прогностической способностью, если одновременно выполнялись следующие четыре условия: 1) численные значения коэффициентов детерминации R , R 0, R '0, Q F1, Q F2 и критерия CCC были близки между собой и стремились к

единице; 2) численное значение критерия R~m > 0.5 при LR2m < 0.2; 3) численное значение средней абсолютной ошибки МАЕ не превышало 10% от интервала активности соединений 1р1С50, моделируемой обучающей выборки ÜB^ 4) величина параметра MAE + 3*S.D. (где S.D. - стандартное отклонение) не превышала 20% от интервала активности соединений JpIC50, моделируемой обучающей выборки ÜBj.

Прогностическая способность модели М,- оценивалась как низкая, если одновременно: 1) ленные значения коэффициентов детерминации R'

чис-

T

RT0,

RT'0, Qtf1, Qtft и критерия CCC не превышали

порогового значения 0.6; 2) численное значение критерия Я2т < 0.5 при - Я2т < 0.2; 3) численное значение средней абсолютной ошибки МАЕ превышало 15% от интервала активности соединений 1р1С50, моделируемой обучающей выборки ОВ1; 4) величина параметра МАЕ + З^.Б. превышала 25% от интервала активности соединений ^Юзо, моде-

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

В ходе выполненных исследований установлено, что все модели QSAR обладали довольно высокой описательной и умеренной прогностической способностью, табл. 1. В качестве примера на рис. 4 (а, б) представлены статистические параметры, оцененные по результатам прогноза параметра pIC50 для ингибиторов dUTPase, содержащихся в тестовых выборках ТВ1 и ТВ2, по консенсус-моделям QSAR М3 и М6 соответственно.

Незначительная разница между численными значениями разных типов коэффициентов детерминации в сочетании с приемлемыми значениями параметров МАЕ и MAE + 3*S.D., представленными в табл. 1 и на рис. 4, свидетельствует о том, что строить валидные модели QSAR, ориентированные на прогноз параметра pIC50 для ингибиторов dUT-Pase среди производных урацила с общей структурной формулой I (рис. 1), [18-22] можно как с использованием одного определенного типа дескрипторов (QNA- или MNA-дескрипторов), так и их сочетания в консенсус-подходе.

ЛИТЕРАТУРА

1. Toth J., Varga B., Kovacs M., Malnasi-Csizmadia A., Vertes-sy B.G. Kinetic Mechanism of Human dUTPase, an Essential Nucleotide Pyrophosphatase Enzyme // J. Biol. Chem. 2007. V. 282. P. 33572-33582.

2. Shlomai J., Kornberg, A. Deoxyuridine triphosphatase of Escherichia coli: purification, properties and use as a reagent to reduce uracil incorporation into DNA // J. Biol. Chem. 1978. V. 25. P. 3305-3312.

3. Gadsden M. H., Mcintosh E. M., Game J. C., Wilson P. J., Haynes R. H. dUTP pyrophosphatase is an essential enzyme in Saccharo-myces cerevisiae // EMBO J. 1993. V. 12. P. 4425-4431.

4. El-Hajj H. H., Zhang H., Weiss B. Lethality of a Dut (Deoxyuridine Triphosphatase) Mutation in Escherichia coli // J. Bacteriol. 1988. V. 170. P. 1069-1075.

5. Traut T. W. Physiological concentrations of purines and pyrimidines // J. Mol. Cell Biochem. 1994. V. 140. P. 1-22.

6. Pearl L. H., Savva R. The problem with pyrimidines // J. Nat. Struct. Biol. 1996 V. 3. P. 485-487.

7. Pugacheva E. N., Ivanov A. V., Kravchenko J. E., Kop-nin B. P., Levine A. J., Chumakov P. M. Novel gain of function activity of p23 mutants: activation of the dUTRPase gene expression leading to resistance to 5-fluorouracil // J. Oncogene 2002. V. 21. P. 4595-4600.

8. Chano T., Mori K., Scotlandi K., Benini S., Lapucci C., Manara M. C., Serra M., Picci P., Okabe H., Baldini N. Differentially expressed genes in multidrug resistant variants of U-2 OS human osteosarcoma cells // J. Oncol. Rep. 2004. V. 11. P. 1257-1263.

9. Ladner R. D., Lynch F. J., Groshen S., Xiong Y. P., Sher-rod A., Caradonna S. J., Stoehlmacher J., Lenz H. J. dUTP nucleotidohydrolase isoform expression in normal and neop-lastic tissues: association with survival and response to 5-fluorouracil in colorectal cancer // J. Cancer Res. 2000. V. 60. P. 3493-3503.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

10. Romeike B. F., Bockeler A., Kremmer E., Sommer P., Krick C., Grasser F. Immunohistochemical detection of dUTPase in intracranial tumors // J. Pathol. Res. Pract. 2005. V. 201. P. 727-732.

11. Vassylyev D.G., Morikawa K. Precluding uracil from DNA // J. Structure 1996. V. 4. P. 1381-1385.

12. Miyakoshi H., Miyahara S., Yokogawa T., Chong Kh.T., Taguchi J., Endoh K., Yano W., Wakasa T., Ueno H., Ta-kao Y., Nomura M., Shuto S., Nagasawa H., Fukuoka M. Synthesis and Discovery of N-Carbonylpyrrolidine- or N-Sulfonylpyrrolidine-Containing Uracil Derivatives as Potent Human Deoxyuridine Triphosphatase Inhibitors // J. Med. Chem. 2012. V. 55. P. 2960-2969.

13. Miyahara S., Miyakoshi H., Yokogawa T., Chong Kh.T., Taguchi J., Muto T., Endoh K., Yano W., Wakasa T., Ueno H., Takao Y., Fujioka A., Hashimoto A., Itou K., Yamamura K., Nomura M., Nagasawa H., Shuto S., Fukuoka M. Discovery of a Novel Class of Potent Human Deoxyuridine Triphos-phatase Inhibitors Remarkably Enhancing the Antitumor Activity of Thymidylate Synthase Inhibitors // J. Med. Chem. 2012. V. 55. P. 2970-2980.

14. Патент №8530490 B2 US, Appl. №12/996,079. Uracil compound or salt thereof having human deoxyuridine triphospha-tase inhibitory activity / Fukuoka M., Yokogawa T., Miyahara S., Miyakoshi H., Yano W., Taguchi J., Takao Y. Khayrullina V.R., Gerchikov A.Ya., Lagunin A.A., Zarudii F.S. Quantitative Analysis of Structure-Activity Relationships of Tetrahy-dro-2H-isoindole Cyclooxygenase-2 Inhibitors V// J. Biokhi-miya. 2015. V. 80. P. 96-110.

15. Хайруллина В. Р., Герчиков А. Я., Зарудий Ф. С. Анализ взаимосвязи «структура-ингибирующая активность циклоок-сигеназы-2» в ряду производных ди-трет-бутилфенола, тиа-золона и оксазолона // Вестник Башкирского государственного университета. 2014. Т. 19. №2. С. 417-422.

16. Zakharov A. V., Lagunin A. A., Filimonov D. A., Poroi-kov V. V. Quantitative prediction of antitarget interaction profiles for chemical compounds // Chemical Research in Toxicology. 2012. V. 25. №11. P. 2378-2385.

17. Хайруллина В. Р., Герчиков А. Я., Лагунин А. А., Зарудий Ф. С. Количественный анализ взаимосвязи «структура-активность» ингибиторов циклооксигеназы-2 среди производных тетрагидро-2Н-изоиндола // Биохимия. 2015. Т. 80. С. 96-110.

18. Lagunin A., Zakharov A., Filimonov D., Poroikov V. QSAR Modelling of Rat Acute Toxicity on the Basis of PASS Prediction // Molecular Informatics. 2011. V. 30. P. 241-250.

19. Filimonov D. A., Zakharov A. V., Lagunin A. A., Poroikov V. V. QNA based «Star Track» QSAR approach // SAR and QSAR in Environmental Research. 2009. V. 20. №7-8. P. 679-709.

20. Masanda V. H., Mahajana D. T., Patil K. N., Dawale N. E., Hadda T. B., Alafeefy A. A., Chinchkhede K. D. General Unrestricted Structure Activity Relationships based evaluation of quinoxaline derivatives as potential influenza NS1A protein inhibitors // Der PharmaChemica. 2011. V. 3. P. 517-525.

21. Roy K., Mitra I., Kar S., Ojha P.K., Das R.N., Kabir H. Comparative Studies on Some Metrics for External Validation of QSPR Models // J. Chem. Inf. Model. 2012. V. 52. P. 396-408.

22. Roy P. P., Paul S., Mitra I., Roy K. On Two Novel Parameters for Validation of Predictive QSAR Models // J. Molecules. 2009. V. 14. P. 1660-1701.

23. Roy K., Das R. N., Ambure P., Aher R. B. Be aware of error measures. Further studies on validation of predictive QSAR models // J. Chemom. Intell. Lab. Syst. 2016. V. 152. P. 18-33.

24. Zakharov A. V., Peach M. L., Sitzmann M., Nicklaus M. C. A New Approach to Radial basis function approximation and Its application to QSAR // J. Chem. Inf. Model. 2014. V. 54. P. 713-719.

25. Xternal Validation Plus, version 1.2, 2016, DTC Lab., Kolka-ta, India. URL: https://sites.google.com/site/dtclabxvplus/

Поступила в редакцию 25.07.2018 г.

99T

ХИМИЯ

QSAR-SIMULATION OF SOME URACIL DERIVATIVES WITH EXPRESSED INHIBITOR ACTIVITY IN DESOXYURIDINPHOSPHYLASE

© V. R. KhairuUina1*, Yu. Z. Akbasheva1, K. R. Jangirova1, A. R. Gimadieva2, A. G. Mustafin1,2

1Bashkir State University 32 Zaki Validi Street, 450076 Ufa, Republic of Bashkortostan, Russia.

2Ufa Institute of Chemistry, Ufa Federal Research Center ofRAS 69 Oktyabrya Avenue, 450054 Ufa, Republic of Bashkortostan, Russia.

*Email: veronika1979@yandex.ru

Due to the widespread prevalence, deoxyuridinetriphosphatase (UTPase) is considered by modern biochemists and physicians as a promising substance for development of drugs with a wide range of activities. The therapeutic effect of such drugs is conditioned by suppression of DNA biosynthesis in various viruses, bacteria, and protozoa. In order to rationalize the search for new dUTPase inhibitors, Russian and foreign researchers are actively using the QSAR methodology at the stage of selection of compounds. However, the practical application of this methodology is impossible without the existence of valid QSAR models. With the use of the GUSAR 2013 software, a quantitative analysis of the relationship between the structure and efficacy of 135 deoxyuridinetriphosphatase inhibitors based on uracil derivatives was performed in the IC50 range of 30 ^ 185 000 nmol/l. Six statistically significant valid consensus models characterized by high descriptive ability and moderate prognostic ability on the structures of training and test samples were constructed. To design valid QSAR models for dUTPase inhibitors, one can use QNA or MNA descriptors and their combinations in a consensus approach.

Keywords: deoxyuridinetriphosphatase inhibitors, QSAR modelling, GUSAR 2013 software, QNA-descriptors, MNA-descriptors.

Published in Russian. Do not hesitate to contact us at bulletin_bsu@mail.ru if you need translation of the article.

REFERENCES

1. Toth J., Varga B., Kovacs M. J. Biol. Chem. 2007. Vol. 282. Pp. 33572-33582.

2. Shlomai J., Kornberg, A. J. Biol. Chem. 1978. Vol. 25. Pp. 3305-3312.

3. Gadsden M. H., Mcintosh E. M., Game J. C., Wilson P. J., Haynes R. H. EMBO J. 1993. Vol. 12. Pp. 4425-4431.

4. El-Hajj H. H., Zhang H., Weiss B. J. Bacteriol. 1988. Vol. 170. Pp. 1069-1075.

5. Traut T. W. J. Mol. Cell Biochem. 1994. Vol. 140. Pp. 1-22.

6. Pearl L. H., Savva R. J. Nat. Struct. Biol. 1996 V. 3. Pp. 485-487.

7. Pugacheva E. N., Ivanov A. V., Kravchenko J. E., Kopnin B. P., Levine A. J., Chumakov P. M. J. Oncogene 2002. Vol. 21. Pp. 45954600.

8. Chano T., Mori K., Scotlandi K., Benini S., Lapucci C., Manara M. C., Serra M., Picci P., Okabe H., Baldini N. J. Oncol. Rep. 2004. Vol. 11. Pp. 1257-1263.

9. Ladner R. D., Lynch F. J., Groshen S., Xiong Y. P., Sherrod A., Caradonna S. J., Stoehlmacher J., Lenz H. J. J. Cancer Res. 2000. Vol. 60. Pp. 3493-3503.

10. Romeike B. F., Bockeler A., Kremmer E., Sommer P., Krick C., Grasser F. J. Pathol. Res. Pract. 2005. Vol. 201. Pp. 727-732.

11. Vassylyev D.G., Morikawa K. J. Structure 1996. Vol. 4. Pp. 1381-1385.

12. Miyakoshi H., Miyahara S., Yokogawa T., Chong Kh.T., Taguchi J., Endoh K., Yano W., Wakasa T., Ueno H., Takao Y., Nomura M., Shuto S., Nagasawa H., Fukuoka M. J. Med. Chem. 2012. Vol. 55. Pp. 2960-2969.

13. Miyahara S., Miyakoshi H., Yokogawa T., Chong Kh.T., Taguchi J., Muto T., Endoh K., Yano W., Wakasa T., Ueno H., Takao Y., Fujioka A., Hashimoto A., Itou K., Yamamura K., Nomura M., Nagasawa H., Shuto S., Fukuoka M. J. Med. Chem. 2012. Vol. 55. Pp. 2970-2980.

14. Patent No. 8530490 B2 US, Appl. No. 12/996,079. Uracil compound or salt thereof having human deoxyuridine triphosphatase inhibitory activity / Fukuoka M., Yokogawa T., Miyahara S., Miyakoshi H., Yano W., Taguchi J., Takao Y. Khayrullina V.R., Gerchikov A.Ya., Lagunin A.A., Zarudii F.S. J. Biokhimiya. 2015. Vol. 80. Pp. 96-110.

15. Khairullina V. R., Gerchikov A. Ya., Zarudii F. S. Vestnik Bashkirskogo gosudarstvennogo universiteta. 2014. Vol. 19. No. 2. Pp. 417422.

16. Zakharov A. V., Lagunin A. A., Filimonov D. A., Poroikov V. V. Chemical Research in Toxicology. 2012. Vol. 25. No. 11. Pp. 23782385.

ISSN 1998-4812

BeciHHK EamKHpcKoro yHHBepcHTeTa. 2018. T. 23. №4

993

17. Khairullina V. R., Gerchikov A. Ya., Lagunin A. A., Zarudii F. S. Biokhimiya. 2015. Vol. 80. Pp. 96-110.

18. Lagunin A., Zakharov A., Filimonov D., Poroikov V. Molecular Informatics. 2011. Vol. 30. Pp. 241-250.

19. Filimonov D. A., Zakharov A. V., Lagunin A. A., Poroikov V. V. SAR and QSAR in Environmental Research. 2009. Vol. 20. No. 7-8. Pp. 679-709.

20. Masanda V. H., Mahajana D. T., Patil K. N., Dawale N. E., Hadda T. B., Alafeefy A. A., Chinchkhede K. D. Der PharmaChemica. 2011. Vol. 3. Pp. 517-525.

21. Roy K., Mitra I., Kar S., Ojha P.K., Das R.N., Kabir H. J. Chem. Inf. Model. 2012. Vol. 52. Pp. 396-408.

22. Roy P. P., Paul S., Mitra I., Roy K. J. Molecules. 2009. Vol. 14. Pp. 1660-1701.

23. Roy K., Das R. N., Ambure P., Aher R. B. J. Chemom. Intell. Lab. Syst. 2016. Vol. 152. Pp. 18-33.

24. Zakharov A. V., Peach M. L., Sitzmann M., Nicklaus M. C. J. Chem. Inf. Model. 2014. Vol. 54. Pp. 713-719.

25. Xternal Validation Plus, version 1.2, 2016, DTC Lab., Kolkata, India. URL: https://sites.google.com/site/dtclabxvplus/

Received 25.07.2018.

i Надоели баннеры? Вы всегда можете отключить рекламу.