Background: Based on the previous studies, antibiotics can have affected biological properties of biomass and fouling properties of mixed liquor in aeration tank. The present study was conducted to explore the fouling mechanisms of membrane bioreactor (MBR) system during the treatment of wastewater containing erythromycin (ERY) antibiotic under several mixed liquor suspended solids (MLSS) concentrations. Methods: A lab-scale two-chamber MBR system equipped with a polypropylene hollow fiber submerged membrane was fed with synthetic wastewater containing different initial concentrations of ERY. MBR system was operated under the constant flux mode and different MLSS concentrations (5.0-13.0 g/L) and the obtained results were evaluated using different individual and combined fouling models. Results: The variation of MLSS concentrations had not significantly affected the kind of best-fitted model. From the individual models, the standard model indicated the best performance for permeate prediction under different MLSS concentrations (R2 adj>0.997). For all studied MLSS concentrations, the R2 adj values of combined fouling models were higher than 0.986 and demonstrated good fitness performance of combined models compared to individual models. Overall, the cake-intermediate model showed the lowest fitness, and cake-complete and complete-standard models were the most successful models in filtrated volume prediction in comparison with other combined fouling models. Conclusion: This study indicated that mechanistic models are suitable for fouling prediction of MBR systems in ERY removal and under a wide range of MLSS concentrations and provide valuable information on fouling mechanisms of full-scale MBR systems.