Introduction Bipolar disorder (BD) patients present an increased risk of suicide attempts. Most current machine learning (ML) studies predicting suicide attempts are cross-sectional, do not employ time-dependent variables, and do not assess more than one modality. Therefore, we aimed to predict 12-month suicide attempts in a sample of BD patients, using clinical and brain imaging data. Methods A sample of 163 BD patients were recruited and followed up for 12 months. Gray matter volumes and cortical thickness were extracted from the T1-weighted images. Based on previous literature, we extracted 56 clinical and demographic features from digital health records. Support Vector Machine was used to differentiate BD subjects who attempted suicide. First, we explored single modality prediction (clinical features, GM, and thickness). Second, we implemented a multimodal stacking-based data fusion framework. Results During the 12 months, 6.13% of patients attempted suicide. The unimodal classifier based on clinical data reached an area under the curve (AUC) of 0.83 and balanced accuracy (BAC) of 72.7%. The model based on GM reached an AUC of 0.86 and BAC of 76.4%. The multimodal classifier (clinical + GM) reached an AUC of 0.88 and BAC of 83.4%, significantly increasing the sensitivity. The most important features were related to suicide attempts history, medications, comorbidities, and depressive polarity. In the GM model, the most relevant features mapped in the frontal, temporal, and cerebellar regions. Conclusions By combining models, we increased the detection of suicide attempts, reaching a sensitivity of 80%. Combining more than one modality proved a valid method to overcome limitations from single-modality models and increasing overall accuracy.