An Investigation on the Vibrational Analysis, Molecular Structure and Interaction With the Cellulase Enzymes of the Cellulose I Alpha and Cellulose I Beta: Dft Calculations, Docking Simulations

Abstract

Molecular docking is one of the most widely used techniques for simulating molecular interactions between molecules and forecasting the mode of binding and affinities between them. Due to the presence of structure-function relationship, in this study firstly, the molecular structures of the cellulose I(alpha) and I(beta) molecules were optimized and their most stable structures were determined by density functional theory (DFT) using B3LYP method with 6-31G(d,p) basis set. The vibrational wavenumbers of 1-ring, 2-ring, 3-ring, and 4-ring structures of cellulose I(alpha) and I(beta) were calculated using the same level of theory. Reliable vibrational assignments were made based on potential energy distribution (PED %) of the vibrational modes. The energy gap (Eg = ELUMO - EHOMO) of the cellulose Iα and cellulose Iβ was found to be 8.286 eV and 7.965 eV, respectively. To identify the molecular interactions between cellulose Iα and Iβ ligands and the cellulase enzymes, molecular docking studies were performed. The molecular docking results revealed the strong interaction of the cellulose Iα and Iβ with Endoglucanase enzyme (-6.4 and -6.3 kcal/mol, respectively), enzyme β-glucosidase (-5.3 and -5.2 kcal/mol, respectively), and Exoglucanase enzyme (-6.1 and -6.2 kcal/mol, respectively).