Supporting the CIF file format of proteins in molecular dynamics simulations

Abstract

Molecular dynamics (MD) simulations can capture the dynamic behavior of proteins in full atomic detail and at very fine temporal resolution, so they have become an important tool in the study of protein dynamics. To date, several MD packages are widely used. An MD simulation starts from an initial structure that is generally taken from the Protein Data Bank (PDB). Until 2014, the PDB format was the standard file format for protein structures. However, there are certain intrinsic limitations in the PDB format, such as the storage of structural information in a fixed-width format, which is an issue for very large protein complexes. Therefore, the CIF (crystallographic information framework) format has been proposed, which is characterized by its superior expansibility. To our knowledge, the current mainstream MD packages support only the PDB format but do not support the CIF format directly. In this study, we modified the source code of one of the MD packages, GROMACS, which enables it to support CIF-formatted structure files as input and subsequently generate molecular topology files. This work simplifies the preprocessing of large protein complexes for MD simulations.