Время жизни синглетного кислорода в столкновительных комплексах o2-co2 Текст научной статьи по специальности «Нанотехнологии»

УДК 577.33; 577.345; 544.147

Щепин А.С., Пешков С.А., Пеш кова Т.В.

Оренбургский государственный университет E-mail: thisiscaster@mail.com

ВРЕМЯ ЖИЗНИ СИНГЛЕТНОГО КИСЛОРОДА В СТОЛКНОВИТЕЛЬНЫХ КОМПЛЕКСАХ O2-CO2

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

Методами самосогласованного поля (SCF/ROHF) и функционала плотности (DFT/B3LYP) проведен конформационный анализ для определения геометрических характеристик комплексов кислорода с углекислым газом и выявления среди них наиболее устойчивой структуры, обладающей минимумом энергии. Мультиконфигурационным методом CASSCF (14, 11) рассчитаны электронные термы комплекса 1,3[O2-CO2]. Методом спинорбитального взаимодействия (СОВ) рассчитаны моменты запрещенных переходов a1 - X3 b1 + - a1 , коэффициенты Эйнштейна (вероятности переходов), а также время жизни синглетного кислорода (a1 g, b1 g+) при комплексообразовании молекулы кислорода с молекулой углекислого газа.

Определена наиболее устойчивая структура комплекса [O2-CO2] - II, обладающая точеной группой симметрии C2. Выявлено, что основное и возбужденные состояния кислорода в комплексе 1,3[O2-CO2] O2(X3 ) + Co^A1 g+), O2(a1 g ; a'1 g) + CO2(A1 g+), O2(b1 g+) + CO2(A1 g+) являются малоустойчивыми, так как обладают низкими энергиями диссоциации. Показана возможность проявления запрещенных индуцированных переходов a1 g - X3g-, b1 g+ - a1 g в молекуле O2 под действием молекулы столкновения углекислого газа, которые практически не проявляются в отдельной молекуле О2.

Рассчитанные значения времени жизни синглетного кислорода в столкновительных комплексах 13[O2-CO2], характеризуют возможность образования комплексов кислорода с углекислым газом, и другими низкомолекулярными соединениями, обладающими мембранным транспортом, в тканях живого организма при облучении низкоинтенсивным лазером, за счет чего, увеличивается время жизни синглетного кислорода O2(a1 g, b1 g+).

Ключевые слова: синглетный кислород, возбужденные состояния, столкновительные комплексы, углекислый газ.

Синглетный кислород - общий термин, применяемый для двух метастабильных состояний (а1 А Ь1^4") молекулярного кислорода (02) с более высокой энергией, чем в триплетном состоянии (Х3^-), являющимся основным. Первоначально синглетный кислород в состоянии 02(Ь1Её+) представлял особый интерес к изучению только в физике верхних слоёв атмосферы, за счет проявляемого ночного свечения, образуемого при переходе Ь1^4 - Х3^-, именуемого позже люминесценцией кислорода [1]. В многочисленных экспериментальных исследованиях обнаружилось, что синглетный кислород имеет очень высокую окислительную активность при взаимодействии с органическими веществами, участвует в биологических процессах живых организмов [2], [3]. Отмечена наибольшая значимость использования синглетного кислорода в медицине в качестве агента, применяемого при лечении многих заболеваний, в том числе и онкологических [4].

Современная онкология, насчитывает большое количество методов, в том числе и комби-

нированных по устранению злокачественных опухолей в живом организме человека, и одним из них перспективно развивающимся в наше время является фотодинамическая терапия. Интерес к ней обусловлен тем, что разрушение опухоли происходит при облучении ее лазерным излучением (низкой интенсивности), исключающим опасность неконтролируемого термического повреждения стенки органа [5]. Механизм проведения процедуры подробно описан в работах [6], [7]. Известно, что при ведении фотодинамической терапии, в живой организм вводятся молекулы фотосенсибилизатора преимущественно органической природы, заданной концентрации, которые избирательно фиксируются на мембранах опухолевых клеток и митохондриях [7]. Затем на эту часть опухолевой ткани, воздействуют лазерным излучением, в результате чего происходит возбуждение фотосенсибилизатора, который передает энергию собственного возбуждения на молекулы кислорода О2. В результате происходит конверсия спина на молекулярной орбитали п молекулы

О2 (изменение спинового состояния) с последующим переходом кислорода в возбужденные синглетные состояния (а:Д , ЫЕ +).

У Б' % '

Учитывая то, что синглетный кислород в растворенном виде (в воде) имеет малое время жизни т(а1Д%) = 10-11 с, т(Ь11%+) = 10-5 с [8], то остается необъяснимым, почему в живом организме, его время жизни увеличивается, так как он успевает уничтожить большое количество раковых клеток.

Механизм увеличения времени жизни может заключаться в том, что получив энергию от сенсибилизатора, синглетный кислород может образовывать комплексы слабой прочности с молекулами газов, свободно проходящими через билипидный слой мембран клеток в тканях. В основном это низкомолекулярные соединения: неполярные (азот, бензол) и полярные (углекислый газ, вода, мочевина) молекулы.

Немногочисленные экспериментальные и теоретические исследования оксокомплексов с газами и 3 ^переходными металлами в различных спиновых состояниях указывают на значительное увеличение жизни таких комплексов (на несколько порядков). Объясняется это тем, что в оксокомплексах под воздействием излучения (ИК - инфракрасного, УФ - ультрафиолетового) происходит генерация долгоживущего синглетного кислорода стабилизированного молекулой столкновения посредством сложных индуцированных запрещенных переходов а^ -Х3Е ЫЕ + - а1 Д , ЫЕ + - Х3Е - и квадрупольных

Б ' Б Б' Б Б

переходов из высоковозбужденных состояний (системы Герцберга, системы Шумана-Рунге) [1], [9]-[14].

Нельзя однозначно утверждать, что ок-сокомплекс сенсибилизатора с молекулой кислорода обладает высокой прочностью. В современном теоретическом исследовании взаимодействия крупной органической молекулы 5-метилрезорцина с молекулами кислорода в возбужденных состояниях говорится о том, что прочность межмолекулярных связей полученных комплексов очень мала и характеризуется значениями энергий связывания 1,06 - 2,03 кДж/моль [10].

Так как, провести теоретические расчеты всех соединений обладающих мембранным транспортом в живых системах с молекулой кислорода, требует значительных затрат време-

ни при расчете высококонфигурационными методами, то для объяснения механизма в качестве модельной системы нами был выбран комплекс кислорода с углекислым газом [O2-CO2].

Методами математического моделирования при использовании современных высокоуровневых методов на качественной модели комплекса [O2-CO2] можно рассчитать геометрические конфигурации, прочность межмолекулярных связей и оценить полученные данные с подобными теоретическими и экспериментальными работами.

Цель работы: выявить наиболее устойчивые конформации комплексов [O2-CO2], определить их геометрические характеристики; оценить значения энергий связей основного и возбужденных состояний молекулы кислорода с молекулой столкновения CO2, объяснить механизм увеличения времени жизни синглетного кислорода при комплексообразовании с молекулой столкновения CO2.

Методы расчёта

Методами самосогласованного поля SCF/ ROHF и функционала плотности DFT/B3LYP были рассчитаны равновесные структуры стол-кновительных комплексов кислорода с углекислым газом в различных конформациях. Для определения основного и возбужденных состояний (синглетные и триплетные состояния молекулы кислорода) в комплексах 1,3[O2-CO2] произведен расчет мультиконфигурационным методом, включающим в себя полное пространство активных молекулярных орбиталей CASSCF (14, 11) - распределение 14 электронов на 11 молекулярных орбиталях. Структура ор-биталей, включенных в активное пространство при расчете комплексов приведена в таблице 1. Все расчеты произведены при использовании базисного набора функций 6-311G++(2p, 2d), содержащего диффузные и поляризационные функции для более полного описания Ван-дер-Вальсовых взаимодействий комплексов молекул газов различного состава [11]. Методом спин-орбитального взаимодействия (SOC) на основе функций CASSCF рассчитаны моменты запрещенных излучательных переходов a:Ag -

Х3Е ЫЕ + - a1 A . На основе значений момен-

g ' g g

тов переходов были рассчитаны коэффициенты Эйнштейна A (вероятности перехода), время

жизни синглетного кислорода x(a1Ag, b1Sg+). Расчеты проведены с помощью программного обеспечения для математического моделирования Firefly 8.1 [15].

Обсуждение результатов

Методами SCF/ROHF и DFT/B3LYP был произведен расчет устойчивости комплексов молекулы кислорода с молекулой углекислого газа для различных конформаций. Значения

энергий, и геометрические структура комплексов 3[O2-CO2] приведены в таблице 2. В ходе оптимизации, все конформации I-VI сохранили свое геометрическое положение в одной плоскости XY. Расстояния между молекулами O2 и CO2, рассчитанные методами DFT и ROHF, отличаются незначительно (0,1-0,3 Á). Среди всех конформаций 3[O2-CO2], следует отметить, что II - обладает минимумом энергии, поэтому именно эта структура взята за основу, при

Таблица 1 - Структура молекулярных орбиталей молекулы кислорода 02 и углекислого газа С02 в комплексе 13[02-С02], метод расчета CASSCF(14, 11)

Молекулярная орбиталь

_(O)_

Структура

Молекулярная орбиталь

(CO,)

Структура

2

u

(O,)

t

1n z

(CO,)

>

3

2

(O,)

1n x (CO,)

2n z

u

(O,)

К

(CO,)

2n y

u

(O2)

J

z

(CO,)

$

2n y

g

(O2)

g

(O,)

3

2

(O,)

%

&

Примечание: в ходе образования комплекса

3[02-С02], произошло искажение молекулярных орбиталей углекислого газа.

моделировании комплекса столкновения син-глетных состояний кислорода О2(а1ДБ, Ь1еб+) с углекислым газом.

На рисунке 1 приведены энергетические термы образования комплекса 1,3[02-С02]. Для основного триплетного и возбужденных син-глетных состояний кислорода, термы - слабосвязанные, имеют минимумы энергий при равновесном расстоянии между 02 и С02 R = 3,67 А. Устойчивость таких состояний: О2(Х3Е ) + СО^ЕД О^а^ ; а'1Д'ё) + СО^ЕД О2(Ь1Её+) + СО2(А1Её+) составляет 0,43, 0,48 и 0,49 кДж/моль соответственно. Значения энергий устойчивости качественно согласуются с расчетами столкновительных комплексов

1,3

[О2-Н2] [12], ^-ед] [13]. Энергии расщепления термов комплекса между состояниями О^ЕДО/а^) и О2(Х3Её-)-О2(Ь1Её+) равны 0,98 и 1,59 эВ соответственно. К примеру, экспериментально определенные значения энергий расщепления в отдельной молекуле О2 между соответствующими состояниями - 0,98 и 1,63 эВ [14].

