Objective: The motion trajectory prediction (MTP) based brain-computer interface (BCI) leverages electroencephalography (EEG) signals to reconstruct the three-dimensional trajectory of upper limb motion, which is pivotal for the advancement of prosthetic devices that can assist motor-disabled individuals. Most research focused on improving the performance of regression models while neglecting the correlation between the implicit information extracted from EEG features across various frequency bands with limb kinematics. Current work aims to identify key channels that capture information related to various motion execution movements from different frequency bands and reconstruct three-dimensional motion trajectories based on EEG features. Methods: We propose an interpretable motion trajectory regression framework that extracts bandpower features from different frequency bands and concatenates them into multi-band fusion features. The extreme gradient boosting regression model with Bayesian optimization and Shapley additive explanation methods are introduced to provide further explanation. Results: The experimental results demonstrate that the proposed method achieves a mean Pearson correlation coefficient (PCC) value of 0.452, outperforming traditional regression models. Conclusion: Our findings reveal that the contralateral side contributes the most to motion trajectory regression than the ipsilateral side which improves the clarity and interpretability of the motion trajectory regression model. Specifically, the feature from channel C5 in the Mu band is crucial for the movement of the right hand, while the feature from channel C3 in the Beta band plays a vital role. Significance: This work provides a novel perspective on the comprehensive study of movement disorders.