This study introduces a novel approach to predict the shear bearing capacity of FRP-concrete interfaces using explainable machine learning. Eight algorithms are employed: three standalone models (Artificial Neural Network, Support Vector Regression, and Decision Tree) and five ensemble learning models (Bagging, Random Forest, Adaptive Boosting, Gradient Boosting, and Extreme Gradient Boosting). Four scenarios with varying input features, including engineered features inspired by mechanics-based bearing capacity equations, are examined. Notably, the inclusion of engineered features such as the stiffness of the FRP strip (Kf) significantly enhanced prediction accuracy and efficiency, although the width correction coefficient (bf/bc) did not yield significant benefits, contrary to findings in mechanics-based models. The Extreme Gradient Boosting (XGBoost) algorithm emerged as the top-performing model, achieving an impressive R-square value of 0.949. In terms of explainability, the study utilized the Shapley Additive Explanation (SHAP) technique to comprehensively elucidate the significance, dependency, and interaction effects of features within the best-performing model. Remarkable, approximately 75 % of the total mean absolute SHAP values were attributed to Kf and bf of the FRP, highlighting their pivotal role in the prediction process. This detailed explanation of the model using SHAP instills trust in the predictions and facilitates its practical implementation in real-world scenarios.