The importance of finding or making thermostable enzymes in different industries have been highlighted. Therefore, it is inevitable to understand the features involving in enzymes' thermostability. Different approaches have been employed to extract or manufacture thermostable enzymes. Here we have looked at features contributing to Endo-1,4,βxylanase (EC 3.2.1.8) thermostability, the key enzyme with possible applications in waste treatment, fuel and chemical production and paper industries. We trained different neural networks with/without feature selection and classification modelling on all available xylanase enzymes amino acids sequences to find features contributing to enzyme thermal stability. Frequency of Met (-0.006) and Lys (-0.010) showed the weakest correlation with xylanase enzymes' optimum temperature; the count of Lie (0.326) and Glu (0.324) showed the strongest direct correlation while the count of oxygen (-0.38) and frequency of Gln (-0.299) reversely correlated to xylanase enzyme thermostability Six modelling methods (Quick, Dynamic, Multiple, Prune, Exhaustive Prune and RBFN) applied on all available xylanase sequences with/without validation set and/or feature selection (24 neural networks); with estimated accuracy between 80% to 90%; the best one (90.638%) in Multiple method of neural network without validation set and without feature selection, exactly in the most complicated neural network. The weakest accuracy (80.560%) found in Dynamic method of neural network without feature selection and with validation set. In 6 out of 24 neural networks generated here, the frequency of Gln was the most important feature contributing to optimum xylanase temperature and in 4 networks count of other charged residues were the most important features. Considering the analytical and performance evaluation of different models examined here, we found Multiple model generated in modelling without feature selection and validation set a good candidate to use for testing thermostability in 7030 virtually generated Bacillus halodurans mutants. We applied this model Manuscript