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

Краткое содержание

Спектроскопия обратного рассеяния нейтронов обеспечивает неразрушающий и свободный от меток доступ к ps-ns динамике белков и их гидратационной воде. Рабочий процесс представлен двумя исследованиями амилоидных белков: динамики лизоцима во время агрегации с временным разрешением и динамики гидратационной воды тау при образовании волокон.

Аннотация

Рассеяние нейтронов дает возможность исследовать динамику в образцах для широкого диапазона энергий неразрушающим способом и без маркировки, кроме дейтерия. В частности, спектроскопия обратного рассеяния нейтронов регистрирует сигналы рассеяния при нескольких углах рассеяния одновременно и хорошо подходит для изучения динамики биологических систем на шкале времени ps-ns. Используя_D2O- и, возможно, дейтерированные буферные компоненты, метод позволяет контролировать как диффузию центра масс, так и движения основной цепи и боковой цепи (внутреннюю динамику) белков в жидком состоянии.

Кроме того, динамика гидратационной воды может быть изучена с использованием порошков пердейтерированных белков, гидратированных с H₂O. В этой статье представлен рабочий процесс, используемый на приборе IN16B в Институте Лауэ-Ланжевена (ILL) для исследования динамики белковой и гидратационной воды. Объясняется приготовление образцов раствора и образцов гидратированного протеинового порошка с использованием парообмена. Описана процедура анализа данных как белковой, так и гидратационной динамики воды для различных типов наборов данных (квазиупругие спектры или сканирование с фиксированным окном), которые могут быть получены на спектрометре обратного рассеяния нейтронов.

Метод проиллюстрирован двумя исследованиями с участием амилоидных белков. Показано, что агрегация лизоцима в сферические агрегаты размером с мкм, обозначенные частицами, происходит в одностадийном процессе в пространственно-временном диапазоне, исследованном на IN16B, в то время как внутренняя динамика остается неизменной. Далее изучали динамику гидратации воды тау на гидратированных порошках пердейтерированного белка. Показано, что поступательные движения воды активируются при образовании амилоидных волокон. Наконец, обсуждаются критические шаги в протоколе относительно того, как рассеяние нейтронов позиционируется по отношению к изучению динамики по отношению к другим экспериментальным биофизическим методам.

Введение

Нейтрон представляет собой массивную частицу без заряда, которая на протяжении многих лет успешно использовалась для зондирования образцов в различных областях от фундаментальной физики до биологии¹. Для биологических применений широко используются малоугловое рассеяние нейтронов, неупругое рассеяние нейтронов, а также нейтронная кристаллография и рефлектометрия ^2,3,4. Неупругое рассеяние нейтронов обеспечивает усредненное по ансамблю измерение динамики, не требующее специфической маркировки как таковой, и качество сигнала, которое не зависит от размера или белка⁵. Измерение может быть выполнено с использованием очень сложной среды для исследуемого белка, которая имитирует внутриклеточную среду, такую как дейтерированный бактериальный лизат или даже in vivo ^3,6,7. Для изучения динамики могут быть использованы различные экспериментальные установки, а именно: i) доступ к динамике суб-ps-ps, дающий время пролета, ii) доступ к динамике ps-ns, дающий обратное рассеяние, и iii) спин-эхо-дающий доступ к динамике от нс до сотен нс. Обратное рассеяние нейтронов использует закон Брэгга 2d sinθ = nλ, где d — расстояние между плоскостями в кристалле, θ — угол рассеяния, n — порядок рассеяния и λ — длина волны. Использование кристаллов для обратного рассеяния по направлению к детекторам позволяет достичь высокого разрешения по энергии, обычно ~ 0,8 мкВ. Для измерения энергообмена используется либо доплеровский привод, несущий кристалл в обратном рассеянии, для определения и настройки входящей длины волны ^{нейтрона 8,9,10} (рис. 1), либо времяпролетная установка может быть использована за счет уменьшения энергетического разрешения ¹¹.

