This study introduces machine learning (ML) models to predict CO₂ solubility in mixtures of a deep eutectic solvent (tetrabutylammonium bromide: ethylene glycol, 1:4 molar ratio) and methyldiethanolamine (MDEA). To ensure robust extrapolation capabilities and prevent potential data leakage, a rigorous grouped validation strategy, known as “Leave-One-Composition- Out,” was implemented, replacing conventional random data splitting. Four ML algorithms, Gradient Boosting (GBoost), Gaussian Process Regression (GPR), AdaBoost-SVR, and Multilayer Perceptron (MLP), were evaluated using temperature, pressure, and MDEA weight% as input features. Results revealed that the GPR model outperformed all others across all statistical metrics, achieving an R2 value of 0.996 and a mean absolute error (MAE) of 0.2839%. Furthermore, SHAP (SHapley Additive exPlanations) analysis identified pressure as the most significant driving factor for CO2 solubility, while temperature exhibited an inverse correlation. Crucially, this perfect alignment with fundamental principles clearly demonstrates that the model successfully learned the underlying thermodynamic relationships rather than merely memorizing numerical data. Finally, leverage validation confirmed that 99.16% of the data points safely reside within the model’s valid applicability domain. This work establishes the developed GPR framework as a highly reliable and physically consistent approach for predicting CO₂ solubility and optimizing complex multi-component capture systems.