Structural and Molecular Dynamics Analysis of Novel Bacterial Esterases: Insights into Substrate Binding and Thermostability for Rational Enzyme Desig

• T. Ndiitwani

  Department of Biochemistry, School of Physical and Chemical Sciences, Faculty of Natural and Agri-cultural Sciences, North-West University, Mafikeng Campus, MMabatho 2790, South Africa. Department of Biotechnology and Food Science, Faculty of Applied Sciences, Durban University of Technology, P.O. Box 1334, Durban, 4000, South Africa.

  Author

• F. Mokoena

  Department of Biochemistry, School of Physical and Chemical Sciences, Faculty of Natural and Agri-cultural Sciences, North-West University, Mafikeng Campus, MMabatho 2790, South Africa.

  Author

• S. Sabiu

  Department of Biotechnology and Food Science, Faculty of Applied Sciences, Durban University of Technology, P.O. Box 1334, Durban, 4000, South Africa.

  Author

• Israa M Shamkh

  Chief Computational Chemistry Department EBO Bio Solution Company, London, EC1V2NX, United Kingdom

  Author

• I Achilonu

  Protein Structure-Function and Research Unit, School of Molecular and Cell Biology, Faculty of Science, University of the Witwatersrand, Braamfontein, Johannesburg 2050, South Africa.

  Author

• MG. Tlou

  Department of Biochemistry, School of Physical and Chemical Sciences, Faculty of Natural and Agri-cultural Sciences, North-West University, Mafikeng Campus, MMabatho 2790, South Africa.

  Author

Bacterial esterases are pivotal biocatalysts with diverse industrial applications, yet their structural and dynamic mechanisms remain underexplored. This study investigates three novel esterases extracted and sequenced from Bacillus pumilus, Bacillus licheniformis, and Geobacillus kaustophilus, which were computationally modeled and functionally characterized for the first time. High-accuracy structural predictions were achieved using SWISS-MODEL (for B. pumilus and G. kaustophilus) and AlphaFold v2 (for B. licheniformis), with sequence identities exceeding 80–95% against their templates. Structural validation via MolProbity, PROCHECK, and QMEANDisCo confirmed stereochemical robustness (Ramachandran favored regions >90%, clash scores <2.0). Active site analysis using CASTp revealed conserved catalytic triads (Ser-His-Glu/Asp) and substrate-binding residues. Molecular docking with 4-nitrophenyl butyrate identified strong binding affinities (−4.93 to −5.21 kcal/mol), with B. pumilus exhibiting the lowest RMSD (0.78 Å). Molecular dynamics simulations (100 ns) highlighted dynamic stability, with progressive ligand burial (SASA reduction: 11–12%) and rigid torsional profiles. Notably, B. pumilus and G. kaustophilus displayed enhanced binding over time due to optimized van der Waals (−27.68 and −27.46 kcal/mol) and lipophilic interactions, while B. licheniformis showed weakened affinity, attributed to elevated ligand strain. Secondary structures remained stable (α-helices: ~70–75%), underscoring preserved enzymatic folds. This study advances the mechanistic understanding of esterase-ligand interactions, emphasizing the role of dynamic simulations in capturing time-dependent binding shifts. The findings position B. pumilus as a prime candidate for thermostability engineering and industrial biocatalysis, while B. licheniformis necessitates structural refinement. These insights bridge computational modeling with biotechnological applications, offering a framework for rational enzyme design.

• PDF

• Computational Science Section

Copyright (c) 2025 T. Ndiitwani, F. Mokoena, S. Sabiu, Israa M Shamkh, I Achilonu, MG. Tlou (Author)

This work is licensed under a Creative Commons Attribution 4.0 International License.

This work is licensed under a Creative Commons Attribution 4.0 International License.