Количественная оценка вторичной структуры методом кругового дихроизма: референсные наборы, нормировка по концентрации и пределы критериев сходства

1 Введение

Среди методов характеристики структуры высшего порядка КД занимает особое место: он внедрён раньше других, проще всего в исполнении и легче всего поддаётся избыточной интерпретации. С КД-спектрометром, кварцевой кюветой и продувкой азотом спектр получают за полчаса; один щелчок «подогнать» в программе даёт четыре процента. Ни один шаг не выдаёт ошибки, поэтому эти четыре числа часто попадают в отчёты как результаты измерений.

Однако по меньшей мере четыре звена этой цепочки вносят существенную неопределённость, и ни одно из них не подаёт сигнала тревоги:

• Референсный набор. Подгонка исходит из того, что «неизвестный спектр = линейная комбинация нескольких спектров известной структуры». Для глобулярных белков это примерно верно, для богатых β-слоями антител — в меньшей степени: в референсных наборах не хватает спектров, близко гомологичных вариабельным и Fc-доменам антител.
• Алгоритм. Разные алгоритмы дают для одного спектра разные результаты, обычно с разницей 5–10 процентных пунктов. Это не ошибка программирования: они по-разному ограничивают «решение».
• Концентрация. Средняя остаточная эллиптичность имеет концентрацию в знаменателе: насколько ошибочна концентрация, настолько же и в том же направлении ошибочна сообщаемая [θ]_MRW. Важно, однако, что подогнанные доли определяются главным образом формой спектра и мало чувствительны к общему масштабированию — поэтому ошибка концентрации не делает четыре числа «абсурдными», но смещает абсолютные значения, сравнения абсолютной интенсивности и сопоставимость между образцами в целом; это один из самых трудно обнаруживаемых видов смещения.
• Рабочий диапазон прибора. Световой поток в дальнем УФ (особенно 190–200 нм) очень мал; при недостаточной продувке азотом или стареющей лампе измеряется, возможно, шум, а не сигнал. Проявляется это тем, что «спектр всё ещё выглядит гладким».

Цель статьи — разобрать эти четыре звена по отдельности и для каждого дать основания для суждения и требования к отчётности.

2 Принципы

2.1 Физическая природа кругового дихроизма

Хиральная молекула по-разному поглощает лево- и правоциркулярно поляризованный свет; разность ΔA = A_L − A_R называется круговым дихроизмом. Прибор измеряет разность интенсивностей двух пучков и пересчитывает её в эллиптичность θ (в милиградусах, mdeg): θ = 32980 × ΔA.

Сигнал КД возникает только у хромофоров в хиральном окружении. В белках это главным образом пептидная связь (переходы в дальнем УФ) и ароматические боковые цепи (переходы в ближнем УФ). Эта физическая природа задаёт две естественные области КД:

• Дальний УФ (190–250 нм): сигнал дают переходы n→π* и π→π* пептидной связи; он отражает упорядоченность остова, то есть вторичную структуру.
• Ближний УФ (250–320 нм): сигнал дают боковые цепи фенилаланина, тирозина и триптофана и дисульфидные связи; он отражает микроокружение этих остатков, то есть третичную структуру.

Информация двух областей совершенно разной природы, и одна не может заменить другую.

2.2 Дальний УФ: линейная суперпозиция вторичных структур

Три основные вторичные структуры имеют характерную форму спектра в дальнем УФ: α-спираль — два минимума близкой интенсивности при 208 и 222 нм; β-слой — минимум при 216–218 нм и переход в положительную область около 195 нм (его форма сильно зависит от параллельного/антипараллельного типа и длины тяжей, это наименее «стандартный» из трёх типов); неупорядоченная структура — сильный минимум около 200 нм и переход в положительную область на более коротких волнах.

Измеренный спектр белка — взвешенная суперпозиция этих компонентов, веса равны долям остатков в каждой структуре. Это физическая основа обращения при подгонке: по набору референсных спектров известной структуры ищут неотрицательные веса, взвешенная сумма с которыми приближает измеренный спектр.

