автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.18, диссертация на тему:Математическая модель процессов взаимодействия диоксида серы со структурными элементами клеточной мембраны

На правах рукописи

2 7 АВ1 2009

ЗОЛОТАРЕВА НАТАЛЬЯ ВАЛЕРЬЕВНА

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПРОЦЕССОВ ВЗАИМОДЕЙСТВИЯ ДИОКСИДА СЕРЫ СО СТРУКТУРНЫМИ ЭЛЕМЕНТАМИ КЛЕТОЧНОЙ МЕМБРАНЫ

Специальность 05.13.18 - Математическое моделирование, численные методы и комплекты программ

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата технических наук

Астрахань-2009

003475723

Работа выполнена па кафедре аналитической и физической химии химического факультета государственного образовательного учреждения высшего профессионального образования «Астраханский государственный университет».

Научный руководитель: доктор химических наук, профессор

Алыков Нариман Мирзаевич

Официальные оппоненты: доктор физико-математических наук,

профессор

Элькин Михаил Давыдович

кандидат физико-математических наук, доцент

Коваленко Илья Борисович

Ведущая организация: Московский государственный

университет экономики, статистики и информатики, г. Москва

Защита состоится «25» сентября 2009 г. в 13:00 часов на заседании диссертационного совета ДМ 212.009.03 при Астраханском государственном университете по адресу: 414056, г. Астрахань, ул. Татищева, 20 а, конференц-зал.

С диссертацией можно ознакомиться в библиотеке Астраханского государственного университета.

Автореферат разослан «^ I

Ученый секретарь диссертационного /

совета, к.т.н. [/ж// ^ О. В. Щербинина

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

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

Цель и задачи исследования. Целыо данной работы является разработка математической модели процессов взаимодействия диоксида серы со структурными элементами клеточной мембраны. Это позволит выявить активные центры на поверхности макромолекулярной системы, подверженной воздействию токсиканта. Для реализации поставленной цели предусмотрено решение следующих задач:

• разработка модели описания взаимодействий в макромолекулярных системах и выбор программного обеспечения для решения поставленной цели;

• разработка алгоритма взаимодействия диоксида серы с компонентами клеточной мембраны на основе математической модели;

• разработка приёмов оптимизации макромолекул и сопоставление результатов расчета с известными данными с целыо получения геометрически оптимальной конфигурации систем;

• создание базы данных и программного модуля с целыо построения молекулярных диаграмм биополимеров клеточной мембраны;

• разработка программного модуля с целью формирования схем взаимодействий диоксида серы с биополимерами мембраны.

Методы исследований. Для создания «первичной» молекулярной модели использована программа Морас, расчеты в которой проведены полуэмпирическим РМЗ методом. Для дальнейшей оптимизации, корректировки г-матрицы и поиска энергетических параметров исходных и моделирующих систем использована программа РС Сатезз, расчёты в которой проведены с применением РМЗ метода и базиса 6-31 в в хартри-фоковском приближении. Задача геометрической оптимизации молекул решена при использовании метода следования собственному вектору.

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

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

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

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

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

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

На защиту выносятся следующие положения:

• математическая модель и алгоритм математического моделирования процессов взаимодействия диоксида серы со структурными элементами клеточной мембраны;

• база данных и программный модуль BíoMolDiagrams для составления молекулярных диаграмм биополимеров мембраны;

• квантово-химические расчеты для установления геометрии молекул и сопоставление результатов с известными данными;

• программный модуль МоскШиегасйопз для моделирования взаимодействий в макромолекулах;

• методика определения активных центров на поверхности структурных компонентов мембраны, максимально подверженных воздействию со стороны диоксида серы.

Практическая значимость. На основании математической модели и алгоритма разработан программный комплекс Molecular Model (исх. № 17-4 от 13.07.2009), позволяющий:

• структурировать расчеты макромолекул в виде молекулярных диаграмм, которые необходимы химикам-органикам для синтеза новых соединений;

• формировать схемы взаимодействий, на основании которых появляется информация о потенциальных «мишенях» в биологических мембранах, необходимая химикам для установления антидотов к токсичным соединениям, а биоэкологам для изучения экологического последствия.

Разработанный алгоритм впервые был успешно внедрен в область моделирования взаимодействий компонентов на установках очистки газового конденсата от кислых газов.

В соавторстве по данному направлению созданы и зарегистрированы базы данных:

• № 2009620009 «Токсическое воздействие на биологические структуры»;

• № 2009620103 «Воздействие блокаторов Н-рецепторов и антигис-таминных препаратов на биологические структуры»;

• № 2009620395 «Молекулярные диаграммы структурированных поверхностей»;

• исх. № 17-4 от 13.07.2009 «Структурные характеристики биополимеров и двухкомпонентных взаимодействующих систем», которые нашли свое применение при создании новой учебной дисциплины для вузов «Квантовая экологическая химия». Этот предмет изучают студенты химического факультета Астраханского государственного университета. В результате, по программе «У.М.Н.И.К. - 2009» получен грант.

Результаты диссертационной работы используются в учебном процессе по дисциплинам: «квантовая механика и квантовая химия», «экология на стыке математики, физики и химии».

Апробация работы. Результаты исследований доложены на различных Международных и Всероссийских научных конференциях, среди которых: «Средства и методы обеспечения экологической безопасности» (Астрахань, 2005) «Фундаментальные и прикладные проблемы современной химии в исследованиях молодых ученых» (Астрахань, 2006); III школа-семинар «Квантовохимические расчеты: структура и реакционная способность органических и неорганических молекул» (Иваново, 2007); «Фундаментальные и прикладные проблемы современной химии» (Астрахань, 2008); «Инновационные технологии в управлении, образовании, промышленности АСТИНТЕХ-2009» (Астрахань, 2009); IV школа-семинар

«Квантово-химические расчеты: структура и реакционная способность органических и неорганических молекул» (Иваново, 2009).

Публикации. По материалам диссертации опубликовано 5 статей в журналах, рекомендованных ВАК и 5 статей в материалах конференций.

Структура п объём диссертации. Диссертация состоит из введения, трёх глав, заключения, приложения и библиографического списка (132 наименования). В приложении имеется 4 акта о внедрении результатов диссертационных исследований в учебный процесс.

СОДЕРЖАНИЕ РАБОТЫ

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

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

К настоящему времени скопился большой арсенал методов моделирования, в основе которых лежит решение уравнения Шредингера для стационарных состояний:

НЧ'((1) = ЕЧЧдг) (1)

где - волновая функция, удовлетворяющая вероятностному уравнению:

}} }|//:ауау<ъ = I (2)

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

