Rapid simulation of protein motion: merging flexibility, rigidity and normal mode analyses

TL;DR

This paper introduces a rapid, combined computational approach using normal mode, rigidity, and geometric simulations to efficiently explore large-scale protein conformational changes on standard desktop computers.

Contribution

It presents an integrated method combining three inexpensive simulation techniques to quickly analyze protein motions, significantly reducing computational costs compared to traditional methods.

Findings

Method identifies specific motion types in diverse proteins.

It determines amplitude limits of conformational changes.

Large-amplitude motions are explored within minutes.

Abstract

Protein function frequently involves conformational changes with large amplitude on timescales which are difficult and computationally expensive to access using molecular dynamics. In this paper, we report on the combination of three computationally inexpensive simulation methods-normal mode analysis using the elastic network model, rigidity analysis using the pebble game algorithm, and geometric simulation of protein motion-to explore conformational change along normal mode eigenvectors. Using a combination of ELNEMO and FIRST/FRODA software, large-amplitude motions in proteins with hundreds or thousands of residues can be rapidly explored within minutes using desktop computing resources. We apply the method to a representative set of six proteins covering a range of sizes and structural characteristics and show that the method identifies specific types of motion in each case and determines their amplitude limits.