Shape-preserving elastic solid models of macromolecules

Mass-spring models have been a standard approach in molecular modeling for the last few decades, such as elastic network models (ENMs) that are widely used for normal mode analysis. In this work, we present a vastly different elastic solid model (ESM) of macromolecules that shares the same simplicity and efficiency as ENMs in producing the equilibrium dynamics and moreover, offers some significant new features that may greatly benefit the research community. ESM is different from ENM in that it treats macromolecules as elastic solids. Our particular version of ESM presented in this work, named αESM, captures the shape of a given biomolecule most economically using alpha shape, a well-established technique from the computational geometry community. Consequently, it can produce most economical coarse-grained models while faithfully preserving the shape and thus makes normal mode computations and visualization of extremely large complexes more manageable. Secondly, as a solid model, ESM's close link to finite element analysis renders it ideally suited for studying mechanical responses of macromolecules under external force. Lastly, we show that ESM can be applied also to structures without atomic coordinates such as those from cryo-electron microscopy. The complete MATLAB code of αESM is provided.