A computational study of the dechlorination of β-hexachlorocyclohexane (β-HCH) catalyzed by the haloalkane dehalogenase LinB

Insights into generalization of the rate-limiting steps of the dehalogenation by LinB and DhaA: A computational approach

LinB and DhaA are well-known haloalkane dehalogenases (HLDs) capable of converting a plethora of halogenated alkanes, also those considered persistent pollutants. The dehalogenation reaction that these two enzymes catalyze has been studied to determine its rate-limiting step (rls) for the last two decades now. As a result, it has been determined that HLDs can show different rate-limiting steps for individual substrates, and at this point we do not have a basis for any generalization in this matter. Therefore, in this work we aimed at gaining insights into the enzymatic dehalogenation of selected dibromo- and bromochloro-ethanes and propanes by LinB and DhaA using computational approach to determine whether defined structural similarities of the substrates result in a unified mechanism and the same rls. By predicting halogen binding isotope effects (BIEs) as well as computing interaction energy for each HLD-ligand complex the nature of the protein-ligand interactions has been characterized. Furthermore, C and Br kinetic isotope effects (KIEs) as well as the minimum free energy paths (MFEPs) were computed to investigate the chemical reaction for the selected systems. Accuracy of the approach and robustness of the computational predictions were validated by measuring KIEs on the selected reactions. Overall results strongly indicate that any generalization with respect to the enzymatic process involving various ligands in the case of DhaA is impossible, even if the considered ligands are structurally very similar as those analyzed in the present study. Moreover, even small structural differences such as changing of one of the (non-leaving) halogen substituents may lead to significant changes in the enzymatic process and result in a different rls in the case of LinB. It has also been demonstrated that KIEs themselves cannot be used as rls indicators in the reactions catalyzed by the studied HLDs.

Removal of organochlorine pesticides and metagenomic analysis by multi-stage constructed wetland treating landfill leachate

Constructed wetlands (CWs) can effectively treat landfill leachate (LL). However, there is limited research on the removal of organochlorine pesticides (OCPs) refractory organics during LL treatment in CWs. In this study, multi-stage subsurface flow CWs was used to treat LL, and the removal fate of hexachlorocyclohexane (HCH) and dichlorodiphenyltrichloroethane (DDT) in CWs was investigated. The structural differences between plant roots and substrate microbial communities were compared and the Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathway of organic matter was analyzed based on metagenomic analysis. The results showed that substrate adsorption (50.55%-72.74%) and microbial degradation (20.38%-27.89%) were the main ways to remove OCPs. The Proteobacteria occupied a dominant position in the CWs system, among which Betaproteobacteria (34.37%-35.90%) were contained in the substrate, and Alphaproteobacteria (21.19%-23.84%) was a more dominant microorganism in plant roots. Formaldehyde assimilation and serine pathway were the main pathways of methane metabolism. This study provides a reference for the removal mechanism of OCPs to promote the application of CWs technology in LL treatment.

Transformation of short-chain chlorinated paraffins by the bacterial haloalkane dehalogenase LinB - Formation of mono- and di-hydroxylated metabolites

Short-chain chlorinated paraffins (SCCPs) are listed as persistent organic pollutants (POPs) under the Stockholm Convention. Such substances are toxic, bioaccumulating, transported over long distances and degrade slowly in the environment. Certain bacterial strains of the Sphingomonadacea family are able to degrade POPs, such as hexachlorocyclohexanes (HCHs) and hexabromocyclododecanes (HBCDs). The haloalkane dehalogenase LinB, expressed in certain Sphingomonadacea, is able to catalyze the transformation of haloalkanes to hydroxylated compounds. Therefore, LinB is a promising candidate for conversion of SCCPs. Hence, a mixture of chlorinated tridecanes was exposed in vitro to LinB, which was obtained through heterologous expression in Escherichia coli. Liquid chromatography mass spectrometry (LC-MS) was used to analyze chlorinated tridecanes and their transformation products. A chloride-enhanced soft ionization method, which favors the formation of chloride adducts [M+Cl]^- without fragmentation, was applied. Mathematical deconvolution was used to distinguish interfering mass spectra of paraffinic, mono-olefinic and di-olefinic compounds. Several mono- and di-hydroxylated products including paraffinic, mono-olefinic and di-olefinic compounds were found after LinB exposure. Mono- (rt = 5.9-6.9 min) and di-hydroxylated (rt = 3.2-4.5 min) compounds were separated from starting material (rt = 7.7-8.5 min) by reversed phase LC. Chlorination degrees of chlorinated tridecanes increased during LinB-exposure from n_Cl = 8.80 to 9.07, indicating a preferential transformation of lower chlorinated (Cl_<9) tridecanes. Thus, LinB indeed catalyzed a dehalohydroxylation of chlorinated tridecanes, tridecenes and tridecadienes. The observed hydroxylated compounds are relevant CP transformation products whose environmental and toxicological effects should be further investigated.

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.

Mechanistic study of the ATP hydrolysis reaction in dynein motor protein