Из линейной суперпозиции прямо следует главное ограничение КД: он чувствителен только к изменению «количества», но не «положения». Если остаток переместится с третьей на четвёртую позицию α-спирали или изменится сайт гликозилирования, спектр КД практически не изменится, пока общие доли вторичной структуры остаются прежними. КД не видит «деталей конформации», только её «состав».

Для α-спирали есть дополнительный критерий — отношение эллиптичности при 222 и 208 нм. Отношение меньше 1 (обычно 0,8–0,9) указывает на одиночные спирали; заметно больше 1 (часто выше 1,1) — на переплетение спиралей по типу coiled-coil. Это отношение полезно при оценке структуры слитых белков и пептидов. Для коллагеноподобных белков есть свой критерий: интактная тройная спираль в дальнем УФ даёт «положительный пик при 220–225 нм + отрицательный пик при 196–198 нм», и отношение интенсивностей положительного и отрицательного пиков R_pn (отношение абсолютных значений θ₂₂₁ и θ₁₉₈) характеризует целостность тройной спирали (около 0,10–0,15 у интактного коллагена, при денатурации падает до 0 или становится отрицательным)¹¹. Это один из немногих случаев, когда КД даёт бинарный ответ «определённая структура есть/нет», и он надёжнее подгонки долей.

2.3 Ближний УФ: окно в третичную структуру через ароматические остатки и дисульфиды

Сигналы КД в ближнем УФ на один-два порядка слабее, чем в дальнем, но информация иного рода: она отражает полярность и подвижность окружения ароматических остатков и их взаимодействие с соседними группами.

Ключевой признак — форма спектра, а не интенсивность. Когда часть белка разворачивается и ароматические остатки выходят из гидрофобного ядра в воду, положение пиков и их относительная интенсивность в ближнем УФ меняются, тогда как дальний УФ может почти не измениться. Поэтому в цепочке доказательств HOS ближний УФ КД и собственная флуоресценция относятся к одному уровню (третичная структура) и могут подтверждать друг друга.

Практическая трудность ближнего УФ КД — отношение сигнал/шум: для пригодного сигнала концентрацию нужно поднять до 0,5–1 мг/мл и использовать длину оптического пути 10 мм, что, в свою очередь, создаёт риск агрегации и эффекта внутреннего фильтра. Поэтому данные ближнего УФ обычно хуже воспроизводятся, чем данные дальнего УФ, и в исследованиях сопоставимости им следует придавать соответственно меньший вес.

2.4 Кривые термической денатурации и температура плавления

При фиксированной длине волны (обычно 222 нм, где α-спираль даёт наибольший сигнал) температуру непрерывно повышают и регистрируют эллиптичность как функцию температуры, получая S-образную кривую денатурации. Точка перегиба или экстремум первой производной — температура плавления T_m.

Важны три момента. Во-первых, одиночная длина волны 222 нм отражает только потерю α-спирали. У белков с преобладанием β-слоя (например, IgG) изменение сигнала при 222 нм мало и точность определения T_m падает; вместо этого можно использовать 208 нм или анализ главных компонент по многим длинам волн. Во-вторых, T_m — термодинамическая величина, а не вся стабильность: высокая T_m означает, что для разворачивания нужно больше энергии, но не означает отсутствия агрегации при температуре хранения; это определяется кинетикой и должно дополняться температурой начала агрегации T_on по DLS. В-третьих, скорость нагрева влияет на кажущуюся T_m: при быстром нагреве система не успевает прийти в равновесие и кажущаяся T_m завышается; скорость должна быть зафиксирована и указана в методе (обычно 1 ℃/мин), а результаты, полученные при разных скоростях, между лабораториями несопоставимы.

3 Методика эксперимента

3.1 Концентрация, длина оптического пути и азот

Примечание: концентрацию и длину пути нужно выбирать совместно, а не по отдельности. Концентрации в таблице соответствуют пути 1 мм; при пути 0,1 мм для сохранения той же оптической плотности концентрацию нужно повысить примерно на порядок (обычно до 1–2 мг/мл, при этом нужно следить за агрегацией). Это соответствие следует внести в метод и экспериментально проверить при квалификации прибора.

