A landslide susceptibility assessment is one of the critical steps in planning for landslide disaster prevention. Advanced machine learning methods can be used as data-driven approaches for landslide susceptibility zonation with several landslide conditioning factors. Despite there being a number of studies on landslide susceptibility assessment, the literature is limited in several contexts, such as parameter optimization, an examination of the factors in detail, and study area. This study addresses these lacks in the literature and aims to develop a landslide susceptibility map of Kentucky, US. Four machine learning methods, namely artificial neural network (ANN), k-nearest neighbor (KNN), support vector machine (SVM), and stochastic gradient boosting (SGB), were used to train the dataset comprising sixteen landslide conditioning factors after pre-processing the data in terms of data encoding, data scaling, and dimension reduction. The hyperparameters of the machine learning methods were optimized using a state-of-the-art artificial bee colony (ABC) algorithm. The permutation importance and Shapley additive explanations (SHAP) methods were employed to reduce the dimension of the dataset and examine the contributions of each landslide conditioning factor to the output variable, respectively. The findings show that the ABC-SGB hybrid model achieved the highest prediction performance. The SHAP summary plot developed using the ABC-SGB model shows that intense precipitation, distance to faults, and slope were the most significant factors affecting landslide susceptibility. The SHAP analysis further underlines that increases in intense precipitation, distance to faults, and slope are associated with an increase in the probability of landslide incidents. The findings attained in this study can be used by decision makers to develop the most effective resource allocation plan for preventing landslides and minimizing related damages.