Trichoderma virens β-glucosidase I (BGLI) gene; expression in Saccharomyces cerevisiae including docking and molecular dynamics studies

Gammadde Hewa Ishan Maduka Wickramasinghe, Pilimathalawe Panditharathna Attanayake Mudiyanselage Samith Indika Rathnayake, Naduviladath Vishvanath Chandrasekharan, Mahindagoda Siril Samantha Weerasinghe, Ravindra Lakshman Chundananda Wijesundera, Wijepurage Sandhya Sulochana Wijesundera

Cellulose, a linear polymer of β 1-4, linked glucose, is the most abundant renewable fraction of plant biomass (lignocellulose). It is synergistically converted to glucose by endoglucanase (EG) cellobiohydrolase (CBH) and β-glucosidase (BGL) of the cellulase complex. BGL plays a major role in the conversion of randomly cleaved cellooligosaccharides into glucose. As it is well known, Saccharomyces cerevisiae can efficiently convert glucose into ethanol under anaerobic conditions. Therefore, S.cerevisiae was genetically modified with the objective of heterologous extracellular expression of the BGLI gene of Trichoderma virens making it capable of utilizing cellobiose to produce ethanol. The cDNA and a genomic sequence of the BGLI gene of Trichoderma virens was cloned in the yeast expression vector pGAPZα and separately transformed to Saccharomyces cerevisiae. The size of the BGLI cDNA clone was 1363 bp and the genomic DNA clone contained an additional 76 bp single intron following the first exon. The gene was 90% similar to the DNA sequence and 99% similar to the deduced amino acid sequence of 1,4-β-D-glucosidase of T. atroviride (AC237343.1). The BGLI activity expressed by the recombinant genomic clone was 3.4 times greater (1.7 x 10(-3) IU ml(-1)) than that observed for the cDNA clone (5 x 10(-4) IU ml(-1)). Furthermore, the activity was similar to the activity of locally isolated Trichoderma virens (1.5 x 10(-3) IU ml(-1)). The estimated size of the protein was 52 kDA. In fermentation studies, the maximum ethanol production by the genomic and the cDNA clones were 0.36 g and 0.06 g /g of cellobiose respectively. Molecular docking results indicated that the bare protein and cellobiose-protein complex behave in a similar manner with considerable stability in aqueous medium. The deduced binding site and the binding affinity of the constructed homology model appeared to be reasonable. Moreover, it was identified that the five hydrogen bonds formed between the amino acid residues of BGLI and cellobiose are mainly involved in the integrity of enzyme-substrate association. The BGLI activity was remarkably higher in the genomic DNA clone compared to the cDNA clone. Cellobiose was successfully fermented into ethanol by the recombinant S.cerevisiae genomic DNA clone. It has the potential to be used in the industrial production of ethanol as it is capable of simultaneous saccharification and fermentation of cellobiose. Homology modeling, docking studies and molecular dynamics simulation studies will provide a realistic model for further studies in the modification of active site residues which could be followed by mutation studies to improve the catalytic action of BGLI.