Catalytic mechanism of the maloalkane dehalogenase LinB from Sphingomonas paucimobilis UT26

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.

Transient binding sites at the surface of haloalkane dehalogenase LinB as locations for fine-tuning enzymatic activity

The haloalkane dehalogenase LinB is a well-known enzyme that contains buried active site and is used for many modelling studies. Using classical molecular dynamics simulations of enzymes and substrates, we searched for transient binding sites on the surface of the LinB protein by calculating maps of enzyme-ligand interactions that were then transformed into sparse matrices. All residues considered as functionally important for enzyme performance (e.g., tunnel entrances) were excluded from the analysis to concentrate rather on non-obvious surface residues. From a set of 130 surface residues, twenty-six were proposed as a promising improvement of enzyme performance. Eventually, based on rational selection and filtering out the potentially unstable mutants, a small library of ten mutants was proposed to validate the possibility of fine-tuning the LinB protein. Nearly half of the predicted mutant structures showed improved activity towards the selected substrates, which demonstrates that the proposed approach could be applied to identify non-obvious yet beneficial mutations for enzyme performance especially when obvious locations have already been explored.

Defluorination Capability of l-2-Haloacid Dehalogenases in the HAD-Like Hydrolase Superfamily Correlates with Active Site Compactness

l-2-Haloacid dehalogenases, industrially and environmentally important enzymes that catalyse cleavage of the carbon-halogen bond in S-2-halocarboxylic acids, were known to hydrolyse chlorinated, brominated and iodinated substrates but no activity towards fluorinated compounds had been reported. A screen for novel dehalogenase activities revealed four l-2-haloacid dehalogenases capable of defluorination. We now report crystal structures for two of these enzymes, Bpro0530 and Rha0230, as well as for the related proteins PA0810 and RSc1362, which hydrolyse chloroacetate but not fluoroacetate, all at ∼2.2 Å resolution. Overall structure and active sites of these enzymes are highly similar. In molecular dynamics (MD) calculations, only the defluorinating enzymes sample more compact conformations, which in turn allow more effective interactions with the small fluorine atom. Structural constraints, based on X-ray structures and MD calculations, correctly predict the defluorination activity of the homologous enzyme ST2570.

Structural Analysis of the Ancestral Haloalkane Dehalogenase AncLinB-DmbA

Haloalkane dehalogenases (EC 3.8.1.5) play an important role in hydrolytic degradation of halogenated compounds, resulting in a halide ion, a proton, and an alcohol. They are used in biocatalysis, bioremediation, and biosensing of environmental pollutants and also for molecular tagging in cell biology. The method of ancestral sequence reconstruction leads to prediction of sequences of ancestral enzymes allowing their experimental characterization. Based on the sequences of modern haloalkane dehalogenases from the subfamily II, the most common ancestor of thoroughly characterized enzymes LinB from Sphingobium japonicum UT26 and DmbA from Mycobacterium bovis 5033/66 was in silico predicted, recombinantly produced and structurally characterized. The ancestral enzyme AncLinB-DmbA was crystallized using the sitting-drop vapor-diffusion method, yielding rod-like crystals that diffracted X-rays to 1.5 Å resolution. Structural comparison of AncLinB-DmbA with their closely related descendants LinB and DmbA revealed some differences in overall structure and tunnel architecture. Newly prepared AncLinB-DmbA has the highest active site cavity volume and the biggest entrance radius on the main tunnel in comparison to descendant enzymes. Ancestral sequence reconstruction is a powerful technique to study molecular evolution and design robust proteins for enzyme technologies.

H2O-Bridged Proton-Transfer Channel in Emitter Species Formation in Obelin Bioluminescence

Bioluminescence of a number of marine organisms is conditioned by Ca²⁺-regulated photoprotein (CaRP) with coelenterazine as the reaction substrate. The reaction product, coelenteramide, at the first singlet excited state (S₁) is the emitter of CaRP. The S₁-state coelenteramide is produced via the decomposition of coelenterazine dioxetanone. Experiments suggested that the neutral S₁-coelenteramide is the primary emitter species. This supposition contradicts with theoretical calculations showing that the anionic S₁-coelenteramide is a primary product of the decomposition of coelenterazine dioxetanone. In this study, applying molecular dynamic (MD) simulations and the hybrid quantum mechanics/molecular mechanics (QM/MM) method, we investigated a proton-transfer (PT) process taking place in CaRP obelin from Obelia longissima for emitter formation. Our calculations demonstrate a concerted PT process with a water molecule as a bridge between anionic S₁-coelenteramide and the nearest histidine residue. The low activation barrier as well as the strong hydrogen-bond network between the proton donor and the proton acceptor suggests a fast PT process comparable with that of the lifetime of excited anionic S₁-coelenteramide. The existence of the PT process eliminates the discrepancy between experimental and theoretical studies. The fast PT process at emitter formation can also take place in other CaRPs.

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.

Description of Transport Tunnel in Haloalkane Dehalogenase Variant LinB D147C+L177C from Sphingobium japonicum