Рисунок 1: Эскиз спектрометра обратного рассеяния нейтронов с доплеровским приводом. Входящий луч попадает на прерыватель⁴² преобразования фазового пространства (PST), который увеличивает поток в положении образца. Затем он обратно рассеивается по направлению к образцу доплеровским приводом, который выбирает энергию E₁ (голубая стрелка). Затем нейтроны рассеиваются образцом (с различными энергиями, представленными цветом стрелок), и анализаторы, изготовленные из кристаллов Si 111, будут рассеивать только нейтроны с определенной энергией E₀ (здесь красные стрелки). Следовательно, передача импульса q получается из обнаруженного положения нейтрона на матрице детекторов, а передача энергии получается из разности E₁ - E₀. Время пролета, ожидаемое для нейтронного импульса, создаваемого PST, используется для отбрасывания сигнала от нейтронов, рассеянных непосредственно к детекторным трубкам. Аббревиатура: PST = преобразование фазового пространства. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Для спектроскопии обратного рассеяния основной вклад в сигнал от образцов, богатых протонами водорода, таких как белки, приходится на некогерентное рассеяние, для которого интенсивность рассеяния S_inc(q, ω) показана уравнением (1)¹²

(1)

Где σ_inc — некогерентное поперечное сечение рассматриваемого элемента, k' — норма рассеянного волнового вектора, k — норма входящего волнового вектора, q (= k - k') — передача импульса, r j(t) — вектор положения атома _j в момент времени t, а ω — частота, соответствующая переносу энергии между входящим нейтроном и системой. Угловые скобки обозначают среднее значение ансамбля. Следовательно, некогерентное рассеяние исследует усредненную по ансамблю самокорреляцию положений атомов во времени и дает самодинамику, усредненную по всем атомам в системе и различным временным началам (среднее значение ансамбля). Функция рассеяния является преобразованием Фурье во времени промежуточной функции рассеяния I(q, t), которое можно рассматривать как преобразование Фурье в пространстве корреляционной функции Ван Хова, показанное уравнением (2):

(2)

Где ρ(r,t) — плотность вероятности нахождения атома в положении r и времени t ¹³.

Для фиковского диффузионного процесса функция самодиффузии возникает (см. Eq (3)) после двойного преобразования Фурье в функции рассеяния, состоящей из лоренциана ширины линии, заданной γ = Dq².

(3)

Были разработаны и признаны полезными более сложные модели, такие как модель скачкообразной диффузии Сингви и Шёландера для внутренней динамики белка^{ps-ns 14} или модель вращения Sears для гидратационной воды^15,16,17.

На приборе обратного рассеяния нейтронов (NBS) IN16B ^8,9 в ILL, Гренобль, Франция (дополнительный рисунок S1), установка, обычно используемая с белками, состоит из кристаллов Si 111 для анализаторов с доплеровским приводом для настройки входящей длины волны (дополнительный рисунок S2A), тем самым предоставляя доступ к диапазону передачи импульса ~ ^0,2 Å-1 < q < ~ 2 ^Å-1 и диапазону передачи энергии ^-30 мкэВ < < 30 мкэВ, что соответствует временным шкалам в диапазоне от нескольких пс до нескольких нс и расстояний в несколько Å. Кроме того, IN16B предлагает возможность выполнять упругое и неупругое сканирование с фиксированным окном (E/IFWS)¹⁰, которое включает сбор данных при фиксированной передаче энергии. Поскольку поток ограничен при работе с нейтронами, E/IFWS позволяет максимизировать поток для одной передачи энергии, тем самым сокращая время сбора данных, необходимое для получения удовлетворительного отношения сигнал/шум. Более поздним вариантом является режим¹¹ спектрометра обратного рассеяния и времяпролетного спектрометра (BATS), который позволяет измерять широкий диапазон переносов энергии (например, -150 мкВ < < 150 мкВ) с более высоким потоком, чем с доплеровским приводом, но за счет более низкого энергетического разрешения (дополнительный рисунок S2B).

