Simulations reveal the key role of Arg15 in the promiscuous activity in the HisA enzyme

Multiscale QM/MM modelling of catalytic systems with ChemShell

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods are a powerful computational tool for the investigation of all forms of catalysis, as they allow for an accurate description of reactions occurring at catalytic sites in the context of a complicated electrostatic environment. The scriptable computational chemistry environment ChemShell is a leading software package for QM/MM calculations, providing a flexible, high performance framework for modelling both biomolecular and materials catalysis. We present an overview of recent applications of ChemShell to problems in catalysis and review new functionality introduced into the redeveloped Python-based version of ChemShell to support catalytic modelling. These include a fully guided workflow for biomolecular QM/MM modelling, starting from an experimental structure, a periodic QM/MM embedding scheme to support modelling of metallic materials, and a comprehensive set of tutorials for biomolecular and materials modelling.

Combined MD and QM/MM Calculations Reveal Allostery-Driven Promiscuity in Dipeptide Epimerases of Enolase Family

The adaptability of the active site to amplify the secondary function is supposed to be the fundamental cause of the promiscuity and the evolution of new functions in enzymes. In most cases, mutations occur close to the active site and/or in the catalytic site to change the active site plasticity to accommodate the non-native substrate. In the present study, using MD simulations and hybrid QM/MM calculations, we have shown a way to enhance the promiscuity, i. e., the allostery-driven promiscuity. Using a case study of the AEE enzyme where the capping loop recognizes the substrate, herein, we show that a single site mutation (D321G) far from the capping loop can induce a large conformational change in the capping loop to recognize different substrates for different functions. The QM/MM calculations for the WT and mutated enzyme provide a first validation of the mechanism of 1,1-proton transfer and dehydration by the AEE enzyme. Since AEE epimerase possesses a highly conserved TIM-barrel fold, we believe that our study provides a crucial lead to understanding the mechanism of emergence of secondary function which can be useful to repurpose ancient enzymes for modern usage.

Complex Loop Dynamics Underpin Activity, Specificity, and Evolvability in the (βα)8 Barrel Enzymes of Histidine and Tryptophan Biosynthesis

Enzymes are conformationally dynamic, and their dynamical properties play an important role in regulating their specificity and evolvability. In this context, substantial attention has been paid to the role of ligand-gated conformational changes in enzyme catalysis; however, such studies have focused on tremendously proficient enzymes such as triosephosphate isomerase and orotidine 5'-monophosphate decarboxylase, where the rapid (μs timescale) motion of a single loop dominates the transition between catalytically inactive and active conformations. In contrast, the (βα)₈-barrels of tryptophan and histidine biosynthesis, such as the specialist isomerase enzymes HisA and TrpF, and the bifunctional isomerase PriA, are decorated by multiple long loops that undergo conformational transitions on the ms (or slower) timescale. Studying the interdependent motions of multiple slow loops, and their role in catalysis, poses a significant computational challenge. This work combines conventional and enhanced molecular dynamics simulations with empirical valence bond simulations to provide rich details of the conformational behavior of the catalytic loops in HisA, PriA, and TrpF, and the role of their plasticity in facilitating bifunctionality in PriA and evolved HisA variants. In addition, we demonstrate that, similar to other enzymes activated by ligand-gated conformational changes, loops 3 and 4 of HisA and PriA act as gripper loops, facilitating the isomerization of the large bulky substrate ProFAR, albeit now on much slower timescales. This hints at convergent evolution on these different (βα)₈-barrel scaffolds. Finally, our work reemphasizes the potential of engineering loop dynamics as a tool to artificially manipulate the catalytic repertoire of TIM-barrel proteins.