Managing Parkinson's disease (PD) through medication can be challenging due to varying symptoms and disease duration. This study aims to demonstrate the potential of sequence-by-sequence algorithms in recommending personalized medication combinations for patients with PD based on their previous visits. Our proposed method employs a gated recurrent unit model to predict accurate combinations of critical medication types for PD based on each patient's motor symptoms and prescribed medication from previous visits. We built a multi-label model with gated recurrent units on two data architectures: (1) personalized input using each patient's previous visits as a sample and (2) non-personalized input treating each visit as an independent sample. The 10-fold cross-validation results showed that the personalized architecture model outperforms the non-personalized model in accuracy (0.92), precision (0.94), recall (0.94), F1-score (0.94), Hamming loss (0.03), and macro average area under the receiver operating characteristic (0.94). To interpret the model's predictions, we employed SHapley Additive exPlanations (SHAP) values, which provide insights into the importance of variables both globally (across the entire model) and at the individual patient level. The results contribute to the sequential-based decision support system potentially enhancing the remote management of PD pharmacologic issues.