The activity of enzymes with active sites buried inside their protein core highly depends on the efficient transport of substrates and products between the active site and the bulk solvent. The engineering of access tunnels in order to increase or decrease catalytic activity and specificity in a rational way is a challenging task. Here, we describe a combined experimental and computational approach to characterize the structural basis of altered activity in the haloalkane dehalogenase LinB D147C+L177C variant. While the overall protein fold is similar to the wild type enzyme and the other LinB variants, the access tunnels have been altered by introduced cysteines that were expected to form a disulfide bond. Surprisingly, the mutations have allowed several conformations of the amino acid chain in their vicinity, interfering with the structural analysis of the mutant by X-ray crystallography. The duration required for the growing of protein crystals changed from days to 1.5 years by introducing the substitutions. The haloalkane dehalogenase LinB D147C+L177C variant crystal structure was solved to 1.15 Å resolution, characterized and deposited to Protein Data Bank under PDB ID 6s06.

A Haloalkane Dehalogenase from Saccharomonospora viridis Strain DSM 43017, a Compost Bacterium with Unusual Catalytic Residues, Unique (S)-Enantiopreference, and High Thermostability

Haloalkane dehalogenases can cleave a carbon-halogen bond in a broad range of halogenated aliphatic compounds. However, a highly conserved catalytic pentad composed of a nucleophile, a catalytic base, a catalytic acid, and two halide-stabilizing residues is required for their catalytic activity. Only a few family members, e.g., DsaA, DmxA, or DmrB, remain catalytically active while employing a single halide-stabilizing residue. Here, we describe a novel haloalkane dehalogenase, DsvA, from a mildly thermophilic bacterium, Saccharomonospora viridis strain DSM 43017, possessing one canonical halide-stabilizing tryptophan (W125). At the position of the second halide-stabilizing residue, DsvA contains the phenylalanine F165, which cannot stabilize the halogen anion released during the enzymatic reaction by a hydrogen bond. Based on the sequence and structural alignments, we identified a putative second halide-stabilizing tryptophan (W162) located on the same α-helix as F165, but on the opposite side of the active site. The potential involvement of this residue in DsvA catalysis was investigated by the construction and biochemical characterization of the three variants, DsvA01 (F165W), DsvA02 (W162F), and DsvA03 (W162F and F165W). Interestingly, DsvA exhibits a preference for the (S)- over the (R)-enantiomers of β-bromoalkanes, which has not been reported before for any characterized haloalkane dehalogenase. Moreover, DsvA shows remarkable operational stability at elevated temperatures. The present study illustrates that protein sequences possessing an unconventional composition of catalytic residues represent a valuable source of novel biocatalysts.IMPORTANCE The present study describes a novel haloalkane dehalogenase, DsvA, originating from a mildly thermophilic bacterium, Saccharomonospora viridis strain DSM 43017. We report its high thermostability, remarkable operational stability at high temperatures, and an (S)-enantiopreference, which makes this enzyme an attractive biocatalyst for practical applications. Sequence analysis revealed that DsvA possesses an unusual composition of halide-stabilizing tryptophan residues in its active site. We constructed and biochemically characterized two single point mutants and one double point mutant and identified the noncanonical halide-stabilizing residue. Our study underlines the importance of searching for noncanonical catalytic residues in protein sequences.

CaverDock: a molecular docking-based tool to analyse ligand transport through protein tunnels and channels

MOTIVATION

Protein tunnels and channels are key transport pathways that allow ligands to pass between proteins' external and internal environments. These functionally important structural features warrant detailed attention. It is difficult to study the ligand binding and unbinding processes experimentally, while molecular dynamics simulations can be time-consuming and computationally demanding.

RESULTS

CaverDock is a new software tool for analysing the ligand passage through the biomolecules. The method uses the optimized docking algorithm of AutoDock Vina for ligand placement docking and implements a parallel heuristic algorithm to search the space of possible trajectories. The duration of the simulations takes from minutes to a few hours. Here we describe the implementation of the method and demonstrate CaverDock's usability by: (i) comparison of the results with other available tools, (ii) determination of the robustness with large ensembles of ligands and (iii) the analysis and comparison of the ligand trajectories in engineered tunnels. Thorough testing confirms that CaverDock is applicable for the fast analysis of ligand binding and unbinding in fundamental enzymology and protein engineering.

AVAILABILITY AND IMPLEMENTATION

User guide and binaries for Ubuntu are freely available for non-commercial use at https://loschmidt.chemi.muni.cz/caverdock/. The web implementation is available at https://loschmidt.chemi.muni.cz/caverweb/. The source code is available upon request.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Fluorescent substrates for haloalkane dehalogenases: Novel probes for mechanistic studies and protein labeling

Haloalkane dehalogenases are enzymes that catalyze the cleavage of carbon-halogen bonds in halogenated compounds. They serve as model enzymes for studying structure-function relationships of >100.000 members of the α/β-hydrolase superfamily. Detailed kinetic analysis of their reaction is crucial for understanding the reaction mechanism and developing novel concepts in protein engineering. Fluorescent substrates, which change their fluorescence properties during a catalytic cycle, may serve as attractive molecular probes for studying the mechanism of enzyme catalysis. In this work, we present the development of the first fluorescent substrates for this enzyme family based on coumarin and BODIPY chromophores. Steady-state and pre-steady-state kinetics with two of the most active haloalkane dehalogenases, DmmA and LinB, revealed that both fluorescent substrates provided specificity constant two orders of magnitude higher (0.14-12.6 μM^-1 s^-1) than previously reported representative substrates for the haloalkane dehalogenase family (0.00005-0.014 μM^-1 s^-1). Stopped-flow fluorescence/FRET analysis enabled for the first time monitoring of all individual reaction steps within a single experiment: (i) substrate binding, (ii-iii) two subsequent chemical steps and (iv) product release. The newly introduced fluorescent molecules are potent probes for fast steady-state kinetic profiling. In combination with rapid mixing techniques, they provide highly valuable information about individual kinetic steps and mechanism of haloalkane dehalogenases. Additionally, these molecules offer high specificity and efficiency for protein labeling and can serve as probes for studying protein hydration and dynamics as well as potential markers for cell imaging.

