Научная статья на тему 'КВАНТОВО-ХИМИЧЕСКИЙ АНАЛИЗ АЛКИЛЬНЫХ ПРОИЗВОДНЫХ 3,4-ДИГИДРОХИНОКСАЛИН-2(1H)-ОНА'

КВАНТОВО-ХИМИЧЕСКИЙ АНАЛИЗ АЛКИЛЬНЫХ ПРОИЗВОДНЫХ 3,4-ДИГИДРОХИНОКСАЛИН-2(1H)-ОНА Текст научной статьи по специальности «Химические науки»

CC BY
60
9
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
КВАНТОВОХИМИЧЕСКИЙ РАСЧЕТ / ПРОИЗВОДНЫЕ 3 / 4-ДИГИДРОХИНОКСАЛИН-2(1Н)-ОНА / СПЕКТРАЛЬНЫЙ АНАЛИЗ / ЭНЕРГИЯ ДИССОЦИАЦИИ

Аннотация научной статьи по химическим наукам, автор научной работы — Пешков С.А., Кунавина Е.А., Левенец Т.В., Казаев К.А.

Рассматриваются проблемы выбора метода и базиса для оптимизации соединений исследуемого ряда 3-[2 -алкил-2 -оксоэтилиден]-3,4-дигидрохиноксалин-2(1Н)-онов с помощью квантово-химических программ, а также приводятся, рассчитанные с помощью выбранного приближения, термохимические данные. Воспроизводимость параметров расчета доказана на основе сравнения практических ЯМР и ИК-спектров с теоретическими. Ввиду отсутствия термохимической информации для синтезируемых веществ, в качестве тестовых систем использовались соединения, принимающие участие в их синтезе. Доказано, что хорошие термодинамические показатели, в сравнении с экспериментом, получаются при использовании приближения PBE/jorge-dzp . Показана величина влияния молекул растворителя на энтальпии образования соединений по двум моделям: в форме континуальной модели и физического учета молекул воды. При использовании любой модели немного повышается теоретическая энтальпия образования, однако использование «физических» молекул воды в расчете приводит к большим значениям, но не сильно отличающимся от континуальной модели IEFPCM-SMD. Предложен механизм растворимости алкильных производных 3,4-дигидрохиноксалин-2(1Н)-она. Рассчитаны энергии диссоциации пяти представителей исследуемого ряда веществ. Предложен способ использования, полученных теоретически, значений Gaq и lnKsp , как относительных показателей показывающих растворимость в ряду гомологов, что очень важно при исследовании механизма связывания с молекулами белков, поскольку такие реакции протекают в водной среде

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

Похожие темы научных работ по химическим наукам , автор научной работы — Пешков С.А., Кунавина Е.А., Левенец Т.В., Казаев К.А.

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

QUANTUM-CHEMICAL ANALYSIS OF ALKYL DERIVATIVES OF 3,4-DIHYDROQUINOXALIN-2(1H)-ONE

The article discusses the problems of choosing a method and basis for the optimization of compounds of the series 3-[2 -alkyl-2 -oxoethylidene]-3,4-dihydroquinoxalin-2(1H)-ones using quantum-chemical programs, approximation of thermo-chemical data. The reproducibility of the calculation parameters has been proven by comparing practical NMR and IR spectra with theoretical ones. Due to the lack of thermochemical information for the synthesized substances, compounds participating in their synthesis were used as test systems. It has been proven that good thermodynamic performance, in comparison with experiment, is obtained using the PBE/jorge-dzp approximation. The magnitude of the influence of solvent molecules on the enthalpies of formation of compounds is shown according to two models: in the form of a continual model and a discrete approach. When using any model, the value of the theoretical enthalpy of formation increases slightly, however, the use of «physical» water molecules in the calculation leads to higher values, but not much different from the results obtained using the IEFPCM-SMD continuum model. Solubility is thought to be related to the dissociation of a proton at the amino group position. The dissociation energies of five representatives of the investigated series of substances are calculated. A method is proposed for using the theoretically obtained values of Gaq and lnKsp as relative indicators showing solubility in a series of homologues, which is very important when studying the mechanism of binding to protein molecules, since such reactions occur in an aqueous medium.

Текст научной работы на тему «КВАНТОВО-ХИМИЧЕСКИЙ АНАЛИЗ АЛКИЛЬНЫХ ПРОИЗВОДНЫХ 3,4-ДИГИДРОХИНОКСАЛИН-2(1H)-ОНА»

Раздел 02.00.17

Математическая и квантовая химия

УДК 544.15 DOI: 10.17122/bcj-2022-1-57-64

С. А. Пешков (к.х.н., доц.), Е. А. Кунавина (к.х.н., доц.), Т. В. Левенец (к.х.н., доц.), К. А. Казаев (асп.)

