УДК 004.9/.942
СТРУКТУРНО-АДСОРБЦИОННАЯ МОДЕЛЬ ВОЗДЕЙСТВИЯ НИЗКОМОЛЕКУЛЯРНЫХ ТОКСИЧНЫХ СОЕДИНЕНИЙ НА БИОПОЛИМЕРЫ КЛЕТОЧНОЙ МЕМБРАНЫ
Н.В. Золотарева, Н.М. Алыков
Предложена модель процесса воздействия низкомолекулярных токсичных соединений на структурные элементы биологической мембраны на основе структурно-адсорбционного представления механизма образования связей. Квантово-механический аппарат позволяет решить задачу поиска минимума энергии, а корреляция с экспериментальными данными по геометрии позволяет оценить предложенную модель.
Ключевые слова: адсорбционная модель, квантово-механический аппарат, биополимеры, геометрическая оптимизация, активные центры.
Key words: model of adsorption, the quantum-mechanical apparatus, biopolymers, geometrical optimization, the active centers.
В последние десятилетия экологические проблемы отличаются остротой наблюдающихся процессов и явлений, вызванных мощным антропогенным давлением на природные системы. В настоящее время к основным источникам токсических выбросов (H2S, SO2, СО и др.) относятся продувка скважин, факельное хозяйство (на примере Астраханского газоперерабатывающего завода - АГПЗ), различного рода утечки газа, площадки розлива жидкой серы, установки Клауса типа «Сульфрен».
Токсическому воздействию диоксида серы подвергаются не только рабочие, но и жители ближайших населенных мест. Общее действие диоксида серы заключается в нарушении углеводного и белкового обмена, угнетении окислительных процессов в головном мозге, печени, мышцах, снижении содержания витаминов Bi и С. Диоксид серы также влияет на кроветворные органы, способствует образованию метгемоглобина, вызывает изменения в эндокринных органах и костной ткани.
Современные исследования и методы борьбы ограничиваются лишь оказанием медицинской помощи в виде лечения конкретного заболевания. Однако если подойти к этому вопросу со стороны прогнозирования очагов воспаления на молекулярном уровне, то можно смоделировать возможные процессы, происходящие на поверхности клеточной мембраны при воздействии на них токсикантов, это позволит установить активные центры в сложных биополимерных структурах. Знание о потенциальных мишенях воздействия в биологических мембранах позволит с помощью современных компьютерных программ изучить свойства веществ (антидотов-блокаторов), при взаимодействии с которыми наш организм не будет «реагировать» на природу токсикантов.
Данная проблема является актуальной, интересной и в то же время сложной с позиций изучения строения и превращения макромолекул, поскольку в процессе образования комплексов «мембрана - токсикант» необходимо учитывать следующие факторы: дальнодейст-
вующие электростатические взаимодействия между компонентом мембраны (например, белком) и токсикантом, геометрическое и энергетическое соответствие областей связывания, водородные связи, ван-дер-ваальсовы взаимодействия и прочее.
К сожалению, в настоящее время не существует универсального метода моделирования процессов образования комплексов «мембрана - токсикант», учитывающего вышеперечисленные факторы и точно прогнозирующего структуру образовавшегося комплекса. Однако возможным является исследование отдельных стадий. Таким образом, определены задачи исследования: поиск способа решения; применение методов для достижения поставленной цели; обработка полученных результатов и визуализация результатов расчетов в виде схем взаимодействий.
В качестве теоретической модели предлагаем рассмотреть ситуацию, когда между компонентом мембраны и токсикантом образуются водородные связи, а расчет характеристических показателей осуществляется с использованием квантово-механического аппарата. Моделирование процесса воздействия одного компонента на другой неразрывно связано с изучением процессов сорбции, в этом случае поверхность мембраны является «сорбционным материалом».
Поскольку мембраны представляют собой сложные системы, состоящие из отдельных молекул белков, углеводов и липидов, то расчеты можно свести к исследованию процессов воздействия низкомолекулярных соединений на отдельные фрагменты мембраны. Дефрагментирование кластера мембраны на отдельные составляющие позволяет упростить расчеты. Модель белка представима в виде последовательности аминокислотных остатков, например, пентапептида С\,N^N50782, модель углевода - в виде олигосахарида С12Н22О11, липида - в виде диглицерида С15Н24О5.
Если использовать теорию графов, то можно построить структуру молекул следующим образом: вершины графа соответствуют атомам (О, 8, Н и др.) молекул, а ребра - связям (0-Н, 8-Н и др.) [3, 5]. В результате, молекулярный граф задается в виде матрицы связей (1). Если молекула токсиканта образована по типу А-В-С, как большинство низкомолекулярных токсичных соединений (80г, Н28, N02), ей соответствует матрица вида 1.
А в с
А У 1 0
В 1 У 1
С 0 1 У
(1)
где у = (а - I'.) />: а - кулоновский интеграл; /? - резонансный интеграл; Е - полная энергия молекулярной орбитали (МО).
Каждая МО представлена в виде линейной комбинации атомных орбиталей (рп\
¥, = са(р1 + сп(р2 + ... + с1к<рк + ... + ст(рп (2)
где Т, - /-я молекулярная орбиталь, сл - орбитальный коэффициент.
Задача сводится к отысканию наилучшего набора значений орбитальных коэффициентов, который соответствует минимальному значению энергии МО и отвечает наиболее стабильному распределению электронной плотности. Этого можно добиться, используя вариационный принцип:
[ %Н^;с!т \\V\dT
где с!т - элементарный объем.
СИСТЕМНЫЙ АНАЛИЗ
И МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ
Процедура минимизации функции по орбитальным коэффициентам позволяет найти оптимальные конфигурации исходных молекул:
дЕ Л
— = 0 (4)
При моделировании процесса воздействия необходимо учитывать различные способы расположения токсиканта по отношению к конкретному типу атома в структуре биополимера. По количеству реакционно-способных атомов в молекуле биополимера можно судить о количестве образующихся адсорбционных систем.
При построении адсорбционной модели воздействия будем опираться на статическую модель Коулсона и Лонге-Хиггинса [4]. За основными изменениями в структуре биополимера можно следить по ближайшему окружению при образовании межмолекулярных связей. Поэтому можно ограничиться анализом трех-пяти атомов, например, если к атому Н в реагирующей молекуле (биополимер) приближается молекула реагента А (диоксид серы) по рис. 1, необходимо изучить изменения, происходящие с ближайшими атомами (5, р, I) в структуре биополимера при приближении и удалении реагента А.
. 5
Н----А
7
Рис. 1. Модель образования водородной связи
Количество взаимодействий зависит от типа и индивидуальных особенностей атомов. В молекуле 802 присутствуют два типа атома - электроположительный атом серы и электроотрицательный атом кислорода. Обозначим Х[ - 8 и х2 - О и составим матрицу образующихся адсорбционных систем (4) при взаимодействии токсиканта с биополимером: 802 -пентапептид, / = 40, где / соответствует количеству атомов водорода (электроположительный) и кислорода (электроотрицательный) в молекуле биополимера; 802 - олигосахарид, у = 33; 802 - диглицерид,у = 29. Характеристическими показателями протекания процесса воздействия являются энергетические параметры Е^.Жп образующихся систем.
802 - пентапептид Ег е2 Е3 е4 е5 Е40
XI 1 1 0 0 0 1
Х2 0 0 1 1 1 0
802 - олигосахарид Ех е2 Е3 е4 е5 Е33
XI 1 1 1 1 1 0
Х2 0 0 0 0 0 1
802 - диглицерид Ех е2 Е3 е4 е5 Е29
XI 1 1 1 0 0 0
Х2 0 0 0 1 1 1
(5)
В таблице приведены некоторые модели адсорбционных систем [2].
Таблица
Модели образования адсорбционных систем при воздействии диоксида серы
на ст руктуру пентапептида
802 - пентапептид Ег Е2 Еъ е4
XI я Vcн~c\ <о О-Н - -
Х2 — - О. О—Н 0-5.0 V \ 1 о <^сн „ /З-н-о^о н2с—сн2 ч
Согласно рис. 1, длина связи между атомом Н биополимера и ближайшим атомом молекулы А (диоксида серы) варьирует от 1,75 до 1,90 А, энергия адсорбционной системы (Ек), в нашем случае может быть представлена в виде:
Ег = (с,2 (Я) •«(£) + с2 (Я) • а (Я)) + 2((с1 (Я) • с2 (О) • ДО.. .О) + (с, (Я) • с2 (О) • ДО.. .О))
' 1?75 ' ' Ц» ' ' 1?75 ' ' Ц» ' (6)
Поиск орбитальных коэффициентов и собственных значений энергии проводится численными методами решения систем линейных уравнений. Начальными условиями являются С/(0) = сД г = 1,2 (электроположительный и электроотрицательный атомы в молекуле токсиканта),/ = 40 (для пентапептида):
40
^(а-Е^)-С;= 0 (7)
7 = 1
Система уравнений (7) нормируется по орбитальным коэффициентам с учетом волновой функции (ру каждой пары молекулярной орбитали:
40
с.1^ + сп(р2 + • •. + с. 40'(р40 = ^ суф] = 1 (8)
7=1
После того, как закончился поиск энергии адсорбционных систем, осуществляется расчет энергии адсорбции процесса АЕадс:
ЛЕадс = Епош {система) - [Епош {компонент) + Епош (НО,)] (9)
Данная процедура позволяет из «-количества адсорбционных систем отобрать только те, которые обладают значительным запасом энергии и являются стабильными конфигурациями.
Минимум энергии адсорбции свидетельствует о наиболее удачном сцеплении между молекулами. В связи с этим задача сводится к анализу образующихся устойчивых конфигураций, т.е. решается вопрос об активных центрах в структуре биополимера: на каком расстоянии токсиканта от мембраны проявляется повышенная склонность к атаке. Если вос-
пользоваться уравнением:
)^’таг=ел (ю)
>0
где Го соответствует минимальному расстоянию взаимодействующих атомов, которое со-ставляет 1,75 А, а г соответствует максимальному удалению атомов в адсорбционной систе-
ме 1,90 А, Ч> - волновая функция, характеризующая общую молекулярную орбиталь в образующейся системе, то можно рассчитать расстояние (г - г о), при котором образуется оптимальная связь между токсикантом и биополимером клеточной мембраны. Это позволяет установить наиболее устойчивые адсорбционные системы и определить участки биополимеров, которые с наибольшей вероятностью подвержены максимальной атаке со стороны токсиканта. В нашем случае процесс моделирования носит многошаговый итерационный характер. В результате на каждой итерации циклическая процедура моделирования повторяется до тех пор, пока не будет достигнута сходимость.
Геометрическая оптимизация процесса взаимодействия
Для решения задачи геометрической оптимизации установим критерии, при которых обеспечивается образование стабильной формы адсорбционных систем:
• минимум полной энергии Ер или равенство градиента функции энергии нулю;
• минимум энергии адсорбции АЕ;.
Основное ограничение накладывается на расстояние между атомами взаимодействующих молекул, которое лежит в интервале 1,75А < у < 1? 90А.
При моделировании процесса реализована градиентная форма оптимизации с использованием метода следования собственному вектору. Алгоритм следования собственному вектору включает несколько последовательных стадий. На первой стадии определяется стартовая точка поиска, в которой вычисляется значение среднеквадратичного градиента (О).
3 п
где суммирование производится по всем п атомам молекулярной модели; хг, уи 2, - декартовы координаты /-ого атома.
Процедуру повторяют до тех пор, пока полная энергия для двух последовательных итераций не будет совпадать с заданной точностью А, где 1\к> - Е^^ < А. Далее производится поиск критической точки в направлении собственного вектора проверкой условия окончания поиска по заданной величине нормы градиента или по критерию максимального числа итераций.
Молекула с начальными координатами не является равновесной, поэтому критерием окончания поиска точки экстремума на первой стадии является установление «завышенной» величины градиента - 10"2 кДж/(моль-А). Это позволяет выявить локальные стационарные точки, которым соответствуют метастабильные конфигурации. При достижении сходимости между итерациями все изменения в геометрии отражаются в матрице связей для молекул (1), для адсорбционных систем (5). Таким образом, на первой стадии осуществляется многочисленный отбор возможных вариантов образования адсорбционных систем.
Большинство расчетов в первом приближении, как правило, имеет разногласие с параметрами, регистрирующими стабильность конфигураций, поэтому задачу оптимизации можно свести ко второй стадии моделирования. Дальнейшее снижение градиента до 10"5 кДж/(моль-А) позволяет скорректировать геометрию найденных адсорбционных систем с целью обнаружения критической точки на поверхности энергии.
Нами проведена комплексная проверка применимости разработанной модели к изучению процессов воздействия на примере таких компонентов, как сероводород [5], диоксид азота и другие. В результате вычислительный эксперимент показал, что необходимо использовать условия выбора оптимальных конфигураций, это позволяет из п-количества образующихся систем установить наиболее стабильные, обозначив соответствующие атомы (О),
и отсеять менее устойчивые, обозначив ( ), в знак того, что атака со стороны токсиканта
возможна.
Рассмотрим критерии оптимального выбора конфигураций и условия их выполнения, предложенные для решения поставленной цели.
1. Если длина связи между атомами взаимодействующих молекул лежит в интервале 1,750-1,820 А, то осуществляется поиск минимального значения энергии формирования АЕадс, который соответствует наиболее выгодной системе. При выполнении условия взаимодействующий атом молекулы биополимера обозначается (О), это свидетельствует о высокой реакционной способности атома; если условие не выполняется, то процесс переходит на проверку п. 3.
2. Если разность б/; / = (ДЛт||1 - ДЛ'; /) < 5 кДж/моль, то взаимодействующий атом био-
полимера в системеу+7 также обозначается (О); если условие не выполняется, то фиксируется ( ).
3. Если длина связи между атомами взаимодействующих молекул лежит в интервале 1,821—
1,900 А, то осуществляется поиск минимального значения Аиз выбранных систем. При выполнении условия, взаимодействующий атом в молекуле биополимера обозначается (О).
4. Если разность 4+; = (АЛт||1 - АЕ]+1) < 3 кДж/моль, то атом в системе (/ + 1) фиксируется (О); если условие не выполняется, то фиксируется ( ), что свидетельствует о сла-
бой активности атомов вступать во взаимодействия.
В случае если на первом или на третьем этапах присутствует только одна система, предусмотрен дополнительный случай, когда энергия адсорбции ^-системы подлежит сравнительному анализу с другими системами, не входящими в указанный диапазон длины связи. Это позволяет адекватно оценить весь спектр рассматриваемых моделей взаимодействия. В результате на каждом этапе отбирается от одной до нескольких наиболее устойчивых адсорбционных систем.
Компьютерная реализация позволяет продемонстрировать схемы формирования адсорбционных систем на поверхности фрагмента клеточной мембраны (пентапептида) под воздействием диоксида серы.
Рис. 2. Модель установления активных центров на поверхности биополимера (пентапептида)
атаки со стороны токсиканта
Разработанная структурно-адсорбционная модель воздействия низкомолекулярных токсикантов на биополимеры клеточной мембраны ориентирована на изучение адсорбционных процессов и установление активных центров в сложных структурах. Это позволяет представить информацию о возможных механизмах направленного воздействия токсикантов
на объекты окружающей среды. В данной работе мы постарались осветить проблему и необходимость разработки модели, а компьютерная реализация позволяет согласованно подходить к решению поставленной задачи.
Библиографический список:
1. Жарких, Л. И. Квантово-химическое кластерное моделирование процесса адсорбции сероводорода на поверхности белковой мембраны / Л. И. Жарких // Вестник МГОУ. - 2006. - № 9. - С. 56-59.
2. Казанцева, Н. В. (Золотарева) Теоретическое обоснование сорбции диоксида серы на структурных элементах клеточных мембран / Н. В. Казанцева (Золотарева) // Экологические системы и приборы. - 2007. - №9. - С. 35-37.
3. Соловьев, М. Е. Компьютерная химия / М. Е. Соловьев, М. М. Соловьев. - М. : СОЛОН -Пресс, 2005. - 536 с.
4. Хигаси, К. Квантовая органическая химия / К. Хигаси, X. Баба, А. Рембаум - М. : Мир, 1967. - 380 с.
5. Химические приложение топологии и теории графов / под ред. Р. Кинга. - М.: Мир, 1987. - 560 с.
УДК 539.194
СТРУКТУРНО-ДИНАМИЧЕСКИЕ МОДЕЛИ МОНОГАЛОИДОЗАМЕЩЕННЫХ БЕНЗОНИТРИЛА М.Д. Элъкин, А.П. Смирнов, И.И. Гордеев
В рамках гибридного функционала плотности DFT/B3LYP построены структурнодинамические модели моногалоидозамещенных бензонитрила. Исследовано влияние ангар-монизма колебаний, выяснен механизм межмолекулярного взаимодействия.
Ключевые слова: колебательные спектры галоидозамещенных бензонитрила, резонансные эффекты, ангармонизм колебаний, межмолекулярное взаимодействие.
Key words: halogen monosubstituted benzonitrile, vibrational spectra, resonance effects, anharmonicity of vibration, intermolecular interaction.
Теоретический и практический интерес к замещенным бензола как продуктам деятельности целого ряда химических предприятий в последнее время только усилился. В первую очередь, это связано с экологическим мониторингом водных ресурсов. Использование методов оптической физики в данном направлении связано с построением структурнодинамических моделей исследуемых молекулярных объектов.
Анализ структурно-динамических моделей замещенных бензонитрила, основанный на использовании спектральных данных по колебательным спектрам соединений, производился неоднократно. В работе [9] была предпринята попытка обобщить имеющиеся экспериментальные данные по спектрам инфракрасного поглощения и комбинационного рассеяния монометил- и моногалоидозамещенных бензонитрила. Теоретической основой для классификации и интерпретации фундаментальных колебаний монозамещенных бензонитрила являлась классическая схема расчетов [3], не выходящая за рамки гармонического приближения теории молекулярных колебаний [6]. В этой схеме гармонические силовые поля переносятся из молекулярных фрагментов, родственных по электронной структуре исследуемым молекулярным объектам, а затем в рамках методов решения обратных задач уточняется ис-