Рентгенівська кристалографія — як ми бачимо атоми

1. Кристалічні ґратки та елементарні комірки

Кристал — це тверде тіло, в якому атоми, іони або молекули розташовані в абсолютно періодичному тривимірному візерунку. Ця періодичність описується ґраткою Браве — нескінченною множиною точок, породженою трьома векторами ґратки a, b, c, де будь-яку точку ґратки можна досягти з початку координат цілочисельною комбінацією n₁a + n₂b + n₃c.

Елементарна комірка — найменша повторювана одиниця ґратки — паралелепіпед, визначений трьома векторами ґратки та кутами між осями α, β, γ. Уся структура кристала закодована у вмісті й геометрії однієї елементарної комірки. Сім кристалічних систем (триклінна, моноклінна, орторомбічна, тетрагональна, тригональна, гексагональна, кубічна) і 14 окремих типів ґраток Браве вичерпують усі можливі трансляційні симетрії у трьох вимірах.

Усередині елементарної комірки атоми займають конкретні дробові координати (x, y, z), кожна з яких коливається від 0 до 1. Повна структура кристала потім генерується застосуванням усіх кристалографічних операцій симетрії (обертань, відбиттів, гвинтових осей, площин ковзання), визначених просторовою групою кристала. Існує рівно 230 окремих просторових груп.

2. Індекси Міллера та кристалографічні площини

Будь-яку сукупність паралельних площин, що проходять через точки ґратки, можна описати трьома цілими числами (h, k, l), які називаються індексами Міллера. Вони визначаються як обернені значення дробових точок перетину площини з осями елементарної комірки a, b, c, зведені до найменших цілих чисел.

Наприклад, площина, що перетинає вісь a у 1/2, вісь b у 1/3 і йде паралельно c (перетин у нескінченності), має дробові перетини (1/2, 1/3, ∞), обернені значення (2, 3, 0), що дає індекси Міллера (230). Відстань d між сусідніми площинами сімейства (hkl) залежить як від індексів Міллера, так і від параметрів ґратки. Для кубічного кристала з параметром ґратки a:

Індекси Міллера забезпечують повний каталог усіх площин у кристалі. Кожне сімейство площин утворює окрему дифракційну пляму (рефлекс), а повний набір спостережених рефлексів становить дифракційну картину, яка кодує структуру кристала.

Обернена ґратка — математичний дуал прямої ґратки: кожна точка в оберненому просторі з координатами (h, k, l) відповідає одному сімейству площин у реальному просторі. Вектори оберненої ґратки a*, b*, c* задовольняють a·a* = 1, a·b* = 0 тощо. Дифракція відбувається, коли вектор розсіювання дорівнює вектору оберненої ґратки — це умова Лауе, еквівалентна закону Брегга.

3. Закон Брегга: 2d·sinθ = nλ

У 1913 році Вільям Генрі Брегг та його син Вільям Лоуренс Брегг вивели надзвичайно просту умову конструктивної інтерференції рентгенівських променів, розсіяних кристалом. Розглядаючи дифракцію як відбиття від кристалічних площин, вони показали, що чітка дифракційна пляма (брегівський рефлекс) спостерігається лише коли:

де:

• d — міжплощинна відстань площин (hkl) (в ангстремах або нанометрах)
• θ — кут між падаючим рентгенівським променем і відбивальною площиною (не нормаллю до площини)
• n — порядок дифракції (додатне ціле число; зазвичай n = 1, оскільки вищі порядки відповідають площинам з відстанню d/n)
• λ — довжина хвилі рентгенівського випромінювання (зазвичай 0,5–2,5 Å для кристалографії; характеристичне випромінювання Kα міді становить 1,5418 Å)