Важным свойством рассеяния нейтронов является то, что некогерентное поперечное сечение σ_{включительно} имеет в 40 раз более высокое значение для водорода, чем для дейтерия, и пренебрежимо мало для других элементов, обычно встречающихся в биологических образцах. Следовательно, динамика белков в жидкой среде может быть изучена с использованием дейтерированного буфера, а состояние порошка позволяет изучать либо внутреннюю динамику белка с гидрогенизированным протеиновым порошком, гидратированным с_D2O, либо исследование гидратационной воды для пердейтерированного протеинового порошка, гидратированного с H₂O. В жидком состоянии обратное рассеяние нейтронов обычно позволяет одновременно получить доступ к самодиффузии белков в центре масс (диффузия фиковского типа) и их внутренней динамике. Последние представляют собой движения основной цепи и боковой цепи, обычно описываемые так называемой скачкообразной диффузионной моделью или другими ^3,18. В гидрогенизированных протеиновых порошках диффузия белка отсутствует и необходимо моделировать только внутреннюю динамику. Для гидратационной воды вклад поступательных и вращательных движений молекул воды представляет различную зависимость от передачи импульса q, что позволяет различать их в процессе анализа данных¹⁷.

Эта работа иллюстрирует метод обратного рассеяния нейтронов с изучением белков, которые, как было установлено, способны разворачиваться, агрегироваться в каноническую форму, состоящую из стопок β-нитей - так называемого перекрестного β паттерна^19,20 - и образовывать удлиненные волокна. Это так называемая амилоидная агрегация, которая широко изучается из-за ее центральной роли в нейродегенеративных расстройствах, таких как болезни Альцгеймера или Паркинсона^21,22. Изучение амилоидных белков также мотивировано функциональной ролью, которую они могут играть ^23,24, или их высоким потенциалом для разработки новых биоматериалов²⁵. Физико-химические детерминанты агрегации амилоида остаются неясными, и нет общей теории агрегации амилоида, несмотря на огромный прогресс за последние^{годы 21,26}.

Агрегация амилоида подразумевает изменения структуры и стабильности белка с течением времени, изучение которых, естественно, подразумевает динамику, связанную со стабильностью конформации белка, функцией белка и энергетическим ландшафтом^{белка 27}. Динамика напрямую связана со стабильностью определенного состояния через энтропийный вклад для самых быстрых движений²⁸, и функция белка может поддерживаться движениями в различных временных масштабах от суб-пс для светочувствительных белков²⁹ до мс для доменных движений, что может быть облегчено пикосекундно-наносекундной динамикой³⁰.

Будут представлены два примера использования спектроскопии обратного рассеяния нейтронов для изучения амилоидных белков: один в жидком состоянии для изучения динамики белка и один в гидратированном порошковом состоянии для изучения динамики гидратационной воды. Первый пример касается агрегации лизоцима в сферы размером мкм (называемые частицами), за которой следуют в режиме^{реального времени 5}, а второй - сравнение динамики воды в нативном и агрегированном состояниях человеческого белка тау-31.

Лизоцим является ферментом, участвующим в иммунной защите и состоит из 129 аминокислотных остатков. Лизоцим может образовывать твердые частицы в дейтерированном буфере при pD 10,5 и температуре 90 °C. С помощью рассеяния нейтронов мы показали, что временная эволюция коэффициента диффузии лизоцима в центре масс следует за одной экспоненциальной кинетикой флуоресценции тиофлавина Т (флуоресцентный зонд, используемый для мониторинга образования амилоидных перекрестных β паттернов³²), что указывает на то, что образование сверхструктур частиц и перекрестных β происходит за один шаг с одинаковой скоростью. Более того, внутренняя динамика оставалась постоянной на протяжении всего процесса агрегации, что можно объяснить либо быстрым конформационным изменением, которое невозможно наблюдать на приборах NBS, либо отсутствием существенного изменения внутренней энергии белка при агрегации.

Человеческий белок тау представляет собой внутренне неупорядоченный белок (IDP), состоящий из 441 аминокислоты для так называемой изоформы 2N4R, которая, в частности, участвует в болезни Альцгеймера³³. Используя обратное рассеяние нейтронов на порошках пердейтерированного белка тау, мы показали, что динамика гидратационной воды увеличивается в состоянии волокна с более высокой популяцией молекул воды, претерпевающих поступательные движения. Результат свидетельствует о том, что увеличение энтропии гидратации воды может привести к фибрилляции амилоида тау.

