Soil texture and particle size fractions (PSFs) are a critical characteristic of soil that influences most physical, chemical, and biological properties of soil; furthermore, reliable spatial predictions of PSFs are crucial for agro-ecological modeling. Here, series of hybridized artificial neural network (ANN) models with bio-inspired metaheuristic optimization algorithms such as a genetic algorithm (GA-ANN), particle swarm optimization (PSO-ANN), bat (BAT-ANN), and monarch butterfly optimization (MBO-ANN) algorithms, were built for predicting PSFs for the Mazandaran Province of northern Iran. In total, 1595 composite surficial soil samples were collected, and 64 environmental covariates derived from terrain, climatic, remotely sensed, and categorical datasets were used as predictors. Models were tested using a repeated 10-fold nested cross-validation approach. The results indicate that the hybridized ANN methods were far superior to the reference approach using ANN with a backpropagation training algorithm (BP-ANN). Furthermore, the MBO-ANN approach was consistently determined to be the best approach and yielded the lowest error and uncertainty. The MBO-ANN model improved the predictions in terms of RMSE by 20% for clay, 10% for silt, and 24% for sand when compared to BP-ANN. The physiographical units, soil types, geology maps, rainfall, and temperature were the most important predictors of PSFs, followed by the terrain and remotely sensed data. This study demonstrates the effectiveness of bio-inspired algorithms for improving ANN models. The outputs of this study will support and inform sustainable soil management practices, agro-ecological modeling, and hydrological modeling for the Mazandaran Province of Iran.