Фізична картина: рентгенівські промені, розсіяні послідовними паралельними площинами, проходять різницю ходу 2d·sinθ. Коли вона дорівнює цілому числу довжин хвиль, хвилі, розсіяні всіма площинами, складаються конструктивно. Коли умова не виконана, хвилі від багатьох тисяч площин реального кристала гасять одна одну деструктивною інтерференцією — розсіяна інтенсивність не спостерігається.

Закон Брегга одразу пояснює, чому для структур атомного масштабу потрібні рентгенівські промені, а не видиме світло: щоб дослідити міжплощинні відстані 1–10 Å, довжина хвилі має бути порівнянною, що ставить необхідне випромінювання твердо в рентгенівський діапазон (0,1–10 нм довжини хвилі). Видиме світло (400–700 нм) не може розрізнити атомні відстані.

4. Структурний фактор F(hkl)

Закон Брегга каже нам, коли виникає рефлекс, але не наскільки він сильний. Інтенсивність кожного рефлексу залежить від розташування атомів усередині елементарної комірки. Це кількісно виражається структурним фактором F(hkl) — комплексним числом, квадрат модуля якого дає інтенсивність:

де сума береться за всіма j атомами в елементарній комірці, кожен з яких має дробові координати (xⱼ, yⱼ, zⱼ), а:

• fⱼ — атомний фактор розсіювання атома j; пропорційний числу електронів (важчі атоми розсіюють сильніше). Також залежить від sinθ/λ через скінченний розмір електронної хмари.
• exp[2πi(hxⱼ + kyⱼ + lzⱼ)] — фазовий множник, що кодує положення атома j; скалярний добуток вектора індексів Міллера (h,k,l) з атомним положенням (xⱼ,yⱼ,zⱼ) визначає внесок у фазу

Виміряна інтенсивність рефлексу (hkl) дорівнює:

На практиці застосовується кілька поправкових факторів: фактор Лоренца (геометрична поправка на час збору даних), фактор поляризації (рентгенівські промені частково поляризовані) та фактор Дебая-Валлера B, що враховує теплові коливання — атоми, що коливаються навколо рівноважних положень, розмивають свою електронну густину, послаблюючи висококутові рефлекси:

5. Електронна густина та перетворення Фур'є

Структурні фактори F(hkl) є коефіцієнтами Фур'є електронної густини ρ(x, y, z) в елементарній комірці. Електронна густина — фізична величина, яку насправді потрібно визначити, — відновлюється оберненим перетворенням Фур'є:

де V — об'єм елементарної комірки, а сума береться за всіма виміряними рефлексами. Пряме перетворення пов'язує реально-просторову густину зі структурними факторами в дифракційному просторі; обернене перетворення повертає назад. Цей взаємний зв'язок центральний для кристалографії.

Карти електронної густини високої роздільної здатності (обчислені за багатьма тисячами рефлексів до великих sinθ/λ) чітко показують піки в атомних положеннях. При роздільній здатності 2 Å окремі атоми з'являються як добре розділені піки. При 1 Å стає видимою електронна густина зв'язків між атомами. Нижче 0,5 Å («надвисока роздільна здатність») можна картографувати неподілені пари й зв'язкові електрони.

На практиці карта електронної густини обчислюється на сітці, контурується на обраному числі стандартних відхилень над середньою густиною і підганяється атомною моделлю. Модель уточнюється ітеративною мінімізацією методом найменших квадратів R-фактора:

Хороші структури малих молекул досягають R ~ 0,03–0,05; білкові структури зазвичай сходяться до R ~ 0,15–0,25, причому вільний R-фактор (обчислений на відкладеному тестовому наборі) використовується для контролю перенавчання.

6. Проблема фаз

Тут криється центральна перешкода кристалографії: дифракційний експеримент вимірює інтенсивності I(hkl) ∝ |F(hkl)|², що дає нам амплітуду |F(hkl)| кожного структурного фактора. Але F(hkl) — комплексне число: воно має і амплітуду, і фазовий кут φ(hkl). Перетворення Фур'є, що відновлює електронну густину, потребує обох.