The impact of tunnel mutations on enzymatic catalysis depends on the tunnel-substrate complementarity and the rate-limiting step

Transport of ligands between bulk solvent and the buried active sites is a critical event in the catalytic cycle of many enzymes. The rational design of transport pathways is far from trivial due to the lack of knowledge about the effect of mutations on ligand transport. The main and an auxiliary tunnel of haloalkane dehalogenase LinB have been previously engineered for improved dehalogenation of 1,2-dibromoethane (DBE). The first chemical step of DBE conversion was enhanced by L177W mutation in the main tunnel, but the rate-limiting product release was slowed down because the mutation blocked the main access tunnel and hindered protein dynamics. Three additional mutations W140A + F143L + I211L opened-up the auxiliary tunnel and enhanced the product release, making this four-point variant the most efficient catalyst with DBE. Here we study the impact of these mutations on the catalysis of bulky aromatic substrates, 4-(bromomethyl)-6,7-dimethoxycoumarin (COU) and 8-chloromethyl-4,4'-difluoro-3,5-dimethyl-4-bora-3a,4a-diaza-s-indacene (BDP). The rate-limiting step of DBE conversion is the product release, whereas the catalysis of COU and BDP is limited by the chemical step. The catalysis of COU is mainly impaired by the mutation L177W, whereas the conversion of BDP is affected primarily by the mutations W140A + F143L + I211L. The combined computational and kinetic analyses explain the differences in activities between the enzyme-substrate pairs. The effect of tunnel mutations on catalysis depends on the rate-limiting step, the complementarity of the tunnels with the substrates and is clearly specific for each enzyme-substrate pair.

Droplet-based optofluidic systems for measuring enzyme kinetics

The study of enzyme kinetics is of high significance in understanding metabolic networks in living cells and using enzymes in industrial applications. To gain insight into the catalytic mechanisms of enzymes, it is necessary to screen an enormous number of reaction conditions, a process that is typically laborious, time-consuming, and costly when using conventional measurement techniques. In recent times, droplet-based microfluidic systems have proved themselves to be of great utility in large-scale biological experimentation, since they consume a minimal sample, operate at high analytical throughput, are characterized by efficient mass and heat transfer, and offer high levels of integration and automation. The primary goal of this review is the introduction of novel microfluidic tools and detection methods for use in high-throughput and sensitive analysis of enzyme kinetics. The first part of this review focuses on introducing basic concepts of enzyme kinetics and describing most common microfluidic approaches, with a particular focus on segmented flow. Herein, the key advantages include accurate control over the flow behavior, efficient mass and heat transfer, multiplexing, and high-level integration with detection modalities. The second part describes the current state-of-the-art platforms for high-throughput and sensitive analysis of enzyme kinetics. In addition to our categorization of recent advances in measuring enzyme kinetics, we have endeavored to critically assess the limitations of each of these detection approaches and propose strategies to improve measurements in droplet-based microfluidics. Graphical abstract.

Assessing Aerobic Biotransformation of Hexachlorocyclohexane Isomers by Compound-Specific Isotope Analysis

Contamination of soils and sediments with the highly persistent hexachlorocyclohexanes (HCHs) continues to be a threat for humans and the environment. Despite the existence of bacteria capable of biodegradation and cometabolic transformation of HCH isomers, such processes occur over time scales of decades and are thus challenging to assess. Here, we explored the use of compound-specific isotope analysis (CSIA) to track the aerobic biodegradation and biotransformation pathways of the most prominent isomers, namely, (-)-α-, (+)-α-, β-, γ-, and δ-HCH, through changes of their C and H isotope composition in assays of LinA2 and LinB enzymes. Dehydrochlorination of (+)-α-, γ-, and δ-HCH catalyzed by LinA2 was subject to substantial C and H isotope fraction with apparent ¹³C- and ²H-kinetic isotope effects (AKIEs) of up to 1.029 ± 0.001 and 6.7 ± 2.9, respectively, which are indicative of bimolecular eliminations. Hydrolytic dechlorination of δ-HCH by LinB exhibited even larger C but substantially smaller H isotope fractionation with ¹³C- and ²H-AKIEs of 1.073 ± 0.006 and 1.41 ± 0.04, respectively, which are typical for nucleophilic substitutions. The systematic evaluation of isomer-specific phenomena showed that, in addition to contaminant uptake limitations, diffusion-limited turnover ((-)-α-HCH), substrate dissolution (β-HCH), and potentially competing reactions catalyzed by constitutively expressed enzymes might bias the assessment of HCH biodegradation by CSIA at contaminated sites.

Recombinant Ca2+-regulated photoproteins of ctenophores: current knowledge and application prospects

