An isolated gene from Neosartorya fischeri NRRL181 encoding a β-glucosidase (BGL) was cloned, and its nucleotide sequence was determined. DNA sequence analysis revealed an open reading frame of 1,467 bp, capable of encoding a polypeptide of 488 amino acid residues. The gene was over-expressed in Escherichia coli, and the protein was purified using nickel-nitrilotriacetic acid chromatography. The purified recombinant BGL showed a high level of catalytic activity, with V (max) of 886 μmol min(-1) mg-protein(-1) and a K (m) of 68 mM for p-nitrophenyl-β-D: -glucopyranoside (pNPG). The optimal temperature for enzyme activity was about 40°C, and the optimal pH was about 6.0. A homology model of N. fischeri BGL1 was constructed based on the X-ray crystal structure of Phanerochaete chrysosporium BGLA. Molecular dynamics simulation studies of the enzyme with the pNPG and cellobiose shed light on the unique substrate specificity of N. fischeri BGL1 only towards pNPG.