Unraveling the potential of uninvestigated thermoalkaliphilic lipases by molecular docking and molecular dynamic simulation: an in silico characterization study.

Thermoalkaliphilic lipase enzymes are mostly favored for use in the detergent industry. While there has been considerable research on Geobacillus lipases, a significant portion of these enzymes remains unexplored or undocumented in the scientific literature. This work performed in silico phylogeny, sequence alignment, structural and enzyme-substrate interaction analyses of the five thermoalkaliphilic lipases belonging to different Geobacillus species (Geobacillus stearothermophilus lipase = GsLip, Geobacillus sp. B4113_201601 lipase = Gb4Lip, Geobacillus kaustophilus HTA426 lipase = GkLip, Geobacillus sp. SP22 lipase = GspLip, Geobacillus sp. NTU 03 lipase = GntLip). For this purpose, unreviewed enzyme sequences of five Geobacillus thermoalkaliphilic lipases were analyzed at sequence and phylogeny levels. 3D homology enzyme models were built, validated, and investigated by different bioinformatics tools. The ligand interactions screening using seven para-nitrophenyl (pNP) esters and enzyme-ligand interactions were analyzed on Gb4Lip:pNP-C12 and BTL2:pNP-C12 by MD simulation. Biophysicochemical characteristic analysis showed that Gb4Lip had a theoretical T m value of above 65 ºC, and a higher aliphatic index indicating greater thermal stability. Sequence alignment showed a hydrophilic threonine in the α6 helix of Gb4Lip, indicating high enzymatic activity. A normalized temperature factor B (B'-factor) analysis showed that the lid domains of five lipases significantly possessed lower B'-factor values, compared to G. thermocatenulatus lipase 2 (BTL2), indicating that they had higher rigidity. Molecular docking results indicated that the five lipases had the highest binding affinity toward pNP-C12. The RMSF investigation revealed that the thermostability of Gb4Lip is influenced by specific molecular elements: D202-S203 within the αB region of the lid domain, and E274-Q275 within the b3 strand, as well as W278 in the b3-b4 loop, and H282 in the b4 strand of the Ca2+-binding region. MD simulation analysis showed that catalytic residue S114 and at least one oxyanion hole residue (F17 and/or Q114) in Gb4Lip frequently formed hydrogen bonds with the pNP-C12 ligand at 343 K and 348 K throughout the simulation process, indicating that Gb4Lip might catalyze relatively long-chain ligand pNP-C12 with high performance. In conclusion, Gb4Lip might be a more suitable candidate as the detergent additive. In addition, this investigation can offer valuable perspectives on Family I.5 lipases such as Gb4Lip for future exploration in the field of protein engineering. Supplementary Information: The online version contains supplementary material available at 10.1007/s13205-024-04023-5.

Revealing the role of the X25 domains through the characterization of truncated variants of amylopullulanase enzyme from Thermoanaerobacter brockii brockii.

To understand the role of the X25 domains of the amylopullulanase enzyme from Thermoanaerobacter brockii brockii (T. brockii brockii), four truncated variants that are TbbApuΔX25-1-SH3 (S130-A1484), TbbApuΔX25-2-SH3 (T235-A1484), TbbApuΔX25-1-CBM20 (S130-P1254), and TbbApuΔX25-2-CBM20 (T235-P1254) were constructed, expressed and characterized together with the SH3 and CBM20 domain truncated variants (TbbApuΔSH3 (V1-A1484) and TbbApuΔCBM20 (V1-P1254). TbbApuΔSH3 showed improved affinity and specificity for both pullulan and soluble starch than full-length TbbApu with lower Km and higher kcat/Km values. It indicates that SH3 is a disposable domain without any effect on the activity and stability of the enzyme. However, TbbApuΔX25-1-SH3, TbbApuΔX25-2-SH3, TbbApuΔX25-1-CBM20, TbbApuΔX25-2-CBM20 (T235-P1254) and TbbApuΔCBM20 showed higher Km and lower kcat/Km values than TbbApuΔSH3 to both soluble starch and pullulan. It specifies that the X25 domains and CBM20 play an important role in both α-amylase and pullulanase activity. Also, it is revealed that while truncation of the CBM20 domain as starch binding domain (SBD) did not affect on raw starch binding ability of the enzyme, truncation of both X25 domains caused almost complete loss of the raw starch binding ability of the enzyme. All these results enlightened the function of the X25 domains that play a more crucial role than CBM20 in the enzyme's binding to raw starch and also play a crucial role in its activity.

Exploring the Structural Insights of Thermostable Geobacillus esterases by Computational Characterization.

This study conducted an in silico analysis of two biochemically characterized thermostable esterases, Est2 and Est3, from Geobacillus strains. To achieve this, the amino acid sequences of Est2 and Est3 were examined to assess their biophysicochemical properties, evolutionary connections, and sequence similarities. Three-dimensional models were constructed and validated through diverse bioinformatics tools. Molecular dynamics (MD) simulation was employed on a pNP-C2 ligand to explore interactions between enzymes and ligand. Biophysicochemical property analysis indicated that aliphatic indices and theoretical T m values of enzymes were between 82-83 and 55-65 °C, respectively. Molecular phylogeny placed Est2 and Est3 within Family XIII, alongside other Geobacillus esterases. DeepMSA2 revealed that Est2, Est3, and homologous sequences shared 12 conserved residues in their core domain (L39, D50, G53, G55, S57, G92, S94, G96, P108, P184, D193, and H223). BANΔIT analysis indicated that Est2 and Est3 had a significantly more rigid cap domain compared to Est30. Salt bridge analysis revealed that E150-R136, E124-K165, E137-R141, and E154-K157 salt bridges made Est2 and Est3 more stable compared to Est30. MD simulation indicated that Est3 exhibited greater fluctuations in the N-terminal region including conserved F25, cap domain, and C-terminal region, notably including H223, suggesting that these regions might influence esterase catalysis. The common residues in the ligand-binding sites of Est2-Est3 were determined as F25 and L167. The analysis of root mean square fluctuation (RMSF) revealed that region 1, encompassing F25 within the ß2-α1 loop of Est3, exhibited higher fluctuations compared to those of Est2. Overall, this study might provide valuable insights for future investigations aimed at improving esterase thermostability and catalytic efficiency, critical industrial traits, through targeted amino acid modifications within the N-terminal region, cap domain, and C-terminal region using rational protein engineering techniques.