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

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

М.А. Кременная, М.А. Солдатов, А.П. Чайников, Ю.С. Подковырина, А.Л. Бугаев, К.А. Ломаченко, А.Н. Кравцова

В настоящее время болезнь Альцгеймера считается неизлечимым заболеванием [1]. В России от нее страдают около полутора миллионов человек, и частота поражения болезнью Альцгеймера будет неуклонно возрастать в связи с общей тенденцией к старению населения [2]. Главную морфологическую особенность болезни Альцгеймера составляют нерастворимые отложения бета-амилоида в ткани головного мозга [3]. Однако современные исследования показали, что бета-амилоид является нормальным белком организма и выполняет важные физиологические функции [4,5]. Не так давно было обнаружено, что в нерастворимых отложениях бета-амилоида - амилоидных бляшках - повышена концентрация меди и при этом ионы меди связаны с белком [6]. Структура центра связывания меди отвечает за основные физико-химические свойства белка и представляет несомненный интерес для понимания биохимических процессов. Установление фундаментальных взаимосвязей структуры центра связывания меди в бета-амилоиде с процессами, протекающими во время образования амилоидных бляшек важно для понимания этиологии этого заболевания и в будущем может способствовать разработке эффективного метода борьбы с болезнью Альцгеймера.

В нашей работе представлено исследование локального окружения иона меди Cu1+ в пептиде бета-амилоида при помощи комбинации методов квантово-химических расчетов и XANES спектроскопии (международный термин - X-ray Absorption Near Edge Structure). Спектроскопия XANES -мощный метод исследования локальной атомной и электронной структуры

вещества, чувствительный к ближайшему окружению поглощающего атома и позволяющий исследовать структуру соединений без дальнего порядка в расположении атомов [7]. В частности ранее с помощью рентгеновской спектроскопии поглощения была исследована динамика атомной и электронной структуры материалов водородной энергетики [8] и разработан новый метод для определения параметров 3D наноразмерной атомной и электронной структуры нанокатализаторов [9].

Известно, что при определённых условиях связывание меди с бета-амилоидом более вероятно в линейной конформации [10,11] (рис. 1А). В данной работе были определены расстояния и углы связи для модели центра связывания иона меди Cu1+ в бета-амилоиде с помощью комбинации методов рентгеновской спектроскопии XANES и квантово-химических расчётов. За основу для моделирования была взята простая структура - ион меди, координированный двумя имидазольными кольцами.

Теоретический анализ XANES спектров проводился на основе метода конечных разностей в полном потенциале (программный комплекс FDMNES [12]). Моделирование структуры центра связывания было выполнено при помощи методов молекулярной механики и теории функционала электронной плотности (программный комплекс ADF - Amsterdam Density Functional [13]). Проводилась геометрическая оптимизация локальной структуры путем минимизации полной энергии системы с использованием обменно-корреляционного потенциала в приближении обобщенных градиентов Becke-Perdew и базисным набором DZ (Double Z). Рассматривались модели с нулевым суммарным зарядом.

На основе схемы (рис. 1А) из статьи [10] была построена модель М0 (рис. 1Б). Эта модель является заведомо неточной, так как имидазольные кольца были расположены случайным образом, но с помощью нее можно наглядно показать чувствительность метода к изменению ближайшего окружения иона меди для данной структуры. Далее, методом молекулярной

механики была получена модель ММ (рис. 1В). Модель DFT (рис. 1Г) была рассчитана в рамках теории функционала плотности.

Б

Модель МО:

Си-Ы=1.98 А N-^-N=161.4^

Г

Модель БРТ:

Си-К=1.87 А М-Си-Ы=178.4:'

Рис. 1. - А - схематическое изображения центра связывания иона меди с бета-амилоидом, Б - расчетная модель М0, В - расчетная модель ММ, Г -расчетная модель DFT

Полученные модели (ММ и DFT) имеют минимальные различия, в частности разница в длине связи Си-Ы равна 0.02 А, что находится в пределах чувствительности метода рентгеновской спектроскопии. Так для модели ММ длина связи Си-Ы=1.89 А, а угол: N-^-N=179.8, а для модели DFT - ^-N=1.87 А, угол: N-^-N=178.4°.

Для приведённых выше моделей были рассчитаны спектры рентгеновского поглощения за К-краем меди (I). Расчеты проводились для разных радиусов сетки 3.5 А (в расчет включено 11 атомов), 4.0 А (11 атомов), 4.5 А (15 атомов) 5.0 А (17 атомов), шаг сетки во всех случаях составлял 0.25 А. На рис. 2 представлено сравнение рассчитанных спектров рентгеновского поглощения для трёх моделей с экспериментальным.

В

Модель ММ:

Си-Ы=1.89 А М-Си-Ы=179.8Г

Экспериментальный спектр рентгеновского поглощения для жидкого раствора бета-амилоида с ионами одновалентной меди был измерен при комнатной температуре [10]. Следует отметить, что сходимость теоретических спектров для всех моделей наблюдается при радиусе кластера 5.0 А.

8970 8930 8990 9000 9010 9020 9030 9040 90 50 9060 9070 8970 8980 8990 9000 9010 9020 9030 9040 9050 9060 9070

Энергия, эВ

Энергия, »В

Рис. 2. - Теоретические спектры рентгеновского поглощения: А - для модели М0, Б - для модели ММ, В - для модели ОБТ, Г - сравнение спектров рентгеновского поглощения для всех рассмотренных моделей На рис. 2А представлены спектры рентгеновского поглощения для модели М0. Спектры, рассчитанные для модели М0 не показывают хорошего согласия с экспериментом: интенсивность предпика а на рассчитанном спектре завышена, но энергетическое положение соответствует эксперименту. Стоит отметить, что расчёты проводились без самосогласования электронной структуры, поэтому интенсивность предпика а может быть оценена не совсем верно. Более важным является тот факт, что энергетическое положение основного максимума Ь не соответствует

