The neurodegenerative progression of Parkinson's disease affects brain structure and function and, concomitantly, alters the topological properties of brain networks. The network alteration accompanied by motor impairment and the duration of the disease has not yet been clearly demonstrated in the disease progression. In this study, we aim to resolve this problem with a modeling approach using the reduced Jansen-Rit model applied to large-scale brain networks derived from cross-sectional MRI data. Optimizing whole-brain simulation models allows us to discover brain networks showing unexplored relationships with clinical variables. We observe that the simulated brain networks exhibit significant differences between healthy controls ( = 51) and patients with Parkinson's disease ( = 60) and strongly correlate with disease severity and disease duration of the patients. Moreover, the modeling results outperform the empirical brain networks in these clinical measures. Consequently, this study demonstrates that utilizing the simulated brain networks provides an enhanced view of network alterations in the progression of motor impairment and identifies potential biomarkers for clinical indices.