Продувка азотом не является факультативной. Кислород заметно поглощает ниже 200 нм (особенно при 185–190 нм); при недостаточной продувке поглощение кислорода примешивается к измеренному «сигналу белка» и проявляется как аномальное отрицательное смещение и рост шума в области 190–200 нм. Критерий — не «гладкий ли спектр» (такая помеха часто не даёт ни пиков, ни скачков и выглядит совершенно нормально), а сравнение спектров одного образца до и после полноценной продувки: если после продувки абсолютная эллиптичность при 190–200 нм уменьшилась и шум снизился, кислородная помеха была.

Длина пути связана с концентрацией: 0,2 мг/мл при пути 1 мм и 0,1 мг/мл при пути 2 мм дают одинаковое поглощение белка и одинаковый сигнал, но во втором случае вдвое возрастает и поглощение буфера и растворённого кислорода, поэтому суммарная оптическая плотность выше, а ошибка от рассеянного света больше. Если сигнала достаточно, предпочтительнее короткий путь и более высокая концентрация; цена — повышенный риск агрегации, который нужно контролировать по мутности или DLS.

3.2 Напряжение HT: критерий пригодности данных дальнего УФ

Это самый важный практический критерий раздела.

Во время сканирования КД-спектрометр синхронно регистрирует высокое напряжение (HT, high tension) фотоумножителя. При достаточном световом потоке HT остаётся низким (обычно 300–600 В; некоторые приборы показывают кВ, их нужно пересчитать перед сравнением); при недостатке света (слабая лампа на коротких волнах, сильное поглощение образца, поглощение кюветы) прибор автоматически повышает HT для поддержания сигнала, и шум растёт вместе с ним¹².

Отсюда правило разработки метода: сначала скан HT, затем нижняя граница длины волны. Если HT уже при 200 нм превышает порог, диапазон подгонки должен быть 200–250 нм, а не механически 190–250 нм.

3.3 Параметры и пригодность системы

Параметры сканирования (время отклика, ширина полосы, шаг, число накоплений) задаются по отраслевому стандарту Китая JY/T 0572—2020 «Общие правила спектроскопии кругового дихроизма»². Ширина полосы и шаг должны соответствовать друг другу: слишком широкая полоса сглаживает узкие пики, слишком большой шаг искажает положение пиков.

Пригодность системы проверяют по стандартному веществу: камфорсульфонат аммония (ACS) имеет определённую удельную эллиптичность при 290,5 нм (литературное значение [θ] ≈ +7 260 deg·cm²·dmol⁻¹) и служит для проверки шкалы эллиптичности прибора; его особенности в дальнем УФ можно использовать и для проверки точности длины волны. Его измеряют при каждом определении и подтверждают, что отклонение находится в допустимых пределах. Конкретное вещество, концентрацию и пределы нужно задать и валидировать в самом методе: отраслевой стандарт даёт лишь общие требования и не называет конкретного вещества.

4 Анализ данных и нормировка

4.1 Расчёт средней остаточной эллиптичности и её зависимость от концентрации

Исходные данные — эллиптичность θ (mdeg), которая линейно зависит от концентрации и длины пути и поэтому не может напрямую сравниваться между образцами. Нормировка переводит её в среднюю остаточную эллиптичность:

[θ]_MRW = θ × MRW / (10 × c × l)

где MRW — средняя масса остатка (молекулярная масса, делённая на число остатков), c — концентрация (мг/мл), l — длина пути (см); рассчитанная так [θ]_MRW выражается в deg·cm²·dmol⁻¹. Если θ вводится в mdeg, число получается в «mdeg·cm²·dmol⁻¹»; в отчёте единицу нужно пересчитать или указать — это один из самых частых источников расхождений при межлабораторных сравнениях.

В этой формуле MRW задаётся последовательностью, l — номинальной длиной кюветы; измерять нужно только c. А c обычно получают как A280 того же образца, делённую на коэффициент экстинкции. Получается цепочка:

спектр КД → [θ]_MRW → c → A280 → коэффициент экстинкции ε → измеренный аминокислотный состав