экспериментальному, что означает несоответствие структурных параметров модели М0 исследуемой структуре.

Далее рассматривая спектры рентгеновского поглощения для моделей ММ (рис. 2Б) и DFT (рис.2В), можно видеть, что спектры дают хорошее согласие с экспериментом при размере атомного кластера 5.0 А. На спектрах, рассчитанных для моделей ММ и DFT, воспроизводится энергетическое положение предпика а, однако эти расчеты так же проводились без самосогласования электронной структуры, и вследствие этого интенсивность особенности а недооценена Энергетическое положение и интенсивность основного максимума Ь, а также поведение спектра рентгеновского поглощения за основным максимумом соответствует экспериментальному спектру, что позволяет говорить о том, что модельные структуры ММ и DFT адекватно описывают исследуемую структуру.

Теоретический анализ спектров рентгеновского поглощения для центра связывания в пептиде бета-амилоида показал, что спектры,

рассчитанные для моделей ММ и DFT (рис. 2Г) дают хорошее согласие с экспериментом, и структуры ММ и DFT подходят для описания центра связывания иона меди ^ (I) в бета-амилоиде.

Расчёты проводились на вычислительных кластерах ЮГИНФО ЮФУ. Исследование выполнено при поддержке Министерства образования и науки Российской Федерации: соглашения 14.A18.21.1940 и 14.A18.21.1070,

стипендия Президента Российской Федерации СП-2373.2013.5.

Литература:

1. Болезнь Альцгеймера //Большая иллюстрированная энциклопедия. В 32 томах. Т1. Б 79 А-АМА - М.: АСТ:Астрель. 2010. -499 ^

2. Белоусов Ю. Б., Зырянов С. К., Белоусов Д. Ю., Бекетов А. С. // Качественная клиническая практика, Спецвыпуск. 2009. № 28.

3. Дэйвис К. Л. Болезнь Альцгеймера: поиск новых способов сохранения функции мозга // Международный Медицинский Журнал. 1999. №9-10.

4. Kamal A., Almenar-Queralt A., LeBlanc J. F., Roberts E. A., Goldstein L.S. B. Kinesin-mediated axonal transport of a membrane compartment containing b-secretase and presenilin-1 requires APP // Nature. 2001. - Vol. 414. P. 603-618.

5. Prangkio P, Yusko EC, Sept D, Yang J, Mayer M (2012) Multivariate Analyses of Amyloid-Beta Oligomer Populations Indicate a Connection between Pore Formation and Cytotoxicity// PLoS ONE. 2012. - Vol.7(10). P. e47261.

6. Kessler H., Bayer T. A., Bach D., Schneider-Axmann T., Supprian T., Herrmann W., Haber M., Multhaup G., Falkai P., Pajonk F.-G. Intake of copper has no effect on cognition in patients with mild Alzheimer’s disease: a pilot phase 2 clinical trial // Journal of Neural Transmission. 2008. Vol. 115. P. 1181-1187.

7. Солдатов А. В. От спектроскопии EXAFS к спектроскопии XANES: новые возможности исследования материи // Соросовский образовательный журнал. 1998. №12.

8. Бугаев А.Л., Гуда А.А., Дмитриев В.П., Ломаченко К.А., Панкин И.А., Смоленцев Н.Ю., Солдатов М.А., Солдатов А.В. Динамика наноразмерной атомной и электронной структуры материалов водородной энергетики при реалистичных технологических условиях [Электронный ресурс] // «Инженерный вестник Дона», 2012, №4. - Режим доступа: http://www.ivdon.ru/magazine/archive/n4p1y2012/1269 (доступ свободный) -Загл. с экрана. - Яз. рус.

9. Смирнова Ю.О., Положенцев О.Е., Леонтьева Д.В., Чайников А.П.,

Сучкова С.А., Гуда А.А., Ломаченко К.А., Смоленцев Н.Ю., Подковырина Ю.С., Солдатов М.А., Кравцова А.Н., Солдатов А.В. Разработка нового комплексного метода определения параметров 3D наноразмерной атомной и электронной структуры материалов на основе методик XAFS, XRD и Raman [Электронный ресурс] // «Инженерный вестник Дона», 2012, №4. - Режим доступа: http://www.ivdon.ru/magazine/archive/n4p1y2012/1268 (доступ

свободный) - Загл. с экрана. - Яз. рус.

10. Hureau C., Balland V., Coppel Y., Solari P. L., Fonda E., Faller P. Importance of dynamical processes in the coordination chemistry and redox

conversion of copper amyloid-b complexes // JBIC Journal of Biological Inorganic Chemistry. 2009. - Vol. 14. P. 995.

11. Shearer J., Szalai V. A. The amyloid-P peptide of Alzheimer's disease binds CuI in a linear bis-His coordination environment: Insight into a possible neuroprotective mechanism for the amyloid-P peptide // J Am Chem Soc. 2008. -Vol. 130(52). P. 17826-17835.

12. Joly Y. X-ray absorption near edge structure calculations beyond the muffin-tin approximation // Phys. Rev. B. 2001. - Vol. 63. 125120.

13. te Velde G., Bickelhaupt F.M., van Gisbergen S.J.A., Fonseca Guerra C., Baerends E.J., Snijders J.G., Ziegler T. // Chemistry with ADF, Journal of Computational Chemistry. 2001. - Vol. 22. P. 931.