The application of the response surface methodology and the central composite design (CCD) technique for modeling and optimization of the influence of some operating variables on copper, molybdenum and rhenium recoveries in a bioleaching process was investigated. Three main bioleaching parameters, namely pH, solid concentration and inoculum percent, were changed during the bioleaching tests based on CCD. The ranges of the bioleaching process variables used in the design were as follows: pH 1.46-2.14, solid concentration 0.95%-11.05%, and inoculum percent 1.59%-18.41%. A total of 20 bioleaching tests were carried out by the CCD method according to software-based designed matrix. Empirical model equations were developed according to the copper, molybdenum and rhenium recoveries obtained with these three parameters. Model equations of responses at the base of parameters were achieved by using statistical software. The model equations were then individually optimized by using quadratic programming to maximize copper, molybdenum and rhenium recoveries individually within the experimental range. The optimum conditions for copper recovery were pH 1.68, solid concentration 0.95% and the inoculum 18.41% (v/v), while molybdenum and rhenium recoveries were 2.18% and 24.41%, respectively. The predicted values for copper, molybdenum and rhenium recoveries were found to be in good agreement with the experimental values. Also jarosite formation during bioleaching tests was also investigated.