Смысл этой цепочки: неопределённость ε передаётся в абсолютное значение [θ]_MRW в соотношении 1:1 — смещение ε на 5% даёт смещение на 5% всех нормированных значений эллиптичности, а форма спектра остаётся неизменной, и на графике КД этого не видно. Следует подчеркнуть, что это смещение не равно пропорциональному смещению подогнанных долей: доли определяются главным образом формой спектра и относительными амплитудами референсных спектров компонентов и мало чувствительны к общему масштабированию спектра (все три алгоритма решают задачу при ограничении «сумма компонентов = 1», поэтому результат больше зависит от формы, чем от абсолютной амплитуды); влияние ошибки концентрации на доли компонентов частичное и непропорциональное, обычно заметно меньше 1:1. Опасность именно в том, что оно «не подаёт сигнала»: искажённые абсолютные значения, искажённые критерии по абсолютной интенсивности (θ₂₂₂/θ₂₀₈, R_pn, пороги эллиптичности в одной точке) и общий сдвиг сопоставимости между сериями и лабораториями не оставляют следа ни в спектре, ни в показателях качества подгонки.

На практике отсюда три требования: КД и УФ должны использовать одно и то же разведение и одно и то же определение концентрации; результаты КД сообщаются вместе с концентрацией и её источником; для исследований сопоставимости с высокими требованиями к точности сначала определяют коэффициент экстинкции по измеренному аминокислотному составу и нормируют по полученной концентрации.

4.2 Выбор алгоритма подгонки и его неопределённость

Для подгонки вторичной структуры существует несколько алгоритмов, и три основных различаются уже по принципу: SELCON3 основан на подгонке методом наименьших квадратов по базисным функциям с проверкой самосогласованности референсного набора; CONTINLL использует гребневую регрессию с регуляризацией, устойчивее к шуму и даёт более гладкие решения; CDSSTR сначала понижает размерность референсного набора сингулярным разложением, а затем подгоняет в пространстве малой размерности, сильнее завися от полноты референсного набора.

Различие между ними — не в точности, а в том, что они по-разному ограничивают «решение»: может ли сумма компонентов отклоняться от 1, допускаются ли необъяснённые остатки, какова сила регуляризации. Поэтому разброс доли α-спирали в 5–10 процентных пунктов для одного спектра, обработанного всеми тремя, находится в пределах нормы¹⁰. На результат влияют также полнота и версия референсного набора; Whitmore и Wallace систематически сравнили работу разных референсных наборов на однотипных белках⁹, а протокол Greenfield даёт полную рабочую схему от спектра до вторичной структуры⁸.

Отсюда два методических требования. Во-первых, метод должен фиксировать алгоритм и референсный набор: нельзя «пересчитать в другой программе и посмотреть» — смена алгоритма означает смену метода. Во-вторых, отчёт должен содержать показатели качества подгонки: остаток подгонки (NRMSD), отклонение суммы компонентов от 1 и название использованного референсного набора. Без этих показателей по четырём процентам читатель не может судить о надёжности результата.

4.3 Количественная оценка сходства

В исследованиях сопоставимости часто накладывают два спектра КД и смотрят, «совпадают ли» они. Это нужно выражать количественно, иначе вывод субъективен. Доступны три вида количественных показателей:

Коэффициент корреляции спектров — коэффициент Пирсона двух спектров, рассчитанный по длинам волн; чем ближе к 1, тем больше сходство. Недостаток — нечувствительность к общей интенсивности: два спектра одинаковой формы, различающиеся по амплитуде на 10%, всё равно дают корреляцию около 1.

Среднеквадратичное значение разностного спектра (RMSD) — поточечное вычитание с последующим взятием среднеквадратичного; отражает различия амплитуды, но чувствительно к единичным выбросам.

Разность подогнанных компонентов — подгонка каждого спектра и сравнение долей компонентов с указанием разностей и их неопределённости. Это показатель, ближе всего отражающий «структурное различие», но он ограничен неопределённостью алгоритма.

На практике их следует использовать вместе: корреляция подтверждает форму, RMSD — амплитуду, разность компонентов — структуру. Любого одного из них недостаточно для вывода.

5 Результаты и обсуждение

5.1 Насколько расходятся алгоритмы

