Эта статья может быть слишком технической для понимания большинства читателей . ( Май 2012 ) |
Молекулярное замещение (MR) [1] — это метод решения фазовой проблемы в рентгеновской кристаллографии . MR опирается на существование ранее решенной структуры белка, которая похожа на нашу неизвестную структуру, из которой получены данные дифракции. Это может быть гомологичный белок или структура ЯМР белка с более низким разрешением того же белка. [2]
Первой целью кристаллографа является получение карты электронной плотности, которая связана с дифрагированной волной следующим образом:
При использовании обычных детекторов измеряется интенсивность , а вся информация о фазе ( ) теряется. Тогда, при отсутствии фаз (Φ), мы не можем завершить показанное преобразование Фурье, связывающее экспериментальные данные рентгеновской кристаллографии (в обратном пространстве ) с электронной плотностью в реальном пространстве, в которое встроена атомная модель. MR пытается найти модель, которая наилучшим образом соответствует экспериментальным интенсивностям среди известных структур.
Мы можем вывести карту Паттерсона для интенсивностей, которая является межатомной векторной картой, созданной путем возведения в квадрат амплитуд структурных факторов и установки всех фаз в ноль. Эта векторная карта содержит пик для каждого атома, связанного с каждым другим атомом, с большим пиком в 0,0,0, где векторы, связывающие атомы с собой, «накапливаются». Такая карта слишком шумная, чтобы вывести какую-либо структурную информацию высокого разрешения — однако, если мы сгенерируем карты Паттерсона для данных, полученных из нашей неизвестной структуры и из структуры ранее решенного гомолога, в правильной ориентации и положении внутри элементарной ячейки , две карты Паттерсона должны быть тесно коррелированы. Этот принцип лежит в основе MR и может позволить нам вывести информацию об ориентации и местоположении неизвестной молекулы с ее элементарной ячейкой.
Из-за исторических ограничений вычислительной мощности поиск MR обычно делится на два этапа: вращение и перемещение .
В функции вращения наша неизвестная карта Паттерсона сравнивается с картами Паттерсона, полученными из нашей известной гомологической структуры в разных ориентациях. Исторически для оценки функции вращения использовались r-факторы и/или коэффициенты корреляции , однако современные программы используют алгоритмы на основе максимального правдоподобия . Самая высокая корреляция (и, следовательно, оценки) получаются, когда две структуры (известная и неизвестная) находятся в схожей ориентации(ях) — затем их можно вывести в углах Эйлера или сферических полярных углах.
В функции перевода теперь правильно ориентированная известная модель может быть правильно позиционирована путем перевода ее в правильные координаты в пределах асимметричной единицы. Это достигается путем перемещения модели, вычисления новой карты Паттерсона и сравнения ее с неизвестной полученной картой Паттерсона. Этот поиск методом грубой силы является вычислительно затратным, и теперь чаще используются быстрые функции перевода. Позиции с высокой корреляцией выводятся в декартовых координатах .
Благодаря улучшению прогнозирования структуры белка de novo многие протоколы, включая MR-Rosetta, QUARK, AWSEM-Suite и I-TASSER-MR, могут генерировать множество ложных структур, подобных нативным, которые полезны для решения фазовой проблемы путем молекулярной замены. [3]
После этого у нас должны быть правильно ориентированные и переведенные модели фазирования, из которых мы можем вывести фазы, которые (надеюсь) достаточно точны для получения карт электронной плотности. Их можно использовать для построения и уточнения атомной модели нашей неизвестной структуры.