Insight into the process of product expulsion in cellobiohydrolase Cel6A from Trichoderma reesei by computational modeling

Glycoside hydrolase cellulase family 6 from Trichoderma reesei (TrCel6A) is an important cellobiohydrolase to hydrolyze cellooligosaccharide into cellobiose. The knowledge of enzymatic mechanisms is critical for improving the conversion efficiency of cellulose into ethanol or other chemicals. However, the process of product expulsion, a key component of enzymatic depolymerization, from TrCel6A has not yet been described in detail. Here, conventional molecular dynamics and steered molecular dynamics (SMD) were applied to study product expulsion from TrCel6A. Tyr103 may be a crucial residue in product expulsion given that it exhibits two different posthydrolytic conformations. In one conformation, Tyr103 rotates to open the −3 subsite. However, Tyr103 does not rotate in the other conformation. Three different routes for product expulsion were proposed on the basis of the two different conformations. The total energy barriers of the three routes were calculated through SMD simulations. The total energy barrier of product expulsion through Route 1, in which Tyr103 does not rotate, was 22.2 kcal·mol⁻¹. The total energy barriers of product expulsion through Routes 2 and 3, in which Tyr103 rotates to open the −3 subsite, were 10.3 and 14.4 kcal·mol⁻¹, respectively. Therefore, Routes 2 and 3 have lower energy barriers than Route 1, and Route 2 is the thermodynamically optimal route for product expulsion. Consequently, the rotation of Tyr103 may be crucial for product release from TrCel6A. Results of this work have potential applications in cellulase engineering.