CC BY
21
В.А. Потёмкин, М.А. Гришина, А.В. Белик, ВН. Чарушин, ОН. Чупахин* Теоретическое исследование реакций синтеза
полициклических фторхиниолонов
Проведено квантовохимическое исследование реакций, приводящих к образованию полициклических фторхинолонов. Обнаружено, что направление реакций циклизации подчиняется орбитальному контролю и зависит как от плотности НСМО на кольцевом атоме углерода, претерпевающем внутримолекулярную нуклеофильную атаку гете-роатомом заместителя, так и от плотности ВЗМО атакующего ге-тероатома. Полученные закономерности использованы для прогноза направления реакций новых соединений.
В современном арсенале фармацевтических препаратов важное место принадлежит производным фторхинолона как высокоактивным антибактериальным и противотуберкулезным средствам. Создание новых производных фторхинолона направлено на поиск еще более эффективных лекарственных средств. Изучение закономерностей при протекании реакций образования производных хинолона позволяет выяснить механизмы превращений и определить наилучшие пути синтеза новых соединений. В данной работе проведено исследование реакций образования полициклических фторхинолонов.
Сопоставлена реакционная способность соединений, приведенных на рис.1 (заместители Я1 которых приведены в табл.1), к реакциям циклизации. Обнаружено, что для соединений с У=Е реакция протекает по направлению 1, а для соединений с У=И - по направлению 2 (рис. 1).
Направление 1
Направление 2
р
Конформер а У=И, Б
Конформер Ь У=И, Б
* Работа выполнена при финансовой поддержке РФФИ (грант № 01-03-96414).
73
Для пары соединений 4: Для пары соединений 4:
У=И, Б У=И, Б
Рис.1. Возможные направления реакции циклизации
Таблица 1
Заместители Щ для пар соединений 1-3
№ пары соединений Заместитель Rl
1 -ох; / ; си'
2 в ^—'
3 УО"
В рамках полуэмпирического квантовохимического метода РМ3 определена оптимальная геометрия, рассчитаны энергетические и зарядовые характеристики конформеров а и Ь (рис.1) исходных соединений и гипотетически возможных продуктов направления реакций 1 и 2. Анализ показал, что разница в полной энергии конформеров а и Ь и барьер конформационного перехода невелики, -следовательно, в реакционной смеси должны присутствовать оба конформера.
Анализ энергетических характеристик продуктов и исходных соединений показал, что термодинамически более выгодным процессом является циклизация с образованием хинолонового цикла для соединений как с У=И, так и с У=Е. Выигрыш в полной энергии при протекании процесса по направлению 1 составляет 11.5 - 19.0 ккал/моль, а по направлению 2-2.0 - 11.3
ккал/моль. Однако результаты расчетов показывают, что термодинамический фактор далеко не единственный при выборе направления реакции.
Для выяснения основных причин направления реакций в рамках неэмпирических квантовохимических методов (базис STO 3G) проведен детальный анализ структур а и Ь, который показал следующее. Обнаружено, что направления реакций 1 и 2 являются конкурирующими процессами. С одной стороны, протекание реакции по направлению 1 является более термодинамически выгодным процессом, но для взаимодействия с нуклеофильной частицей соединениям в конформации а нужно преодолеть достаточно большой энергетический барьер, поскольку они обладают большей величиной энергии НСМО по сравнению с конформерами Ь (табл.2). Протекание же реакции по направлению 2 определяется кинетическими факторами, а именно распределением (коэффициентом волновой функции) НСМО на атоме N(2).
Таблица 2
Энергетические характеристики конформеров a и b исходных молекул, рассчитанные ab initio базис функций STO-3G
№ пары соединений Y Конформер Энергия Евзмо (Эв) Энергия Енсмо (эВ) &Е—Енсмо-Евзмо (эВ)
1 F а -6.677 5.445 12.122
1 F b -5.867 5.354 11.221
1 H b -6.506 -3.724 2.782
1 H a -6.333 5.889 12.222
2 F a -5.793 5.593 11.386
2 F b -5.368 5.240 10.608
2 H b -6.722 -2.816 3.906
2 H a -5.776 5.766 11.542
3 F a -7.224 5.212 12.436
3 F b -7.166 5.168 12.334
3 H b -7.540 -4.386 3.154
3 H a -7.245 5.301 12.546
4 F a -7.013 5.469 12.482
4 F b -6.755 5.329 12.084
4 H b -6.696 5.560 12.256
4 H a -7.161 5.751 12.912
Если плотность орбитали на этом атоме становится достаточной для взаимодействия с избыточной электронной плотностью кислорода, то про-
75
цесс протекает по направлению 2. Из данных табл.3 видно, что величина распределения НСМО на атоме N(2) для конформеров Ь в соединениях с У=Н больше, чем для соединений с У=Е, что и определяет направление их циклизации по пути 2. Более того, соединения с У=Н в конформации Ь чрезвычайно реакционно способны для взаимодействия с нуклеофильной частицей (исключение составляет соединение 4 с У=Н), поскольку обладают низкой величиной энергии НСМО (ЕНСМО) и энергии электронного перехода (ЛЕ=ЕНСМО-Евзмо) (табл.2).
Отмечено также, что из пар соединений (У=Е и У=Н) молекулы с У=Е, реагирующие по направлению 1, характеризуются большим распределением НСМО на кольцевом атоме С(2), а соединения с У=Н обладают большим распределением НСМО на атоме Н при N(2) (рНСМО) (по сравнению с соединениями с У=Е). Можно предположить, что отсутствие вакансии НСМО на атоме С(2) не способствует протеканию нуклеофильной атаки по направлению 1. Поскольку, как было показано ранее, в реакционной смеси присутствуют конформеры а и Ь, то для иллюстрации использованы средние значения распределения НСМО на атомах (табл.3).
Таблица 3
Величины распределения НСМО на атомах для конформеров а и Ь
№ пары соединений У р И<2> Унсмо конформе-ра а рНСМОМ(22 конформе-ра Ь _ С2 рнсмо для конформеров а и Ь Рнсмо атома Н при N(2) для конфор-меров а и Ь
1 Р 0.008 0.102 0.251 0.008
1 И 0.139 0.383 0.044 0.060
2 Р 0.008 0.109 0.233 0.006
2 И 0.014 0.194 0.155 0.025
3 Р 0.171 0.178 0.022 0.042
3 И 0.171 0.452 0.010 0.052
4 Р 0.064 0.136 0.304 0.040
4 И 0.042 0.184 0.259 0.047
Часто направление 1 реакции циклизации, приведенной на рис.1, не останавливается на образовании фторхинолона, и продукты в зависимости от природы радикала R претерпевают дальнейшую внутримолекулярную циклизацию с образованием конденсированной трехкольцевой системы.
Поэтому нами изучены факторы, определяющие дальнейший процесс циклизации (рис.2).
I II
Рис.2. Внутримолекулярная циклизация производных фторхинолона
Сопоставлена реакционная способность соединения I в зависимости от природы радикала R и гетероатома X. Экспериментальные данные приведены в табл.4.
Таблица 4
Примеры конкурентных реакций внутримолекулярного замещения атома F и циклизации с участием карбонильной группы (жесткие условия - продолжительное нагревание в толуоле в присутствии К2С03, мягкие условия - нагревание в бензоле или толуоле без добавления каких-либо реагентов, связывающих HF)
№ соеди- Заместители Протекание циклизации
нения
1 X=NH, R=Ph Не протекает
2 X=O, R=CH2CN Не протекает
3 X=O, R=CHз Не протекает
4 X=O, R=3-NO2C6H4 Протекает в жестких условиях
5 X=O, R=1-этил-6-фтор-7- Протекает в жестких условиях
циклоалкилимино -1,4-дигидро -4-
оксохинолин-3 -ил
6 X=O, R=Ph Протекает в жестких условиях
7 X=O, R=пиридин-3 -ил Протекает в жестких условиях
8 X=O, R=пиридин-4-ил Протекает в жестких условиях
9 X=S, R=циклогексиламино Протекает в мягких условиях
10 X=S, R=морфолин-1-ил Протекает в мягких условиях
11 X=S, R=пиперидин-1-ил Протекает в мягких условиях
12 X=S, R=тиоморфолин-1-ил Протекает в мягких условиях
13 X=S, R=пирролидин-1-ил Протекает в мягких условиях
Обнаружено, что процесс циклизации, по всей видимости, протекает из таутомерной формы 1б (рис.3). Для легко циклизующихся соединений она является наиболее выгодной и имеет достаточно большой выигрыш в полной энергии 11-13 ккал/моль. Остальные соединения в большинстве случаев преимущественно образуют форму 1а. Исключение составляет молекула с Х=0, Я= 1 -этил-6-фтор-7-циклоалкил-1,4-дигидро-оксохинолин-3-ил. Интересно, что и нециклизующаяся молекула с X = NH, Я = РЬ также преимущественно образует форму 1б. Однако разность энергий таутомерных форм для этих двух молекул минимальна и составляет соответственно 3 и 5 ккал/моль. Это говорит о несущественном преимуществе формы 1б по сравнению с формой 1а для этих соединений, что приводит к понижению количества формы 1б в реакционной смеси в сравнении с быстрореагирующими соединениями. Следует отметить, что важной особенностью соединений с высокой реакционной способностью является высокое напряжение валентных углов (ЕБто,табл.5), это снижает их термодинамическую стабильность и, соответственно, приводит к повышению их реакционной способности.
Таблица 5
Энергия напряжения валентных углов Евшо, выигрыш в энергии при протекании реакции циклизации Е1, орбитальные характеристики форм
1б
№ соединения ЕБЕ№>, ккал/моль Ей ккал/моль С р НСМО рВЗМО
1 12 23 0.1698 0.1099
2 11 3 0.3349 0.0677
3 11 3 0.3493 0.0813
4 49 3 0.0089 0.0758
6 51 3 0.0333 0.0770
7 61 3 0.0181 0.0779
8 67 3 0.0177 0.0767
9 15 16 0.3358 0.1489
10 19 16 0.3246 0.1418
11 17 17 0.3515 0.1179
12 19 15 0.3225 0.1445
13 21 17 0.3605 0.1209
Рассмотрено предположение, что образование циклического продукта может протекать через образование анионной формы (1в) ( рис.3).
O
O
n
F
F
+ Н+
COOEt
R
Рис.3. Образование анионной формы соединения I.
Однако энергетические, орбитальные и зарядовые характеристики анионов как циклизующихся, так и нециклизующихся соединений имеют близкие значения.
Из последнего может быть сделано заключение о том, что протекание реакции определяется образованием таутомерной формы 1б, диссоциация которой не лимитирует протекание процесса. Обнаружено, что выигрыш в энергии при образовании циклического продукта составляет:
1) для легко циклизующихся соединений (с Х=8) - 15+17 ккал/моль;
2) для нециклизующихся и трудно циклизующихся соединений - 1+3
Исключение составляет молекула 1 (X=N, R=Ph), которая не вступает в процесс циклизации в данных условиях, хотя расчеты показывают, что образование циклического продукта возможно. Близки по значениям и выигрыши в энергии циклизующихся и нециклизующихся соединений 2 - 8.
Все говорит о том, что термодинамическая выгодность циклических продуктов реакции является не единственным условием ее протекания.
Проведенный более детальный анализ характеристик таутомерной формы 1б (ab initio, базис STO-3G) показал, что наиболее информативным фактором, отличающим циклизующиеся и нециклизующиеся соединения, является плотность распределения ВЗМО (рВЗМО) на гетероатоме X и НСМО на атоме С8 хинолонового цикла (рНСМО). Последнее вполне согласуется с выводами, полученными при рассмотрении реакций, представленных на рис.1. Характеристики рВЗМО и рНСМО приведены в табл.5. Наибольшей заселенностью ВЗМО обладают легко циклизующиеся соединения, что может характеризовать их как соединения с более ярко выраженными нуклеофиль-ными свойствами верхней орбитали гетероатома Х. При этом распределение НСМО на атоме С8 также должно иметь достаточно высокое значение, то есть избыточная электронная плотность верхней орбитали атома Х должна взаимодействовать с вакансией НСМО атома С8. В случае недостаточной "нуклеофильности" атома Х избыточная вакансия НСМО атома С8 играет
ккал/моль (см. табл.5).
негативную роль, поскольку С8 в этом случае ассоциирует с каким-либо другим атомом данной или соседней молекулы или с растворителем. Введение щелочи в такую систему является малоэффективным, поскольку ассоциация атома С8 в этом случае происходит со значительно более электроноизбыточ-ным анионом СО32-, который стерически блокирует центр циклизации. Этот эффект, вероятно, наблюдается в молекулах 2, 3 (отчасти, может быть, 1). А вот в случае малого распределения НСМО на атоме С8 циклизацию можно "вынудить" добавлением К2СО3, поскольку последний не ассоциирует с С8, а усиливает диссоциацию связи Х-Н, что увеличивает распределение ВЗМО на атоме Х (усиливая его "нуклеофильность"). Кроме того, введение карбоната калия приводит к нейтрализации HF.
Особенно ярко орбитальный запрет выражен у соединения 1 , которое преимущественно находится в таутомерной форме 1б, продукт циклизации -термодинамически выгоден, однако заселенность ВЗМО атома Х ниже, чем у соединений 9 - 13 при относительно большой вакансии НСМО на атоме С8.
Таким образом, реакция циклизации протекает из таутомерной формы 1б. Эффективное протекание реакции обусловливается термодинамической стабильностью данной формы, а также продуктов циклизации. Дополнительный запрет на протекание циклизации определяется орбитальным контролем (высокая заселенность ВЗМО атома Х и ее соотношение с НСМО атома С8).
Полученные результаты реакции циклизации производных фторхино-лона (рис. 2) использованы для предсказания направления процессов, приведенных на рис. 4. Расчет характеристик молекул показал следующее:
1) Величина энергии напряжения валентных углов соединений III (Ebend=28.77 ккал/моль) и V (Ebend=22.17) выше, чем величина Ebend 1б-форм легко циклизующихся соединений (табл.5).
2) Величины выигрыша в энергии при протекании реакций, приведенных на рис. 4, составляют соответственно E1=ETOTAL(ucx.)+ETOTAL(HF)-ETOTAL(np.)= - 3.5 и -1 ккал/моль. Поэтому можно сделать вывод, что эти реакции термодинамически выгодны.
3) Обнаружено, что соединение III обладает плотностью ВЗМО на ге-тероатоме X, равной рВЗМО=0.1003 (здесь X=N), и плотностью НСМО на атоме С(8) хинолонового кольца, равной рС8НСМО=0.3369. Для соединения V рвЗМО=0.03955 (здесь X=O) и рС8НСМО=0.25629. Сравнивая эти данные с данными табл.5, видно, что протекание реакций, приведенных на рис. 4, попадает под орбитальный запрет, так как величины плотности рВЗМО этих соединений недостаточны для протекания процесса. Плотность рС8НСМО достаточно велика для обоих соединений и может играть здесь негативную роль, поскольку С8 в этом случае может ассоциировать с каким-либо другим атомом данной или соседней молекулы или с растворителем. Введение щелочи в такую систему является малоэффективным, поскольку ассоциация атома С8 в этом случае может произойти со значительно более электроноизбыточным анионом СО32-, который стерически блокирует центр циклизации.
V VI
Рис.4. Предполагаемые реакции с образованием новых фторхинолонов
Исходя из полученных результатов, можно сделать вывод, что циклизация соединений III и V с образованием соединений IV и VI соответственно протекать не будет.
Таким образом, проведен детальный анализ реакционной способности большой группы фторсодержащих ароматических соединений в реакциях циклизации. Обнаружено, что циклизация с образованием хинолонового цикла является термодинамически выгодным процессом, однако влияние кинетических факторов может привести к протеканию конкурирующего процесса, а именно к циклизации по атому С карбонильной группы заместителя. Выяснено, что основным кинетическим фактором при протекании конкурирующего процесса является высокое распределение НСМО на атоме азота, участвующем в циклизации. Недостаточное распределение НСМО на атоме углерода кольца также мешает протеканию процесса циклизации с образованием хинолонового цикла.
Показано, что протекание процессов циклизации производных фторхи-нолона также определяется величиной распределения НСМО на атоме углерода хинолонового цикла. Помимо этого одним из основных факторов, определяющих возможность превращения, является высокое распределение ВЗМО на атакующем атом С(8) цикла гетероатоме заместителя.
Полученные данные использованы для предсказания возможности протекания ряда реакций.
