КОМПЬЮТЕРНОЕ ПРЕДСКАЗАНИЕ ПАТОГЕННЫХ СВОЙСТВ МУТАНТОВ
БЕЛКА SOD1, АССОЦИИРОВАННЫХ С БОКОВЫМ АМИОТРОФИЧЕСКИМ
СКЛЕРОЗОМ
Н.А. Алемасов, мл. науч. сотр. Н.В. Иванисенко, мл. науч. сотр. В.А. Иванисенко, вед. науч. сотр.
Федеральный исследовательский центр Институт цитологии и генетики Сибирского отделения Российской академии наук (Россия, г. Новосибирск)
DOI: 10.24411/2500-1000-2018-10053
Работа выполнена при поддержке гранта РНФ №14-44-00011 «Программируемая клеточная гибель, индуцируемая через рецепторы смерти: идентификация молекулярных механизмов инициации апоптоза с помощью молекулярного моделирования»
Аннотация. Внутриклеточная агрегация белков является одной из причин многих нейродегенеративных заболеваний. В частности, мутации в гене SOD1 связаны с наследственной формой бокового амиотрофического склероза (БАС). Ранее нами была предложена регрессионная модель, предсказывающая дожитие пациентов с БАС, которые являются носителями мутаций в белке SOD1. Эта модель была построена на основе анализа стабильности водородных связей, возникающих в мутантах SOD1 при моделировании методом молекулярной динамики (МД). В настоящей работе данная регрессионная модель была улучшена с помощью учёта новой характеристики водородных связей, отражающей показатель консервативности водородной связи в пространстве конформационных состояний белка. Расчёт конформационной консервативности водородных связей, осуществлённый с помощью эластичных сетевых (ЭС) моделей позволил увеличить коэффициент корреляции между предсказанными с помощью новой модели и наблюдаемыми значениями дожития пациентов с мутациями в SOD1 ^ = 0.92).
Ключевые слова: боковой амиотрофический склероз, SOD1, водородные связи, регрессионные модели, молекулярная динамика, эластичные сетевые модели.
Боковой амиотрофический склероз (БАС) - неизлечимое нейродегенеративное заболевание. Выделяется две его формы: наследственная и спорадическая [1]. Известно, что второй по распространённости и наиболее изученной причиной наследственной формы БАС являются мутации в гене SOD1, кодирующем фермент супер-оксиддисмутазу-1 [2]. Одной из гипотез относительно механизма возникновения заболевания на молекулярном уровне является агрегация белка SOD1, вызванная его неверной укладкой в результате мутаций [2].
Известно, что водородные связи играют важную роль, в частности, в поддержании
структуры белков [3]. В одном из т silico исследований структуры и динамики SOD1 уже анализировались водородные связи и сеть, которую они образуют в этом белке [4]. Ранее нами с помощью метода молекулярной динамики (МД) было исследовано 39 мутаций в белке SOD1 [5]. Были обнаружены водородные связи, доля времени существования которых в траектории МД (или стабильность) достоверно коррелирует с дожитием пациентов с БАС. На основе найденных водородных связей была построена регрессионная модель, позволяющая предсказать дожитие пациентов с мутациями в SOD1.
Настоящее исследование направлено на расширение ранее созданной модели с учётом крупномасштабной динамики белка SOD1. Для этого проведено моделирование динамики структуры белка с использованием эластичных сетевых моделей (ЭС) [6]. На основе результатов анализа структуры и динамики мутантов 8ОБ1 предложена комплексная регрессионная модель, учитывающая как конформации, полученные с помощью МД, так и кон-формации, полученных в данной работе с помощью ЭС. Для построения регрессионной модели проводился расчёт конфор-мационной консервативности водородных связей, как отношение числа конформа-ций, в которых связь присутствовала, к общему числу полученных в ЭС конфор-маций. Анализ конформационной консервативности водородных связей, позволил улучшить точность расширенной модели по сравнению с ранее созданной моделью. Коэффициент корреляции между предсказанными и наблюдаемыми значениями дожития пациентов, носителей мутаций в SOD1, составил R = 0.92.
Материалы и методы
В работе исследована структура белка SOD1 дикого типа и 35 его мутантных форм, взятые из базы данных ALSOD (http://alsod.iop.kcl.ac.uk/). Получение пространственных структур мутантных белков SOD1 производилось с помощью моделирования мутаций в структуре PDBID 2V0A с использованием FoldX [7]. Проверка и исправление полученных структур мутантов, включая удаление клэшей и исправление неоптимальных конформаций также проводилась с использованием БоЫХ.
Построение ЭС для белка дикого типа и его мутантов проводилось с помощью системы ЕШето с параметрами по-
умолчанию [8]. В рамках каждой такой модели был получен ансамбль из 11 кон-формаций для каждой из 5 нетривиальных мод колебаний с наивысшим показателем коллективности. Таким образом, для белка дикого типа и каждого из изучаемых мутантов были рассчитаны 55 конформаций, представляющих крупномасштабные флуктуации структуры этих молекул. С помощью пакета программ AMBER [9] структура, соответствующая каждой кон-формации впоследствии минимизировалась с целью избежать физически недостижимых торсионных углов, длин связей и других параметров SOD1.
Водородные связи в каждой из 55 кон-формаций белка находились с помощью утилиты cpptraj из пакета программ AMBER. Статистический анализ стабильности водородных связей проводился с помощью библиотек numpy, scipy и sklearn для языка Python.
Результаты и обсуждение
Следует заметить, что характер движений атомов молекул в моделях ЭС и МД различается. Различия состоят в том, что в моделировании ЭС все конформации молекулы отличаются и соответствуют позициям атомов вдоль выбранной моды колебаний в некоторый промежуток времени. В МД атомы движутся неподалёку от некоторого локального минимума и могут принимать одну и ту же конформацию не один раз. По этой причине доля существования водородной связи среди конформаций белка в ЭС названа конформационной
консервативностью, в то время как в нашей предыдущей работе [5] доля конформаций белка в траектории МД, в которых связь была сформирована называлась стабильностью.
рормационная консервативность,
Рис. 1. Корреляция между функцией конформационной консервативности водородных связей, найденных в рамках моделирования эластичных сетей для мутантов SOD1 и дожитием пациентов с БАС - носителей соответствующих мутациями в белке.
Показатель, характеризующий
патогенность водородных связей в МД рассчитывался с помощью суммарной стабильности От,К, описанной нами ранее в работе [5]. Для оценки патогенности водородных связей в ЭС применялась суммарная конформационная консервативность Си,^ которая рассчитывается на основе
конформационной консервативности
отдельных водородных связей аналогично От,К. Регрессионная модель включала два независимых фактора - От,К и Си,К
Коэффициент корреляции между функцией конформационной
консервативности Си^ и дожитием пациентов с мутациями в SOD1 составил R = 0.84 (р = 2.8*10-9) (см. рис. 1).
20
h
Ф 15
X 10
о -
0) о Е
0)
та § 5
ч
щ
та
% Мутанты SO Dl
---y=0.4xi+0.7x2-0.6 (R= -0.92, p = 5.77e-14) • у * • ✓ • / ✓ у ✓
У •
% ч / • V
•A
•
0.0 2,5 5.0 7.5 10.0 12.5 15.0 17.5 20.0 Предсказанное дожитие пациентов, лет
Рис. 2. Корреляция между предсказанным временем жизни пациентов с БАС с мутациями в SOD1, полученным с помощью двухфакторной линейной регрессии на основе конформационной консервативности патологичных водородных связей в ЭС и стабильностей патологичных водородных связей в МД
Конформационная консервативность и Заключение
стабильность водородных связей, Учёт показателя консервативности
соответственно, в ЭС и МД являются водородных связей в пространстве различными по своей природе. Поэтому на конформационных состояний белка SOD1 основе данных двух независимых позволил улучшить созданную ранее характеристик была построена модель, которая была построена на множественная регрессия. Коэффициент использовании информации о
корреляции для данной регрессии оказался стабильности водородных связей. выше (R = 0.92, p = 5.8*10-14), чем для Коэффициент корреляции между регрессий на основе каждого из предсказанными с помощью новой модели параметров в отдельности (см. рис. 2). и наблюдаемыми значениями дожития Таким образом, комбинация результатов пациентов с мутациями в SOD1 составил R моделирования ЭС и МД позволила = 0.92, что позволяет объяснить более 84% получить более точную регрессионную (R2 = 0,846) дисперсии в дожитии модель для предсказания дожития пациентов, носителей рассматриваемых пациентов с мутациями в SOD1, мутаций в SOD1. страдающих БАС.
Библиографический список
1. KiernanM.C. et al. Amyotrophic lateral sclerosis // Lancet. 2011. Vol. 377. №9769. P. 942955.
2. Tokuda E., Furukawa Y. Copper Homeostasis as a Therapeutic Target in Amyotrophic Lateral Sclerosis with SOD1 Mutations // Int. J. Mol. Sci. Multidisciplinary Digital Publishing Institute, 2016. Vol. 17, № 5. P. 636.
3. Fernández A., Scheraga H.A. Insufficiently dehydrated hydrogen bonds as determinants of protein interactions. // Proc. Natl. Acad. Sci. U. S. A. 2003. Vol. 100, № 1. P. 113-118.
4. Khare S.D., Dokholyan N. V. Common dynamical signatures of familial amyotrophic lateral sclerosis-associated structurally diverse Cu, Zn superoxide dismutase mutants. // Proc. Natl. Acad. Sci. U. S. A. 2006. Vol. 103, № 9. P. 3147-3152.
5. Alemasov N.A. et al. Dynamic properties of SOD1 mutants can predict survival time of patients carrying familial amyotrophic lateral sclerosis // J. Biomol. Struct. Dyn. Taylor & Francis, 2017. Vol. 35, № 3. P. 645-656.
6. Tirion M. Large Amplitude Elastic Motions in Proteins from a Single-Parameter, Atomic Analysis. // Phys. Rev. Lett. American Physical Society, 1996. Vol. 77, № 9. P. 1905-1908.
7. Guerois R., Nielsen J.E., Serrano L. Predicting Changes in the Stability of Proteins and Protein Complexes: A Study of More Than 1000 Mutations // J. Mol. Biol. 2002. Vol. 320, № 2. P. 369-387.
8. Suhre K., Sanejouand Y.-H. ElNemo: a normal mode web server for protein movement analysis and the generation of templates for molecular replacement. // Nucleic Acids Res. Oxford University Press, 2004. Vol. 32, № Web Server issue. P. W610-4.
9. Salomon-Ferrer R., Case D.A., Walker R.C. An overview of the Amber biomolecular simulation package // Wiley Interdiscip. Rev. Comput. Mol. Sci. 2013. Vol. 3, № 2. P. 198-210.
COMPUTATIONAL PREDICTION OF PATHOGENIC PROPERTIES OF SOD1 MUTANTS WHICH ARE ASSOCIATED WITH AMYOTROPHIC LATERAL
SCLEROSIS
N.A. Alemasov, junior researcher N.V. Ivanisenko, junior researcher V.A. Ivanisenko, senior research
Federal research center institute of cytology and genetics, Siberian branch of the Russian academy of sciences (Russia, Novosibirsk)
Abstract. Intracellular protein aggregation is one of the causes of many neurodegenerative diseases. In particular, mutations in the SOD1 gene are associated with the familial form of amyotrophic lateral sclerosis (ALS). We have previously proposed a regression model that predicts the survival of patients with ALS who are carriers of mutations in the SOD1 protein. This model was built on the basis of the analysis of the stability of hydrogen bonds forming in SOD1 mutants when simulated by the molecular dynamics (MD). In the present work, this regression model has been improved by taking into account the new characteristic of hydrogen bonds, reflecting the conservativeness index of the hydrogen bond in the space of conformational states of the protein. The calculation of conformational conservatism of hydrogen bonds, carried out using elastic network (EN) models, made it possible to increase the correlation coefficient between the predicted using the new model and the observed survival values of patients with mutations in SOD1 (R = 0.92).
Keywords: amyotrophic lateral sclerosis, SOD1, hydrogen bonds, regression models, molecular dynamics, elastic network models.