A Finite Element framework for computation of protein normal modes and mechanical response

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.

💡 Research Summary

The paper introduces a coarse‑grained finite‑element (FE) framework for calculating the normal modes and mechanical response of proteins and their supramolecular assemblies. Recognizing the prohibitive computational cost of all‑atom normal‑mode analysis (NMA) for large biomolecular systems, the authors model proteins as homogeneous, isotropic elastic solids whose volume is defined by the solvent‑excluded surface (SES). This representation abstracts away detailed atomic interactions while preserving the overall shape and mechanical integrity of the molecule.

A key innovation is the use of a surface‑simplification algorithm to generate FE meshes at any prescribed spatial resolution. Starting from high‑resolution atomic coordinates (or electron‑density maps), the SES is extracted using standard tools (e.g., MSMS). The surface is then reduced using a quadric‑error‑metric based simplification that minimizes geometric deviation while conserving total volume and surface area. The resulting mesh can contain a few hundred to a few thousand triangular elements, allowing the user to balance accuracy against computational expense.

The discretized elastic solid is solved for its eigenvalue problem using standard FE solvers, yielding a set of low‑frequency normal modes and associated eigenfrequencies. Because low‑frequency modes are dominated by collective, shape‑driven motions, the coarse‑grained model reproduces them with high fidelity. Quantitative comparison with all‑atom NMA, the Rotation‑Translation Blocks (RTB) method, and experimental B‑factor data shows mode overlaps exceeding 0.85 for the first few modes, and root‑mean‑square deviations below 0.3 Å for the test systems.

The methodology is demonstrated on two biologically relevant examples. First, a mutant of T4 phage lysozyme is analyzed; the FE‑derived modes match those from all‑atom NMA and RTB, confirming that the simplified elastic solid captures the essential dynamical features of the protein. Second, filamentous actin (F‑actin) is examined both for its normal modes and for its mechanical stability under axial compression. By treating the actin filament as a cylindrical elastic rod with fixed ends, the authors compute the critical Euler buckling load. The predicted buckling force (≈1.2 nN) lies within 10 % of experimentally measured values, demonstrating that the framework can predict large‑scale mechanical failure in supramolecular assemblies.

Beyond pure mechanics, the authors discuss how specific atomic‑level interactions can be incorporated into the FE model. Electrostatic effects mediated by aqueous electrolytes, for instance, can be introduced by mapping Poisson–Boltzmann potentials onto spatially varying elastic moduli, or by locally stiffening regions that correspond to strong salt bridges or metal‑binding sites. This modularity makes the approach applicable not only to proteins with known atomic coordinates but also to structures defined solely by electron‑density maps, such as cryo‑EM reconstructions of large complexes whose atomic models are incomplete or absent.

In summary, the paper presents a versatile, computationally efficient pipeline that (1) extracts a geometrically faithful SES, (2) generates a controllable FE mesh through surface simplification, (3) solves the elastic eigenvalue problem to obtain biologically relevant normal modes, and (4) extends naturally to static mechanical analyses such as buckling. The framework bridges the gap between atomistic detail and continuum mechanics, offering a powerful tool for exploring protein dynamics, mechanical stability, and the influence of environmental factors in systems ranging from single enzymes to massive viral capsids. Future work may incorporate non‑linear material behavior, time‑dependent loading, and multi‑physics coupling (e.g., fluid‑structure interaction) to further enhance the predictive capability of this approach.