Parkinson’s disease (PD) exhibits heterogeneity in terms of symptoms and prognosis, likely due to diverse neuroanatomical alterations. This study employs a contrastive deep learning approach to analyze Magnetic Resonance Imaging (MRI) data from 932 PD patients and 366 controls, aiming to disentangle PD-specific neuroanatomical alterations. The results reveal that these neuroanatomical alterations in PD are correlated with individual differences in dopamine transporter binding deficit, neurodegeneration biomarkers, and clinical severity and progression. The correlation with clinical severity is verified in an external cohort. Notably, certain proteins in the cerebrospinal fluid are strongly associated with PD-specific features, particularly those involved in the immune function. The most notable neuroanatomical alterations are observed in both subcortical and temporal regions. Our findings provide deeper insights into the patterns of brain atrophy in PD and potential underlying molecular mechanisms, paving the way for earlier patient stratification and the development of treatments to slow down neurodegeneration.