КВАНТОВО-ХИМИЧЕСКИЙ АНАЛИЗ АЛКИЛЬНЫХ ПРОИЗВОДНЫХ 3,4-ДИГИДРОХИНОКСАЛИН-2(1Н)-ОНА

Оренбургский государственный университет, кафедра химии 460018, г. Оренбург, пр. Победы 13; e-mail: peshkov.sergey@internet.ru

S. A. Peshkov, E. A. Kunavina, T. V. Levenets, K. A. Kazaev

QUANTUM-CHEMICAL ANALYSIS OF ALKYL DERIVATIVES OF 3,4-DIHYDROQUINOXALIN-2(1H)-ONE

Orenburg State University 13, Prospekt Pobedy Str, Orenburg, 460018, Russia; e-mail: peshkov.sergey@internet.ru

Рассматриваются проблемы выбора метода и базиса для оптимизации соединений исследуемого ряда 3-[2'-алкил-2'-оксоэтилиден]-3,4-дигидро-хиноксалин-2(1Н)-онов с помощью квантово-химических программ, а также приводятся, рассчитанные с помощью выбранного приближения, термохимические данные. Воспроизводимость параметров расчета доказана на основе сравнения практических ЯМР и ИК-спектров с теоретическими. Ввиду отсутствия термохимической информации для синтезируемых веществ, в качестве тестовых систем использовались соединения, принимающие участие в их синтезе. Доказано, что хорошие термодинамические показатели, в сравнении с экспериментом, получаются при использовании приближения РВЕ/ ]отдв-йгр. Показана величина влияния молекул растворителя на энтальпии образования соединений по двум моделям: в форме континуальной модели и физического учета молекул воды. При использовании любой модели немного повышается теоретическая энтальпия образования, однако использование «физических» молекул воды в расчете приводит к большим значениям, но не сильно отличающимся от континуальной модели ШБРСМ^МБ. Предложен механизм растворимости алкильных производных 3,4-дигидрохи-ноксалин-2(1Н)-она. Рассчитаны энергии диссоциации пяти представителей исследуемого ряда веществ. Предложен способ использования, полученных теоретически, значений АОад и ¡пКзр, как относительных показателей показывающих растворимость в ряду гомологов, что очень важно при исследовании механизма связывания с молекулами белков, поскольку такие реакции протекают в водной среде.

The article discusses the problems of choosing a method and basis for the optimization of compounds of the series 3-[2'-alkyl-2'-oxoethylidene]-3,4-dihydroquinoxalin-2(1H)-ones using quantum-chemical programs, approximation of thermo-chemical data. The reproducibility of the calculation parameters has been proven by comparing practical NMR and IR spectra with theoretical ones. Due to the lack of thermochemical information for the synthesized substances, compounds participating in their synthesis were used as test systems. It has been proven that good thermodynamic performance, in comparison with experiment, is obtained using the PBE/jorge-dzp approximation. The magnitude of the influence of solvent molecules on the enthalpies of formation of compounds is shown according to two models: in the form of a continual model and a discrete approach. When using any model, the value of the theoretical enthalpy of formation increases slightly, however, the use of «physical» water molecules in the calculation leads to higher values, but not much different from the results obtained using the IEFPCM-SMD continuum model. Solubility is thought to be related to the dissociation of a proton at the amino group position. The dissociation energies of five representatives of the investigated series of substances are calculated. A method is proposed for using the theoretically obtained values of AGaq and lnKSp as relative indicators showing solubility in a series of homologues, which is very important when studying the mechanism of binding to protein molecules, since such reactions occur in an aqueous medium.

Ключевые слова: квантово-химический расчет; производные 3,4-дигидрохиноксалин-2( 1 Н)-она; спектральный анализ; энергия диссоциации.

Key words: derivatives of 3,4-dihydro-quinoxalin-2 (1H) -one; dissociation energy; quantum-chemical calculation; spectral analysis.

Дата поступления 13.09.21

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

Ранее был осуществлен синтез 3-[2-оксо-пропилиден]-3,4-дигидрохиноксалин-2(1Н)-она двухстадийным и тандемным способами. Структура соединения подтверждена спектральными методами анализа 5. При нормальных условиях образующиеся соединение находится в твердом состоянии.

В данной статье рассматриваются проблемы выбора метода и базиса для оптимизации свойств соединений исследуемого ряда 3-[2'-алкил-2'-оксоэтилиден]-3,4-дигидрохинокса-лин-2(1Н)-онов. Теоретические данные, полученные в выбранном приближении, будут использованы для идентификации получаемых соединений и последующего молекулярного докинга. Расчет взаимодействия молекул белков или их упрощенных аналогов с 3-[2'-ал-кил-2'-оксоэтилиден]-3,4-дигидрохиноксалин-2(1Н)-онами по модели ОМ/ММ+ 6 или ОШОМ позволит определить силу и механизм связывания их с белками и оценить биологическую активность. Предварительный прогноз биологической активности может освободить исследователей от нерационального синтеза заведомо неэффективных противомикробных препаратов, что значительно ускоряет и удешевляет их поиск.

