Ученые разработали вычислительную модель, точно описывающую поведение положительно заряженных частиц — протонов — в молекулах с прочными водородными связями. Авторы установили, что при расчетах важно учитывать не только тип растворителя, в котором находится исследуемая молекула, но и самопроизвольное движение протона и «неопределенность» его расположения в молекуле. Использование этих параметров значительно повысит точность компьютерных моделей для анализа химических соединений, используемых при разработке новых лекарств и материалов. Результаты исследования, поддержанного грантом Российского научного фонда, опубликованы в Journal of Chemical Information and Modeling.
Водородные связи — взаимодействие между атомом водорода и двумя электроотрицательными атомами других элементов, например кислорода, — присутствуют во многих биологических соединениях и играют важную роль в их взаимодействии между собой, с молекулами воды и другими веществами. Поэтому изучение структуры таких связей — важная задача современной химии. Один из способов, которым пользуются ученые, — метод спектроскопии ядерного магнитного резонанса. В рамках этого подхода вещество помещают в магнитное поле и отслеживают, как в таких условиях ведут себя атомы. В результате ученые получают спектры — наборы «пиков», по которым можно исследовать особенности строения химических соединений, так как спектр ядерного магнитного резонанса представляет собой «отпечаток пальцев» исследуемого соединения.
Под воздействием внешнего магнитного поля собственные магнитные поля атомов тоже меняются, и эти изменения зависят от типа и расположения в пространстве атомов-соседей. Это явление называется химическим сдвигом. Их экспериментальное наблюдение либо компьютерное моделирование позволяет установить химическую структуру молекул и прочность химических связей. Обычно для таких измерений к исследуемому веществу добавляют специальный растворитель и таким образом оценивают, как окружающая среда влияет на поведение протонов и показатель химического сдвига. Однако при проведении компьютерного эксперимента по вычислению сдвигов оказывается не всегда достаточным просто добавить некоторое количество молекул растворителя к исследуемой системе, так как такой подход не учитывает особенности динамики частиц — их движение согласно законам механики под влиянием внешних факторов. Точность же в таких расчетах очень важна — она позволяет ученым лучше понимать тонкие процессы, происходящие на атомарном уровне, такие как электронное строение и принцип работы разных типов веществ, например ферментов — «ускорителей» химических реакций.
Исследователи из Санкт-Петербургского государственного университета впервые при моделировании учли и количественно оценили влияние трех основных факторов на химические сдвиги атомов водорода: присутствия молекул растворителя, ядерной динамики протона и квантовой делокализации (отсутствия «закрепленности» протона в конкретной области пространства, возможность его нахождения в разных точках одновременно).
Чтобы оценить влияние этих факторов, химики использовали несколько методов: квантово-химические расчеты, компьютерное моделирование методом молекулярной динамики, а также серию решений уравнения Шредингера. Уравнение Шредингера широко применяется в квантовой химии для вычисления вероятности обнаружения делокализованной частицы в заданной точке пространства. Все вычисления ученые провели на суперкомпьютерах.
Авторы смоделировали химические сдвиги в атоме водорода на примере трех соединений с сильными водородными связями: аниона бифторида (соединение из одного атома водорода и двух атомов фтора), катиона Цунделя (соединение, в котором дополнительный протон «зажат» между двумя молекулами воды) и катиона пиридина-пиридиния, состоящего из атомов водорода, азота и углерода. Эти вещества хорошо изучены химиками, что позволило авторам сравнить полученные результаты с измерениями других ученых и лучше оценить их точность.
Расчеты показали, что включение в модель факторов движения атомов и делокализации протона существенно меняет рассчитанные значения химического сдвига. Например, если при растворении статичного катиона пиридина-пиридиния во фреоне — соединении, содержащем углерод, водород и атом галогена, — показатель сдвига составил 20,8 миллионной доли (единиц, в которых принято оценивать химический сдвиг), то с учетом динамики протона при растворении в том же веществе — уже 22 миллионные единицы, что ближе к данным, полученным ранее в ходе экспериментов. Результаты подчеркивают, что комплексный анализ факторов, влияющих на химический сдвиг, повышает точность его вычисления.
«Результаты нашего исследования важны для специалистов в области моделирования систем с водородными связями, так как дают готовые "рецепты" улучшения согласованности расчетных данных с экспериментальными. Кроме того, наше исследование помогает лучше понять феномен водородной связи, что важно для разработки новых лекарств, материалов и катализаторов, поскольку связи этого типа играют ключевую роль в биологических процессах и современных технологиях. Более точные методы моделирования позволят лучше понять, как работают ферменты, и разработать новые материалы для энергетики и медицины», — рассказывает руководитель проекта, поддержанного грантом РНФ, Елена Тупикина, кандидат физико-математических наук, доцент кафедры физической органической химии Института химии Санкт-Петербургского государственного университета.
Подписывайтесь на InScience.News в социальных сетях: ВКонтакте, Telegram, Одноклассники.
