Active natural compounds perturb the melanoma risk-gene network

Cutaneous melanoma is an aggressive type of skin cancer with a complex genetic landscape caused by the malignant transformation of melanocytes. This study aimed at providing an in silico network model based on the systematic profiling of the melanoma-associated genes considering germline mutations, somatic mutations, and genome-wide association study signals accounting for a total of 232 unique melanoma risk genes. A protein-protein interaction network was constructed using the melanoma risk genes as seeds and evaluated to describe the functional landscape in which the melanoma genes operate within the cellular milieu. Not only were the majority of the melanoma risk genes able to interact with each other at the protein level within the core of the network, but this showed significant enrichment for genes whose expression is altered in human melanoma specimens. Functional annotation showed the melanoma risk network to be significantly associated with processes related to DNA metabolism and telomeres, DNA damage and repair, cellular ageing, and response to radiation. We further explored whether the melanoma risk network could be used as an in silico tool to predict the efficacy of anti-melanoma phytochemicals, that are considered active molecules with potentially less systemic toxicity than classical cytotoxic drugs. A significant portion of the melanoma risk network showed differential expression when SK-MEL-28 human melanoma cells were exposed to the phytochemicals harmine and berberine chloride. This reinforced our hypothesis that the network modeling approach not only provides an alternative way to identify molecular pathways relevant to disease but it may also represent an alternative screening approach to prioritize potentially active compounds.