While deep brain stimulation (DBS) remains an effective therapy for Parkinson's disease (PD), sources of variance in patient outcomes are still not fully understood, underscoring a need for better prognostic criteria. Here we leveraged routinely collected T1-weighted (T1-w) magnetic resonance imaging (MRI) data to derive patient-specific measures of brain structure and evaluate their usefulness in predicting changes in PD medications in response to DBS. Preoperative T1-w MRI data from 231 patients with PD were used to extract regional measures of fractal dimension (FD), sensitive to the structural complexities of cortical and subcortical areas. FD was validated as a biomarker of Parkinson's disease (PD) progression through comparison of patients with PD and healthy controls (HCs). This analysis revealed significant group differences in FD across nine brain regions which supports its utility as a marker of PD. We evaluated the impact of adding imaging features (FD) to a clinical model that included demographics and clinical parameters-age, sex, total number and location of DBS electrodes, and preoperative motor response to levodopa. This model aimed to explain variance and predict changes in medication following DBS. Regression analysis revealed that inclusion of the FD of distributed brain areas correlated with post-DBS reductions in medication burden, explaining an additional 13.6% of outcome variance (R 2 =0.388) compared to clinical features alone (R 2 =0.252). Hypergraph-based classification learning tasks achieved an area under the receiver operating characteristic curve of 0.64 when predicting with clinical features alone, versus 0.76 when combining clinical and imaging features. These findings demonstrate that PD effects on brain morphology linked to disease progression influence DBS outcomes. The work also highlights FD as a potentially useful imaging biomarker to enhance DBS candidate selection criteria for optimized treatment planning.