Leveraging protein dynamics to identify cancer mutational hotspots in 3D-structures

Abstract

Large-scale exome sequencing of tumors has enabled the identification of cancer drivers using recurrence and clustering-based approaches. Some of these methods also employ three-dimensional protein structures to identify mutational hotspots in cancer-associated genes. In determining such mutational clusters in structures, existing approaches overlook protein dynamics, despite the essential role of dynamics in protein functionality. In this work, we present a framework to identify driver genes using a dynamics-based search of mutational hotspot communities. After partitioning 3D structures into distinct communities of residues using anisotropic network models, we map variants onto the partitioned structures. We then search for signals of positive selection among these residue communities to identify putative drivers. We applied our method using the TCGA pan-cancer atlas missense mutation catalog. Overall, our analyses predict one or more mutational hotspots within the resolved structures of 434 genes. Ontological and pathway enrichment analyses implicate genes with predicted hotspots to be enriched in biological processes associated with tumor progression. Additionally, a comparison between our approach and existing hotspot detection methods that use structural data suggests that the inclusion of dynamics significantly increases the sensitivity of driver detection.