Seeking the Source of Catalytic Efficiency of Lindane Dehydrochlorinase, LinA

Herein we present the results of an in-depth simulation study of LinA and its two variants. In our analysis, we combined the exploration of protein conformational dynamics with and without bound substrates (hexachlorocyclohexane (HCH) isomers) performed using molecular dynamics simulation followed by the extraction of the most frequently visited conformations and their characteristics with a detailed description of the interactions taking place in the active site between the respective HCH molecule and the first shell residues by using symmetry-adapted perturbation theory (SAPT) calculations. A detailed investigation of the conformational space of LinA substates has been accompanied by description of enzymatic catalytic steps carried out using a hybrid quantum mechanics/molecular mechanics (QM/MM) potential along with the computation of the potential of mean force (PMF) to estimate the free energy barriers for the studied transformations: dehydrochlorination of γ-, (-)-α-, and (+)-α-HCH by LinA-type I and -type II variants. The applied combination of computational techniques allowed us not only to characterize two LinA types but also to point to the most important differences between them and link their features to catalytic efficiency each of them possesses toward the respective ligand. More importantly it has been demonstrated that type I protein is more mobile, its active site has a larger volume, and the dehydrochlorination products are stabilized more strongly than in the case of type II enzyme, due to differences in the residues present in the active sites. Additionally, interaction energy calculations revealed very interesting patterns not predicted before but having the potential to be utilized in any attempts of improving LinA catalytic efficiency. On the basis of all these observations, LinA-type I protein seems to be more preorganized for the dehydrochlorination reaction it catalyzes than the type II variant.