Оскільки детектори вимірюють інтенсивності (потужність), інформація про фазу втрачається. Без фаз обернене перетворення Фур'є обчислити неможливо, і структуру визначити не можна. Це проблема фаз — фундаментальний бар'єр кристалографії, і її розв'язання є мистецтвом цієї дисципліни.

Три основні стратегії розв'язують проблему фаз:

• Прямі методи — використовують статистичні співвідношення між фазами структурних факторів, накладені вимогою, що електронна густина має бути додатною скрізь. Працюють для малих молекул (до ~200 атомів); лежать в основі пакету програм SHELX. Джером Карле та Герберт Гауптман отримали Нобелівську премію з хімії 1985 року за розробку прямих методів.
• Ізоморфне заміщення — вимочування кристала в розчині важкого атома (платини, ртуті, золота). Важкі атоми сильно розсіюють і передбачувано змінюють інтенсивності. Порівнюючи нативну й похідну дифракційні картини, можна обчислити фази. Розалінд Франклін та інші використовували цей підхід.
• Молекулярне заміщення — якщо гомологічна структура вже відома, її можна розташувати (обертанням і зсувом) в елементарній комірці нового кристала, щоб отримати початкові фази. Це найпоширеніший сьогодні метод, оскільки PDB тепер містить так багато еталонних структур.

7. Функції Паттерсона та методи фазування

До розробки прямих методів функція Паттерсона надавала безфазовий спосіб отримати інформацію про міжатомні вектори. Артур Ліндо Паттерсон (1935) показав, що перетворення Фур'є від даних інтенсивності — а не від самих структурних факторів — дає корисну величину:

Функція Паттерсона має пік у (u, v, w) щоразу, коли два атоми в структурі розділені вектором (u, v, w). Це автокореляція електронної густини. Для структури з N атомами карта Паттерсона має N² − N піків (без урахування початку координат). Для малих структур ці піки можна інтерпретувати безпосередньо; для більших структур карта стає надто заповненою.

Підхід Паттерсона особливо потужний для локалізації важких атомів: атом ртуті має 80 електронів проти ~8 у середньому атомі білка, тому вектори Hg-Hg дають дуже виразні піки Паттерсона. Після знаходження положень важких атомів їхній внесок у фази можна обчислити й використати для фазування всього набору даних методом множинного ізоморфного заміщення (MIR).

Аномальна дисперсія (SAD/MAD-фазування) використовує той факт, що на енергіях рентгенівського випромінювання поблизу краю поглинання важкого елемента атомний фактор розсіювання набуває уявної складової. Це порушує закон Фріделя (зазвичай |F(hkl)| = |F(−h,−k,−l)|) і надає інформацію про фазу з різниць між парами Фріделя. Селен — включений у білки як селенометіонін — має свій край поглинання зручно в діапазоні синхротронних пучків, що робить Se-SAD домінантним методом фазування на сучасних синхротронах.

8. Синхротронне випромінювання та сучасні пучки

Лабораторні джерела рентгенівського випромінювання (генератори з обертовим анодом) випромінюють характеристичне випромінювання на фіксованих довжинах хвиль. Синхротронне випромінювання перевершує їх кардинально: електрони або позитрони, що циркулюють з релятивістськими швидкостями в накопичувальному кільці, випромінюють надзвичайно інтенсивні, високо колімовані, налаштовувані пучки рентгенівського випромінювання як побічний продукт доцентрового прискорення.

Ключові переваги синхротронних джерел над лабораторними:

• Потік — у 10⁸–10¹² разів більше фотонів за секунду. Набір даних, що вимагає днів на лабораторному джерелі, займає секунди на синхротроні.
• Налаштовуваність — довжина хвилі неперервно регульована в рентгенівському діапазоні, що дозволяє фазування аномальною дисперсією на специфічних для елементів краях поглинання
• Колімація — пучки діаметром до 1 мкм, що дозволяє мікрокристалографію крихітних кристалів і серійну кристалографію
• Часова роздільна здатність — експерименти типу «накачка-зондування» можуть відстежувати хімічні реакції в кристалах у діапазоні від мікросекунд до фемтосекунд

