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

Development of molecular dynamics and quantum chemical approaches to predict and optimize linear and nonlinear photophysical properties of fluorescent biomarkers

Источник финансирования НИР

грант РНФ

Этапы НИР

# Сроки Название
1 28 июля 2021 г.-30 июня 2022 г. Развитие молекулярно-динамических и квантово-химических подходов для предсказания и оптимизации линейных и нелинейных фотофизических свойств флуоресцентных биомаркеров
Результаты этапа: Разработан комбинированный молекулярно-динамический и квантово-химический подход для моделирования формы и положения полос в спектрах однофотонного поглощения хромофоров в реальном окружении с учетом однородного и неоднородного уширения. Для работы с большим числом конфигураций при усреднении вдоль молекулярно-динамической траектории реализована серия автоматизированных протоколов, написанных на языке Python. Разработанный теоретический подход объединяет методы молекулярной динамики, вибронные модели, методы квантовой химии высокого уровня точности, такие как расширенная многоконфигурационная квазивырожденная теория возмущений второго порядка XMCQDPT2 и комбинированные методы квантовой и молекулярной механики КМ/ММ на ее основе. Предложенная методика моделирования спектров поглощения комплексообразующих красителей в водном растворе в присутствии различных ионов металлов включает в себя следующие этапы: (1) молекулярно-динамическое моделирование в NPT ансамбле при 300 К с целью выбора конфигураций для последующих этапов моделирования и усреднения; (2) оптимизация геометрических параметров выбранных структур с помощью метода КМ/ММ на основе теории функционала электронной плотности; (3) расчёт энергий вертикального возбуждения, вероятности переходов и параметров вибронной модели с помощью метода XMCQDPT2 с учетом влияния окружения в рамках метода потенциала эффективных фрагментов EFP; (4) моделирование вибронной структуры спектра при гомогенном уширении на основе расчета перекрывания колебательных волновых функций различных электронных состояний в модели линейной связи в приближении Кондона; (5) расчёт интегрального электронно-колебательного спектра путем усреднения вдоль молекулярно-динамической траектории. С использованием разработанного подхода установлены основные факторы, влияющие на форму и положение полос в спектрах поглощения кальциевого индикатора Fura-2 в водном растворе. Установлено, что основной вклад в спектр поглощения свободного красителя вносит негомогенное уширение - наличие большого числа конформаций в основном электронном состоянии, связанных с внутримолекулярным вращением и влияющих на сопряжение в пи-электронной системе молекулы красителя. Это позволяет объяснить большую ширину спектра поглощения красителя Fura-2 в свободном виде. Хелатный комплекс красителя с ионами Ca2+ конформационно значительно менее подвижен, при этом сопряжение в системе сильно нарушается за счет связывания карбоксилатных групп с ионами металла, что приводит к сдвигу максимума поглощения в более коротковолновую область и значительному уменьшению спектральной ширины. Анализ характера низколежащего возбужденного состояния красителя Fura-2 свидетельствует о наличии внутримолекулярного переноса заряда. Перенос заряда при возбуждении приводит к увеличению дипольного момента в возбужденном состоянии и объясняет сильную зависимость фотофизических свойств красителя от окружения, в первую очередь, от полярности среды. Сделан вывод о том, что спектр поглощения красителя Fura-2 в свободном виде наиболее подвержен влиянию среды. По результатам проекта подготовлена и принята в печать 1 статья, а также представлено 3 доклада на международных и всероссийских конференциях.
2 1 июля 2022 г.-30 июня 2023 г. Развитие молекулярно-динамических и квантово-химических подходов для предсказания и оптимизации линейных и нелинейных фотофизических свойств флуоресцентных биомаркеров
Результаты этапа: 1) Разработан метод расчёта спектров двухфотонного поглощения конформационно нежестких хромофоров в растворе и белковом окружении с учётом однородного и неоднородного уширения. 2) С помощью разработанного метода проведено моделирование спектров поглощения кальциевого индикатора Fura-2 в водном растворе. Установлено, что основной вклад в спектр поглощения свободного красителя как при однофотонном, так и двухфотонном возбуждении вносит негомогенное уширение – наличие большого числа конформаций в основном электронном состоянии, связанных с внутримолекулярным вращением и приводящего к значительному спектральному уширению, что подтверждается экспериментальными данными. 3) Анализ характера низколежащего электронно-возбужденного состояния в молекуле Fura-2 свидетельствует о наличии внутримолекулярного переноса заряда, что позволяет объяснить сильную зависимость фотофизических свойств красителя от среды, в первую очередь, от её полярности. 4) Показано, что протонирование нехелатообразующей карбоксильной группы свободного красителя не влияет на спектр его поглощения. 5) Для флуоресцентных белков сечение двухфотонного поглощения, как правило, составляет от 1 до 100 ГМ. Ратиометрический краситель Fura-2, согласно экспериментальным данным, имеет сечения двухфотонного поглощения порядка 10-15 ГМ, что является не самым высоким показателем. В связи с сильной зависимостью сечения двухфотонного поглощения от разницы средних дипольных моментов в начальном и конечном состоянии, а также дипольного момента перехода, сечение двухфотонного поглощения можно увеличивать, изменяя локальное окружение красителя, а также вводя дополнительные функциональные группы для увеличения переноса заряда в системе. 6) Построена полная атомистическая модель модифицированного белка EGFP T203I, оптимизирована ее геометрия, рассчитано сечение двухфотонного поглощения и проведен анализ факторов, влияющих на него. С помощью методов молекулярной динамики проведено сэмплирование конфигурационного пространства активного центра мутированного белка. Показано, что рассчитанные значения сечения двухфотонного поглощения варьируются в широких пределах и сильно зависят от структуры активного центра и сетки водородных связей около хромофорной группы. 7) Показано, что увеличение сечения двухфотонного поглощения EGFP может быть достигнуто путем модификации белкового окружения хромофора при введении точечных мутаций, влияющих на электростатическое поле белка вблизи хромофора. Увеличению сечения отвечает усиление напряженности электрического поля вдоль направления от фенольного кольца хромофора к имидазолоновому (противоположному направлению вектора изменения постоянного дипольного момента при возбуждении), которое может быть достигнуто, например, изменением заряда окружения фенольного кольца хромофора на более положительный, а имидазолонового – на более отрицательный, что подтверждается экспериментальными данными.

Прикрепленные к НИР результаты

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