Concentrations of glucose and total reducing sugars obtained by chemical hydrolysis of three different lignocellulosic feedstocks were maximized. Two response surface methodologies were applied to model the amount of sugars produced: (1) classical quadratic least-squares fit (QLS), and (2) artificial neural networks based on radial basis functions (RBF). The results obtained by applying RBF were more reliable and better statistical parameters were obtained. Depending on the type of biomass, different results were obtained. Improvements in fit between 35% and 55% were obtained when comparing the coefficients of determination (R2) computed for both QLS and RBF methods. Coupling the obtained RBF models with particle swarm optimization to calculate the global desirability function, allowed to perform multiple response optimization.