При образовании комплекса 1,3[О2-СО2] в молекуле кислорода инициируются запрещенные электродипольные переходы а1ДБ - X3Еg", Ь1Еб+ - а1ДБ. По значениям моментов переходов (таблица 3), были рассчитаны: коэффициенты Эйнштейна, время жизни синглетного кислорода. Выяснилось, что значения х(а1Дё), т(Ь1Е%+)

Таблица 2 - Значения энергий и геометрические структура комплексов 3[O2-CO2] в различных конформациях I-VI. Метод расчета SCF/ROHF и DFT/B3LYP. В круглых скобках приведены значения длин связей по методу DFT/B3LYP

Комплекс

Геометрическая структура

Энергия комплекса, а.е.м.

Методы расчета

SCF/ROHF DFT/B3LYP

Конформация I

-337,348579

-339,024583

Конформация II

-337,348741

-339,024606

Конформация III

-337,348701

-339,024634

Конформация IV

° ° t!j Иг-" ° с °

-337,348637

-339,024601

Конформация V

-337,348562

-339,024588

Конформация VI

-337,348625

-339,024536

Таблица 3 -Значения моментов переходов М (еа0), коэффициентов Эйнштейна А (с-1), времени жизни т (с) состояний кислорода в комплексе 1,3[02-С02]. Метод расчета СА88СЕ(14, 11).

