The evolutionary variability of a protein's residues is highly dependent on protein region and protein function. Solvent-exposed residues, excluding those at interaction interfaces, are more variable than buried residues. Active site residues are considered to be conserved as they ensure an enzyme's activity and selectivity. The abovementioned rules apply also to α/β-hydrolase fold proteins - an example of enzymes with buried active sites equipped with tunnels linking the reaction site with the exterior. We hypothesised two scenarios: (1) tunnels are lined by mostly variable residues, allowing adaptation to the evolutionary pressures of a changeable environment; or (2) tunnels are lined by mostly conserved amino acids, and are equipped with a number of specific variable residues that are able to respond to evolutionary pressure. We also wanted to check if evolutionary analysis can help distinguish functional and non-functional tunnels. Soluble epoxide hydrolases (sEHs) represent a good case study for the analysis of the evolution of tunnels in an α/β-hydrolase fold family due to their size and architecture. Here, we propose methods for the comparison of tunnels detected in both crystal structures and molecular dynamics simulations, as well as the assignment of tunnel functionality, and we identify critical steps for careful tunnel inspection. We also compare the entropy values of the tunnel-lining residues and system-specific compartments in seven selected sEHs from different clades. We present three different cases of entropy distribution among tunnel-lining residues. As a result, we propose a 'perforation' model for tunnel evolution via the merging of internal cavities or surface perforations. We also report an approach for the identification of highly variable tunnel-lining residues as potential targets to be used for the fine-tuning of selected enzymes.