Протокол

1. Подготовьте дейтерированный буфер для белков в жидком состоянии

1. Растворите все компоненты буфера в чистом виде D₂O.
2. Если pH-электрод был откалиброван в H₂O, отрегулируйте pD по формуле pD = pH + 0,4 с помощью NaOD или DCl³⁴.
   ПРИМЕЧАНИЕ: Использование D 2 O вместо H₂O может повлиять на растворимость белка, и, возможно, потребуется адаптировать буферные условия (например, путем небольшого изменения концентрации соли).

2. Приготовьте H₂O-гидратированные порошки пердейтерированного белка.

1. Подготовьте держатель образца.
   1. Тщательно очистите плоский алюминиевый держатель образца с уплотнением из индиевой проволоки и винтами водой и этанолом и дайте ему высохнуть.
      ПРИМЕЧАНИЕ: Плоский держатель образца используется таким образом, чтобы порошок мог быть равномерно распределен по поверхности. Количество порошка должно быть достаточным, чтобы его можно было удерживать между стенками и не падать при вертикальном размещении держателя образца.
   2. Взвесьте различные части держателя образца - дно, крышку и индиевую проволоку - отдельно на прецизионных весах.
   3. Поместите уплотнение из индиевой проволоки диаметром 1 мм в канавку нижней части держателя образца, оставив небольшое перекрытие в месте соединения двух концов (рис. 2A).
   4. Поместите соответствующее количество лиофилизированного белка (обычно ~ 100 мг белка) так, чтобы он заполнил внутреннюю поверхность нижней части держателя образца.
2. Увлажните протеиновый порошок.
   1. Поместите держатель образца в эксикатор с чашкой Петри, содержащей порошок P₂O₅ , на 24 часа, чтобы полностью высушить протеиновый порошок³⁵ (рис. 2B). Взвесьте сухую нижнюю часть держателя образца, содержащую индиевую печать и сухой порошок, чтобы получить_{m сухой}.
      ВНИМАНИЕ: порошок P₂O₅ очень агрессивен.
   2. Извлеките P 2 O₅ из эксикатора и поставьте чашку Петри с D₂O внутри. Регулярно контролируйте массу порошка, чтобы проверить уровень гидратации h = m hyd / m dry, где m_hyd и m_dry - это масса гидратированного порошка и сухого порошка соответственно.
      ПРИМЕЧАНИЕ: Для высокогидрофобных белков, таких как инсулин, может потребоваться повысить температуру внутри эксикатора, чтобы получить более высокое давление пара и достичь желаемого уровня гидратации h.
   3. Повторите шаги 2.2.1 и 2.2.2 не менее трех раз, чтобы правильно преобразовать все обменные атомы водорода в дейтроны.
      ПРИМЕЧАНИЕ: В качестве альтернативы для лучшего обмена H/D можно использовать циклы сублимационной сушки и растворения в чистом D₂O при условии, что он не влияет на белок.
   4. Гидратируйте порошок до уровня чуть выше желаемого, оставьте нижнюю часть держателя образца с индиевой проволокой и гидратированным порошком на прецизионных весах и подождите, пока масса медленно уменьшится до желаемого значения, чтобы получить целевой h (обычно 0,2-0,4, если глобулярный белок среднего размера должен быть покрыт одним полным гидратационным слоем).
   5. Быстро наденьте крышку на нижнюю часть и сначала закройте держатель образца четырьмя винтами, чтобы остановить парообмен (дополнительный рисунок S3A).
   6. Установите и затяните все оставшиеся винты до тех пор, пока между нижней частью и крышкой не останется зазора (дополнительный рисунок S3B).
   7. Взвесьте герметичный держатель образца, чтобы проверить возможные потери гидратации из-за утечек после нейтронного эксперимента.

3. Провести эксперимент по некогерентному рассеянию нейтронов