Материалы и методы

Основной проблемой квантово-химических расчетов является учет корреляционной энергии. Прямые расчеты почти всегда дают неудовлетворительные результаты. Более точным можно назвать гомодесмический подход 7 8. В нашей работе из-за сложности выбранных систем мы не можем использовать такой способ расчета. В целом точность получаемых значений энер-

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

Fk /СН, Оч ОС2Н

бензол, 1:1:1

СК /R

Схема 1. Синтез 3,4-дигидрохиноксалин-2(1Н)-онов

Для проверки было выбрано пять базисных наборов: jorge-dzp 9, cc-pvdz, aug-cc-pvdz, 6-311++G(2d,2p), 6-31++G**. Все базисы, за исключением cc-pvdz, aug-cc-pvdz, являются валентными и хорошо оптимизированы для расчетов по методу Хартри-Фока. Наборы были выбраны таким образом, чтобы в них постепенно увеличивалось количество базисных функций, и, как следствие, росли сложность и точность расчета. Расчет проводили по формулам 8:

H298 = Etotal + ZPE + TC>

где полная энергия системы;

ZPE — энергия нулевых колебаний;

TC — термическая поправка.

H

О

298

= X H

продуктов

X H

исходных'

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

Выбор базиса. Экспериментальных данных энтальпии образования диметилового эфира этандиовой кислоты реакцией этерификации метанола и щавелевой кислоты в жидкой фазе найти не удалось. В газовой фазе энергия образования равна (—709 ± 1) кДж/моль, по разным источникам 10-12. Однако есть термодинамические данные реакции гидролиза диме-тилоксалата (схема 2), в которой изменение энтальпии равно 17.428 ± 0.065 кДж/моль 12. Эти значения и были взяты для сравнения реакции, представленной на схеме 1.

о о-сьц о он

\\ / 3 \\ /

2Н20+ с-с -2НХ-ОН + с-с

^ / \\ J / \\

н3с-о о но о

Схема 2. Гидролиз диметилоксалата

Наиболее близкое к экспериментальному значение энтальпии образования было получено при использовании корреляционно-согласованного базисного набора cc-pvdz. Однако базисный набор, разработанный группой ученых Francisco, Elias, Jorge (jorge-dzp), для текущих расчетов почти не уступает в точности, но значительно превосходит по скорости

cc-pvdz. Расчеты с использованием базисных наборов, разработанных Поплом, дают немного завышенные результаты (табл. 1). Скорость расчета и сходимость в задачах оптимизации молекул для этих базисов превосходят все используемые, но энергия образования диметилоксалата, относительно экспериментальных значений, по результатам расчета, отличается на 90%, поэтому базисные наборы 6-311++0 (2d,2p), 6-31++0** для оптимизации производных 3,4-дигидрохиноксалин-2( 1 Н)-она в дальнейшем не применялись.

Оптимизация молекул производилась условно в газовой фазе, но реакция гидролиза идет в водной среде. Очевидно, что неучет молекул растворителя может повлиять на конечное значение энтальпии образования. Чтобы проверить данный факт, был сделан расчет той же реакции в приближении НГ/]отдв^гр с учетом континуальной модели ¡ЕБРСМ-БМО (табл. 2). В качестве молекул растворителя выбрана вода. Также известно, что расчет с использованием теории функционала плотности (ОБТ) демонстрирует более быструю сходимость. Данные базисные наборы разрабатывались не для метода ОБТ, но успешно применя-

Таблица 1

Абсолютные и относительные значения энтальпий молекул, участвующих в реакции гидролиза диметилоксалата, рассчитанные по методу Хартри-Фока

Базисный набор Вода Метанол Щавелев. кислота Метил-оксалат Диметил-оксалат AH1/AH2 (кДж/моль)

jorge-dzp -76.015151 -115.006900 -376.401896 -415.399436 -454.394931 -15.20/-25.04 (-7.61)

cc-pvdz -76.000260 -114.989728 -376.348341 -415.341897 -454.335203 -10.73/-20.81 (-3.38)

aug-cc-pvdz -76.015032 -115.002571 -376.377797 -415.371842 -454.365236 -17.08/-32.45 (-15.03)

6-311++G (2d,2p) -76.029917 -115.025965 -376.442218 -415.444660 -454.446288 -16.79/-31.44 (-14.01)

6-31++G** -76.004137 -114.992142 -376.330948 -415.325003 -454.319598 -15.88/-33.19 (-15.76)

Примечание. В скобках указана разница в сравнении с экспериментальным значением энтальпии. АН1 — энтальпия образования метилоксалата, АН2 — диметилоксалата.

Таблица 2

Абсолютные и относительные значения энергий молекул, участвующих в реакции гидролиза диметилоксалата, рассчитанные по методу ОРТ

