A robust Quantitative Structure-Property Relationship (QSPR) model was presented to predict the surface tension property of psychoanaleptic (psychostimulant and antidepressant) drugs. A dataset of 112 molecules was utilized, and three feature selection methods were applied: genetic algorithm combined with Ordinary Least Squares (GA-OLS), Partial Least Squares (GA-PLS), and Support Vector Machines (GA-SVM), each identifying ten pertinent AlvaDesc descriptors. The models were constructed using the Dragonfly Algorithm combined with the Support Vector Regressor (DA-SVR), with the GA-SVM-based model emerging as the top performer. Rigorous statistical metrics validate its superior predictive accuracy (R2 = 0.98142, Q2LOO = 0.98142, RMSE = 1.12836, AARD = 0.78746). Furthermore, an external test set of ten compounds was employed for model validation and extrapolation, along with assessing the applicability domain, further underscoring the model's reliability. The selected descriptors (X0Av, VE1sign_B(e), ATSC1e, MATS6v, P_VSA_ppp_A, TDB01u, E1s, R2m+, N-067, SssO) collectively elucidate the key structural factors influencing surface tension in the studied drugs. The model provides excellent predictions and can be used to determine the surface tension of new psychoanaleptic drugs. Its outcomes will guide the design of novel medications with targeted surface tension properties.