[논문 리뷰] Ab initio structure solution by charge flipping

G. Oszlányi, András Sütö|arXiv (Cornell University)|2003. 08. 07.

Enzyme Structure and Function인용 수 16

한 줄 요약

이 논문은 고해상도 X선 회절 데이터를 위한 단순한 아비지노 구조 해법 방법인 전하 뒤집기 전략을 소개한다. 이 방법은 푸리에 변환을 통해 실공간과 역공간을 번갈아 가며 진행된다. 임계값 δ 이하의 전하 밀도에 부호를 뒤집고, 가중치 없이 관측된 구조 인자 크기 |Fobs|를 강제로 적용함으로써, 사전에 위상 또는 대칭 정보 없이도 소분자 구조를 높은 성공률로 해결할 수 있으며, 고정점이 아닌 한계 사이클로 수렴한다.

ABSTRACT

In this paper we present an extremely simple structure solution method termed charge flipping. It works ab initio on high resolution x-ray diffraction data in the manner of Fourier recycling. The real space modification simply changes the sign of charge density below a threshold, while in reciprocal space the modification is the Fobs map without any weighting. We test the algorithm using synthetic data for a wide range of structures, analyse the solution statistics and check the quality of reconstruction. Finally, we reconsider mathematical aspects of the algorithm in detail, showing that in this chaotic iteration process the solution is a limit cycle and not a fixed point.

연구 동기 및 목표

• 사전 화학적 또는 위상 정보 없이도 작동하는 최소한의 개념적으로 단순한 아비지노 구조 해법 방법을 개발하는 것.
• 이dealized 고해상도 합성 데이터를 사용하여 중심대칭 및 비중앙대칭 소분자 구조의 광범위한 셋에 대해 전하 뒤집기 알고리즘의 타당성과 강건성을 테스트하는 것.
• 반복 과정의 수학적 행동, 특히 수렴 성질과 해의 성격을 분석하는 것.
• 분리된 정련 프로그램 없이도 모든 원자를 회복할 수 있음을 보여주는 것.

제안 방법

• 알고리즘은 프리델의 법칙을 만족하는 무작위 위상 세트로 시작하여, 역푸리에 변환을 통해 초기 실공간 전하 밀도를 생성한다.
• 각 반복 단계에서, 임계값 δ 이하의 전하 밀도 픽셀의 부호가 뒤집히고, 그 이상의 값은 그대로 유지된다.
• 수정된 실공간 밀도는 FFT를 통해 다시 역공간으로 변환되어 일시적인 구조 인자 G(k)를 계산한다.
• G(k)의 위상은 유지되지만, 크기는 관측된 구조 인자 진폭 Fobs(k)로 대체되며, F(0)는 자유롭게 변할 수 있다.
• 결과적으로 얻어진 F(k)는 역푸리에 변환을 통해 새로운 실공간 전하 밀도를 생성함으로써, 푸리에 재활용의 한 사이클을 완료한다.
• 이 과정은 최적화 없이 반복되며, 목적 함수로는 R-요구를 사용하지 않고, 단지 감시용으로만 사용된다.

실험 결과

연구 질문

• RQ1 실공간에서의 단순한 부호 뒤집기 규칙과 역공간에서의 가중치 없는 Fobs 제약 조건을 조합하여 아비지노 결정 구조를 성공적으로 해결할 수 있는가?
• RQ2 해의 수학적 성격은 무엇인가—특히 알고리즘이 고정점이 아니라 한계 사이클로 수렴하는가?
• RQ3 임계값 δ의 선택이 수렴, 총 전하 및 해의 품질에 어떻게 영향을 미치는가?
• RQ4 이상적인 데이터 조건 하에서 다양한 소분자 구조, 비중앙대칭 구조를 포함하여 이 방법이 어느 정도 성공을 거두는가?
• RQ5 유한한 해상도와 노이즈 조건에서 알고리즘이 어떻게 행동하며, 데이터의 완벽하지 못한 점에 대해 얼마나 강건한가?

주요 결과

• 전하 뒤집기 알고리즘은 0.8 Å 해상도의 합성 데이터를 사용하여 캠브리지 구조 데이터베이스에 포함된 모든 테스트 소분자 구조를 고해상도 성공률로 해결하였다.
• 모든 원자가 별도의 정련 단계 없이 회복되었으며, 이는 아비지노 구조 해법에서의 방법의 효과성을 입증한다.
• 총 전하 F(0)는 위상 무작위화로 인해 초도에 과도한 오버슈트를 보였고, 이후 급격히 감소하고 장기간 정체된 후 진짜 값의 약 30% 정도 낮은 값으로 수렴하였다.
• 반복 과정은 혼돈스러운 성질을 지녔으며, 고정점이 아닌 한계 사이클로 수렴하는 것으로 확인되었고, 위상 변화 및 R-요구의 행동으로부터 뒷받침되었다.
• 알고리즘의 성능은 포isson 노이즈 및 절대 척도와 등방성 온도 계수의 오차와 같은 데이터 오류에 대해 강건하였지만, 주요 분석에는 이상적인 데이터를 사용하였다.
• 수렴 후 δ를 줄이면 진동을 억제하고 총 전하를 증가시켜 해의 품질을 향상시켰지만, δ = 0일 경우 영역이 사라지고 전하가 과다해지는 현상이 발생하였다.