Parki so 's disease PD has ee reported to i ol e postga glio i s patheti failure and a wide spectrum of autonomic dysfunctions including cadiovascular, sexual, bladder, gastrointestinal and sudo-motor abnormalities. While these symptoms may have a significant impact on daily activities, as well as quality of life, the evaluation of autonomic nervous system (ANS) dysfunctions relies on a large and expensive battery of autonomic tests only accessible in highly specialized laboratories. In this paper we aim to devise a comprehensive computational assessment of disease-related heartbeat dynamics based on instantaneous, time-varying estimates of spontaneous (resting state) cardiovascular oscillations in PD. To this end, we combine standard ANS-related heart rate variability (HRV) metrics with measures of instantaneous complexity (dominant Lyapunov exponent and entropy) and higher-order statistics (bispectra). Such measures are computed over 600-s recordings acquired at rest in 29 healthy subjects and 30 PD patients. The only significant group-wise differences were found in the variability of the dominant Lyapunov exponent. Also, the best PD vs. healthy controls classification performance (balanced accuracy: 73.47%) was achieved only when retaining the time-varying, non-stationary structure of the dynamical features, whereas classification performance dropped significantly (balanced accuracy: 61.91%) when excluding variability-related features. Additionally, both linear and nonlinear model features correlated with both clinical and neuropsychological assessments of the considered patient population.Our results demonstrate the added value and potential of instantaneous measures of heartbeat dynamics and its variability in characterizing PD-related disabilities in motor and cognitive domains.3