1. Обсудите и перепроверьте конфигурацию прибора, необходимого для эксперимента, с локальным контактом за несколько недель до назначенного времени луча.
2. Подготовьте образец жидкого состояния.
   1. Растворите белок в дейтерированном буфере.
   2. Определите подходящий объем жидкости, который будет помещен в держатель для образца, используя воду (убедитесь, что при закрытом держателе для образца нет перелива; Рисунок 2В).
      ПРИМЕЧАНИЕ: Следующие этапы (3.3 и 3.4) описывают эксперимент, проведенный на спектрометре NBS IN16B на ILL ^8,9 с использованием криопечи в качестве среды для образцов. Система управления прибором будет меняться от одного прибора к другому, но принципы работы останутся прежними.
3. Вставьте образец.
   1. Тщательно высушите стержень для образца (рис. 2D) и удалите предыдущий образец, если таковой имеется, убедившись, что доза ионизирующего излучения ниже 100 мкЗв/ч, прежде чем приступать к работе с любым материалом (в ILL).
   2. Поместите образец, проверьте правильность центрирования относительно центра луча (дополнительный рисунок S4) и вставьте палочку для образца в криопечь (рис. 2D). Включите вакуумный насос до температуры менее ^10-3 бар и промойте воздух внутри криопечи, повторив следующие три раза: заполните криопечь газообразным гелием до достижения атмосферного давления и снова удалите газ с помощью вакуумного насоса.
      ПРИМЕЧАНИЕ: В случае плоского держателя образца держатель образца должен быть ориентирован под углом 45° относительно входящего луча. Полезный диапазон передачи импульса может быть уменьшен из-за поглощения и рассеяния ячейкой. Сильный поглотитель нейтронов, такой как кадмий, может быть использован для маскировки определенных частей держателя образца (например, винтов, толстых деталей).
   3. Введите немного газообразного гелия в криопечь так, чтобы давление составляло ~ 0,05 бар.
4. Сбор данных (например, при использовании NOMAD на IN16B в ILL предполагается, что пользователь предпочитает температуру 200 К перед получением спектра квазиупругих нейтронов (QENS), затем E/IFWS во время изменения температуры до 310 K при 0,5 K в минуту и, наконец, QENS при 310 K).
   1. С помощью NOMAD на вкладке «Выполнение» перетащите контроллер FurnaceCryostat на стартовую площадку. Установите температуру 200 К. Используйте быстрый режим и тайм-аут в 30 минут, чтобы температура успела стабилизироваться. Нажмите на значок вращающихся стрелок, чтобы запустить его в фоновом режиме, чтобы данные можно было получить во время снижения температуры.
   2. Перетащите контроллер IN16DopplerSettings , установите профиль скорости на Точная скорость , установленная Max ΔE, значение 0,00 мкВ и 128 каналов , чтобы получить конфигурацию EFWS.
   3. Перетащите контроллер подсчета , заполните поле «Субтитры » именем, позволяющим легко идентифицировать данные, и установите 60 повторений 30-секундных сканирований (дополнительный рисунок S5A).
   4. Перетащите контроллер IN16DopplerSettings , установите профиль скорости на синусоиду , установленную скоростью , со значением 4.5 м/с и 2,048 каналами , чтобы получить конфигурацию QENS.
   5. Перетащите контроллер Count с 4 повторениями по 30 минут сканирования (дополнительный рисунок S5B).
   6. Для изменения температуры перетащите контроллер FurnaceCryostat , установите температуру на 310 К, установите для параметра Рампа значение SetPoint с Δ = 0,05 К и 6 с. Используйте тайм-аут 220 минут (дополнительный рисунок S6A).
   7. Используйте цикл for с 65 повторениями. Внутрь вставьте контроллер IN16DopplerSettings , как показано на шаге 3.4.2, а затем отсчет 30 с. Затем вставьте IN16DopplerSettings, как описано ранее, но с использованием смещения энергии 1.5 мкВ и 1,024 канала с последующим однократным отсчетом в 3 минуты (дополнительный рисунок S6B).
   8. Чтобы получить последний QENS с разрешением 310 K, перетащите контроллеры IN16DopplerSettings и Count , настроенные так, как описано в шагах 3.4.4 и 3.4.5 соответственно.
   9. Нажмите кнопку «Пуск » (прямоугольный треугольник в нижней части окна), чтобы запустить сценарий.
      ПРИМЕЧАНИЕ: Каждый эксперимент потребует сбора калибровочных данных; то есть пустая ячейка для вычитания или поправок на поглощение, только буфер при различных температурах, используемых для моделирования фона, и измерение ванадия (или, что то же самое, образца при температуре 10 К или ниже) для получения функции разрешения прибора.