Bright bioluminescence of ctenophores is conditioned by Ca²⁺-regulated photoproteins. Although they share many properties characteristic of hydromedusan Ca²⁺-regulated photoproteins responsible for light emission of marine animals belonging to phylum Cnidaria, a substantial distinction still exists. The ctenophore photoproteins appeared to be extremely sensitive to light-they lose the ability for bioluminescence on exposure to light over the entire absorption spectrum. Inactivation is irreversible because keeping the inactivated photoprotein in the dark does not recover its activity. The capability to emit light can be restored only by incubation of inactivated photoprotein with coelenterazine in the dark at alkaline pH in the presence of oxygen. Although these photoproteins were discovered many years ago, only the cloning of cDNAs encoding these unique bioluminescent proteins in the early 2000s has provided a new impetus for their studies. To date, cDNAs encoding Ca²⁺-regulated photoproteins from four different species of luminous ctenophores have been cloned. The amino acid sequences of ctenophore photoproteins turned out to completely differ from those of hydromedusan photoproteins (identity less than 29%) though also similar to them having three EF-hand Ca²⁺-binding sites. At the same time, these photoproteins reveal the same two-domain scaffold characteristic of hydromedusan photoproteins. This review is an attempt to systemize and critically evaluate the data scattered through various articles regarding the structural features of recombinant light-sensitive Ca²⁺-regulated photoproteins of ctenophores and their bioluminescent and physicochemical properties as well as to compare them with those of hydromedusan photoproteins. In addition, we also discuss the prospects of their biotechnology applications.

Computational evidence for the degradation mechanism of haloalkane dehalogenase LinB and mutants of Leu248 to 1-chlorobutane

The catalytic degradation ability of the haloalkane dehalogenase LinB toward 1-chlorobutane (1-CB) was studied using a combined quantum mechanics/molecular mechanics (QM/MM) approach. Two major processes are involved in the LinB-catalyzed removal of halogens: dechlorination and hydrolyzation. The present study confirmed the experimentally proposed reaction path at the molecular level. Moreover, based on nucleophilic substitution mechanism (SN2 reaction), dechlorination was found to be the rate-determining step of the entire reaction process. In this study, the Boltzmann-weighted average barrier for dechlorination was determined to be 17.0 kcal mol-1, which is fairly close to the experimental value (17.4 kcal mol-1). The state of His107 and the influence of Leu248 on the dechlorination process were also explored. In addition, an intriguing phenomenon was discovered: the potential energy barrier decreased by 7.5 kcal mol-1 when the Leu248 residue was mutated into Phe248. This discovery might be of great help to design new mutant enzymes or novel biocatalysts.

In silico design of potentially functional artificial metallo-haloalkane dehalogenase containing catalytic zinc

Artificial metalloenzymes are unique as they combine the good features of homogeneous and enzymatic catalysts, and they can potentially improve some difficult catalytic assays. This study reports a method that can be used to create an artificial metal-binding site prior to proving it to be functional in a wet lab. Haloalkane dehalogenase was grafted into a metal-binding site to form an artificial metallo-haloalkane dehalogenase and was studied for its potential functionalities in silico. Computational protocols regarding dynamic metal docking were studied using native metalloenzymes and functional artificial metalloenzymes. Using YASARA Structure, a simulation box covering template structure was created to be filled with water molecules followed by one mutated water molecule closest to the metal-binding site to metal ion. A simple energy minimization step was subsequently run using an AMBER force field to allow the metal ion to interact with the metal-binding residues. Long molecular dynamic simulation using YASARA Structure was performed to analyze the stability of the metal-binding site and the distance between metal-binding residues. Metal ions fluctuating around 2.0 Å across a 20 ns simulation indicated a stable metal-binding site. Metal-binding energies were predicted using FoldX, with a native metalloenzyme (carbonic anhydrase) scoring 18.0 kcal/mol and the best mutant model (C1a) scoring 16.4 kcal/mol. Analysis of the metal-binding site geometry was performed using CheckMyMetal, and all scores for the metalloenzymes and mutant models were in an acceptable range. Like native metalloenzymes, the metal-binding site of C1a was supported by residues in the second coordination shell to maintain a more coordinated metal-binding site. Short-chain multihalogenated alkanes (1,2-dibromoethane and 1,2,3-trichloropropane) were able to dock in the active site of C1a. The halides of the substrate were in contact with both the metal and halide-stabilizing residues, thus indicating a better stabilization of the substrate. The simple catalytic mechanism proposed is that the metal ion interacted with halogen and polarized the carbon-halogen bond, thus making the alpha carbon susceptible to attack by nucleophilic hydroxide. The interaction between halogen in the metal ion and halide-stabilizing residues may help to improve the stabilization of the substrate-enzyme complex and reduce the activation energy. This study reports a modified dynamic metal-docking protocol and validation tests to verify the metal-binding site. These approaches can be applied to design different kinds of artificial metalloenzymes or metal-binding sites.

Complete genome sequence of Bacillus sp. HBCD-sjtu, an efficient HBCD-degrading bacterium

Environmental pollution caused by the release of industrial chemicals is currently one of the most important environmental harms. Manufacturing chemicals can be biodegraded, and valuable intermediates can be used as pharmacophores in drug targeting and have several other useful purposes. Hexabromocyclododecane (HBCD), a non-aromatic brominated flame retardant, is a toxic compound that consists of a cycloaliphatic ring of 12 carbon atoms to which six bromine atoms are attached. It is formed by bromination of cis-trans-trans-1,5,9-cyclododecatriene, but its use is now restricted in several countries, because it is an environmental pollutant. Little is known about whether bacteria can degrade HBCD. A bacterial strain that degrades HBCD was recently isolated using enrichment culture techniques. Based on morphological, biochemical and phylogenetic analysis this isolate was categorized as Bacillus cereus and named strain HBCD-sjtu. Maximum growth and HBCD-degrading activity were observed when this strain was grown at 30 °C, pH 7.0 and 200 RPM in mineral salt medium containing 0.5 mm HBCD. The genome of strain HBCD-sjtu, which consists of only one circular chromosome, was sequenced. This whole genome sequence will be crucial for illuminating the molecular mechanisms of HBCD degradation.

