How to complete the tautomerization and substrate-assisted activation prior to C–C bond fission by meta-cleavage product hydrolase LigY?

Two feasible binding modes could complete the C–C bond fission of the substrate. One is the bidentate mode and five-coordination, and the other is the monodentate mode and five-coordination.

Finding Reactive Configurations: A Machine Learning Approach for Estimating Energy Barriers Applied to Sirtuin 5

Sirtuin 5 is a class III histone deacetylase that, unlike its classification, mainly catalyzes desuccinylation and demanoylation reactions. It is an interesting drug target that we use here to test new ideas for calculating reaction pathways of large molecular systems such as enzymes. A major issue with most schemes (e.g., adiabatic mapping) is that the resulting activation barrier height heavily depends on the chosen educt conformation. This makes the selection of the initial structure decisive for the success of the characterization. Here, we apply machine learning to a large number of molecular dynamics frames and potential energy barriers obtained by quantum mechanics/molecular mechanics calculations in order to identify (1) suitable start-conformations for reaction path calculations and (2) structural features relevant for the first step of the desuccinylation reaction catalyzed by Sirtuin 5. The latter generally aids the understanding of reaction mechanisms and important interactions in active centers. Using our novel approach, we found eleven key features that govern the reactivity. We were able to estimate reaction barriers with a mean absolute error of 3.6 kcal/mol and identified reactive configurations.

Identification and characterization of a meta-cleavage product hydrolase involved in biphenyl degradation from Arthrobacter sp. YC-RL1

Polychlorinated biphenyls (PCBs) are a group of persistent organic pollutants (POPs) widely existing in the environment. Arthrobacter sp. YC-RL1 is a biphenyl-degrading bacterium that shows metabolic versatility towards aromatic compounds. A 2-hydroxy-6-oxo-6-phenylhexa-2, 4-dienoate (HOPDA) hydrolase (BphD) gene involved in the biodegradation of biphenyl was cloned from strain YC-RL1 and heterologously expressed in Escherichia coli BL21 (DE3). The recombinant BphD_YC-RL1 was purified and characterized. BphD_YC-RL1 showed the highest activity at 45 °C and pH 7. It was stable under a wide range of temperature (20-50 °C). The enzyme had a K_m value of 0.14 mM, K_cat of 11.61 s^-1, and V_max of 0.027 U/mg. Temperature dependence catalysis exhibited a biphasic Arrhenius Plot with a transition at 20 °C. BphD_YC-RL1 was inactivated by SDS, Tween 20, Tween 80, Trition X-100, DTT, CHAPS, NBS, PMSF, and DEPC, but insensitive to EDTA. Site-directed mutagenesis of the active-site residues revealed that the catalytic triad residues (Ser115, His275, and Asp247) of BphD_YC-RL1 were necessary for its activity. The investigation of BphD_YC-RL1 not only provides new potential enzyme resource for the biodegradation of biphenyl but also helps deepen our understanding on the catalytic process and mechanism.

How Many Conformations of Enzymes Should Be Sampled for DFT/MM Calculations? A Case Study of Fluoroacetate Dehalogenase

The quantum mechanics/molecular mechanics (QM/MM) method (e.g., density functional theory (DFT)/MM) is important in elucidating enzymatic mechanisms. It is indispensable to study "multiple" conformations of enzymes to get unbiased energetic and structural results. One challenging problem, however, is to determine the minimum number of conformations for DFT/MM calculations. Here, we propose two convergence criteria, namely the Boltzmann-weighted average barrier and the disproportionate effect, to tentatively address this issue. The criteria were tested by defluorination reaction catalyzed by fluoroacetate dehalogenase. The results suggest that at least 20 conformations of enzymatic residues are required for convergence using DFT/MM calculations. We also tested the correlation of energy barriers between small QM regions and big QM regions. A roughly positive correlation was found. This kind of correlation has not been reported in the literature. The correlation inspires us to propose a protocol for more efficient sampling. This saves 50% of the computational cost in our current case.