4. Анализ данных - QENS

1. Импортируйте набор данных с помощью метода 'IN16B_QENS.process()' в программном обеспечении Python nPDyn v3.x³⁶
   >>>> из nPDyn.dataParsers import IN16B_QENS
   
   >>> выборка = IN16B_QENS(
   ... <путь к файлам данных>
   ... [detGroup=<целочисленный файл или файл группировки детекторов в XML
   ... формат>]
   ... ). process()
   
   >>> выборка = sample.get_q_range(0,3, 1,8)
2. Выполните коррекцию данных (необязательно) с помощью следующих команд (см. документацию nPDyn для получения дополнительной информации, рис. 3):
   #it предполагается, что данные для пустой ячейки, ванадия и буфера
   # были импортированы уже в набор данных под названием 'empty_cell', 'vanadium',
   # и 'buffer', соответственно.
   
   # для вычитания пустых ячеек с коэффициентом масштабирования
   # (ошибки распространяются автоматически)
   >>> выборка = выборка - 0,95 * empty_cell
   
   # для коррекции с использованием коэффициента Паалмана-Пинга
   # (взаимоисключающий с приведенным выше примером)
   >>> образец = sample.absorptionCorrection(empty_cell)
   
   # для нормализации
   >>> образец = sample.normalize(ванадий)
   
   # для биннинга вдоль наблюдаемой оси
   # наблюдаемым является время агрегации здесь
   >>> образец = образец.bin(3, ось=0)
3. Подгонка калибровочных данных. Набор данных - образцы, пустая ячейка, дейтерированный буфер (при необходимости) и ванадий - могут быть установлены с помощью встроенных моделей или пользовательской модели (см. документацию nPDyn):
   >>> из импорта nPDyn.models.builtins (
   ... модельPVoigt,
   ... модельВода,
   ... modelCalibratedD2O,
   ... )
   
   # Встроенные модели используют столбцовый вектор импульса
   # Передача значений Q
   >>> q = ванадий.q[:, Нет]
   
   # ванадий установлен с использованием псевдо-профиля Фойгта
   >>> vanadium.fit(modelPVoigt(q))
4. Используйте встроенную модель для гидратации воды под названием «modelWater». Эта модель читается так, как показано на Eq (4)¹⁷
   (4)
   Где₀,_r и_t являются скалярами, учитывающими относительный вклад упругого сигнала, вращательных движений и поступательных движений соответственно; j₁(qd) — сферическая функция Бесселя^l-го порядка, где q — передача импульса; d расстояние O-H в молекуле воды; δ(ω) — дельта Дирака, которая здесь умножается на EISF; N — высший порядок используемой сферической функции Бесселя (обычно ~5); и являются вращательными и поступательными движениями Лоренца соответственно; b(q) — плоский фоновый термин. Сферические функции Бесселя дают относительный вклад каждого состояния углового момента молекул воды, а число N определяется на основе q-диапазона передачи импульса. В случае типичного NBS-спектрометра термины до N = 4 почти полностью объясняют сигнал (дополнительный рисунок S7).
   # Здесь уравнение 2 используется для гидратации воды
   # свертка с функцией разрешения и добавлением
   # Фон D2O выполняется автоматически с помощью кнопки
   # предоставленные аргументы
   >>> sample.fit(modelWater(q),
   ... res=ванадий,
   ... bkgd=буфер,
   ... volume_fraction_bkgd=0,95
   ... )
   ПРИМЕЧАНИЕ: Вклад вращательных и поступательных движений должен быть запутанным, чтобы быть совершенно строгим. Успех аддитивной модели следует объяснить наличием различных популяций воды на поверхности белка и ограниченным энергетическим диапазоном.
5. Для построения графика данных используйте следующее (рис. 4):
   >>> из графика импорта nPDyn.plot
   >>> сюжет (образец)

