Investigation of an “alternate water supply system” in enzymatic hydrolysis in the processive endocellulase Cel7A from Rasamsonia emersonii by molecular dynamics simulation

Cel7A from Rasamsonia emersonii is one of the processive endocellulases classified under family 7 glycoside hydrolase. Molecular dynamics simulations were carried out to obtain the optimized sliding and hydrolyzing conformations, in which the reducing ends of sugar chains are located on different sites. Hydrogen bonds are investigated to clarify the interactions between protein and substrate in either conformation. Nine hydrogen bonding interactions are identified in the sliding conformation, and six similar interactions are also found correspondingly in the hydrolyzing conformation. In addition, four strong hydrophobic interactions are also determined. The domain cross‐correlation map analysis shows movement correlation of protein including autocorrelation between residues. The root mean square fluctuations analysis represents the various flexibilities of different fragment in the two conformations. Comparing the two conformations reveals the water‐supply mechanism of selective hydrolysis of cellulose in Cel7A. The mechanism can be described as follow. When the reducing end of substrate slides from the unhydrolyzing site (sliding conformation) to the hydrolyzing site (hydrolyzing conformation), His225 is pushed down and rotated, the rotation leads to the movement of Glu209 with the interstrand hydrogen bonding in β‐sheet. It further makes Asp211 close to the hydrolysis center and provides a water molecule bounding on its carboxyl in the previous unhydrolyzing site. After the hydrolysis takes place and the product is excluded from the enzyme, the Asp211 comes back to its initial position. In summary, Asp211 acts as an elevator to transport outer water molecules into the hydrolysis site for every other glycosidic bond.