SARS-CoV-2 mRNA vaccines are highly immunogenic in people living with HIV (PLWH) on effective antiretroviral therapy (ART). However, whether viro-immunologic parameters or other factors affect immune responses to vaccination is debated. This study aimed to develop a machine learning-based model able to predict the humoral response to mRNA vaccines in PLWH and to assess the impact of demographic and clinical variables on antibody production over time. Different machine learning algorithms have been compared in the setting of a longitudinal observational study involving 497 PLWH, after primary and booster SARS-CoV-2 mRNA vaccination. Both Generalized Linear Models and non-linear Models (Tree Regression and Random Forest) were trained and tested. Non-linear algorithms showed better ability to predict vaccine-elicited humoral responses. The best-performing Random Forest model identified a few variables as more influential, within 39 clinical, demographic, and immunological factors. In particular, previous SARS-CoV-2 infection, BMI, CD4 T-cell count and CD4/CD8 ratio were positively associated with the primary cycle immunogenicity, yet their predictive value diminished with the administration of booster doses. In the present work we have built a non-linear Random Forest model capable of accurately predicting humoral responses to SARS-CoV-2 mRNA vaccination, and identifying relevant factors that influence the vaccine response in PLWH. In clinical contexts, the application of this model provides promising opportunities for predicting individual vaccine responses, thus facilitating the development of vaccination strategies tailored for PLWH.