[논문 리뷰] A Finite Element framework for computation of protein normal modes and mechanical response

Mark Bathe|ArXiv.org|2007. 04. 04.

Protein Structure and Dynamics참고 문헌 72인용 수 5

한 줄 요약

이 논문은 단백질을 용매 노출 표면을 기반으로 균일한 등방성 탄성 고체로 모델링하는 유한요소(FE) 프레임워크를 제시한다. 이는 다양한 공간 해상도에서 정규 모드와 기계적 반응을 계산할 수 있도록 하며, 원자 수준의 상호작용(예: 전기적 효과)을 통합할 수 있다. 이 방법은 전자기적 정규 모드 분석 및 실험 결과와 비교해도 정확한 결과를 도출하며, 원자 수준의 상호작용을 고려할 수 있다.

ABSTRACT

A coarse-grained computational procedure based on the Finite Element Method is proposed to calculate the normal modes and mechanical response of proteins and their supramolecular assemblies. Motivated by the elastic network model, proteins are modeled as homogeneous isotropic elastic solids with volume defined by their solvent-excluded surface. The discretized Finite Element representation is obtained using a surface simplification algorithm that facilitates the generation of models of arbitrary prescribed spatial resolution. The procedure is applied to compute the normal modes of a mutant of T4 phage lysozyme and of filamentous actin, as well as the critical Euler buckling load of the latter when subject to axial compression. Results compare favorably with all-atom normal mode analysis, the Rotation Translation Blocks procedure, and experiment. The proposed methodology establishes a computational framework for the calculation of protein mechanical response that facilitates the incorporation of specific atomic-level interactions into the model, including aqueous-electrolyte-mediated electrostatic effects. The procedure is equally applicable to proteins with known atomic coordinates as it is to electron density maps of proteins, protein complexes, and supramolecular assemblies of unknown atomic structure.

연구 동기 및 목표

• 단백질 정규 모드와 기계적 거동을 예측하기 위한 군집화된 계산 프레임워크를 개발하는 것.
• 단백질의 원자 구조가 알려져 있지 않은 경우에도 그들의 모델링을 가능하게 하는 것.
• 수용성-염류 매개 전기적 효과를 포함한 특정 원자 수준의 상호작용을 기계적 반응 예측에 통합하는 것.
• 표면 단순화와 FE 메esh 이산화를 통해 해상도 적응형으로 단백질 기계적 거동을 연구할 수 있는 유연한 방법을 제공하는 것.
• 모든 원자 수준의 정규 모드 분석, Rotation Translation Blocks 방법 및 실험 데이터와의 비교를 통해 접근법을 검증하는 것.

제안 방법

• 단백질은 용매 노출 표면에 의해 정의된 균일한 등방성 탄성 고체로 모델링된다.
• 표면 단순화 알고리즘이 임의의 공간 해상도에서 이산화된 유한요소 메쉬를 생성한다.
• 유한요소 수식은 축방향 압축 하중에 대한 진동 정규 모드와 기계적 붕괴 하중을 구한다.
• 이 방법은 T4 박테리오파지 리소자임(변종)과 필라멘터스 아크틴에 적용되었으며, 모든 원자 수준 및 실험 기준과 비교되었다.
• 전기적 효과는 FE 프레임워크 내에서 연속체 전기역학 모델을 통해 통합된다.
• 이 방법은 원자 좌표 기반 모델과 원자 구조가 알려지지 않은 저해상도 전자 밀도도에 모두 적용 가능하다.

실험 결과

연구 질문

• RQ1 유한요소 프레임워크는 군집화된 탄성 모델을 사용하여 단백질 정규 모드를 정확하게 계산할 수 있는가?
• RQ2 축방향 압축 하중에 대한 필라멘터스 아크틴의 기계적 반응은 실험 및 모든 원자 수준 예측과 어떻게 비교되는가?
• RQ3 전기적 효과와 같은 원자 수준의 상호작용은 단백질의 연속체 탄성 모델에 어느 정도 통합될 수 있는가?
• RQ4 원자 좌표가 이용 불가능할 경우, 전자 밀도도로부터 기계적 거동을 신뢰성 있게 예측할 수 있는가?
• RQ5 FE 메쉬의 해상도가 정규 모드 및 붕괴 하중 예측 정확도에 어떤 영향을 미치는가?

주요 결과

• 유한요소 프레임워크는 T4 박테리오파지 리소자임 변종의 정규 모드를 잘 계산하여 모든 원자 수준의 정규 모드 분석 결과와 양호한 일치를 보였다.
• 필라멘터스 아크틴의 예측된 비틀림 붕괴 하중이 실험 측정값과 모든 원자 수준 결과와 일치하였다.
• 이 방법은 두 단백질 체계에 대해 Rotation Translation Blocks 절차와도 유사한 비교 결과를 도출하였다.
• 표면 단순화를 통해 사용자가 정의한 공간 해상도에서 메쉬 생성이 가능하여 다스케일 모델링을 지원한다.
• 연속체 모델을 통해 직접 전기적 효과를 통합할 수 있어 물리적 현실성 향상이 가능하다.
• 원자 좌표 기반 모델과 원자 구조가 알려지지 않은 전자 밀도도에 대해 모두 강력한 성능을 보였다.