This study used machine learning models and SHAP explanations to identify tear proteins associated with the severity of meibomian gland dysfunction (MGD) in human patients. The central result is that several proteins, including S100-A8, PRP4, and lysozyme C, were identified as important features for classifying MGD levels, with some showing statistically significant differences between severity groups.