Dynein, a large and complex motor protein, harnesses energy from adenosine triphosphate (ATP) hydrolysis to regulate essential cellular activities. The ATP hydrolysis mechanism for the dynein motor is still shrouded in mystery. Herein, molecular dynamics simulations of a dynein motor disclosed that two water molecules are present close to the γ-phosphate of ATP and Glu1742 at the AAA1 site of dynein. We have proposed three possible mechanisms for the ATP hydrolysis. We divulge by using a quantum mechanics/molecular mechanics (QM/MM) study that two water molecules and Glu1742 are crucial for facilitating the ATP hydrolysis reaction in dynein. Moreover, the ATP hydrolysis step is initiated by the activation of lytic water (W1) by Glu1742 through relay proton transfers with the help of auxiliary water (W2) yielding HPO₄^2- and ADP, as a product. In the next step, a proton is shifted back from Glu1742 to generate inorganic phosphate (H₂PO₄^-) via another relay proton transfer event. The overall activation barrier for the Glu1742 assisted ATP hydrolysis is found to be the most favourable pathway compared to other plausible pathways. We also unearthed that ATP hydrolysis in dynein follows a so-called associative-like pathway in its rate-limiting step. Our study ascertained the important indirect roles of the two amino acids (such as Arg2109, Asn1792) and Mg²⁺ ion in the ATP hydrolysis of dynein. Additionally, multiple sequence alignment of the different organisms of dynein motors has conveyed the evolutionary importance of the Glu1742, Asn1742, and Arg2109 residues, respectively. As similar mechanisms are also prevalent in other motors, and GTPase and ATPase enzymes, the present finding spells out the definitive requirement of a proton relay process through an extended water-chain as one of the key components in an enzymatic ATP hydrolysis reaction.

Measurement and Prediction of Chlorine Kinetic Isotope Effects in Enzymatic Systems

Approaches to determine chlorine kinetic isotope effects (Cl-KIEs) on enzymatic dehalogenations are discussed and illustrated by representative examples. Three aspects are considered. First methodology for experimental measurement of Cl-KIEs, with stress being on FAB-IRMS technique developed in our laboratory, is described. Subsequently, we concentrate our discussion on the consequences of reaction complexity in the interpretation of experimental values, a problem especially important in cases of polychlorinated reactants. The most fruitful studies of enzymatic dehalogenations by Cl-KIEs require their theoretical evaluation, hence the computational focus of the second part of this chapter.

Organochlorine pesticides in soil, air, and vegetation at and around a contaminated site in southwestern China: Concentration, transmission, and risk evaluation

Remediation and management of contaminated sites have become a prevalent problem under the current situation in China. The present study was conducted to investigate the concentration, transmission, and health risk of organochlorine pesticides (OCPs) in soil, air, and vegetation at and around a typical pesticide-contaminated site located in southwestern China. Exchange flux between soil and air was calculated to investigate the transmission of OCPs. Hexachlorocyclohexane (ΣHCHs, top soil: 19.1 mg/kg d.w., air: 52.3 ng/m³, vegetation: 0.17 mg/kg d.w.) and dichlorodiphenyltrichloroethane (ΣDDTs, top soil: 200 mg/kg d.w., air: 130 ng/m³, vegetation: 0.78 mg/kg d.w.) were the dominant pollutants at the contaminated site. Around the site, the soil, air, and vegetation samples had higher OCP concentrations than those in the normal areas, which were found to be dominated by ΣHCHs (top soil: 129 ng/g d.w., air: 5.09 ng/m³, vegetation: 81.8 ng/g d.w.) and ΣDDTs (top soil: 360 ng/g, air: 7.47 ng/m³, vegetation: 189 ng/g d.w.). The fugacity fractions of OCPs (>0.7) showed a net volatilization from soil into air, indicating that soil is an obvious pollution source for the atmosphere, especially at the site. Through human health risk evaluation, serious carcinogenic risk (CR_n) and hazardous index (HI_n) were found at the site (CR_n: 7.4 × 10^-6-1.04 × 10^-4, HI_n: 0.02 to 8.97) and the nearby areas (CR_n: 3.37 × 10^-6 for adults and 1.68 × 10^-6 for children on average).

Theoretical study on the hydrolytic step in the biotransformation of β-hexachlorocyclohexane degraded by haloalkane dehalogenase LinB

The hydrolytic process of LinB-catalyzed biotransformation of a notorious contaminant β-HCH was investigated in atomistic detail with a combined quantum mechanics/molecular mechanics approach. The Boltzmann-weighted averaging method amended by disproportionate effect analysis was showed to capture the fluctuation of a single molecule enzyme reaction. With the potential barriers of 18.7 and 2.6 kcal/mol, two elementary steps that refer to formation and decomposition of a tetrahedral intermediate are involved in the hydrolytic reaction, respectively. Polarized by Glu132, His272 serves as a proton carrier along the whole hydrolysis reaction. The electrostatic influence analysis highlighted residue Leu248 as a possible mutation target for rational design of LinB in enzyme modification. Further spatial location analysis provided explanation for the opposite effect of Asn38 toward the two elementary steps. Getting insight into the catalytic details and the structure and function of LinB can enrich the knowledge of it and promote its application in bioremediation of chlorinated hydrocarbon pollutants.