5. Анализ данных - температурный скачок, эластичное сканирование с фиксированным окном (EFWS)

1. Используйте процедуру, аналогичную разделу 4, для нормализации данных о переходе температуры по сигналу при самой низкой температуре (обычно 10 К):
   >>> из nPDyn.dataParsers import IN16B_FWS
   
   >>> выборка = IN16B_FWS(
   ... <путь к файлам данных>,
   ... detGroup=[detGroup=<файл группировки целых чисел или детекторов в формате XML>]
   ... ). process()
   
   # нормализация с 5 первыми точками на наблюдаемом
   # оси, которые соответствуют температуре
   >>> пример /= образец[:5].mean(0)
   
   # Используемый здесь диапазон добротности передачи импульса меньше
   # так как используемая модель действительна только для низкого q
   >>> выборка = sample.get_q_range(0,2, 0,8)
2. Для начала используйте простую гауссову модель, ширина которой определяется так называемым средним квадратичным смещением (MSD). Постройте и подгоните модель с помощью следующих команд:
   >>> Импорт numpy как NP
   >>> из импорта nPDyn.models Параметры, Модель, Компонент
   
   # a - коэффициент масштабирования
   >>> params = Параметры(
   ... a={'value': 1, 'bounds': (0, np.inf)},
   ... msd={'value': 1, 'bounds': (0, np.inf)}
   ... )
   
   >>> model = Model(params)
   >>> model.addComponent(Component(
   ... 'гауссовский',
   ... лямбда x, a, msd: a * np.exp(-x ** 2 * msd / 6)
   ... ))
   
   >>> sample.fit(model, x=sample.q[:, None])
   
   >>> сюжет (образец)
   ПРИМЕЧАНИЕ: Гауссовское приближение всегда справедливо для q²MSD << 1, но для относительного сравнения между образцами можно использовать более широкий диапазон передачи импульса. Более сложные модели, выходящие за рамки гауссовского приближения, были разработаны^37,38,39.

6. Анализ данных - эластичное и неэластичное сканирование с фиксированным окном (E/IFWS)

1. Как и на шаге 4, импортируйте набор данных, но с помощью класса 'IN16B_FWS':
   >>> из nPDyn.dataParsers import IN16B_FWS
   
   >>> выборка = IN16B_FWS(
   ... <путь к файлам данных>
   ... [detGroup=<целочисленный файл или файл группировки детекторов в формате XML>]
   ... ). process()
   
   >>> выборка = sample.get_q_range(0,3, 1,8)
