The study develops an updated genome-scale metabolic model of tomato leaf (CBZ_iSL3433) incorporating carbamazepine (CBZ) transformation reactions to simulate the pollutant's stress impact under phototrophic conditions. Constraint-based analysis predicts significant reductions in biomass and altered fluxes in nutrient assimilation and secondary metabolism, which are largely ameliorated by biostimulants such as proline, spermine, glycerol, and ethanol. This framework demonstrates a systems-level approach for screening strategies to mitigate pharmaceutical pollutant stress in crops.