Dehalogenases: From Improved Performance to Potential Microbial Dehalogenation Applications

The variety of halogenated substances and their derivatives widely used as pesticides, herbicides and other industrial products is of great concern due to the hazardous nature of these compounds owing to their toxicity, and persistent environmental pollution. Therefore, from the viewpoint of environmental technology, the need for environmentally relevant enzymes involved in biodegradation of these pollutants has received a great boost. One result of this great deal of attention has been the identification of environmentally relevant bacteria that produce hydrolytic dehalogenases—key enzymes which are considered cost-effective and eco-friendly in the removal and detoxification of these pollutants. These group of enzymes catalyzing the cleavage of the carbon-halogen bond of organohalogen compounds have potential applications in the chemical industry and bioremediation. The dehalogenases make use of fundamentally different strategies with a common mechanism to cleave carbon-halogen bonds whereby, an active-site carboxylate group attacks the substrate C atom bound to the halogen atom to form an ester intermediate and a halide ion with subsequent hydrolysis of the intermediate. Structurally, these dehalogenases have been characterized and shown to use substitution mechanisms that proceed via a covalent aspartyl intermediate. More so, the widest dehalogenation spectrum of electron acceptors tested with bacterial strains which could dehalogenate recalcitrant organohalides has further proven the versatility of bacterial dehalogenators to be considered when determining the fate of halogenated organics at contaminated sites. In this review, the general features of most widely studied bacterial dehalogenases, their structural properties, basis of the degradation of organohalides and their derivatives and how they have been improved for various applications is discussed.

Catalytic Cycle of Haloalkane Dehalogenases Toward Unnatural Substrates Explored by Computational Modeling

The anthropogenic toxic compound 1,2,3-trichloropropane is poorly degradable by natural enzymes. We have previously constructed the haloalkane dehalogenase DhaA31 by focused directed evolution ( Pavlova, M. et al. Nat. Chem. Biol. 2009 , 5 , 727 - 733 ), which is 32 times more active than the wild-type enzyme and is currently the most active variant known against that substrate. Recent evidence has shown that the structural basis responsible for the higher activity of DhaA31 was poorly understood. Here we have undertaken a comprehensive computational study of the main steps involved in the biocatalytic hydrolysis of 1,2,3-trichloropropane to decipher the structural basis for such enhancements. Using molecular dynamics and quantum mechanics approaches we have surveyed (i) the substrate binding, (ii) the formation of the reactive complex, (iii) the chemical step, and (iv) the release of the products. We showed that the binding of the substrate and its transport through the molecular tunnel to the active site is a relatively fast process. The cleavage of the carbon-halogen bond was previously identified as the rate-limiting step in the wild-type. Here we demonstrate that this step was enhanced in DhaA31 due to a significantly higher number of reactive configurations of the substrate and a decrease of the energy barrier to the S_N2 reaction. C176Y and V245F were identified as the key mutations responsible for most of those improvements. The release of the alcohol product was found to be the rate-limiting step in DhaA31 primarily due to the C176Y mutation. Mutational dissection of DhaA31 and kinetic analysis of the intermediate mutants confirmed the theoretical observations. Overall, our comprehensive computational approach has unveiled mechanistic details of the catalytic cycle which will enable a balanced design of more efficient enzymes. This approach is applicable to deepen the biochemical knowledge of a large number of other systems and may contribute to robust strategies in the development of new biocatalysts.

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.

Enantioselectivity of haloalkane dehalogenase LinB on the degradation of 1,2-dichloropropane: A QM/MM study

The hydrolysis dechlorination mechanism of a chiral organochlorinepollutant, 1,2-dichloropropane (DCP), catalyzed by haloalkane dehalogenase LinB has been investigated by using a combined quantum mechanics/molecular mechanics method. LinB was confirmed to be enantioselective towards the catabolism of the racemic mixture. Based on the S_N2 nucleophilic substitution mechanism, the dechlorination process was identified as the rate-determining step in LinB-catalyzed degradation of 1,2-dichloropropane, the Boltzmann-weighted average potential barrier of which is 18.8kcal/mol for the (R)-isomer and 24.0kcal/mol for the (S)-isomer. A particular water molecule near (S)-DCP in the reaction system can strongly disturb the dechlorination process, which can account for the enantioselectivity of LinB. Further electrostatic influence analysis indicates that proper mutation of Gly37 may improve the catalytic efficiency of LinB towards DCP.

Regio- and Enantioselective Sequential Dehalogenation of rac-1,3-Dibromobutane by Haloalkane Dehalogenase LinB

The hydrolytic dehalogenation of rac-1,3-dibromobutane catalyzed by the haloalkane dehalogenase LinB from Sphingobium japonicum UT26 proceeds in a sequential fashion: initial formation of intermediate haloalcohols followed by a second hydrolytic step to produce the final diol. Detailed investigation of the course of the reaction revealed favored nucleophilic displacement of the sec-halogen in the first hydrolytic event with pronounced R enantioselectivity. The second hydrolysis step proceeded with a regioselectivity switch at the primary position, with preference for the S enantiomer. Because of complex competition between all eight possible reactions, intermediate haloalcohols formed with moderate to good ee ((S)-4-bromobutan-2-ol: up to 87 %). Similarly, (S)-butane-1,3-diol was formed at a maximum ee of 35 % before full hydrolysis furnished the racemic diol product.

