Mutations in the interface of membrane protein (MP) complexes are key contributors to a broad spectrum of human diseases, primarily due to changes in their binding affinities. While various methods exist for predicting the mutation-induced changes in binding affinity (ΔΔG) in protein–protein complexes, none are specific to MP complexes. This study proposes a novel strategy for ΔΔG prediction in MP complexes, which combines linear and nonlinear models, to obtain a more robust model with improved prediction accuracy. We used multiple linear regression to extract informative features that influence the binding affinity in MP complexes, which included changes in the stability of the complex, conservation score, electrostatic interaction, relatively accessible surface area, and interface contacts. Further, using gradient boosting regressor on the selected features, we developed MPA-MutPred, a novel method specific for predicting the ΔΔG of membrane protein–protein complexes, and it is freely accessible at https://web.iitm.ac.in/bioinfo2/MPA-MutPred/. Our method achieved a correlation of 0.75 and a mean absolute error (MAE) of 0.73 kcal/mol in the jack-knife test conducted on a dataset of 770 mutants. We further validated the method using a blind test set of 86 mutations, obtaining a correlation of 0.85 and an MAE of 0.77 kcal/mol. We anticipate that this method can be used for large-scale studies to understand the influence of binding affinity change on disease-causing mutations in MP complexes, thereby aiding in the understanding of disease mechanisms and the identification of potential therapeutic targets.