It is becoming increasingly clear that brain network organization shapes the course and expression of neurodegenerative diseases. Parkinson disease (PD) is marked by progressive spread of atrophy from the midbrain to subcortical structures and, eventually, to the cerebral cortex. Recent discoveries suggest that the neurodegenerative process involves the misfolding and prion-like propagation of endogenous α-synuclein via axonal projections. However, the mechanisms that translate local "synucleinopathy" to large-scale network dysfunction and atrophy remain unknown. Here, we use an agent-based epidemic spreading model to integrate structural connectivity, functional connectivity, and gene expression and to predict sequential volume loss due to neurodegeneration. The dynamic model replicates the spatial and temporal patterning of empirical atrophy in PD and implicates the substantia nigra as the disease epicenter. We reveal a significant role for both connectome topology and geometry in shaping the distribution of atrophy. The model also demonstrates that SNCA and GBA transcription influence α-synuclein concentration and local regional vulnerability. Functional coactivation further amplifies the course set by connectome architecture and gene expression. Altogether, these results support the theory that the progression of PD is a multifactorial process that depends on both cell-to-cell spreading of misfolded proteins and regional vulnerability.