Geo-Measures: A PyMOL plugin for protein structure ensembles analysis

Abstract

Although molecular dynamics encompasses several applications, studies focusing on biomolecular systems are central issues of this research area. Such simulations require the generation of trajectory files, which provide a path for the analysis and interpretation of results with biological significance. However, although several programs have been developed in Python language for the analyses of molecular dynamics (MD) trajectories, they usually require some knowledge of programming languages in order to write or run the scripts using command lines, which certainly hinders the access of MD simulations to many scientists with the necessary biological background to interpret their results. To ease the access to Python packages focusing on MD trajectory analyses, we built a user-friendly and easy-to-install graphical PyMOL interface. Geo-Measures integrates the PyMOL functionalities with MDTraj, a powerful library of trajectory analyses, allowing the users to access up to 14 different types of analyses. Two sample cases are reported here to demonstrate the use of Geo-Measures. In the first example, which involves the use a MD trajectory file of hemoglobin from the MoDEL MD bank, we exemplified the analyses of the following variables: root mean square deviation, radius of gyration, free energy landscape and principal component analysis. In the second case, we built a trajectory file for the ecto-5′-nucleotidase using the LiGRO program to study the carbon alpha pincer angles, to define the secondary structure of the proteins and to analyze the Modevectors. This user-friendly graphical PyMOL plugin, which can be used to generate several descriptive analyses for protein structures, is open source and can be downloaded at: https://pymolwiki.org/index.php/Geo_Measures_Plugin.