PP-MAPS: dynamic pharmacophore signatures of protein-peptide interfaces from molecular dynamics trajectories

Protein-peptide interactions underlie many cellular signaling and regulatory processes and are increasingly exploited in drug discovery. Characterizing such interfaces often requires the analysis of ensembles of conformations obtained by molecular modeling or molecular dynamics (MD) simulations, where transient contacts and alternative binding modes can be critical. Pharmacophore models provide an intuitive, transferable representation of molecular interactions. Dynophore i.e. "dynamic pharmacophore" approaches have been developed for small-molecule ligands with MD information. We present PP-MAPS (Protein-Peptide Molecular dynamics Assisted Pharmacophore Signatures), an open-source workflow that extracts and aggregates pharmacophore interactions along MD trajectories of protein-peptide complexes. PP-MAPS produces per-residue interaction frequencies and pharmacophore heatmaps that facilitate comparison of peptides, binding sites and receptor variants. PP-MAPS is implemented in Python and is available under an open-source license at https://github.com/camilledepenveiller/PP-MAPS. The workflow relies on GROMACS for trajectory processing and can use either LigandScout or the Chemical Data Processing Toolkit (CDPKit) for pharmacophore feature detection.