Enzymatic conversion of lignocellulosic biomass into biofuels and chemicals constitutes a potential route for sustainable development. Cellobiohydrolases are key enzymes used in industrial cocktails for depolymerization of crystalline Cellulose, and their mechanism of action has been intensely studied in the past several years. Provided with a tunnel-like substrate-binding cavity, cellobiohydrolases possess the ability to processively hydrolyze glycosidic bonds of crystalline cellulose, yielding one molecule of cellobiose per catalytic cycle. As such, cellobiose expulsion from the product binding site is a necessary step in order to allow for the processive hydrolysis mechanism. However, the high-affinity binding of cellobiose to the enzyme impairs the process and causes activity inhibition due to reaction products. Here, we use molecular dynamics simulations to study the binding of cellobiose to the Trichoderma reesei Cel7A (TrCel7A) cellobiohydrolase and the effects of mutations that reduce cellobiose binding, without affecting the structural and dynamical integrities of the enzyme. We observe that the product binding site exhibits an intrinsic flexibility that can sterically hinder cellobiose release. Several point mutations in the product binding site reduce cellobiose-enzyme interactions, but not all modifications are able to maintain the structural integrity of the enzyme. In particular, imitation of charged residues in the TrCel7A product binding site causes perturbations that affect the structure of the loops that form the substrate-binding tunnel of the enzyme and, hence, may affect TrCel7A function in other steps of the hydrolysis mechanism. Our results suggest there is a trade-off between product inhibition and catalytic efficiency, and they provide directions for cellulases engineering.