Abstract
Beyond comparative genomics, we identified 85 sugar transporter genes in Cordyceps militaris, clustering into nine subfamilies as sequence- and phylogenetic-based functional classification, presuming the versatile capability of the fungal growths on a range of sugars. Further analysis of the global gene expression patterns of C. militaris showed 123 genes were significantly expressed across the sucrose, glucose, and xylose cultures. The sugar transporters specific for pentose were then identified by gene-set enrichment analysis. Of them, the putative pentose transporter, CCM_06358 gene, was highest expressed in the xylose culture, and its functional role in xylose transport was discovered by the analysis of conserved structural motifs. In addition, a battery of molecular modeling methods, including homology modeling, transport pathway analysis, residue interaction network combined with molecular mechanics Poisson-Boltzmann surface area simulation (MM-PBSA), was implemented for probing the structure and function of the selected pentose transporter (CCM_06358) as a representative of sugar transportome in C. militaris. Considering the network bottlenecks and structural organizations, we further identified key amino acids (Phe38 and Trp441) and their interactions with other residues, contributing the xylose transport function, as verified by binding free energy calculation. The strategy used herein generated remarkably valuable biological information, which is applicable for the study of sugar transportome and the structure engineering of targeted transporter proteins that might link to the production of bioactive compounds derived from xylose metabolism, such as cordycepin.