The final step in lignocellulose enzymatic saccharification is the cellobiose conversion to glucose by beta-glucosidases (BG). In this work, a valid kinetic model to describe cellobiose degradation for an industrial mixture of BG enzymes present in Aspergillus fumigatus is selected. Firstly, the enzyme mixture was characterised in terms of protein content and enzymatic activity on p-NPG (1326 U mL(preparation)(-1)), determining the molecular weight of the only BG activity band observed in zymograms by SDS-PAGE and MALDI-TOF: 94 kDa. Subsequently, to select the correct kinetic model for the enzymatic hydrolysis of cellobiose, a combined strategy was performed: Firstly, non-linear regressions were applied to initial hydrolysis rate data for different enzyme concentrations and initial substrate and product concentrations, observing inhibition by cellobiose and glucose. Secondly, the optimal kinetic model was discriminated by a coupled non-linear regression-DOE numerical integration approach, by fitting several possible kinetic models involving different product inhibition mechanisms to progress curve data from runs at various initial substrate concentrations and temperatures. The best kinetic model involves non-competitive substrate inhibition and product competitive inhibition with two binding sites for glucose. (C) 2018 Institution of Chemical Engineers. Published by Elsevier B.V. All rights reserved.