Если подогнать типичный спектр IgG в дальнем УФ тремя алгоритмами, доли α-спирали обычно различаются на 5–10 процентных пунктов, доли β-слоя — меньше (у IgG преобладает β-слой, хорошо покрытый референсными наборами), а распределение между поворотами и неупорядоченной структурой часто различается так же сильно или сильнее: эти два типа трудно чётко различить в референсных спектрах, и у алгоритмов большая свобода в их распределении. Их сумма (повороты + неупорядоченная структура) обычно стабильнее каждого компонента по отдельности, поэтому некоторые авторы предлагают сообщать только три класса: «α-спираль / β-слой / прочее».

Методический смысл этого факта: у числа «18% α-спирали» только одна значащая цифра. Указывать «18,3%» значит выдавать неопределённость алгоритма за точность измерения. Правильно — указывать интервал или «около 18%» и называть в методе алгоритм и референсный набор.

5.2 Почему спектры КД разных IgG «выглядят одинаково»

Это самый частый пример неверного использования КД при характеристике антител. Во вторичной структуре IgG преобладает β-слой (подгонка КД обычно даёт 40%–45%, в зависимости от алгоритма и референсного набора, см. 4.2), и этот состав довольно постоянен у разных IgG. Гликозилирование меняет локальное окружение вблизи двух доменов C_H2 Fc-области и затрагивает несколько остатков; их вклад в спектр дальнего УФ тонет в общем сигнале β-слоя.

В итоге два спектра отлично накладываются, но это не доказывает идентичность конформации. Различия могут быть в ближнем УФ (микроокружение ароматических остатков), в термостабильности (разница T_m на 1–2 ℃) или в локальной подвижности (скорость водородно-дейтериевого обмена), и ко всему этому дальний УФ КД нечувствителен.

5.3 Ортогональное подтверждение с FTIR

КД и FTIR оба дают доли вторичной структуры, но на разной физической основе: КД — на хиральных различиях электронных переходов, FTIR — на колебательных диполях амидной связи. Различаются и чувствительные полосы: КД очень чувствителен к признаку α-спирали при 222 нм, FTIR — к β-слою при 1620–1640 cm⁻¹.

Поэтому согласие двух методов особенно ценно: когда КД и FTIR дают одинаковую долю β-слоя, вывод сильнее, чем при любом методе в отдельности. Напротив, явное расхождение указывает на помеху в одном из них (КД может страдать от нормировки по концентрации, FTIR — от водяного пара); нужно искать причину, а не усреднять. Доли компонентов, полученные одним методом, по сути являются оценками, поэтому многие лаборатории принимают внутреннее методическое правило «вывод о вторичной структуре вносится в отчёт только при согласии двух методов»; это стандарт качества, установленный компанией, а не нормативное требование, и в документах это следует указывать.

5.4 Место данных КД в досье и правила формулировок

В регистрационном досье данные КД появляются в трёх местах с разными требованиями к формулировкам. Подтверждение структуры относится к ICH Q6B, исследования сопоставимости — к ICH Q5E (сопоставимость биотехнологических/биологических продуктов при изменении производства) и вопросам-ответам по её применению, а спецификации — к тесту подлинности в частной статье на продукт:

Самая частая ошибка — перенос формулировки теста подлинности на сопоставимость: «спектры соответствуют» превращается в «структуры идентичны». Доказательная сила этих утверждений различается на порядок. Другая частая ошибка — писать «различий нет» вместо «различий не обнаружено»: второе отражает лишь разрешающую способность при использованном диапазоне и отношении сигнал/шум, и вместе с ним следует указывать разрешающую способность метода (RMSD между повторными спектрами, нижнюю границу коэффициента корреляции); иначе «не обнаружено» лишено критерия.

5.5 Валидация метода: по «характеристикам пригодности», а не по «полному набору пунктов»

Подход ICH Q2(R2) строит валидацию по характеристикам пригодности, а не по фиксированному перечню: назначение метода определяет, что валидировать, а не требует выполнять все пункты⁴. КД хорошо подходит для такого подхода: это и качественный метод (подлинность, сравнение спектров), и метод с полуколичественным результатом подгонки, а для этих двух задач нужны разные доказательства.

