Temporal analysis of melanogenesis identifies fatty acid metabolism as key skin pigment regulator.

Therapeutic methods to modulate skin pigmentation has important implications for skin cancer prevention and for treating cutaneous hyperpigmentary conditions. Towards defining new potential targets, we followed temporal dynamics of melanogenesis using a cell-autonomous pigmentation model. Our study elucidates 3 dominant phases of synchronized metabolic and transcriptional reprogramming. The melanogenic trigger is associated with high MITF levels along with rapid uptake of glucose. The transition to pigmented state is accompanied by increased glucose channelisation to anabolic pathways that support melanosome biogenesis. SREBF1-mediated up-regulation of fatty acid synthesis results in a transient accumulation of lipid droplets and enhancement of fatty acids oxidation through mitochondrial respiration. While this heightened bioenergetic activity is important to sustain melanogenesis, it impairs mitochondria lately, shifting the metabolism towards glycolysis. This recovery phase is accompanied by activation of the NRF2 detoxication pathway. Finally, we show that inhibitors of lipid metabolism can resolve hyperpigmentary conditions in a guinea pig UV-tanning model. Our study reveals rewiring of the metabolic circuit during melanogenesis, and fatty acid metabolism as a potential therapeutic target in a variety of cutaneous diseases manifesting hyperpigmentary phenotype.

Boolean model for melanogenesis.

Melanogenesis is a highly regulated process through which the pigment melanin is produced in skin cells. Irregularities in the molecular events that govern the process of skin pigmentation can cause disorders like vitiligo. In order to understand the biology of disease progression, it is important to have an in depth understanding of intracellular events. Mathematical models provide an integrated view of intracellular signalling. There are very few models to date that incorporate intracellular processes relevant to melanogenesis and only one to our knowledge that simulates the dynamics of response to varying levels of input. Here, we report the formulation of the largest Boolean model (265 nodes) for melanogenesis to date. The model was built on the basis of a detailed interaction network graph published by Raghunath et al. Through additional manual curation of the reported interactions, we converted the graph into a set of Boolean rules, following the procedure of the first Boolean model (62 nodes) for melanogenesis published by Lee et al. Simulations show that the predicted response to varying UV levels for most of the nodes is similar to the predictions of the existing model. The greater complexity allows investigation of the sensitivity of melanin to additional nodes. We carried out perturbation analysis of the network through node deletion and constitutive activation to identify sensitivity of outcomes, and compared the nodes identified as sensitive to previous reports.

Method for identification of sensitive nodes in Boolean models of biological networks.

Biological systems are often represented as Boolean networks and analysed to identify sensitive nodes which on perturbation disproportionately change a predefined output. There exist different kinds of perturbation methods: perturbation of function, perturbation of state and perturbation in update scheme. Nodes may have defects in interpretation of the inputs from other nodes and calculation of the node output. To simulate these defects and systematically assess their effect on the system output, two new function perturbations, referred to as 'not of function' and 'function of not', are introduced. In the former, the inputs are assumed to be correctly interpreted but the output of the update rule is perturbed; and in the latter, each input is perturbed but the correct update rule is applied. These and previously used perturbation methods were applied to two existing Boolean models, namely the human melanogenesis signalling network and the fly segment polarity network. Through mathematical simulations, it was found that these methods successfully identified nodes earlier found to be sensitive using other methods, and were also able to identify sensitive nodes which were previously unreported.