Cutibacterium acnes, one of the most abundant skin microbes found in the sebaceous gland, is known to contribute to the development of acne vulgaris when its strains become imbalanced. The current limitations of acne treatment using antibiotics have caused an urgent need to develop a systematic strategy for selectively targeting C. acnes, which can be achieved by characterizing their cellular behaviors under various skin environments. To this end, we developed a genome-scale metabolic model (GEM) of virulent C. acnes, iCA843, based on the genome information of a relevant strain from ribotype 5 to comprehensively understand the pathogenic traits of C. acnes in the skin environment. We validated the model qualitatively by demonstrating its accuracy prediction of propionate and acetate production patterns, which were consistent with experimental observations. Additionally, we identified unique biosynthetic pathways for short-chain fatty acids in C. acnes compared to other GEMs of acne-inducing skin pathogens. By conducting constraint-based flux analysis under endogenous carbon sources in human skin, we discovered that the Wood-Werkman cycle is highly activated under acnes-associated skin condition for the regeneration of NAD, resulting in enhanced propionate production. Finally, we proposed potential anti-C. acnes targets by using the model-guided systematic framework based on gene essentiality analysis and protein sequence similarity search with abundant skin microbiome taxa.