In vitro and in silico characterization of alkaline serine protease from Bacillus subtilis D9 recovered from Saudi Arabia

In this study, we have isolated and characterized proteolytic soil bacteria and their alkaline protease. Based on 16S rRNA sequence analysis, 12 isolates with the highest protease activity were classified as B. subtilis and B. cereus groups. B. subtilis D9 isolate showing the highest protease activity was selected for in vitro and in silico analysis for its ِِAKD9 protease. The enzyme has a molecular mass of 48 kDa, exhibiting optimal activity at 50 °C pH 9.5, and showed high stability till 65 °C and pH 8–11 for 1 h. Fe(3+)‏ stimulated, but Zn(2+) and Hg(2+) strongly inhibited the protease activity. Also, the maximum inhibition with PMSF indicated serine protease-type of AKD9 protease. AkD9 alkaline serine protease gene showed high sequence similarity and close phylogenetic relationship with AprX serine protease of B. subtilis isolates. Functional prediction of AKD9 resulted in the detection of subtilase domain, peptidase_S8 family, and subtilase active sites. Moreover, prediction of physicochemical properties indicated that AKD9 serine protease is hydrophilic, thermostable, and alkali-halo stable. Secondary structure prediction revealed the dominance of the coils enhances AKD9 activity and stability under saline and alkaline conditions. Based on molecular docking, AKD9 showed very promising binding affinities towards casein substrate with expected intrinsic proteolytic activities matching our obtained in vitro results. In conclusion, AKD9 alkaline serine protease seems to be a significant candidate for industrial applications because of its stability, hydrophilicity, enhanced thermostability, and alkali-halo stability.