Xylanases are glycosyl hydrolases that catalyze the hydrolysis of internal β-1,4 glycosidic bonds of xylan, the major hemi-cellulose component of the plant cell wall. Enzymes such as xylanases are used considerably in industries. Their industrial usage is especially attractive since they can replace some of the environmental pollutants. We have earlier isolated a family 11-xylanase gene from Bacillus subtilis-AK1, which is active at high temperature as well as at alkaline pH. In order to understand the factors liable for the adaptation of this enzyme, three dimensional model of B. subtilis-AK1 xylanase has now been obtained by homology modeling. Modeling was carried out using Molecular Operating Environment (MOE) software developed by Chemical Computing Group Inc., running on Pentium IV workstation. The model showed that B. subtilis-AK1 xylanase having molecular weight around 20 kDa contains in its fold an α-helix and two β-sheets packed against each other forming a β-sandwich. The conserved active site amino acids E78R, Y80L were mutated in this novel B. subtilis-AK1 strain, but the protein folding and the function was maintained with high thermal stability. Several minor modifications appeared to be responsible for the increased thermo stability of AK1. Docking studies of the substrate xylan with -AK1 shows the possibility of the Arg 78 acting as the nucleophile instead of Glu 78.