Параметр Численное значение

M(a1Ag ^ X3!g-) 1-10-5

M(b1!g+ ^ a1Ag) 7Д-10"4

A(a1Ag ^ X3!g-) 1-10-4

A(b1!g+ ^ a1Ag) 2,15

T(a1Ag) 9980

T(b1!g+) 0,47

гораздо выше в комплексе 1,3[02-С02] нежели в отдельной молекуле синглетного кислорода растворенного в воде [8] на несколько порядков.

Заключение

Электронные термы взаимодействия кислорода в основном триплетном (Х3^-) и син-глетном (а1 А Ь1^4) состояниях с молекулой столкновения С02 являются слабосвязанными. Время жизни синглетного кислорода в столкно-вительных комплексах 1,3[02-С02] на несколько порядков выше, чем в отдельных молекулах 02 в синглетном возбужденном состоянии. Все это характеризует возможность образования комплексов кислорода с углекислым газом, и другими низкомолекулярными соединениями, обладающими мембранным транспортом, в тканях живого организма при их облучении низкоинтенсивным лазером, за счет чего, уве-

Рисунок 1 - Энергетические термы образования комплекса 1,3[02-С02], в конформации 2, метод расчета CASSCF(14, 11)

личивается время жизни синглетного кислорода

O2(a1A , Ь1! +).

