Computational fluid dynamics (CFD) is a powerful tool to estimate accurately the aerodynamic loads on wind turbine blades at the expense of high requirements like the duration of computation. Such requirements grow in the case of blade shape optimization in which several analyses are needed. A fast and reliable way to mimic the CFD solutions is to use surrogate models. In this study, a machine learning technique, the support vector regression (SVR) method based on a set of CFD solutions, is used as the surrogate model. CFD solutions are calculated by solving the Reynolds-averaged Navier-Stokes equation with the k-epsilon turbulence model using a commercial solver. The support vector regression model is then trained to give a functional relationship between the spanwise twist distribution and the generated torque. The smooth twist distribution is defined using a three-node cubic spline with four parameters in total. The optimum twist is determined for two baseline blade cases: the National Renewable Energy Laboratory (NREL) Phase II and Phase VI rotor blades. In the optimization process, extremum points that give the maximum torque are easily determined since the SVR gives an analytical model. Results show that it is possible to increase the torque generated by the NREL VI blade more than 10% just by redistributing the spanwise twist without carrying out a full geometry optimization of the blade shape with many shape-defining parameters. The increase in torque for the NREL II case is much higher.