Наиболее ранняя модель, описывающая реакционную способность молекул, представлена в статическом методе [Коулсон и Лонге-Хиггинс, 1947 г.]. Расчеты сводятся к определению параметров тг-электронной плотности (лу) атомов, степени л--электронного связывания (/>„), индексов свободной валентности (/•',) и значений самополяризуемости (Пгг) атомов:

т..........V,,.......

= (3) рп = -с„-си (4)

|«| 1=1

............(5) (6)

' / £1 е>

где !/>„„„„„ - занятая молекулярная орбиталь, и, — число электронов на /-ой орбитапи, с„ - орбитальный коэффициент, /Утах - максимальное число связей атома /•, /V, - сумма порядков связей атома г, е, - энергия электрона /-ой орбитали.

Согласно методу, взаимодействию подвергается тот атом, который обладает максимумом кг или максимумом //,,. (для злектрофильных агентов), минимумом л> или максимумом //„ (для нуклеофилЕ.ных агентов), максимумом Гг (для свободных радикалов).

Данная модель ограничена анализом одной молекулы, тем самым, неадекватно оценивается ситуация в среде реагента. Аппарат не может быть использован в тех ситуациях, когда сопряженная молекула содержит гете-роатомы или нарушается ароматичность, когда возможны коиформацнон-ные переходы. Основное достоинство метода - это простота математического аппарата.

Рассмотренная модель является базой, принципы которой заложены в модель, предложенную нами для описания взаимодействий в макромолекулах. На рисунке 1 приведена структура алгоритма моделирования взаимодействия диоксида серы с компонентами мембраны.

начало

| Выбор компонентов мембраны |

1 1

Углевод Ли л ид | ( Белок |

" ; ' 1 ' }

нет

[_МыГщр I nKL'ilк;ш M___1

[_,11N)Î.4JI,I(4-|ÎJ.Î ]

Сраммител аналт и ныГюр метода расчета

Математическое моделирование шаимодейстнин

X

-N групп (компонент + SO2) |

,- î: ^ ------------

Расчет структурных характеристик онтимишроклиных моделей

- п:= п+1 ' I! **'.V ,. нет ^_

Выбор оптималь- да -ш.|\ поли-ксюш

Формирование с\ем|.

конец

Рис. 1. Структура алгоритма математического моделирования нчаимодейстипя диоксида оеры со структурными элементами клеточной мембраны

На этапе оптимизации исходных компонентов предлагается разбить процесс на две стадии. На первой стадии, для расчета полной энергии мак-

ромолекулы, можно использовать РМЗ метод. Большинство расчетов в первом приближении, как правило, имеют разногласия с экспериментальными данными, поэтому задачу можно свести ко второй стадии моделирования. Для дальнейшей обработки г-матрицы и уточнения геометрии предлагается использовать расширенный базис 6-310 в приближении самосогласованного поля Хартри-Фока.

Выделяются следующие критерии оптимизации:

• критерий минимума полной энергии Е или равенство градиента функции энергии нулю;

• критерии геометрических параметров (длина связи, валентный угол);

• критерий минимума энергии адсорбции АЕшк..

Важной характеристикой регистрации стабильной формы является задание значения градиента функции энергии, что соответствует норме изменения /-ой координаты.

где суммирование производится по всем п атомам; у„ - декартовы координаты /' - ого атома.

Задача геометрической оптимизации процесса сводится к поиску наиболее устойчивых структур, которым соответствует минимум энергии, поэтому она может быть решена методом следования собственному вектору. Выбор метода основан на рациональном использовании ресурсов ЭВМ и времени расчетов. Критерием окончания поиска точки экстремума на первой стадии является установление «завышенной» величины градиента 10'" кДж/(моль-А), это позволяет выявить локальные стационарные точки, которым соответствуют метастабильные конфигурации. Дальнейшее снижение градиента до 10"? кДжДмоль-А) позволяет скорректировать геометрию найденных конформаци-онных систем с целью обнаружения критической точки. Рассчитанные оптимизированные геометрические показатели (длина связи, валентный угол) должны быть сопоставимы с экспериментальными или справочными данными.

Процесс продолжается до тех пор, пока не будет оптимизирована молекула и выбран подходящий метод расчёта. При выполнении условия, процесс переходит на этап моделирования взаимодействия, здесь производится построение и оптимизация возможных комбинаций взаимодействия диоксида серы с компонентами мембраны.

На этапе расчёта структурных характеристик оптимизированных моделей следят за условными изменениями в электронной конфигурации по значениям зарядов на атомах (с/„ е), определяющих интегральное значение электронной плотности вблизи каждого атома. Об изменениях электронной плотности, при образовании системы из отдельных подсистем, судят по величине переноса заряда Д(/ с биополимера на диоксид серы:

НМХ С,гаЛеп1=(ЪпУ' £

/

(7)

Д<7 = </(0,) + </(52) + <7(0,)

(8)

Тепловые эффекты (Д/w, кДж/моль) формирования предполагаемых систем рассчитывали как:

Д£„а — Ект(систелш) — [кП1т{компонешн)-ь /'-",„,,„ (МЛ)] (9)

По соотношениям величины ДЕм. и длины связи (/•,,) между атомами, взаимодействующих молекул, судят о стабильности образующихся систем. Перед формированием схем взаимодействий проводится анализ п-количества систем, из которых осуществляется отбор. Рассмотрим критерии оптимального выбора положений и условия их выполнения.

1. Если /•/, в системах лежит в интервале 1,750-1,820 А, то осуществляется поиск минимального значения энергии адсорбции, который соответствует наиболее выгодной конфигурации. При выполнении условия, взаимодействующий атом в молекуле биополимера обозначается

Если условие не выполняется то процесс переходит на проверку п. 3.

2. Если Д'■./ = (àEmi„ - АЕ,./} < 5 кДж/моль, то атом в (/1 /)-снстеме также обозначается ( ' ).

Если условие не выполняется, то над атомом фиксируется ( —).

3. Если Гц (1,821 - 1,900 Â), осуществляется поиск минимального значения ДЕ^.. При выполнении условия, взаимодействующий атом в молекуле биополимера обозначается ( ).

4. Если Д'/./ = (Д£"тш - /) < 3 кДж/моль, то атом в системе (iU) фиксируется ('■').

Если условие не выполняется, то над атомом фиксируется ( —»- ).

Итогом всего процесса является формирование схем взаимодействий для структурных компонентов клеточной мембраны.

Во второй главе представлена реализация модели на ЭВМ, которая позволяет описать структуру и свойства макромолекул. Комплексы программ Морас и PC Gamess реализуют только численные методы оптимизации геометрии молекул, а разработанная автоматизированная система Molecular Model позволяет генерировать результаты и составлять молекулярные диаграммы биополимеров (модуль BioMolDiagrams). Пользовательский интерфейс программы реализован в среде Microsoft Visual Studio 2008 на языке программирования С#.

После оптимизации макромолекул, из выходного текстового файла PC Gamess проводится выборка данных и внесение их в базу данных. Для построения использовалась система управления базами данных FireBird, версии 2.1. На рисунке 2 приведена схема связей данных.

Структура базы данных представляет собой трехуровневую схему, это позволяет хранить информацию о структуре молекул, о функциональных группах и о свойствах атомов. Для формирования схем молекулярных диаграмм база данных содержит таблицы с информацией о базовых структурах (basic_struct), об узлах-атомах (basicjiode) и о длинах связей (basicjine).

Сне. 2. Схема спичей данных

Наименование исследуемых структур хранится в таблице typel, а данные о методе расчета в таблице method_analysis. Таблица structure_one содержит информацию об энергетических параметрах молекулы.

Величина теплоты образования и общая энергия автоматически заполняют базу данных только в том случае, если достигнута сходимость и найдена оптимальная конфигурация молекулы. При удачном окончании процесса оптимизации в программе Gamess, текстовый документ содержит запись "Equilibrium Geometry Located". Если структура не достигла оптимального положения, то процесс запускается повторно с новыми координатами.

Для характеристики функциональных групп в таблице group_one приведены поля их графического расположения в системе координат ху. Информация о параметрах зарядов на атомах (с/,) и значения электронной плотности (тг,) хранится в таблице atom_one. Длина химических связей (/",,) между атомами в исследуемой структуре хранится в таблице connection.

Разработанный программный модуль BioMolDiagrams реализует возможность графического изображения молекул, но поскольку изучаемые системы - клеточные мембраны, представляют громоздкие конструкции, то возникает необходимость уменьшать размеры объектов. Поэтому жид-костно-мозаичная структура мембраны дифференцируется на небольшие составляющие фрагменты - белки, липиды и углеводы, а для упрощения составления моделей полипептидов, в программном модуле реализована настройка выбора аминокислот из списка, это позволяет сэкономить время при работе с макросистемами. Хранение всего списка осуществляется в таблицах «базовые параметры» (рис. 2).

Входными параметрами модуля BioMolDiagrams являются:

• вычисленные геометрические параметры молекулы в виде z-матрицы для конкретного метода расчета по программе Морас;

• вычисленные энергетические и зарядовые характеристики молекулы в программе Морас;

• квантово-химический метод расчета;

• название, формула и тип соединения.

• Выходными данными модуля являются:

• таблица оптимизированных геометрических параметров молекулы;

• минимизированные энергетические параметры молекулы;

• зарядовые характеристики молекулы;

• двухмерная молекулярная модель.

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

В качестве белкового компонента выбрана молекула пентапептида, содержащая различные группировки (рис. 3). "он

кэн .5

7 I М-г' "Ч*

,3сн, I язн

• ! „ !

|!СНг й-р сн,

I с — сн—мн--°с-"сн-мн---с—'ш-'^н "с 21сн г1Ут-2'с ""'сн "мн,

но» 4 4 4 4

Рис. 3. Модель пентапептида

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

Таблица 1

Сравнение результатов расчета геометрии индивидуальных молекул и модели пентапептида методами РМЗ и ПГ/6-31С*.

Параметры Справочные результаты Результаты расчеши 1шпш1пд\;ш.ны\ молекул

РМЗ | 11176-31 С'

Ы-ацетилгмщин »0 ,чс-2сн; 'Мн '-"с "сн, НО;, 3,"

ЦС'-Сг). А 1,510±0,04 1.513 1,505

/(С'=0"), А 1,190±0,03 1,219 1,209

/(С'ЧГ), А 1,310±0,03 1,351 1.3-19

/(С'у-ю, А 1,320±0,005 1,436 1,364

124,00±3 116,26 122.75

|20,00±3 ¡20,21 121.6-1

п-крезо.п 4>н в4— ' 1

/(С„-Н), А 1,110±0,02 1,095 1.072

/(С"-СГ), А 1,360±0,01 1,369 1.381

/(С„-С ,„), А 1,392±0,005 1,395 1.3,48

УС С С ° 120,00 120,00 120.00

метшмеркаптан сн3—эн

/(С'-'-Б"), А 1,810±0,001 1,801 1,882

/(5Н-Н), А 1,329±0,004 1,306 ¡.334

Из таблицы видно, что РМЗ метод хорошо воспроизводит параметры С-С, Сар-О, С=0, С-Н и С-Б связей. Намного сложнее метод описывает связи с гетероатомами, сопряженными с другими функциональными группами, например, длина связи N-0 в пептидной группе завышена на 0,1 А, а валентный угол, включающий атом азота или атом кислорода, как правило, оказывается, занижен на 5-9°. Расширенный базис в хартри-фоковском приближении намного лучше воспроизводит длину связей и валентные углы с гетероатомами, но напротив не рассчитан на описание серосодержащих и ароматических соединений (табл. 1).

Аналогично проведены расчеты для липидной составляющей - структуры моноглицсрида, диглицерида, триглицерида и углеводной составляющей - олигосахарид, лактоза. Сопоставление результатов, позволило установить правомерность использования РМЗ метода и расширенного базиса 6-31 в в хартри-фоковском приближении. Однако полуэмпирический метод, в отличие от расширенного базиса, справляется с поставленной задачей в течение нескольких минут. В таблице 2 приведены данные по затрате времени для установления равновесной конфигурации биополимеров мембраны разными подходами.

Таблица 2

№ Молекула Время расчетов, ч

РМЗ HF/6-31G'

1 Пснтамептнл 0.16 528

2 3 Моиопшцсрид Дпглнцсрмд 0,03 0,66

0.06 4,50

4 Тршлицерид 0,10 11.33

5 Олигосахарид 0.03 2,33

Из таблицы 2 следует, что использование результатов расчета (6-31G ) для формирования схем молекулярных диаграмм может быть ограничено только по времени. В результате, для всех биополимеров в программном модуле Bio Mol Diagrams составлены молекулярные диаграммы, которые могут быть представлены в виде наглядных иллюстраций (рис. 4). но—— "

И,

лС — CH — NH — с-

-111)1 -IUH ÍI2III-

но( I Д-

"н\ ,,„ УГп,

SH

СН2 с-н, у сн2

CH-NH-C-ril-NH-C-CH-NH-C-CH-NH, Ii II II

ООО

диаграмма пентапептида

сн2—сн— SI 15 он о I 1.211 ^ t o=c

(1..10 1-41,5.1 I

\ m

»v5H

CH

51

1.15— H, С

0

1

c=o ({hJ2

CH

II

CH I

CH,

диаграмма диглицерида

¿¿a

H,C-OH

диаграмма о.ччгосахарида

Рис. 4. Молекулярные диаграммы пштансптнда, диглицерида и олигосахарида

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

В третьей главе представлена реализация предложенного алгоритма моделирования взаимодействия диоксида серы со структурными компонентами мембраны. Созданная база данных хранит информацию о параметрах взаимодействующих структур. Приведённая на рисунке 5 схема иллюстрирует связывание данных в таблицах.

Рис. 5. Схема свяжи данных взаимодействующих молекул

Энергетические параметры образующейся адсорбционной системы приведены в таблице interaction. Параметры электронной плотности, зарядовые характеристики атомов и информация об их реакционной способности хранится в таблице atom two. Все параметры получены из текстового файла программы Gamess. Таблица connection_two содержит информацию о длине образующейся связи при взаимодействии молекул.

Принцип подготовки входных данных при работе с адсорбционными системами остается таким же, как и при работе с одной структурой. Программный модуль Modellnteractions реализует расчёт энергии формирования данных систем (ДЕидс) и величины переноса зарядов (Aq) на диоксиде серы, а анализ параметров позволяет оценить активность атомов в макро-

молекулах и графически сформировать схемы взаимодействия. В результате, статическая модель включает специальные обозначения активности атомов (О, —►).

Входными параметрами модуля МоскПтегасйопБ являются:

• геометрические параметры системы в виде г-матрицы для конкретного метода расчета по программе Морас;

• вычисленные энергетические и зарядовые характеристики системы в программе Морас;

• квантово-химический метод расчета;

• название и формула соединения;

• длина связи между атомами взаимодействующих молекул. Выходными данными модуля являются:

• таблица оптимизированных геометрических параметров системы;

• минимизированные энергетические параметры системы;

• зарядовые характеристики атомов в системе;

• характеристика активности атомов;

• двухмерная молекулярная модель взаимодействия.

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

Рассмотрим модели адсорбционных систем и схемы взаимодействий пентапептида ^истегшил-фепчлапаичл-ачанш-гачог^стеитш-тирозин). Для уменьшения загромождения в системы введены обозначения Я, Я, обозначающие остаток углеродного скелета.

СООН 6 <0 1

но-V/ ч V. 45 "о=з2

Система 1 Система 2 Система 3

В таблице 3 приведена сравнительная характеристика геометрических и энергетических параметров изучаемых адсорбционных систем 1-3, по-

лученных РМЗ методом.

Таблица 3

Параметры адсорбционных систем «пентапептид~802»

Системы Г, А Дд, с АЕ„дс, кДж/моль

1 0З7...52/Н48...0, 1,832/1,784 -0,210 -37,12

2 н,5...о, 1,803 0,022 -32,84

3 Н54...0, 1,834 0.010 -28,08

Полученные результаты позволяют оценить энергию взаимодействия, а так же выяснить наиболее активный центр пентапептида, способный легко подвергаться направленному влиянию со стороны S02. В первой системе атака происходит по пептидной группе, устойчивость определяется за счет образования шестичленного цикла. Минимальная величина переноса заряда в системе / (-0,210 е) может быть вызвана сильной поляризацией связи >S=0, за счет увеличения на кислороде О1 электронной плотности д,(-0,6474)/д(-0,6979). В системе 2 гидроксогруппа открыта для взаимодействия, поэтому наблюдается перенос заряда на Z/45 </(0,2013)/д(0,2483). Система 3 образована за счет открытого метиленового радикала.

При исследовании 35-ти моделей адсорбционных систем, было отмечено, что в среднем энергия адсорбции диоксида серы на пентапептиде может составить -21,00 кДж/моль. В результате расчетов и выполнения условий выбора оптимальных положений для пентапептида в программном модуле Modellnteractions составлена схема взаимодействий (рис. б).

"Л»

\VH \ г

нт у w

о, СМ I 1 " у

С N С и

Л

I

%

/"' т

ж^-Ги»-

Рис. 6. Схема межмолекулярпых взаимодействий ментанситнда

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

Ниже приведены модели систем, образованных одной из форм липида (|дигексенат глицерина).

0=8, \

■ЮЬ „ ,.1*

0=*-'-</. Чг'нТ N^4,

■1

Система 4 Система 5 Система 6

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

Таблица 4

Параметры адсорбционных систем «липид-БОг»_

Системы г, h Дг/. е 4/w, кДж/моль

4 On...S2 1,864 -0,185 -24,15

5 S2...06/03...H26 1,868/ 1,788 -0.112 -15,60

6 О3...Н3» 1,861 ^ 0.006 -11,98

Полученные результаты позволяют оценить энергию взаимодействия, а так же выяснить наиболее активные центры в структуре диглицерида. Система б образована за счет протона концевой метальной группы. Устойчивость системы 5 определяется возможностью образования четырехчленного цикла. Межмолекулярная связь в системе 4 образована за счет серы, устойчивость данной структуры определяется за счет того, что сера способна принять электронную пару от кислорода на d-орбиталь, на это указывает увеличение длины связи 0"=С8 (межъядерное расстояние в молекуле диглицерида до взаимодействия составляет 1,21 Л; а после 1,25Â).

Анализ 35 возможных систем позволил установить, что в реальной ситуации сорбируемость диоксида серы на диглицериде низкая и составляет -16 кДж/моль. На рисунке 7 графически представлена схема взаимодействий молекулы диглицерида с диоксидом серы.

\

Il CEI ^

о Г--Г11, ( ау

^ \ / \ Y/S /

Y-о о-'с нс-^г

п/ \ 7 V

"'С !!,(■■ ai V

N \

СИ СИ IK

\<••: 4

Рис. 7. Схема межмолскуляриых взаимодействий диглицерида

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

Аналогично рассматриваем системы, образованные за счёт взаимодействия диоксида серы с поверхностью лактозы. Ниже приведены равновесные геометрии наиболее выгодных положений в системе «диоксид серы -лактоза».

, . о—л

R .н" о ..о Г

И / '

° J н

v \Ан

3 \ Система 7 Система 8 Система 9

В таблице 5 приведена сравнительная характеристика геометрических и энергетических параметров оптимизированных систем.

Таблица 5

Па эаметры адсорбционных систем «лактоза-502»

Системы г, А е ЛЕаОс, кДж/моль

7 0,Н...0| 1,809 0,017 -33,05

8 ОцН...О| 1,806 0,020 -32,26

9 С,Н...О| 1.838 0,011 -30,44

Расчеты показывают, что системы 7 и 8 занимают энергетически более выгодные положения, так как атака происходит по гидроксильной группе, которая является открытой для взаимодействия. Расстояние между взаимодействующими атомами в системе 9 увеличивается, о чем свидетельствует значение г в таблице 5.

При анализе 22 систем установлено, что сорбируемость диоксида серы на лактозе высокая, в среднем АЕшк составляет -27,20 кДж/моль. На рисунке 8 показаны активные центры в молекуле олигосахарида.

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

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

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

Разработанный подход математического моделирования может быть использован при изучении свойств антидотов, с целью расчета энергетических параметров процессов дезактивации токсикантов с поверхности клеточной мембраны и блокирование их активных центров. Разница в энергетических характеристиках до и после введения антидота позволит ответить на вопрос применимости выбранного препарата и судить об эффективности его использования в медицине. Необходимость развития данного направления обусловлена применением разработанного подхода в «горячих точках» химии, экологии и медицины.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ

1. Впервые разработана математическая модель взаимодействия диоксида серы с макромолекулами клеточной мембраны, которая ориентирована на изучение адсорбционных процессов и установление активных центров в сложных биополимерных структурах, это позволяет предоставить информацию химикам и экологам о возможных механизмах направленного воздействия токсикантов на объекты окружающей среды.

2. На основании математической модели разработан алгоритм, который может быть использован для моделирования процессов взаимодействия других токсичных соединений с компонентами мембраны.

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

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

5. Разработана база данных и автоматизированная система Molecular Model, состоящая из нескольких модулей. Программный модуль Bio-MolDiagrams позволяет наглядно иллюстрировать биополимеры клеточной мембраны с помощью молекулярных диаграмм. Модуль Modellnteractions предназначен для формирования схем взаимодействий, в результате чего появляется возможность получения информации о потенциальных «мишенях» воздействия токсичного вещества на биологические системы и изучения последствий экологического риска.

6. Разработана методика определения активных центров на поверхности структурных компонентов мембраны, что позволяет устанавливать очередность воздействия диоксида серы, и может быть применена для других токсикантов.

ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

Стятьи в ведущих рецензируемых научных журналах и изданиях, включенных в список ВАК РФ

1. Казанцева, Н. В. (Золотарева) Теоретическое обоснование сорбции диоксида серы на структурных элементах клеточных мембран [Текст] // Экологические системы и приборы. - 2007. -№ 9. - С. 35-37.

2. Казанцева, Н. В. (Золотарева) Квантово-химическое кластерное моделирование адсорбции фенола на поверхности алюмосиликатов [Текст] / Н. В. Казанцева, Л. И. Головацкая, Н. М. Алыков // Экологические системы и приборы. - 2005. - № 9. - С. 27-29.

3. Казанцева, Н. В. (Золотарева) Квантово-химическое моделирование хемоеорбции диоксида серы на структурных элементах клеточных мембран [Текст] / Н. В. Казанцева, Н. Н. Алыков // Изв. вузов. Химия и хим. технология. - 2007. - Т. 50, № 12.-С. 132-133.

4. Казанцева, Н. В. (Золотарева) Квантово-химическое моделирование хемоеорбции диоксида серы на поверхности энергетического субстрата клетки - олигосахарида / Н. В. Казанцева, Н. Н. Алыков // Вестник МГОУ. Сер. Естественные науки. - 2006. - № 2 (24). - С. 179.

5. Алыков, Н. М. Изучение соединений, содержащих диэтаноламин, ингибиторы кислотной коррозии и продукты их разрушения [Текст] / Н. М. Алыков, М. Н. Котельникова, Н. В. Золотарева, С. Н. Фидурова // Технологии нефти и газа. 2008 г. № 1. С. 23-29.

Статьи в материалах международных и всероссийских конференций

6. Алыков, Н. М. Квантовая экологическая химия - новая учебная дисциплина в системе университетского образования [Текст] / Н. М. Алыков, Л. И. Жарких, Н. В. Золотарева // Инновационные технологии в управлении, образовании, промышленности «АСТИНТЕХ-2009» : материалы Международной научной конференции (11-14 мая 2009 г.). -Астрахань, 2009.-С. 125-127.

7. Золотарева, Н. В. Моделирование биополимеров клеточной мембраны и квантово-химический расчет параметров реакционной способности [Текст] // IV школа-семинар: Квантовохимические расчеты: структура и реакционная способность органических и неорганических молекул (20-22 мая 2009 г.). - Иваново, 2009. - С. 72-76.

8. Золотарева, Н. В. Квантово-химическое моделирование сорбции диоксида серы. Математический аппарат образования межмолекулярных взаимодействий [Текст] / Н. В. Золотарева, Е. В. Шитоева, В. И. Юртаева // Фундаментальные и прикладные проблемы современной химии : материалы II Международной конференции (15-17 апреля 2008 г.) -Астрахань, 2008.-С. 177-181.

9. Золотарева, Н. В. Математическое моделирование процессов взаимодействия диоксида серы со структурными элементами клеточной мембраны [Текст] // Фундаментальные и прикладные проблемы современной химии : материалы II Международной конференции (15-17 апреля 2008 г.). - Астрахань, 2008. - С. 167-177.

10. Казанцева, Н. В. (Золотарева) Квантово-химические расчеты энергий взаимодействия диоксида серы со структурными элементами клеточных мембран [Текст] // III школа-семинар: Квантовохимические расчеты: структура и реакционная способность органических и неорганических молекул (14 марта 2007 г.). - Иваново, 2007. - С. 86-90.

Регистрация интеллектуальной собственности

11. Свидетельство о регистрации базы данных. № 2009620009 Токсическое воздействие на биологические структуры [Текст] / Л. И. Жарких, Н. М. Алыков, Н. В. Золотарёва: заявитель и патентообладатель Астраханский государственный университет. - № 200862034! ; заяв. 27.10.08 ; опубл. 11.01.09.

12.Свидетельство о регистрации базы данных. № 2009620103 Воздействие блокаторов Н-рецепторов и антигистаминных препаратов на биологические структуры [Текст] / А. И. Юсупов, Н, В. Золотарёва, Н. М. Алыков, Л. И. Жарких: заявитель и патентообладатель Астраханский государственный университет. -№2008620452 ; заяв. 31.12.08 ; опубл. 27.02.09.

13. Свидетельство о регистрации базы данных. № 2009620395 Молекулярные диаграммы структурированных поверхностей [Текст] / Н. В. Золотарёва, Н. М. Алыков : заявитель и патентообладатель Астраханский государственный университет. - № 2009620306 ; заяв. 27.05.09 ; опубл. 24.07.09.

Заказ № 1869. Тираж ЮОэкз.

_Уч.-изд. л. 1,3. Усл. печ. л. 1,2._

Оттиражировано в Издательском доме «Астраханский университет» 414056, г. Астрахань, ул.Татищева, 20 факс (8512) 25-17-18, тел. (8512) 54-01-87, 54-01-89; E-mail: asiipress@yandcx.ru

ВВЕДЕНИЕ.

ГЛАВА I. ОБЗОР ИНФОРМАЦИИ И ПОИСК ОПТИМАЛЬНЫХ МЕТОДОВ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ ПРОЦЕССОВ МЕЖМОЛЕКУЛЯРНОГО ВЗАИМОДЕЙСТВИЯ.,.

1.1. Математическое описание молекулярных моделей.

IЛ. 1. Анализ неэмпирических методов.

I.1.2. Анализ полуэмпирических методов.

IЛ .3. Анализ точности квантово-химических методов.

1.2. Методы геометрической оптимизации.

1.3. Математическое описание модели взаимодействия.

1.4. Предлагаемая математическая модель.

1.5. Компьютерная реализация квантово-химических методов.

Выводы первой главы.

ГЛАВА И. МОДЕЛИРОВАНИЕ И ОПТИМИЗАЦИЯ СТРУКТУРНЫХ

КОМПОНЕНТОВ РЕАКЦИОННОЙ СРЕДЫ.

И. 1. Структура базы данных однокомпонентных систем.

И.2. Реализация пользовательского интерфейса Molecular Model.

II.3. Моделирование и оптимизация исходных структур.

II.3.1. Моделирование и оптимизация белкового компонента.

11.3.2. Моделирование и оптимизация липидных компонентов.

11.3.3. Моделирование и оптимизация углеводного компонента.

И.З.4. Моделирование и оптимизация диоксида серы.

Выводы второй главы.

ГЛАВА III. МОДЕЛИРОВАНИЕ ПРОЦЕССОВ МЕЖМОЛЕКУЛЯРНОГО ВЗАИМОДЕЙСТВИЯ ДИОКСИДА СЕРЫ С БИОПОЛИМЕРАМИ

КЛЕТОЧНОЙ МЕМБРАНЫ.

III. 1. Структура базы данных двухкомпонентных систем.

111.2. Программный модуль Modellnteractions.

111.3. Моделирование межмолекулярных взаимодействий.

Ш.3.1. Моделирование в системе «диоксид серы - пентапептид».

111.3.2. Моделирование в системе «диоксид серы - липид».

111.3.3. Моделирование в системе «диоксид серы — олигосахарид».

Выводы третьей главы.

Актуальность темы

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

В последние десятилетия экологические проблемы отличаются остротой наблюдающихся процессов и явлений, вызванных мощным антропогенным давлением на природные системы. Ни для кого не секрет, что в настоящее время основным источником токсичных выбросов (SO2, H2S, СО и др.) является нефтедобывающая отрасль. Зонами повышенного внимания являются продувка скважин, факельное хозяйство, площадки розлива жидкой серы, установки Клауса типа "Сульфрен" и различного рода технологические утечки газа.

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

К настоящему времени сложились отдельные математические модели и методы, описывающие взаимодействие в основном низкомолекулярных соединёний [5, 10, 19, 25, 26, 44, 82, 56], в которых существующий математический аппарат позволяет: смоделировать поведение атомов в системе, теоретически исследовать геометрию, электронную структуру и энергетические характеристики стабильных систем. Но по мере усложнения структуры, за счет увеличения количества атомов в системе, данная задача становится практически невыполнимой при использовании персональных компьютеров. Поэтому, при решении обширных химико-экологических задач, связанных с моделированием процессов взаимодействия в макросистемах, приходится искать новые подходы исследования, отталкиваясь от существующих моделей с привлечением численных методов и комплексов квантово-химических программ. Цель и задачи исследования

Целью данной работы является разработка математической модели процессов взаимодействия диоксида серы со структурными элементами клеточной мембраны. Это позволит выявить активные центры на поверхности макромолекулярной системы, подверженной воздействию токсиканта.

Для реализации поставленной цели предусмотрено решение следующих задач: разработка модели описания взаимодействий в макромолекулярных системах и выбор программного обеспечения для решения поставленной цели; разработка алгоритма взаимодействия диоксида серы с компонентами клеточной мембраны на основе математической модели; разработка приёмов оптимизации макромолекул и сопоставление результатов расчета с известными данными с целью получения геометрически оптимальной конфигурации систем; создание базы данных и программного модуля с целью построения; молекулярных диаграмм биополимеров клеточной мембраны: разработка программного модуля с целью формирования; схем взаимодействий диоксида серы с биополимерами мембраны.

Методы исследований!

Для созданиякспервичной» молекулярной; модели-использованашрограмма Морас, расчеты в которой проведены полуэмпирическим РМЗ методом; Для дальнейшей оптимизации, корректировки z-матрпцы и. поиска энергетических параметров исходных и моделирующих систем: использована: программа PG Gamess, расчёты в которой проведены с применением РМЗ? метода и расширенного базиса 6-31G* в хартри-фоковском приближении:. Задача геометрической оптимизации молекул; решена, при5 использовании, метода следования собственному вектору. !

Научная новизна;

Впервые разработана? математическая модель процессов' взаимодействия . диоксида серы с макромолекулярными структурами, которая-позволяет изучать адсорбционные процессы, происходящие на поверхности клеточной мембраны. V

Разработан алгоритм, который включает проведение: пошаговой оптимизации геометрии макромолекул до и после взаимодействия;

Для? макромолекул установлены; критерии; оптимизации, фиксирование которых позволяет получить геометрические параметры, сопоставимые с экспериментальными; и справочными данными. Впервые, на этапе моделирования взаимодействий, установлены критерии выбора оптимальных положений, что позволяет в- структуре, биополимера определить активность взаимодействующих атомов. ■ ' .

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

Разработана концептуальная, модель базы данных, которая позволяет структурировать результаты квантово-химических расчетов молекулы диоксида серы, компонентов • клеточной мембраны — полипептидов, липидов, полисахаридов, а также образующихся в результате взаимодействия адсорбционных систем.

Разработана методика, которая позволяет устанавливать очередность воздействия диоксида серы на структурные компоненты мембраны. Практическая значимость

На основании математической модели и алгоритма разработан программный комплекс Molecular Model (исх. №17-4 от 13.07.2009), позволяющий:

- структурировать расчеты макромолекул в виде молекулярных диаграмм, которые необходимы химикам-органикам длясинтеза новых соединений;

- формировать схемы взаимодействий, на. основании которых появляется информация о потенциальных "мишенях" в биологических мембранах, необходимая химикам для установления антидотов к токсичным соединениям, а биоэкологам для изучения экологического последствия.

Разработанный алгоритм впервые был успешно внедрен в область моделирования взаимодействий» компонентов на установках очистки газового конденсата от кислых газов [1].

В соавторстве по данному направлению созданы и зарегистрированы* базы данных:

-№2009620009 «Токсическое воздействие на биологические структуры»; -№2009620103 «Воздействие блокаторов Н-рецепторов и антигистаминных препаратов на биологические структуры»;

-№2009620395 «Молекулярные диаграммы структурированных поверхностей»; -исх. №17-4 от 13.07.2009 «Структурные характеристики биополимеров и двухкомпонентных взаимодействующих систем», которые нашли свое применение при создании новой учебной дисциплины для ВУЗов «Квантовая экологическая химия» [3]. Этот предмет изучают студенты химического факультета Астраханского государственного университета, специальности -510500. В результате, по программе «У.М.Н.И.К. - 2009» получен грант.

Результаты диссертационной работы используются в учебном процессе по дисциплинам: «квантовая механика и квантовая химия», «экология на стыке математики, физики-и химии». Достоверность полученных результатов

Достоверность полученных результатов обеспечивается адекватным использованием численных методов, что подтверждено сравнительным анализом теоретических результатов расчета с экспериментальными и справочными данными. Апробация работы

Результаты исследований' доложены на различных Международных и Всероссийских научных конференциях, среди которых: «Средства и методы обеспечения экологической безопасности» (Астрахань, 2005) «Фундаментальные и1 прикладные проблемы современной химии в исследованиях молодых ученых» (Астрахань, 2006); III школа-семинар «Квантовохимические расчеты: структура и реакционная способность органических и неорганических молекул» (Иваново, 2007); «Фундаментальные и прикладные проблемы современной химии» (Астрахань, 2008); «Инновационные технологии в управлении, образовании, промышленности АСТИНТЕХ-2009» (Астрахань/ 2009); IV школа-семинар «Квантовохимические расчеты: структура и реакционная способность органических и неорганических молекул» (Иваново, 2009). Публикации

По материалам диссертаций опубликовано 5 статей в журналах, рекомендованных ВАК и 5 статей в материалах конференций. Личный вклад соискателя

Результаты расчётов представленные в диссертации получены лично автором. В работах, выполненных с соавторами, соискателю принадлежит: участие в выборе направления исследований, разработка математической модели процессов взаимодействия в макромолекулярных системах, создание алгоритмов моделирования и способов интерпретации полученных результатов для биополимеров. Структура и объём диссертации

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

Выводы третьей главы

Рассмотрена структура базы данных и пользовательский интерфейс программного модуля 2 — Modellnteractions, который позволяет проводить анализ адсорбционных систем и устанавливать активность взаимодействующих атомов при фиксированных значениях длин связей. В результате исследования

153 конфигураций установлен критерий и условия оптимального выбора стабильных положений, на основании чего выстроена методика по определению активных центров на поверхности «неподвижной» фазы и установлению очередности воздействия диоксидом серы на мембрану.

Реакционными центрами в структуре пентапептида являются пептидная, гидроксильная, метильная и карбоксильная группы. В липидах наиболее стабильные продукты взаимодействия можно ожидать для диглицерида, а реакционными центрами являются концевые группировки; для моноглицерида и диглицерида - гидроксогруппа и метильная группа, системы образующие кольцо и карбонильная группа. Важно отметить, что максимально стабильными будут системы, образующие межмолекулярную связь за счет серы и в меньшей степени будут стабильны системы, образующие водородные связи через кислород SO2. Для олигосахаридов практически все формы существования систем будут стабильны, реакционными центрами выступают гидроксильные группы. В результате расчетов установлено, что при длительном воздействии диоксида серы на клеточную мембрану, частичная деструкция и многочисленные изменения в структуре возникают вначале на поверхности олигосахарида, затем подвергается белковая и в меньшей степени липидная составляющая.

В результате, разработанный подход математического моделирования может быть использован при изучении свойств антидотов (противоядие к токсикантам), с целью расчета энергетических параметров процессов дезактивации токсикантов с поверхности клеточной мембраны и блокирование их активных центров. Разница в энергетических характеристиках до и после введения антидота позволит ответить на вопрос применимости выбранного препарата и судить об эффективности его использования в медицине. Необходимость развития данного направления обусловлена применением разработанного подхода в «горячих точках» химии, экологии и медицины.

ЗАКЛЮЧЕНИЕ

Перечислим основные результаты, полученные в диссертации:

1. Впервые разработана математическая модель взаимодействия диоксида серы с макромолекулами клеточной мембраны, которая ориентирована на изучение адсорбционных процессов и установление активных центров в сложных биополимерных структурах, это позволяет предоставить информацию химикам и экологам о возможных механизмах направленного воздействия токсикантов на объекты окружающей среды.

2. На основании математической- модели разработан алгоритм, включающий проведение пошаговой оптимизации геометрии макромолекул до и после взаимодействия. Данный алгоритм вполне может быть использован для •моделирования процессов взаимодействия других токсичных соединений с компонентами клеточной мембраны, что позволит накапливать теоретический материал для проведения сравнительного анализа по степени воздействия.

3. Для макромолекул установлены критерии оптимизации, фиксирование которых позволяет получить геометрические параметры, а на этапе моделирования взаимодействий установлены критерии отбора оптимальных конфигураций, что позволяет выбрать из «-количества образующихся адсорбционных систем, те формы, которые являются наиболее устойчивыми.

4. На основании квантово-химических расчетов получены результаты, которые сопоставимы с экспериментальными и справочными данными и могут быть использованы при формировании молекулярных диаграмм биополимеров клеточной мембраны.

5. Разработана* база данных и автоматизированная система Molecular Model, состоящая из нескольких модулей. Программный модуль BioMolDiagrams позволяет наглядно иллюстрировать биополимеры клеточной мембраны с помощью молекулярных диаграмм. Модуль Modellnteractions предназначен для формирования схем взаимодействий, в результате чего появляется возможность получения информации о потенциальных "мишенях" воздействия токсичного вещества на биологические системы и изучения последствий экологического риска.

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

Благодарности

Автор выражает глубокую благодарность научному руководителю, гениальному профессору Нариману Мирзаевичу Алыкову, и великолепному преподавателю, к.т.н. Лесе Ивановне Жарких за постоянное внимание и помощь в процессе работы над диссертацией.

ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

Статьи в ведущих рецензируемых научных журналах и изданиях, включенных в список ВАК РФ

1. Алыков Н.М. Изучение соединений, содержащих диэтаноламин, ингибиторы кислотной коррозии и продукты их разрушения [Текст] / Н. М. Алыков, М. Н. Котельникова, Н. В. Золотарева, С. Н. Фидурова // Технологии нефти и газа. 2008 г. № 1. С. 23-29.

2. Казанцева Н.В. (Золотарева) Теоретическое обоснование сорбции диоксида серы на структурных элементах клеточных мембран [Текст] // Экологические системы и приборы. - 2007. — №9. — С. 35-37.

3. Казанцева Н.В. (Золотарева) Квантово-химическое кластерное моделирование адсорбции фенола на поверхности алюмосиликатов [Текст], / Н. В. Казанцева, JI. И. Жарких, Н. М. Алыков // Экологические системы и приборы. - 2005г. -№9. - С.27-29.

4. Казанцева Н.В. (Золотарева) Квантово-химическое моделирование хемосорбции диоксида серы на структурных элементах клеточных мембран [Текст] / Н. В. Казанцева, Н. Н. Алыков // Изв. ВУЗов. Химия и хим. технология. 2007. Т. 50. -№12. - С. 132-133.

5. Казанцева Н. В. (Золотарева) Квантово-химическое моделирование хемосорбции диоксида серы на поверхности энергетического субстрата клетки — олигосахарида / Н. В. Казанцева, Н. Н. Алыков // Вестник МГОУ. Сер. Естественные науки. - 2006. — № 2(24). - С. 179.

Статьи в материалах международных и всероссийских конференций

6. Алыков Н. М. Квантовая экологическая химия — новая учебная дисциплина в системе университетского образования [Текст] / Н. М. Алыков, Л. И. Жарких, Н. В. Золотарева // „Инновационные технологии в управлении, образовании, промышленности «АСТИНТЕХ-2009»: материалы международной научной-конференции. 11-14 мая 2009 г. - Иваново, 2009. -С. 125-127.

7. Золотарева Н. В. Моделирование биополимеров клеточной мембраны- и квантово-химичеекий расчет параметров реакционной способности [Текст] // IV школа-семинар: Квантовохимические расчеты: структура и реакционная способность органических и неорганических молекул: 20-22 мая 2009 г. - Иваново, 2009. - С. 72-76.

8. Золотарева Н. В. Квантово-химическое моделирование сорбции диоксида серы. Математический аппарат образования межмолекулярных взаимодействий'[Текст] / Н. В. Золотарева; Е. В. Шитоева, В. И. Юртаева // Фундаментальные и прикладные проблемы современной химии: Материалы II Международной конференции. 15-17 апр. 2008 г. — Астрахань, 2008. - С. 177-181.

9. Золотарева Н. В. Математическое моделирование процессов взаимодействия диоксида серы со структурными элементами клеточной мембраны [Текст] // Фундаментальные и прикладные проблемы современной химии: Материалы II Международной конференции. 15-17 апр. 2008 г. - Астрахань, 2008. - С. 167-177.

10. Казанцева Н. В. (Золотарева) Квантово-химические расчеты энергий взаимодействия диоксида серы со структурными элементами клеточных мембран [Текст] // III школа-семинар: Квантовохимические расчеты: структура и реакционная способность органических и неорганических молекул: 14. 03. 2007 г. - Иваново, 2007. - С. 86-90.

Регистрация интеллектуальной собственности

11 Свидетельство о регистрации базы данных. №2009620009 Токсическое воздействие на биологические структуры [Текст] / JI. И. Жарких, Н. М. Алыков, Н. В. Золотарёва: заявитель и патентообладатель Астраханский государственный университет. - № 2008620341; заяв. 27.10.08; опубл. 11.01.09.

12 Свидетельство о регистрации базы данных. №2009620103 Воздействие блокаторов Н-рецепторов и антигистаминных препаратов на биологические структуры [Текст] / А. И. Юсупов, Н. В. Золотарёва, Н. М. Алыков, JI. И. Жарких: заявитель и патентообладатель Астраханский государственный университет. - № 2008620452; заяв. 31.12.08; опубл. 27.02.09.

13 Свидетельство о регистрации базы данных. №2009620395 Молекулярные диаграммы структурированных поверхностей [Текст] / Н.В. Золотарёва, Н.М. Алыков: заявитель и патентообладатель Астраханский государственный университет. - № 2009620306; заяв. 27.05.09; опубл. 24.07.09.

1. Алыков Н.М. Изучение соединений, содержащих диэтаноламин, ингибиторы кислотной коррозии и продукты их разрушения / Н.М. Алыков, М.Н. Котельникова, Н.В. Золотарева, С.Н. Фидурова // Технологии нефти и газа. 2008 г. № 1. С. 23-29.

2. Аскадский А.А. Компьютерное материаловедение полимеров, т.1. Атомно-молекулярный уровень. / А.А. Аскадский, В.И. Кондращенко. М.: Научный мир, 1999. - 544 с.

3. Алыков Н.М. Метод молекулярных орбиталей в приближении Хюккеля. — Астрахань: ЦНТИ, 1989 42 с.

4. Апостолова Е.С. Квантово-химическое описание реакций / Е.С. Апостолова, А.И. Михайлюк, В.Г. Цирельсон. М.: РХТУ, 1999. - 61 с.

5. Атавин Е.Г. Построение ЗБ-моделей нециклических молекул в естественных переменных / Е.Г. Атавин, В.О. Тихоненко // Вестник Омского университета. 1998. №2. С.35-37.

6. Барановский В.И. Квантовая механика и квантовая химия. — М.: Академия, 2008.-384 с.

7. Бейдер Р. Атомы в молекулах: Квантовая теория. Пер. с англ. — М.: Мир, 2001.-532с.

8. Беркинблит М.Б. Математическое описание поведения мембран. -Биофизика, 1968. 13. №4. С. 712-717.

9. Бурштейн ' К.Я: Квантово-химические расчеты в органической химии и молекулярной спектроскопии / К.Я. Бурштейн, П.П. Шорыгин. — М.: Наука,1989.- 104 с.

10. П.Голованов И.Б. Элементарное введение в квантовую биохимию / И. Б. Голованов, А.К. Пискунов, Н.М. Сергеев М.: Наука, 1969; 236 с.

11. Грибов В.Д. Квантовая химия / В.Д. Грибов, С.П. Муштахова. М.: Гардарики, 1999.-387 с. . ,

12. Губанов В.А. Полуэмпирические методы молекулярных орбиталей в квантовой химии / В.А. Губанов, В.П; Жуков, А.О. Литинский — М.: Наука, 1976.-124 с.

13. Давтян O.K. Квантовая химия. -М.: Высшая школа, 1962. — 784 с.

14. Дероум Э. Современные методы ЯМР для; химических исследований: -М.\ Мир, 1992. 403 с. , .

15. Джонсон К. Численные методы в химии. М.: Мир, 1983. 238 с.

16. Дыоар М. Теория молекулярных орбиталей в органической химии. М.: Мир, 1972.-114 с.

17. Жарких JI. И: Квантово-химическое кластерное моделирование процесса адсорбции сероводорода на поверхности белковой мембраны // Вестник МГОУ. Естественные науки; Вып. Химия и химическая экология. 2006. -№9. - С. 56-59. ^

18. Жидомиров Г.М. Кластерное приближение в квантово-химических исследованиях хемосорбции и поверхностных структур / Г.М. Жидомиров, И.Д. Михейкин // Итоги науки. Строение молекул и хим; связь. 1984. Т. 9. -С.1-21.

19. Заградник Р. Основы квантовой химии / Р. Заградник, Р. Полак. — М.: Мир,1990. -412 с. ,

20. Золотарева Н.В. Квантово-химическое моделирование сорбции диоксида серы. Математический аппарат образования • межмолекулярных взаимодействий Текст. / Н.В. Золотарева, Е.В. Шитоева, В.И. Юргаева //

21. Фундаментальные и прикладные проблемы современной химии: Материалы II Международной конференции. 15-17 апр. 2008 г. Астрахань, 2008. -С.177-181.

22. Иванов В.В. Методы вычислений на ЭВМ: Справочное пособие. Киев: Наук, думка, 1986. - 584 с.

23. Иванов В.В., Слета JI.A. Расчетные методы прогноза биологической активности'органических соединений. Харьков: ХНУ, 2003. - 71 с.

24. Инструментальные методы анализа функциональных групп органических соединений / Под ред. С. Сиггиа. М.: Мир, 1974. - 464 с.

25. Казанцева Н.В. (Золотарева) Квантово-химическое кластерное моделирование адсорбции фенола на поверхности алюмосиликатов / Н. В. Казанцева, JI. И. Жарких, Н. М. Алыков // Экологические системы и приборы. 2005г. -№9. - С. 27-29.

26. Казанцева Н.В. (Золотарева) Квантово-химическое моделирование хемосорбции диоксида серы на поверхности энергетического субстратаклетки олигосахарида / Н. Bi Казанцева, Н. Н. Алыков // Вестник МГОУ. Сер. Естественные науки. - 2006. - № 2(24). - С. 179.

27. Казанцева Н.В. (Золотарева)- Квантово-химическое моделирование хемосорбции диоксида серы на структурных элементах клеточных мембран Текст. / Н. В. Казанцева, Н. Н. Алыков // Изв. ВУЗов. Химия и хим.технология; 2007.TV 50:-№12. С. 132-133:• • ' ' '

28. Казанцева Н.В: (Золотарева) Теоретическое обоснование сорбции диоксида серы на структурных элементах клеточных, мембран // Экологические, системы и приборы. -2007. №9. — С. 35-37. •

29. Калашников Ю.-Ж Основы молекулярной; биологической информатики. М., 2004:-66 с. '

30. Калашников Ю.Я. Ферменты и белки это молекулярные биологические автоматы с программным управлением. - Mi, 2002. — 25 с.

31. Кларк Т. Компьютерная химия: Пер. с англ. -Mi: Мир, 1990. 383 с.

32. Кобзев Б. И. Применение неэмпирических и. полуэмпирических методов в квантово-химических расчетах — Оренбург: ГОУ ОГУ, 2004: — 150 с.

33. Компьютерный анализ структуры, функции и эволюции генетических макромолекул: Проблемы интеллектуализации. Сб. трудов — Новосибирск: 1989.

34. Кон В. Электронная структура вещества волновые функции и функционалы плотности,- УФН, 2002. - Т. 172, № 3. - С.336-348:

35. Конформационный анализ углеводов и их производных / А.Н. Верещагин, В.Е. Катаев, А.А. Бредихин и др.М.: Наука. 1990.- 296 с.40