Сучасні пучки макромолекулярної кристалографії (MX) на таких установках, як Diamond Light Source (Велика Британія), ESRF (Франція), APS (США) та SPring-8 (Японія), високо автоматизовані: роботизовані змінники зразків міняють кристали, дифрактометри самостійно центрують зразки, а конвеєри обробки даних працюють у реальному часі. Досвідчений кристалограф може зібрати повні набори даних з десятків кристалів за одну зміну.

Рентгенівські лазери на вільних електронах (XFEL) — у LCLS (Стенфорд) та European XFEL (Гамбург) — просувають межі ще далі: імпульси тривалістю ~10 фс і надзвичайна пікова яскравість уможливлюють «дифракцію перед руйнуванням», збираючи дифракційну картину з одного нанокристала до того, як рентгенівський імпульс його зруйнує. Серійна фемтосекундна кристалографія (SFX) може розв'язувати структури чутливих до випромінювання зразків і вловлювати ферментативні проміжні стани, надто швидкоплинні для звичайних методів.

9. Білкова кристалографія та PDB

Білки — молекулярні машини життя: ферменти, мотори, транспортери, рецептори. Розуміння їхньої функції вимагає знання їхньої тривимірної форми з атомною роздільною здатністю. Білкова кристалографія вперше дала такі уявлення: міоглобін (Джон Кендрю, 1958), гемоглобін (Макс Перуц, 1960) і лізоцим (Девід Філліпс, 1965). Кендрю та Перуц поділили Нобелівську премію з хімії 1962 року.

Білкова кристалографія має виклики, не властиві малим молекулам:

• Вирощування кристалів — білки потрібно кристалізувати з розчину, зазвичай повільним додаванням осаджувачів (сульфату амонію, ПЕГ). Якість кристала надзвичайно чутлива до pH, температури, добавок і чистоти білка. Вирощування кристалів дифракційної якості залишається лімітуючим етапом для багатьох мішеней.
• Розмір — типовий білок зі 300 амінокислот має ~2 400 неводневих атомів. Асиметрична одиниця може містити кілька копій білка, молекули ліків та впорядковані молекули води — загалом десятки тисяч атомів.
• Радіаційне пошкодження — інтенсивні рентгенівські пучки генерують вільні радикали, що пошкоджують кристали під час збору даних. Кріоохолодження кристалів до 100 К (температура рідкого азоту) різко зменшує радіаційне пошкодження і тепер є стандартною практикою.
• Роздільна здатність — білкові кристали дифрагують з нижчою роздільною здатністю, ніж кристали малих молекул, зазвичай 1,5–3,5 Å, оскільки білкові кристали містять ~50% розчинника, що зменшує кількість кристалів з ідеальною упаковкою і обмежує якість дифракції.

При роздільній здатності 2 Å бічні ланцюги амінокислот добре розрізняються, видно молекули води, а геометрію зв'язування ліганда можна точно визначити — достатньо для розробки ліків на основі структури. При 1,2 Å (атомна роздільна здатність) стають видимими атоми водню, і можна картографувати розподіли заряду. При поточній практичній межі ~0,8 Å у сприятливих випадках довжини й кути зв'язків наближаються до точності квантово-хімічних розрахунків.

Банк даних білків (PDB), заснований у 1971 році з сімома структурами, тепер містить понад 220 000 макромолекулярних структур (станом на 2026 рік). Приблизно 85% визначено рентгенівською кристалографією; решту — кріоелектронною мікроскопією та ЯМР. Кожна структурна стаття, опублікована у провідних журналах, вимагає депонування координат і експериментальних даних у PDB, що робить його одним з найцінніших відкритих наукових ресурсів, що існують.