Значения в приближении Вода Метанол Щавелевая кислота Диметил-оксалат AHAG (кДж/моль)

HF (142) Без учета PCM -76.015151 -115.006900 -376.401896 -454.394931 -25.04|-21.79

-76.029093 -115.016210 -376.425119 -454.409427 -26.45|-25.76

B3LYP (197) Без учета PCM -76.412910 -115.683697 -378.350074 -456.909194 -21.10|-20.09

-76.425694 -115.690628 -378.369549 -456.909194 -25.671-23.97

PBE (164) Без учета PCM -76.330059 -115.550251 -377.954432 -456.403416 -22.58|-21.31

-76.343057 -115.558528 -377.972508 -456.412811 -24.58|-22.98

Примечание. В скобках указано время в секундах, затраченное на оптимизацию и расчет гессиана для молекулы диметилоксалата.

ются и для него. Поэтому нами было решено проанализировать на точность и скорость еще два способа расчета, в приближении БЗЬУР/ ]отде-йгр и РБЕ/]отде-йгр.

Значения энергий АИ\АО (табл. 2), подтверждают, что с учетом молекул растворителя по континуальной модели 1ЕГРСМ-БМО незначительно повышается теоретическая энтальпия образования диметилоксалата — в среднем на АИ±1—2 кДж/моль, АС±2—4 кДж/моль. Для квантово-химических методов исследования

13

это не является существенным отклонением . В целом все способы расчета дают хорошие результаты, близкие к эксперименту. Более точные значения АИ\АО диметилоксалата получаются при использовании приближения РБЕ/]отде-йгр, однако при такой разнице в энергиях (1—4 кДж/моль) для корректного поиска энтальпий образования возможно применять любой метод.

Отсутствие в окружении «физических» молекул воды при оптимизации также может влиять на значение АИ\АО. Чтобы проверить это, были сделаны расчеты диметилоксалата, окруженного двумя и четырьмя молекулами воды (рис. 1). Оптимизацию проводили в двух приближениях — ИГ/]'отде-йгр и РБЕ/]'отде-йгр с учетом континуальной модели 1ЕЕРСМ-БМО, а вычисление энтальпии проводили по реакции, указанной на схеме 1, с поправкой на диметилоксалат. Таким образом, образование конечного соединения означало бы сольватацию диметилоксалата, которую можно вычислить по уравнению:

Hсольватации = Haqua — (Ндиметилоксалат + HH гО )•

Как и следовало ожидать, «физические» молекулы воды могут еще немного повысить теоретическое значение энтальпии образования диметилоксалата за счет учета энтальпии сольватации (табл. 3).

Таблица 3

Абсолютные и относительные значения энтальпий сольватации аквакомплексов диметилоксалата с учетом континуальной модели IЕРРСМ-БМО

Приближение 2Н2О 4Н2О АН сольватации

-606.469108 -758.529629 -3.93/-10.06

МР2//о/де-йгр -608.123689 -760.591479 -23.46/-45.31

РВЕ//о/де-с/гр -609.108867 -761.802302 -26.10/-45.32

РВЕ/ 6-31++С* -609.026685 -761.706621 -15.18/-27.97

Вычисленные значения энергии относительно хорошо согласуются с экспериментальными данными (—37.6±0.5 кДж/моль) 14, однако метод Хартри-Фока при расчете дает меньшие величины, чем теория функционала плотности. Молекулы воды после оптимизации в двух различных методах располагаются аналогичным образом (рис. 1). Проблема заключается в неполноте учета корреляционной энергии. Доказать это возможно, воспользовавшись теорией возмущений и проведя расчет с поправкой МР2 (табл. 3).

и

Рис. 1. Геометрическое строение аквакомплексов диметилоксалата

В целом, хорошо себя показали два базисных набора ее-рюйг и ]отде-йгр. Так как отличия в энтальпиях между ее-рюйг и ]отде-йгр не столь существенны, но время, необходимое на расчет во втором случае гораздо меньше, то выбор был сделан в пользу второго базисного набора. Обобщая проделанные вычисления реакции получения диметилоксалата, было решено использовать для дальнейших расчетов приближение РБЕ/]отде-йгр с учетом молекул растворителя по континуальной модели IЕРРСМ-БМБ.

Проверка выбранного приближения на 3-[2-оксопропилиден]-3,4-дигидрохинокса-лин-2-оне.

3-[2-Оксопропилиден]-3,4-дигидрохинок-салин-2-он — первый представитель исследуемого ряда хиноксалинов (рис. 2). Его строение подтверждено спектральными методами 5.

Примечание. Через дробь показаны энтальпии сольватации с двумя и четырьмя молекулами воды. Абсолютное значение энтальпии, рассчитанное в приближении МР2/]отде-йгр, для воды Ио, Хартри = —76.229734, для диметилоксалата Ио, Хартри = -455.655287.