Рабочее разделение: для подлинности оцениваются специфичность (образцы после принудительной деградации или термической денатурации должны давать различимые спектры; отрицательный контроль и буфер не должны давать сигнала) и робастность (время продувки, длина пути, концентрация ±10%, скорость нагрева, разные партии кювет); для полуколичественного назначения добавляются повторяемость (6 последовательных измерений одного образца с RSD θ и каждой доли компонентов) и промежуточная прецизионность (разные дни, разные операторы). С линейностью нужна осторожность: эллиптичность действительно пропорциональна концентрации (закон Бугера–Ламберта–Бера), но результат метода — подогнанные доли, которые от концентрации не зависят; включение линейности в валидацию не соответствует фактическому назначению метода и вызывает вопросы об обосновании при проверке. Что касается правильности, у долей вторичной структуры нет общепринятого истинного значения; её можно подтвердить лишь косвенно — референсными белками или сравнением с FTIR и кристаллическими структурами, — и навязывать степень извлечения не следует.

Пригодность системы проверяется каждый раз: ACS — шкала эллиптичности, кривая HT — нижняя граница длины волны, шум базовой линии (RMS) — не выше заданного в методе значения. Все три проверки дешёвы и быстры и способны «подать сигнал» независимо от результата — это гораздо полезнее, чем рассматривать форму спектра постфактум.

6 Фармакопеи и стандарты

Ни в одной из трёх фармакопей нет отдельной статьи по электронному КД^1,6,7. Стоит отметить два момента. Во-первых, в USP 2025 есть статья 〈782〉 о колебательном круговом дихроизме (VCD), но VCD измеряет колебательные переходы (полезен для различения стереоизомеров и уточнения конформации) и не является тем же методом, что КД в дальнем УФ, основанный на электронных переходах; ссылаться одно на другое нельзя. Во-вторых, EP 2.2.42 — это плотность твёрдых веществ, к КД она отношения не имеет; при поиске этот номер легко ошибочно принять за статью по КД.

Практическое следствие: для метода КД нет фармакопейной статьи, на которую можно сослаться как на основание; можно сослаться только на отраслевой стандарт или собственный метод, а полную систему критериев нужно выстроить при валидации. Это принципиально отличается от практики указывать «ChP 0401» как основание для КД в досье: 0401 — абсорбционная спектрофотометрия в УФ и видимой области, оптически не связанная с циркулярно поляризованным светом; она может служить опорой лишь для общих элементов — калибровки прибора, линейности оптической плотности, подбора кювет, — но не специфическим основанием для КД. Ссылка на 0401 как на основание метода КД — формальный недостаток, на который эксперты легко указывают.

7 Выводы

Ценность КД в цепочке доказательств HOS очевидна, но его результат — подогнанная, а не измеренная величина, и это должно отражаться в разработке метода и в формулировках отчёта.

Во-первых, результат определяют референсный набор и алгоритм. Для одного спектра доля α-спирали может различаться на 5–10 процентных пунктов между алгоритмами, поэтому алгоритм и референсный набор нужно фиксировать и сообщать вместе с результатом; у «18% α-спирали» только одна значащая цифра.

Во-вторых, концентрация — самый скрытый источник ошибки. Средняя остаточная эллиптичность имеет концентрацию в знаменателе, а концентрация через A280 и коэффициент экстинкции выводится из аминокислотного состава, поэтому абсолютное значение [θ]_MRW наследует неопределённость УФ и аминокислотного анализа примерно в соотношении 1:1. Подогнанные доли мало чувствительны к общему масштабированию спектра, поэтому это смещение не подаёт сигнала ни в долях, ни в форме спектра; оно искажает критерии по абсолютной интенсивности (θ₂₂₂/θ₂₀₈, R_pn, пороги эллиптичности в одной точке) и сопоставимость между сериями и лабораториями.

В-третьих, объективный критерий пригодности данных дальнего УФ — напряжение HT, а не гладкость спектра. Длины волн, где HT превышает порог, нельзя использовать для количественной подгонки, а граничную длину волны следует определять по кривой HT.

В-четвёртых, сходство спектров не означает идентичности конформации. Спектры КД разных IgG очень похожи, а локальные изменения вроде гликозилирования лежат в пределах погрешности; наложение спектров КД подтверждает лишь «значимых различий не обнаружено», а выводы о HOS следует строить на сочетании методов — FTIR, ближнего УФ КД, собственной флуоресценции, DSC и др.

Литература