Enzyme‐powered micro/nanomotors have myriads of potential applications in various areas. To efficiently reach those applications, it is necessary and critical to understand the fundamental aspects affecting the motion dynamics. Herein, we explored the impact of enzyme orientation on the performance of lipase‐powered nanomotors by tuning the lipase immobilization strategies. The influence of the lipase orientation and lid conformation on substrate binding and catalysis was analyzed using molecular dynamics simulations. Besides, the motion performance indicates that the hydrophobic binding (via OTES) represents the best orienting strategy, providing 48.4 % and 95.4 % increase in diffusion coefficient compared to hydrophilic binding (via APTES) and Brownian motion (no fuel), respectively (with C[triacetin] of 100 mm). This work provides vital evidence for the importance of immobilization strategy and corresponding enzyme orientation for the catalytic activity and in turn, the motion performance of nanomotors, and is thus helpful to future applications.

Enzyme‐powered micro/nanomotors have myriads of potential applications in various areas. To efficiently reach those applications, it is necessary and critical to understand the fundamental aspects affecting the motion dynamics. Herein, we explored the impact of enzyme orientation on the performance of lipase‐powered nanomotors by tuning the lipase immobilization strategies. The influence of the lipase orientation and lid conformation on substrate binding and catalysis was analyzed using molecular dynamics simulations. Besides, the motion performance indicates that the hydrophobic binding (via OTES) represents the best orienting strategy, providing 48.4 % and 95.4 % increase in diffusion coefficient compared to hydrophilic binding (via APTES) and Brownian motion (no fuel), respectively (with C[triacetin] of 100 mm). This work provides vital evidence for the importance of immobilization strategy and corresponding enzyme orientation for the catalytic activity and in turn, the motion performance of nanomotors, and is thus helpful to future applications.

The three‐dimensional structure of the enzymes provides very relevant information on the arrangement of the catalytic machinery and structural elements gating the active site pocket. The recent success of the neural network Alphafold2 in predicting the folded structure of proteins from the primary sequence with high levels of accuracy has revolutionized the protein design field. However, the application of Alphafold2 for understanding and engineering function directly from the obtained single static picture is not straightforward. Indeed, understanding enzymatic function requires the exploration of the ensemble of thermally accessible conformations that enzymes adopt in solution. In the present study, we evaluate the potential of Alphafold2 in assessing the effect of the mutations on the conformational landscape of the beta subunit of tryptophan synthase (TrpB). Specifically, we develop a template‐based Alphafold2 approach for estimating the conformational heterogeneity of several TrpB enzymes, which is needed for enhanced stand‐alone activity. Our results show the potential of Alphafold2, especially if combined with molecular dynamics simulations, for elucidating the changes induced by mutation in the conformational landscapes at a rather reduced computational cost, thus revealing its plausible application in computational enzyme design.

Halohydrin dehalogenases (HHDH) are industrially relevant biocatalysts exhibiting a promiscuous epoxide-ring opening reactivity in the presence of small nucleophiles, thus giving access to novel carbon–carbon, carbon–oxygen, carbon–nitrogen, and carbon–sulfur bonds. Recently, the repertoire of HHDH has been expanded, providing access to some novel HHDH subclasses exhibiting a broader epoxide substrate scope. In this work, we develop a computational approach based on the application of linear and non-linear dimensionality reduction techniques to long time-scale Molecular Dynamics (MD) simulations to study the HHDH conformational landscapes. We couple the analysis of the conformational landscapes to CAVER calculations to assess their impact on the active site tunnels and potential ability towards bulky epoxide ring opening reaction. Our study indicates that the analyzed HHDHs subclasses share a common breathing motion of the halide binding pocket, but present large deviations in the loops adjacent to the active site pocket and N-terminal regions. Such conformational differences affect the available tunnels for epoxide binding to the active site. The superior activity of the HHDH G subclass towards bulkier substrates is explained by the additional structural elements delimiting the active site region, its rich conformational heterogeneity, and the substantially wider and frequently observed active site tunnels. This study therefore provides key information for HHDH promiscuity and engineering.

Halohydrin dehalogenases (HHDHs) are promising enzymes for application in biocatalysis due to their promiscuous epoxide ring-opening activity with various anionic nucleophiles. So far, seven different HHDH subtypes A to G have been reported with subtype D containing the by far largest number of enzymes. Moreover, several characterized members of subtype D have been reported to display outstanding characteristics such as high catalytic activity, broad substrate spectra or remarkable thermal stability. Yet, no structure of a D-type HHDH has been reported to date that could be used to investigate and understand those features on a molecular level. We therefore solved the crystal structure of HheD2 from gamma proteobacterium HTCC2207 at 1.6Å resolution and used it as a starting point for targeted mutagenesis in combination with molecular dynamics (MD) simulation, in order to study the low thermal stability of HheD2 in comparison with other members of subtype D. This revealed a hydrogen bond between conserved residues Q160 and D198 to be connected with a high catalytic activity of this enzyme. Moreover, a flexible surface region containing two α-helices was identified to impact thermal stability of HheD2. Exchange of this surface region by residues of HheD3 yielded a variant with 10°C higher melting temperature and reaction temperature optimum. Overall, our results provide important insights into the structure-function relationship of HheD2 and presumably for other D-type HHDHs.