Рис. 2. Структуры соединений исследуемого ряда 3-[2'-алкил-2'-оксоэтилиден]-3,4-дигидрохинокса-лин-2(1Н)-онов

1

Для установления наиболее вероятной структуры соединения 1 были построены и оптимизированы его различные конформации и конфигурации (табл. 4). Исходя из величины Е, можно сказать, что наиболее вероятной конформацией будет 1а.

Таблица 4 Энергии образования различных конформаций 3-[2-оксопропилиден]-3,4-дигидрохиноксалин-2-она 1

Конфор-мация Хартри Д£°, кДж/моль Структура

1а -684.5263189 0 н

1Ь -684.5163734 26.11 н

1с -684.5061 893 52.85 н ахг н

1а -684.5054562 54.77 оёа н

Сравнение экспериментально полученных ИК- и ЯМР-спектров с теоретическими, говорит о хорошей геометрической воспроизводимости структуры молекулы 1а, указанной на рис. 2. Масштабирующий фактор для коррекции теоретических частот не применялся. Получение спектров ЯМР проводили в среде растворителя диметилсульфоксида, поэтому при

оптимизации соединения 1а для расчета ЯМР-спектра использовали значение е = 46.826.

По результатам экспериментальных исследований ИК спектр, V, см-1 (тв.) 5: 1678 V (С=0), 1614 V (С=О...НЮ, 1575, 1557, 1500 V (САг—САг), 931-642 8(САг—Н, САг—Н). Сравнение теоретического и экспериментального ИК-спектров представлено на рис. 3. Наблюдается небольшое смещение полос в расчетном спектре относительно экспериментального, но результаты показывают относительное совпадение характеристических частот.

Экспериментальный спектр ЯМР !Н, 8, м.д. (ДМС0-^6) 5: 2.17 с (3Н, СН3), 6.03 с (1Н, =СН-), 6.35 с (1Н, КН), 6.47 с (1Н, КН), 7.04-7.14 м (4Н, С6Н4 в составе хиноксалино-вого ядра). Скорректированный теоретический ЯМР *Н-спектр, 8, м.д. (ДМС0-^6): 1.93, 1.94, 1.23 с (3Н, СН3), 5.74 с (1Н, =СН-), 6.01 с (1Н, КН-0=С), 10.45 с (1Н, КН), 6.22, 6.49, 6.62, 6.71 с (4Н, С6Н4 в составе хинокса-линового ядра). Среднее значение константы экранирования для тетраметилсилана, рассчитанное в приближение РБЕ/]отде-йгр, составляет 31.152 (м.д.). Таким образом, ошибка относительно эксперимента составляет -0.048 (м.д.).

Растворимость 3-[2-оксопропилиден]-3,4-дигидрохиноксалин-2-она составляет 0.501 г на 100 г в воде, 3.965 г на 100 г в этиловом спирте, он также хорошо растворим в ДМСО и диоксане, трудно — в бензоле и толуоле. Механизм сольватации предположительно связан с взаимодействием молекулы хиноксалина и молекулы растворителя в положении аминогрупп (рис. 4). Процесс диссоциации исследуемых хиноксалинов представлен на схеме 3.

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

1.672 446

500 1 000 1 500 2 000 2 500 3 000 3 500

Рис. 3. Сравнение теоретического (снизу) и экспериментального (сверху) ИК-спектров 3-[2-оксопропили-ден]-3,4-дигидрохиноксалин-2-она

Рис. 4. Расположение молекул воды относительно иминогрупп 3- [ 2-оксопропилиден]-3,4-дигидрохи-ноксалин-2-она

+ 2Н,0+

N"

Схема 3. Диссоциация 3-[2'-алкил-2'-оксоэтили-ден]-3,4-дигидрохиноксалин-2(1Н)-онов

Для расчета значения энергии диссоциации соединения 1a был использован термодинамический цикл I (схема 4), из которого следует формула для расчета AGaq:

AGaq = [(Н + ) + ^ (A~ ) " Ggas (HA)] +

+[n^Gsol (H +) + AGsol (A-) - №sol (m)];

AGsol = Gaq - Ggas-

lnKsp =

([( (H +) + AGgas (A-) - AGgas (HA)) + RTln 24.46] +j

|+[nAGsol (H +) + AGsoi I

/RT

A -)- Gso

-n272.18

Схема 4. Термодинамический цикл I

Квантово-химическими расчетами невозможно получить энергию для протона. Для ЛGgas(H+) использовали значение —6.28 ккал/ моль, рассчитанное в рамках теории Сакура-Тетроде, а ЛGsol(H+) заменили экспериментально полученной величиной энергии сольватации протона —265.9 ккал/моль 15. Данный подход уже был апробирован в другой статье и давал неплохие результаты 16. Для приведения ЛGgas298 при 1 атм к 1М, которая используется при расчетах ЛGsol, применяли выражение: (1М) = (1атм) + КТ\п 24 . 46

В итоге уравнение для расчета произведения растворимости выглядит так:

^ (НА)]- п

При расчете Ksp для соединения 1а в этаноле был использован дополнительный термодинамический цикл, аналогичный первому (схема 4), по которому высчитывалось значе-

ние ЛGsol(H+) протона в спирте, так как экспериментальных данных сольватации найти не удалось. ЛGaq этанола в воде рассчитывали из экспериментально полученной константы диссоциации спирта (96%) при 20 оС. Полученное значение ЛGsol(H+) для протона в спирте также подставляли в схему 3 для хиноксалина в спирте.

При таком подходе > 1 для соединения 1а, из чего следует, что для корректного расчета энергии диссоциации необходим учет неидеальности растворов. В целом, согласно теоретическому расчету, можно считать, что соединение 1а должно быть хорошо растворимо в этаноле, что коррелирует с экспериментальными данными. Значения энергий Гиббса ионов и нейтральных молекул исследуемых соединений приведены в табл. 5.

Мы напрямую не учитывали энтальпию сольватации, но известно, что данный подход дает существенно завышенные показатели ЛHsol и, следовательно, рассчитанные с их использованием не совсем корректные значения энергий диссоциации 17' 18. Полученные данным способом теоретические значения ЛGaq и lnKsp можно использовать лишь как относительный показатель растворимости в ряду гомологов. Влияние растворителя на физико-химические характеристики в ряду производных хиноксалина не существенно 19. Опираясь на экспериментальный результат первого представителя и совмещая с теоретическими значениями, можно предсказать относительное экспериментальное значение растворимости, но точный результат можно узнать только в эксперименте. Такой метод будет полезен при очень большом числе гомологов или небольшом выходе продуктов реакции.

Таким образом, наблюдаемая хорошая корреляция рассчитанных ИК и ЯМР-спектров с практическими является доказательством правильной воспроизводимости структур молекул при использовании приближения PBE1/jorдe-dzp. Теоретические значения свободной энергии диссоциации получаются завышенными, но в ряду гомологов хорошо показывают изменение растворимости синтезируемых соединений.

Полученные результаты позволяют с уверенностью говорить о том, что использование базиса Jorge-DZP и метода PBE1 для расчета алкильных производных 3,4-дигидрохинокса-лин-2-она дает возможность провести достоверную оптимизацию в рамках последующего

Таблица 5

Абсолютные значения энергий Гиббса ионов и нейтральных молекул 3,4-дигидрохиноксалин-2-1(Н)-онов в воде и энергии их диссоциации

Соединение 1 2 3 4 5

Ggas(HA) -684.354379 -723.5921 1 -762.834073 -802.076339 -841.318825

Gsoi(HA) -684.370121 -723.610852 -762.852846 -802.094293 -841.335125

Gaas(AZ~) -683.126429 -722.370643 -761.61115 -800.853496 -840.097088

Gsoi(A2-) -683.423352 -722.66296 -761.90491 -801.145767 -840.385841

InKsp 1.82E-38 5.57 E-39 5.31 E-39 2.85 E-39 1.28 E-39

AGaq (кдж/моль) 216.06 219.01 219.13 220.68 222.67

В/А 1 3.273 1.048 1.865 2.226

X/А 1 3.273 3.429 6.394 14.235

б 0.501 0.153 0.146 0.0783 0.0352

Примечание. Экспериментальное рассчитанное значение произведения растворимости соединения 1а в воде без учета активности ионов равно 1.502Е-05; В/А — отношение 1пК$р следующего соединения к предыдущему; X/А — отношение 1пКр следующего соединения к первому представителю; д — предсказанная растворимость будущих соединений.

молекулярного докинга. Энергетические и геометрические параметры рассчитываются с малой погрешностью, что необходимо для корректного ОМ/МИ+ или ОМ1ОМ исследования.

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

Литература

1. Pereira J.A., Pessoa A.M., Cordeiro M.N.D.S., Fernandes R., Prudencio C., Noronha J.P., Vieira M. Quinoxaline, its derivatives and applications: A State of the Art review // European Journal of Medicinal Chemistry.- 2015.- T.97.- C.664-672.

2. Писаренко E.A. Фармакологическая характеристика 2,3-бис-(гидроксиметил)-хиноксалина 1,4-ди-^оксида // Международный журнал прикладных и фундаментальных исследований.- 2014.- №3-1.- C.100-101.

3. Gobouri A. Синтез и биологическая оценка некоторых N -замещенных производных хинокса-лина в качестве противоопухолевых агентов // Биоорганическая химия.- 2020.- T.46, №3.-C.327-327.

4. Носова Э.В., Трашахова Т.В., Устюгов В.С., Мочульская Н.Н., Валова М.С., Липунова Г.Н., Чарушин В.Н. Фторсодержащие стирилп-роизводные хинолина и хиноксалина: синтез и фотофизические свойства // Известия Академии наук. Сер. хим.- 2011.- №5.- C.920-925.

5. Казаев К.А., Левенец Т.В., Кунавина E.A. Тан-демный синтез 3-[2-оксопропилиден]-3, 4-ди-гидрохиноксалин-2(1Н)-она // Университетский комплекс как региональный центр образования, науки и культуры.- 2021.— C.2709-2711.

6. Peshkov S.A., Khursan S.L. Complexation of the Zn, Co, Cd, and Pb ions by metallothioneins: A QM/MM simulation // Computational and Theoretical Chemistry.- 2017.- T.1106.- C.1-6.

7. Хурсан С. Л. Сопоставительный анализ теоретических методов определения термохимических характеристик органических соединений // Вестник Башкирского университета.-2014.- T.19, №2.- C.395-402.

References

1. Pereira J.A., Pessoa A.M., Cordeiro M.N.D.S., Fernandes R., Prudencio C., Noronha J. P., Vieira M. [Quinoxaline, its derivatives and applications: A State of the Art review]. European Journal of Medicinal Chemistry, 2015, vol.97, pp.664-672.

2. Pisarenko E.A. Farmakologicheskaya kharak-teristika 2,3-bis-( gidroksimetil )-hinoksalina-1,4-di-N-oksida [Pharmacological characteristics of 2,3-bis-(hydroxymethyl)-quinoxaline 1,4-di-N-oxide]. Mezhdunarodnyi zhurnal prikladnykh i fundamental'nykh issledovaniy [International Journal of Applied and Fundamental Research], 2014, no.3-1, pp. 100-101.

3. Gobouri A. [Synthesis and Biological Evaluation of Some N-Substituted Quinoxaline Derivatives as Antitumor Agents]. Russian Journal of Bioorganic Chemistry, 2020, vol.46, pp.409-416.

4. Nosova E.V., Trashakhova T.V., Ustyugov V.S., Mochul'skaya N.N., Valova M.S., Lipunova G.N., Charushin V.N. [Fluorine-containing quinoline and quinoxaline styryl derivatives: Synthesis and photophysical properties]. Russian Chemical Bulletin, 2011, vol.60, no.5, pp.942-947.

5. Kazaev K.A., Levenets T.V., Kunavina E.A. Tan-demnyi sintez 3-[2-oksopropiliden]-3,4-digid-rokhinoksalin-2(1H)-ona [Tandem synthesis of 3-[2-oxopropylidene]-3,4-dihydroquinoxalin-2(1H)-one]. Universitetskiy kompleks kak regional'nyi tsentr obrazovaniya, nauki i kul'tury [University complex as a regional center of education, science and culture], 2021, pp.2709-2711.

6. Peshkov S.A., Khursan S.L. [Complexation of the Zn, Co, Cd, and Pb ions by metallothioneins: A QM/MM simulation]. Computational and Theoretical Chemistry, 2017, vol.1106, pp.1-6.

7. Khursan S.L. Sopostavitel'nyy analiz teoreti-cheskikh metodov opredeleniya termokhi-

8.

9.

10.

11.

12.

13.

14.

15.

16.

17.

18.

19.

64

Хурсан С.Л. Гомодесмический метод определения энергии диссоциации связей O—H в фенолах // Кинетика и катализ.— 2016.— T.57, №2.- C.164-175.

Jorge F.E., Neto A.C., Camiletti G.G., Machado S.F. Contracted Gaussian basis sets for Douglas-Kroll-Hess calculations: Estimating scalar relativistic effects of some atomic and molecular properties // The Journal of Chemical Physics.-2009.- T.130, №6.- C.064108.

Anthoney M.E., Carson A.S., Laye P.G. The thermochemistry of formates and oxalates // Conf. Int. Thermodyn. Chim. C.R. 4th.- 1975.-C.99-103.

Anthoney M.E., Carson A.S., Laye P.G., Yrekli M. The enthalpy of formation of dimethyl oxalate // J. Chem. Thermodyn.- 1976.- №8.- C.1009-1010. Anthoney M.E., Finch A., Stephens M. The solution of dimethyl oxalate and oxalic acid in water and alkali. Standard enthalpy of formation of dimethyl oxalate // Thermochim. Acta.-1975.- №12.- C.427-431.

Шеповалов К.М., Маслова О. А., Безносюк С.А., Жуковский М.С., Жуковская Т.М. Оптимальность и точность компьютерных вычислений свободной энергии Гиббса гидратации молекул в континуальных моделях сольватации // Известия Алтайского государственного университета.- 2019.- №1(105).- С.56-59.

Samatov A.A., Nagrimanov R.N., Miroshnichenko E.A., Solomonov B.N. Vaporization/sublimation enthalpies of mono-and dimethyl-esters estimated by solution calorimetry method // Thermochimica Acta.-2020.- T.685.- C.178529.

Camaioni D.M., Schwerdtfeger C.A. Comment on «Accurate experimental values for the free energies of hydration of H+, OH-, and H3O+» // The Journal of Physical Chemistry A.- 2005.-V.109, №47.- Pp.10795-10797.

Liptak M.D., Gross K.C., Seybold P.G., Feldgus S., Shields G.C. Absolute pKa determinations for substituted phenols // Journal of the American Chemical Society.- 2002.- V.124, №22.-Pp.6421-6427.

Meloun M., Ferencikova Z. Enthalpy-entropy compensation for some drugs dissociation in aqueous solutions // Fluid Phase Equilibria.-2012.- V.328.- C.31-41.

Dragan A. I., Read C. M., Crane-Robinson C. Enthalpy-entropy compensation: the role of solvation // European Biophysics Journal.-2017.- V.46, №4.- Pp.301-308. Щербаков B.B., Курбатова С.В., Земцова М.Н. Влияние природы растворителя на физико-химические характеристики производных хиноли-на // Бутлеровские сообщения.- 2019.- T.60, №11.- C.31-39.

micheskikh kharakteristik organicheskikh soyedineniy [Comparative analysis of theoretical methods for the determination of thermochemical parameters of organic compounds]. Vestnik Bashkirskogo universiteta [Bulletin of Bashkir University], 2014, vol.19, no.2, pp.395-402.

8. Khursan, S.L. [ Homodesmotic method of determining the O—H bond dissociation energies in phenols]. Kinetics and Catalysis, 2016, vol.57, no.2. pp.159-169.

9. Jorge F.E., Neto A.C., Camiletti G.G., Machado S.F. Contracted Gaussian basis sets for Douglas-Kroll-Hess calculations: Estimating scalar relativistic effects of some atomic and molecular properties // The Journal of Chemical Physics. 2009. Vol. 130. no.6. pp. 064108.

10. Anthoney M.E., Carson A.S., Laye P.G. [The thermochemistry of formates and oxalates]. Conf. Int. Thermodyn. Chim. C.R. 4th, 1975, pp.99-103.

11. Anthoney M.E., Carson A.S., Laye P.G., Yrekli M. [The enthalpy of formation of dimethyl oxalate]. J. Chem. Thermodyn., 1976, no.8. pp.1009-1010.

12. Anthoney M.E., Finch A., Stephens M. [The solution of dimethyl oxalate and oxalic acid in water and alkali. Standard enthalpy of formation of dimethyl oxalate]. Thermochim. Acta, 1975, no.12, pp.427-431.

13. Shepovalov K.M., Maslova O.A., Beznosyuk S.A., Zhukovsky M.S., Zhukovsky T.M. Opti-mal'nost' i tochnost' komp'yuternykh vychisle-niy svobodnoy energii Gibbsa gidratatsii molekul v kontinual'nykh modelyakh sol'vata-tsii [The Optimality and Accuracy of Computer Calculations of the Gibbs Free Energy of Hydration of Molecules in the Continuum Models of Solvation]. Izvestiya Altayskogo gosudar-stvennogo universiteta [Proceedings of the Altai State University], 2019, no.1(105), pp.56-59.

14. Samatov A.A., Nagrimanov R.N., Miroshnichenko E.A., Solomonov B.N. [Vaporization/ sublimation enthalpies of mono- and dimethyl-esters estimated by solution calorimetry method]. Thermochimica Acta, 2020, vol.685, p.178529.

15. Camaioni D.M., Schwerdtfeger C.A. [Comment on «Accurate experimental values for the free energies of hydration of H+, OH-, and H3O+»]. The Journal of Physical Chemistry A., 2005, vol.109, no.47, pp.10795-10797.

16. Liptak M.D., Gross K.C., Seybold P.G., Feldgus S., Shields G.C. [Absolute pKa determinations for substituted phenols]. Journal of the American Chemical Society, 2002, vol.124, no.22, pp.6421-6427.

17. Meloun M., Ferencikova Z. [Enthalpy-entropy compensation for some drugs dissociation in aqueous solutions]. Fluid Phase Equilibria, 2012, vol.328, pp.31-41.

18. Dragan A.I., Read C.M., Crane-Robinson C. [Enthalpy-entropy compensation: the role of solvation]. European Biophysics Journal, 2017, vol.46, no.4, pp.301-308.

19. Shcherbakov V.V., Kurbatova S.V., Zemtsova M.N. Vliyaniye prirody rastvoritelya na fiziko-khimicheskiye kharakteristiki proizvodnykh khinolina [The effect of the solvent nature on the quinoline derivatives physicochemical characteristics]. Butlerovskiye soobshcheniya [Butlerov Communications], 2019, vol.60, no.11, pp.31-39.

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