2. Подгонка калибровочных данных и данных образца.
   1. Проанализируйте данные E/IFWS с помощью обобщенного MSD⁴⁰ или рассматривая их как грубые спектры QENS (имея только несколько точек данных на оси энергии). Когда E/IFWS рассматривается как грубая QENS, модели, используемые для QENS, используются для одновременной подгонки всего набора данных E/IFWS (глобальное соответствие передач энергии и передач импульса).
      ПРИМЕЧАНИЕ: Последнее решение, использующее модели для QENS на данных E / IFWS, используется здесь, где накладывается зависимость передачи импульса от диффузии центра масс и внутренней динамики белка.
   2. Моделирование динамики белков в жидкостях с использованием следующего эквалайзера (5) ('modelProteinJumpDiff' в nPDyn):
      (5)
      Где R(q,ω) — функция разрешения; β скаляр, независимый для каждой передачи импульса q; a₀ - коэффициент упругой некогерентной структуры (EISF); лоренцев, учитывающий диффузию центра масс с шириной, заданной уравнением (6); — лоренциан, включающий диффузию центра масс и вклад в соответствии с моделью скачкообразной диффузии¹⁴ , учитывающей внутреннюю динамику (Eq (7); являющийся подогнанным сигналом от D₂O, перемасштабированным его объемной долей в образце.
      γ = D_sq² (6)
      D_s — коэффициент самодиффузии.
      (7)
      D_i — кажущийся коэффициент диффузии для внутренней динамики, а τ — время релаксации для диффузионных движений.
      >>> из импорта nPDyn.models.builtins (
      ... модельPVoigt,
      ... modelProteinJumpDiff,
      ... modelCalibratedD2O,
      ... )
      
      # Встроенные модели используют столбцовый вектор импульса
      # Передача значений Q
      >>> q = ванадий.q[:, Нет]
      
      # ванадий установлен с использованием псевдо-профиля Фойгта
      >>> vanadium.fit(modelPVoigt(q))
      
      # для чистого D2O, модели с калиброванной шириной линии
      # для разных температур включен в nPDyn
      >>> buffer.fit(modelCalibratedD2O(q, temp=363))
      
      # Здесь уравнение 3 используется для жидких образцов
      # свертка с функцией разрешения и добавлением
      # Фон D2O выполняется автоматически с предоставленными аргументами #
      >>> sample.fit(modelProteinJumpDiff(q),
      ... res=ванадий,
      ... bkgd=буфер,
      ... volume_fraction_bkgd=0,95
      ... )
3. Постройте график подогнанных данных, используя:
   >>> из графика импорта nPDyn.plot
   >>> сюжет (образец)

Результаты

Агрегацию лизоцима в твердые частицы проводили при 90 ° C с концентрацией белка 50 мг / мл в дейтерированном буфере (0,1 М NaCl при pD 10,5). Образование твердых частиц вызывается повышением температуры до 90 °C и происходит в течение 6 часов (дополнительный рисунок S8). Сбор данных проводился на IN16B, как описано в протоколе выше (данные постоянно курируются ILL и доступны по адресу http://dx.doi.org/10.5291/ILL-DATA.8-04-811).

Спектр QENS был получен при 7 ° C, чтобы полнос...

Обсуждение

Нейтронная спектроскопия является единственным методом, позволяющим зондировать усредненную по ансамблю ps-ns динамику образцов белка независимо от размера белка или сложности раствора при использовании дейтерации⁶. В частности, исследуя самодиффузию белковых сборок в растворе, можно однозначно определить гидродинамический размер таких сборок. Тем не менее, метод обычно ограничен низким потоком нейтронов, что подразумевает длительное время сбора и требование большого количества образца (обычно 100...

Материалы

Список материалов, использованных в этой статье

Имя	Компания	Каталожный номер	Комментарии
Алюминиевый держатель			Отсутствует в продаже. Либо местный контактный лицо по прибору может предоставить их, либо они могут быть изготовлены на основе технического чертежа, который может быть предоставлен местным контактным лицом.
Хлорид дейтерия, 35 мас.% в D2O, ≥ 99 атом % D	Sigma-Aldrich	543047
Оксид дейтерия (D, 99,9%)	Eurisotop	DLM-4DR-PK
Dow Corning высоковакуумная силиконовая смазка	Sigma-Aldrich	Z273554-1EA
Этанол 96%, EMSURE Reag. Ph Eur	Sigma-Aldrich	1.5901
Стеклянный десикатор	VWR	75871-660
Стеклянная десикаторная пластина, 140 мм	VWR	89038-068
Индийская проволока, диаметр 1,0 мм (0,04 дюйма), Puratronic, 99,999%	Alfa Aesar	00470.G1
Лизоцим из белка куриного яйца диализированный, лиофилизированный, порошок, ~100 000 Ед/мг	Sigma-Aldrich	62970
nPDyn		v3.x	см. github.com/kpounot/nPDyn, функции модели для фитинга также включены в программное обеспечение
OHAUS AX324 Adventurer весы, внутренняя калибровка	Dutscher	92641
Pentoxide Phosphorus, ReagentPlus, 99%	Sigma-Aldrich	214701
Pipette ErgoOne 0.5-10 μ L	Starlab	S7100-0510
Пипетка ErgoOne 100-1,000 L	Starlab	S7100-1000
Пипетка ErgoOne 20-200 μ L		Starlab S7100-2200
Наконечник для пипетки TipOne 1,000 μ L	Starlab	S1111-6001
Наконечник для пипетки TipOne 10 μ L	Starlab	S1111-3200
Наконечник для пипетки TipOne 200 μ L	Starlab	S1111-0206
Раствор дейтероксида натрия, 40 мас.% в D2O, 99,5 атом % D	Sigma-Aldrich	372072