Discovery of Novel Haloalkane Dehalogenase Inhibitors

Haloalkane dehalogenases (HLDs) have recently been discovered in a number of bacteria, including symbionts and pathogens of both plants and humans. However, the biological roles of HLDs in these organisms are unclear. The development of efficient HLD inhibitors serving as molecular probes to explore their function would represent an important step toward a better understanding of these interesting enzymes. Here we report the identification of inhibitors for this enzyme family using two different approaches. The first builds on the structures of the enzymes' known substrates and led to the discovery of less potent nonspecific HLD inhibitors. The second approach involved the virtual screening of 150,000 potential inhibitors against the crystal structure of an HLD from the human pathogen Mycobacterium tuberculosis H37Rv. The best inhibitor exhibited high specificity for the target structure, with an inhibition constant of 3 μM and a molecular architecture that clearly differs from those of all known HLD substrates. The new inhibitors will be used to study the natural functions of HLDs in bacteria, to probe their mechanisms, and to achieve their stabilization.

Biochemical and biophysical characterisation of haloalkane dehalogenases DmrA and DmrB in Mycobacterium strain JS60 and their role in growth on haloalkanes

Haloalkane dehalogenases (HLDs) catalyse the hydrolysis of haloalkanes to alcohols, offering a biological solution for toxic haloalkane industrial wastes. Hundreds of putative HLD genes have been identified in bacterial genomes, but relatively few enzymes have been characterised. We identified two novel HLDs in the genome of Mycobacterium rhodesiae strain JS60, an isolate from an organochlorine-contaminated site: DmrA and DmrB. Both recombinant enzymes were active against C2-C6 haloalkanes, with a preference for brominated linear substrates. However, DmrA had higher activity against a wider range of substrates. The kinetic parameters of DmrA with 4-bromobutyronitrile as a substrate were Km  = 1.9 ± 0.2 mM, kcat  = 3.1 ± 0.2 s(-1) . DmrB showed the highest activity against 1-bromohexane. DmrA is monomeric, whereas DmrB is tetrameric. We determined the crystal structure of selenomethionyl DmrA to 1.7 Å resolution. A spacious active site and alternate conformations of a methionine side-chain in the slot access tunnel may contribute to the broad substrate activity of DmrA. We show that M. rhodesiae JS60 can utilise 1-iodopropane, 1-iodobutane and 1-bromobutane as sole carbon and energy sources. This ability appears to be conferred predominantly through DmrA, which shows significantly higher levels of upregulation in response to haloalkanes than DmrB.

Isolation and characterization of two novel strains capable of using cyclohexane as carbon source

Two strains capable of degrading cyclohexane were isolated from the soil and sludge of the wastewater treatment plant of the University of Stuttgart and a biotrickling filter system. The strains were classified as gram negative and identified as Acidovorax sp. CHX100 and Chelatococcus sp. CHX1100. Both strains have demonstrated the capability to degrade cycloalkanes (C5-C8), while only strain CHX1100 used as well short linear n-alkanes (C5-C8) as the sole source of carbon and energy. The growth of Acidovorax sp. CHX100 using cyclohexane was much faster compared to Chelatococcus sp. CHX1100. Degenerated primers were optimized from a set sequences of cyclohexanol dehydrogenase genes (chnA) as well as cyclohexanone monooxygenases (chnB) and used to amplify the gene cluster, which encodes the conversion of cyclohexanol to caprolactone. Phylogenetic analysis has indicated that the two gene clusters belong to different groups. The cyclohexane monooxygenase-induced activity which oxidizes also indole to 5-hydroxyindole has indicated the presence of a CYP-type system monooxygenase involved in the transformation of cyclohexane to cyclohexanol.

Stepwise enhancement of catalytic performance of haloalkane dehalogenase LinB towards β-hexachlorocyclohexane

Two haloalkane dehalogenases, LinB_UT and LinB_MI, each with 296 amino acid residues, exhibit only seven amino acid residue differences between them, but LinB_MI’s catalytic performance towards β-hexachlorocyclohexane (β-HCH) is considerably higher than LinB_UT’s. To elucidate the molecular basis governing this difference, intermediate mutants between LinB_UT and LinB_MI were constructed and kinetically characterized. The activities of LinB_UT-based mutants gradually increased by cumulative mutations into LinB_UT, and the effects of the individual amino acid substitutions depended on combination with other mutations. These results indicated that LinB_UT’s β-HCH degradation activity can be enhanced in a stepwise manner by the accumulation of point mutations.

Computationally efficient and accurate enantioselectivity modeling by clusters of molecular dynamics simulations

Computational approaches could decrease the need for the laborious high-throughput experimental screening that is often required to improve enzymes by mutagenesis. Here, we report that using multiple short molecular dynamics (MD) simulations makes it possible to accurately model enantioselectivity for large numbers of enzyme-substrate combinations at low computational costs. We chose four different haloalkane dehalogenases as model systems because of the availability of a large set of experimental data on the enantioselective conversion of 45 different substrates. To model the enantioselectivity, we quantified the frequency of occurrence of catalytically productive conformations (near attack conformations) for pairs of enantiomers during MD simulations. We found that the angle of nucleophilic attack that leads to carbon-halogen bond cleavage was a critical variable that limited the occurrence of productive conformations; enantiomers for which this angle reached values close to 180° were preferentially converted. A cluster of 20-40 very short (10 ps) MD simulations allowed adequate conformational sampling and resulted in much better agreement to experimental enantioselectivities than single long MD simulations (22 ns), while the computational costs were 50-100 fold lower. With single long MD simulations, the dynamics of enzyme-substrate complexes remained confined to a conformational subspace that rarely changed significantly, whereas with multiple short MD simulations a larger diversity of conformations of enzyme-substrate complexes was observed.

