Binding Process and Free Energy Characteristics of Cellulose Chain into the Catalytic Domain of Cellobiohydrolase TrCel7A

It was observed in experiments that the catalytic domain (CD) of Trichoderma reesei Cel7A (TrCel7A) hydrolyzes crystalline cellulose in a processive manner, but the underlying binding mechanism is still unknown. Here, through replica-exchange molecular dynamics simulations, we find that the loading and sucking-in process of the cellulose chain into CD is entropy-driven and enthalpy-unfavorable, which firmly relate to the desolvation of the binding channel of CD. During the loading process, hydrophobic interactions play a dominant role because several aromatic residues have been identified to guide the cellulose chain processing. At the active site, a transition from enthalpy- to entropy-driven is detected for the driving force. Such a finding reveals the indispensability of the catalytic reaction of the glycosidic bond to provide the energy to drive the movements of the cellulose chain. Our study reveals the interaction pictures between the cellulose chain and TrCel7A at the atomic level, which helps better understand the catalytic mechanism of TrCel7A.