Systematic Analysis of Biomolecular Conformational Ensembles with PENSA

TL;DR

PENSA is an open-source Python software package that systematically analyzes biomolecular conformational ensembles, enabling unbiased feature representation, comparison of molecular states, and identification of allosteric pathways to understand molecular mechanisms.

Contribution

This work introduces PENSA, a novel software tool that provides comprehensive, unbiased analysis and comparison of biomolecular conformational ensembles, including a new method for quantifying information flow.

Findings

PENSA effectively identifies significant conformational differences.

It enables visualization and interpretation of biomolecular data.

Demonstrated usefulness in real-world biomolecular studies.

Abstract

Atomic-level simulations are widely used to study biomolecules and their dynamics. A common goal in such studies is to compare simulations of a molecular system under several conditions -- for example, with various mutations or bound ligands -- in order to identify differences between the molecular conformations adopted under these conditions. However, the large amount of data produced by simulations of ever larger and more complex systems often renders it difficult to identify the structural features that are relevant for a particular biochemical phenomenon. We present a flexible software package named PENSA that enables a comprehensive and thorough investigation into biomolecular conformational ensembles. It provides featurizations and feature transformations that allow for a complete representation of biomolecules like proteins and nucleic acids, including water and ion binding sites, thus avoiding bias that would come with manual feature selection. PENSA implements methods to systematically compare the distributions of molecular features across ensembles to find the significant differences between them and identify regions of interest. It also includes a novel approach to quantify the state-specific information between two regions of a biomolecule, which allows, e.g., tracing information flow to identify allosteric pathways. PENSA also comes with convenient tools for loading data and visualizing results, making them quick to process and easy to interpret. PENSA is an open-source Python library maintained at https://github.com/drorlab/pensa along with an example workflow and a tutorial. We demonstrate its usefulness in real-world examples by showing how it helps to determine molecular mechanisms efficiently.