[论文解读] A Finite Element framework for computation of protein normal modes and mechanical response

Mark Bathe|ArXiv.org|Apr 4, 2007

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.

研究动机与目标

• 开发一种粗粒度计算框架，用于预测蛋白质正常模态和机械行为。
• 实现无论是否已知其原子结构，均可对蛋白质和超分子组装体进行建模。
• 将特定原子级相互作用（包括水溶液电解质介导的静电效应）整合到机械响应预测中。
• 提供一种灵活、分辨率自适应的方法，通过表面简化和有限元离散化研究蛋白质力学。
• 通过与全原子正常模态分析、旋转平动块方法及实验数据对比，验证该方法的有效性。

提出的方法

• 蛋白质被建模为由其溶剂可及表面定义的均匀各向同性弹性固体。
• 通过表面简化算法在任意空间分辨率下生成离散化的有限元网格。
• 有限元公式求解轴向压缩下的振动正常模态和机械屈曲载荷。
• 该方法应用于T4噬菌体溶菌酶（突变体）和丝状肌动蛋白，结果与全原子和实验基准进行对比。
• 通过有限元框架内的连续介质静电模型引入静电效应。
• 该方法适用于基于高分辨率原子坐标和低分辨率电子密度图的建模。

实验结果

研究问题

• RQ1 有限元框架能否利用粗粒度弹性模型准确计算蛋白质正常模态？
• RQ2 丝状肌动蛋白在轴向压缩下的机械响应与实验和全原子预测结果相比如何？
• RQ3 在蛋白质连续弹性模型中，原子级相互作用（如静电作用）能在多大程度上被整合？
• RQ4 当缺乏原子坐标时，该方法能否可靠地从电子密度图预测机械行为？
• RQ5 有限元网格的分辨率如何影响正常模态和屈曲载荷预测的准确性？

主要发现

• 有限元框架成功计算了T4噬菌体溶菌酶突变体的正常模态，与全原子正常模态分析结果高度一致。
• 预测的丝状肌动蛋白临界欧拉屈曲载荷与实验测量值和全原子结果相符。
• 该方法在两种蛋白质体系中与旋转平动块方法的比较结果均表现良好。
• 表面简化技术实现了用户定义空间分辨率下的网格生成，支持多尺度建模。
• 该框架可通过连续介质模型直接整合静电效应，提升物理真实性。
• 该方法在基于原子坐标建模和未知原子结构的电子密度图上均表现出鲁棒性。