The aim of the present investigation was to optimize the variables of organic solvent-stable protease fermentation by response surface methodology. Plackett-Burman design used in the first step was an efficient approach to screenfsignificant variables. Then the method of steepest ascent was employed to approach or near the experimental design space. Box-Behnken design and response surface analysis used in the last steps were useful to determine