ВЗАИМОДЕЙСТВИЕ ДИНУКЛЕОТИДОВ С ВОДНЫМ ОКРУЖЕНИЕМ: МОДЕЛИРОВАНИЕ МЕТОДОМ МОНТЕ КАРЛО

Шестопалова А. В.
Институт радиофизики и электроники НАН Украины, Харьков
Для изучения взаимодействия динуклеотидов с водным окружением проведено компьютерное моделирование систем, содержащих динукелеотиды (ApA, ApT, TpT, TpU, UpT, UpU, GpG, GpC, CpC и др.) и водные кластеры различных размеров: 20, 40, 80, 100, 200 и 400 молекул воды. Такие модельные системы воспроизводят взаимодействие компонентов нуклеиновых кислот с водой при разных уровнях относительной влажности. Расчеты выполнены методом Монте Карло в каноническом (N,V,T) ансамбле. Молекулы воды в процессе моделирования смещались согласно алгоритму Метрополиса. В расчетах использовалось кластерное приближение; системы помещались в центр сферы с жесткими отражательными стенками так, что центр масс "раствора" совпадал с центром сферы. Радиусы сфер выбирались в зависимости от размеров водного кластера от 16А до 26А. Для расчетов энергии межмолекулярного взаимодействия вода-динуклеотид и вода-вода использовались полуэмпирические атом-атомные потенциальные функции с параметрами Полтева; торсионные углы динуклеотидов изменялись согласно алгоритму Метрополиса Стартовая геометрия динуклеотидов соответствовала А- или В- форме с оптимизацией в вакууме по алгоритму Метрополиса. Энергетически выгодная в вакууме конфигурация динуклеотида помещалась в водный кластер и выполнялось моделирование гидратации динуклеотидов. Получены следующие энергетические характеристики изученных систем: средняя полная энергия взаимодействия в системе, включающая в себя взаимодействие вода-вода, вода - динуклеотид и конформационную энергию динуклеотида; средняя энергия взаимодействия вода-вода; средняя энергия взаимодействия динуклеотид - вода; удельная энергия взаимодействия вода-вода, приходящаяся на одну молекулу воды; удельная энергия взаимодействия динуклеотид - вода, приходящаяся на одну молекулу воды; избыточная энергия гидратации DE и вычислено t = t0. exp (-DE/RT), (t0 = 0,9. 10-11 сек,) - время дипольной релаксации молекулы воды в связанном состоянии. Значения избыточной энергии гидратации и времен дипольной релаксации молекул воды в связанном состоянии коррелируют с экспериментальными результатами, полученными калориметрически для полинуклеотидов в процессе их дегидратации. Проведен анализ энергетических и структурных характеристик изученных систем. Показано значительное отличие во взаимодействии динуклеотидов АА, АТ, и GG, GС с водным окружением для всех размеров водных кластеров как по уровню гидратации, так и по влиянию на структуру растворителя. Анализ структурных характеристик водных кластеров (число водородных связей, образуемых молекулами воды в кластере, и число водородных связей, приходящихся на одну молекулу воды) и мгновенных конфигураций систем позволяет определить сетки из молекул воды вблизи гидратноактивных центров динуклеотидов. Проанализировано изменение конформационных параметров динуклеотидов в зависимости от размеров водных кластеров.