Structural and functional analysis of a novel haloalkane dehalogenase with two halide-binding sites

The crystal structure of the novel haloalkane dehalogenase DbeA fromBradyrhizobium elkaniiUSDA94 revealed the presence of two chloride ions buried in the protein interior. The first halide-binding site is involved in substrate binding and is present in all structurally characterized haloalkane dehalogenases. The second halide-binding site is unique to DbeA. To elucidate the role of the second halide-binding site in enzyme functionality, a two-point mutant lacking this site was constructed and characterized. These substitutions resulted in a shift in the substrate-specificity class and were accompanied by a decrease in enzyme activity, stability and the elimination of substrate inhibition. The changes in enzyme catalytic activity were attributed to deceleration of the rate-limiting hydrolytic step mediated by the lower basicity of the catalytic histidine.

Computational library design for increasing haloalkane dehalogenase stability

We explored the use of a computational design framework for the stabilization of the haloalkane dehalogenase LinB. Energy calculations, disulfide bond design, molecular dynamics simulations, and rational inspection of mutant structures predicted many stabilizing mutations. Screening of these in small mutant libraries led to the discovery of seventeen point mutations and one disulfide bond that enhanced thermostability. Mutations located in or contacting flexible regions of the protein had a larger stabilizing effect than mutations outside such regions. The combined introduction of twelve stabilizing mutations resulted in a LinB mutant with a 23 °C increase in apparent melting temperature (Tm,app , 72.5 °C) and an over 200-fold longer half-life at 60 °C. The most stable LinB variants also displayed increased compatibility with co-solvents, thus allowing substrate conversion and kinetic resolution at much higher concentrations than with the wild-type enzyme.

Release of halide ions from the buried active site of the haloalkane dehalogenase LinB revealed by stopped-flow fluorescence analysis and free energy calculations

Release of halide ions is an essential step of the catalytic cycle of haloalkane dehalogenases. Here we describe experimentally and computationally the process of release of a halide anion from the buried active site of the haloalkane dehalogenase LinB. Using stopped-flow fluorescence analysis and umbrella sampling free energy calculations, we show that the anion binding is ion-specific and follows the ordering I(-) > Br(-) > Cl(-). We also address the issue of the protonation state of the catalytic His272 residue and its effect on the process of halide release. While deprotonation of His272 increases binding of anions in the access tunnel, we show that the anionic ordering does not change with the switch of the protonation state. We also demonstrate that a sodium cation could relatively easily enter the active site, provided the His272 residue is singly protonated, and replace thus the missing proton. In contrast, Na(+) is strongly repelled from the active site containing the doubly protonated His272 residue. Our study contributes toward understanding of the reaction mechanism of haloalkane dehalogenase enzyme family. Determination of the protonation state of the catalytic histidine throughout the catalytic cycle remains a challenge for future studies.

The effect of a unique halide-stabilizing residue on the catalytic properties of haloalkane dehalogenase DatA from Agrobacterium tumefaciens C58

Haloalkane dehalogenases catalyze the hydrolysis of carbon-halogen bonds in various chlorinated, brominated and iodinated compounds. These enzymes have a conserved pair of halide-stabilizing residues that are important in substrate binding and stabilization of the transition state and the halide ion product via hydrogen bonding. In all previously known haloalkane dehalogenases, these residues are either a pair of tryptophans or a tryptophan-asparagine pair. The newly-isolated haloalkane dehalogenase DatA from Agrobacterium tumefaciens C58 (EC 3.8.1.5) possesses a unique halide-stabilizing tyrosine residue, Y109, in place of the conventional tryptophan. A variant of DatA with the Y109W mutation was created and the effects of this mutation on the structure and catalytic properties of the enzyme were studied using spectroscopy and pre-steady-state kinetic experiments. Quantum mechanical and molecular dynamics calculations were used to obtain a detailed analysis of the hydrogen-bonding patterns within the active sites of the wild-type and the mutant, as well as of the stabilization of the ligands as the reaction proceeds. Fluorescence quenching experiments suggested that replacing the tyrosine with tryptophan improves halide binding by 3.7-fold, presumably as a result of the introduction of an additional hydrogen bond. Kinetic analysis revealed that the mutation affected the substrate specificity of the enzyme and reduced its K(0.5) for selected halogenated substrates by a factor of 2-4, without impacting the rate-determining hydrolytic step. We conclude that DatA is the first natural haloalkane dehalogenase that stabilizes its substrate in the active site using only a single hydrogen bond, which is a new paradigm in catalysis by this enzyme family.

Crystal structure and site-directed mutagenesis analyses of haloalkane dehalogenase LinB from Sphingobium sp. strain MI1205

