Three-dimensional models of exoinulinase from Bacillus stearothermophilus and endoinulinase from Aspergillus niger were built up by means of homology modeling. The crystal structure of exoinulinase from Aspergillus awamori was used as a template, which is the sole structure of inulinase resolved so far. Docking and molecular dynamics simulations were performed to investigate the differences between the two inulinases in terms of substrate selectivity. The analysis of the structural differences between the two inulinases provided the basis for the explanation of their different regio-selectivity and for the understanding of enzyme-substrate interactions. Surface analysis was performed to point out structural features that can affect the efficiency of enzymes also after immobilization. The computational analysis of the three-dimensional models proved to be an effective tool for acquiring information and allowed to formulate an optimal immobilized biocatalyst even more active that the native one, thus enabling the full exploitation of the catalytic potential of these enzymes.