2Л р" g '

16.02.2016

Список литературы:

1. Кобзев Г.И. Зависимость люминесценции молекулярного кислорода от сорта и числа атомов, входящих в состав комплекса, и числа молекул окружения кислорода // Вестник Оренбургского государственного университета. - 2005. - №1. - С. 150 - 156.

2. Briviba K., Klorz I.-O., Sics H. Toxic and signaling effects of photochemically or chemically generated singlet oxygen in biological systems.// Biol. Chem., 1997, № 378, P. 1259-1265.

3. Landry M.P. Characterization of photoactivated singlet oxygen damage in single-molecule optical trap experiments.// Biophys. J., 2009, № 97, V. 8, P. 2128-2136.

4. Trabanco A.A. A seco-porphyrazine: Superb sensitizer for singlet oxygen generation and endoperoxide synthesis.// Synlett., 2000, № 7, P. 1010-1012.

5. Гельфонд М.Л. Фотодинамическая терапия в онкологии // Практическая онкология, Т. 8, № 4. - 2007, C. 204-210.

6. Henderson B.W., Dougherty T.J. How does photodynamic therapy work? // Photochem Photobiol., 1992, V. 55, P.145-157.

7. Jori G., Reddi E. The role of lipoproteins in the delivery of tumour-targeting photo-sensitizers // Int. J. Biochem., 1993, V. 25, P. 1369-1375.

8. Шинкаренко Н. В., Алесковский В. Б. Синглетныи кислород, методы получения и обнаружения// Успехи химии, 1981, Том L, Вып. 3, С. 406-428.

9. Minaev B.F., Kobzev G.I. Response calculations of electronic and vibrational transitions in molecular oxygen induced by interaction with noble gases // Spectrochim. Acta Part A, 2003, V. 59, P. 3387-3410.

10. Кобзев Г.И., Казаева А.Н., Футтеркнехт З. Влияние 5-метилрезорцина на спиновую конверсию в молекулярном кислороде // Журнал структурной химии 2014, Т. 55, №1, С. 14-19.

11. Урваев Д.Г., Кобзев Г.И. Интермедиаты газофазной реакции Sc + O2: Квантовохимическое исследование // Вестник ОГУ №10 (159) - 2013 - С. 359-361.

12. Minaev B.F., Lunell S., Kobzev G.I. Collision-Induced Intensity of the ЫЕ + -а'Д Transition in Molecular Oxygen: Model Calculations for the Collision Complex O2 + H2 // Intern. J. of Quant. Chem., 1994, V.g50, P. 279-292.

13. Minaev B. F., Kukueva V. V. Configuration Interaction Study of the O2-C2H4 Exciplex : Collision-induced Probabilities of Spin-forbidden Radiative and Non-radiative Transitions // J. Chem. Soc. Faraday Trans., 1994, № 11, V. 90, P. 1479-1486.

14. Schweitzer C., Schmidt R. Physical mechanisms of generation and deactivation of singlet oxygen. // J. Chem. Rev., 2003, V. 103, P. 1685-1757.

15. Schmidt M.W., Baldridge K.K., Boatz J.A., Elbert S.T., Gordon M.S., Jensen J.H., Koseki S., Matsunaga N., Nguyen K.A., Su S.J., Windus T.L., Dupuis M., Montgomery J.A. General Atomic and Molecular Electronic Structure System / J. Comput. Chem., 1993, V. 14, № 11, P. 1347-1363.

Сведения об авторах:

Щепин Антон Сергеевич, студент Оренбургского государственного университета, 460018, Оренбург, пр-т Победы, 13, ауд. 3431, e-mail: thisiscaster@gmail.com

Пешков Сергей Алексеевич, ведущий инженер кафедры химии Оренбургского государственного университета, 460018, Оренбург, пр-т Победы, 13, ауд. 3431, e-mail: darvin156@mail.ru

Пешкова Татьяна Владимировна, ведущий инженер кафедры химии Оренбургского государственного университета, 460018, Оренбург, пр-т Победы, 13, ауд. 3431, e-mail: radislavy@mail.ru