BackgroundOur study was designed to identify the differential attractor modules related with hypertrophic cardiomyopathy (HCM) by integrating clustering-based on maximal cliques algorithm and Attract method. MethodsWe firstly recruited the HCM-related microarray data from ArrayExpress database. Next, protein–protein interaction (PPI) networks of normal and HCM were constructed and re-weighted using spearman correlation coefficient (SCC). Then, maximal cliques were found from the PPI networks through the clustering-based on maximal cliques approach. Afterwards, highly overlapped cliques were eliminated or merged according to the interconnectivity, and then modules were obtained. Subsequently, we used Attract method to identify differential attractor modules, following by the pathway enrichment analyses for genes in differential attractor modules. ResultsAfter removing the cliques with nodes less than or equal to 4, 926 and 1118 maximal cliques in normal and HCM PPI networks were obtained for module analysis. Then, we obtained 32 and 55 modules from the PPI networks of normal and HCM, respectively. By comparing with normal condition, there were 5 module pairs with the same or similar gene composition. Significantly, based on attract method, we found that these 5 modules were differential attractors. Pathway enrichment analyses indicated that proteasome, ribosome and oxidative phosphorylation were the significant pathways. ConclusionsProteasome, ribosome and oxidative phosphorylation might play pathophysiological roles in HCM.