Hundreds of NPC1 variants cause highly heterogeneous phenotypes. This study aims to explore the genotype–phenotype correlation of NPC1, especially for missense variants. In a well‐characterized cohort, phenotypes are graded into three clinical forms: mild, intermediate, and severe. Missense residue structural location was stratified into three categories: surface, partially, and fully buried. The association of phenotypes with the topography of the amino acid substitution in the protein structure was investigated in our cohort and validated in two reported cohorts. One hundred six unrelated NPC1 patients were enrolled. A significant correlation of genotype–phenotype was found in 81 classified individuals with two or one (the second was null variant) missense variant (p < 0.001): of 25 patients with at least one missense variant of surface (group A), 19 (76%) mild, six (24%) intermediate, and none severe; of 31 cases with at least one missense variant of partially buried without surface variants (group B), 11 (35%) mild, 16 (52%) intermediate, and four (13%) severe; of the remaining 25 patients with two or one buried missense variants (group C), eight (32%) mild, nine (36%) intermediate, and eight (32%) severe. Additionally, 7‐ketocholesterol, the biomarker, was lower in group A than in group B (p = 0.024) and group C (p = 0.029). A model was proposed that accurately predicted phenotypes of 72 of 90 (80%), 73 of85 (86%), and 64 of 69 (93%) patients in our cohort, Italian, and UK cohort, respectively. This study proposed a novel genotype–phenotype correlation in NPC1, linking the underlying molecular pathophysiology with clinical phenotype and aiding genetic counseling and evaluation in clinical practice.