Problems associated with Organochlorine pesticide (OCP)-contaminated soils have received wide attention. To understand the anaerobic biodegradation process constraints, innovative mathematical analysis methods are effective. Response surface methodology (RSM) and Tenax TA extraction method combined with the first-three-compartment model were employed to systematically investigate the role of nitrate concentration and bioaccessibility enhancer (methyl-β-cyclodextrin, MCD) in the anaerobic biodegradation of OCPs in contaminated soil. The sole addition of either KNO3 or MCD could facilitate the anaerobic biodegradation of OCPs. The highest biodegradation for total OCPs, hexachlorocyclohexanes, endosulfans, and chlordanes were 71.6, 82.1, 68.3, and 55.6 %, respectively, when 20 mM KNO3 and 3.0 % (w/w) MCD were applied simultaneously. As predicted by RSM, the theoretical maximum biodegradation for total OCPs ranged from 60 to 80 % when 20 to 25 mM KNO3 and >2.5 % (w/w) MCD were applied simultaneously. Tenax TA extraction method demonstrated the enhancement of OCP bioaccessibility caused by MCD addition. Changes in the soil microbial activities also suggested the positive effects of adding suitable amounts of KNO3 as a cosubstrate to facilitate the anaerobic biodegradation of OCPs. The amount of KNO3 and MCD are crucial in influencing OCP biodegradation. RSM was demonstrated to be a powerful tool to estimate and predicting the optimal OCP biodegradation under KNO3 and MCD application simultaneously.