Machine learning has been proposed to utilize D-amino acid oxidase (DAO) and DAO activator (DAOA [or pLG72]) protein levels to ascertain disease status in schizophrenia. However, it remains unclear whether machine learning can effectively evaluate clinical features in relation to DAO and DAOA in schizophrenia patients. We employed an interpretable machine learning (IML) framework including linear regression, least absolute shrinkage and selection operator (Lasso) models, and generalized additive models (GAMs) to analyze DAO/DAOA levels using 380 Taiwanese schizophrenia patients. Additionally, we incorporated 27 parameters encompassing demographic variables, clinical assessments, functional outcomes, and cognitive function as features. The IML framework facilitated linear and non-linear relationships between features and DAO/DAOA. DAO levels demonstrated significant associations with the 17-item Hamilton Depression Rating Scale (HAMD17) based on linear regression. The Lasso model identified four features-HAMD17, age, working memory, and overall cognitive function (OCF)-and highlighted HAMD17 as the most significant feature, using DAO from chronically stable patients. Utilizing DAOA from acutely exacerbated patients, the Lasso model also identified four features-OCF, Scale for the Assessment of Negative Symptoms 20-item, quality of life scale (QLS), and category fluency-and emphasized OCF as the most significant feature. Furthermore, GAMs revealed a non-linear relationship between category fluency and DAO in chronically stable patients, as well as between QLS and DAOA in acutely exacerbated patients. The study suggests that an IML framework holds promise for assessing linear and non-linear relationships between DAO/DAOA and various features in clinical assessments, functional outcomes, and cognitive function in patients with schizophrenia.