The enzymes LinB(UT) and LinB(MI) (LinB from Sphingobium japonicum UT26 and Sphingobium sp. MI1205, respectively) catalyze the hydrolytic dechlorination of β-hexachlorocyclohexane (β-HCH) and yield different products, 2,3,4,5,6-pentachlorocyclohexanol (PCHL) and 2,3,5,6-tetrachlorocyclohexane-1,4-diol (TCDL), respectively, despite their 98% identity in amino acid sequence. To reveal the structural basis of their different enzymatic properties, we performed site-directed mutagenesis and X-ray crystallographic studies of LinB(MI) and its seven point mutants. The mutation analysis revealed that the seven amino acid residues uniquely found in LinB(MI) were categorized into three groups based on the efficiency of the first-step (from β-HCH to PCHL) and second-step (from PCHL to TCDL) conversions. Crystal structure analyses of wild-type LinB(MI) and its seven point mutants indicated how each mutated residue contributed to the first- and second-step conversions by LinB(MI). The dynamics simulation analyses of wild-type LinB(MI) and LinB(UT) revealed that the entrance of the substrate access tunnel of LinB(UT) was more flexible than that of LinB(MI), which could lead to the different efficiencies of dehalogenation activity between these dehalogenases.

Oxygen activation of apo-obelin-coelenterazine complex

Ca(2+) -regulated photoproteins use a noncovalently bound 2-hydroperoxycoelenterazine ligand to emit light in response to Ca(2+) binding. To better understand the mechanism of formation of active photoprotein from apoprotein, coelenterazine and molecular oxygen, we investigated the spectral properties of the anaerobic apo-obelin-coelenterazine complex and the kinetics of its conversion into active photoprotein after exposure to air. Our studies suggest that coelenterazine bound within the anaerobic complex might be a mixture of N7-protonated and C2(-) anionic forms, and that oxygen shifts the equilibrium in favor of the C2(-) anion as a result of peroxy anion formation. Proton removal from N7 and further protonation of peroxy anion and the resulting formation of 2-hydroperoxycoelenterazine in obelin might occur with the assistance of His175. It is proposed that this conserved His residue might play a key role both in formation of active photoprotein and in Ca(2+) -triggering of the bioluminescence reaction.

A single mutation in a tunnel to the active site changes the mechanism and kinetics of product release in haloalkane dehalogenase LinB

Many enzymes have buried active sites. The properties of the tunnels connecting the active site with bulk solvent affect ligand binding and unbinding and also the catalytic properties. Here, we investigate ligand passage in the haloalkane dehalogenase enzyme LinB and the effect of replacing leucine by a bulky tryptophan at a tunnel-lining position. Transient kinetic experiments show that the mutation significantly slows down the rate of product release. Moreover, the mechanism of bromide ion release is changed from a one-step process in the wild type enzyme to a two-step process in the mutant. The rate constant of bromide ion release corresponds to the overall steady-state turnover rate constant, suggesting that product release became the rate-limiting step of catalysis in the mutant. We explain the experimental findings by investigating the molecular details of the process computationally. Analysis of trajectories from molecular dynamics simulations with a tunnel detection software reveals differences in the tunnels available for ligand egress. Corresponding differences are seen in simulations of product egress using a specialized enhanced sampling technique. The differences in the free energy barriers for egress of a bromide ion obtained using potential of mean force calculations are in good agreement with the differences in rates obtained from the transient kinetic experiments. Interactions of the bromide ion with the introduced tryptophan are shown to affect the free energy barrier for its passage. The study demonstrates how the mechanism of an enzymatic catalytic cycle and reaction kinetics can be engineered by modification of protein tunnels.

Comparative bioremediation potential of four rhizospheric microbial species against lindane

Four microbial species (Kocuria rhizophila, Microbacterium resistens, Staphylococcus equorum and Staphylococcus cohnii subspecies urealyticus) were isolated from the rhizospheric zone of selected plants growing in a lindane contaminated environment and acclimatized in lindane spiked media (5-100 μg mL⁻¹). The isolated species were inoculated with soil containing 5, 50 and 100 mg kg⁻¹ of lindane and incubated at room temperature. Soil samples were collected periodically to evaluate the microbial dissipation kinetics, dissipation rate, residual lindane concentration and microbial biomass carbon (MBC). There was a marked difference (p < 0.05) in the MBC content and lindane dissipation rate of microbial isolates cultured in three different lindane concentrations. Further, the dissipation rate tended to decrease with increasing lindane concentrations. After 45 d, the residual lindane concentrations in three different spiked soils were reduced to 0%, 41% and 33%, respectively. Among the four species, S. cohnii subspecies urealyticus exhibited maximum dissipation (41.65 mg kg⁻¹) and can be exploited for the in situ remediation of low to medium level lindane contaminated soils.

Redesigning dehalogenase access tunnels as a strategy for degrading an anthropogenic substrate

Engineering enzymes to degrade anthropogenic compounds efficiently is challenging. We obtained Rhodococcus rhodochrous haloalkane dehalogenase mutants with up to 32-fold higher activity than wild type toward the toxic, recalcitrant anthropogenic compound 1,2,3-trichloropropane (TCP) using a new strategy. We identified key residues in access tunnels connecting the buried active site with bulk solvent by rational design and randomized them by directed evolution. The most active mutant has large aromatic residues at two out of three randomized positions and two positions modified by site-directed mutagenesis. These changes apparently enhance activity with TCP by decreasing accessibility of the active site for water molecules, thereby promoting activated complex formation. Kinetic analyses confirmed that the mutations improved carbon-halogen bond cleavage and shifted the rate-limiting step to the release of products. Engineering access tunnels by combining computer-assisted protein design with directed evolution may be a valuable strategy for refining catalytic properties of enzymes with buried active sites.