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

Haloalkane dehalogenases other than LinB can contribute to the γ-hexachlorocyclohexane utilization

Haloalkane dehalogenases (HLDs) convert halogenated compounds to corresponding alcohols by a simple hydrolytic mechanism. Although many bacterial strains possess HLDs or HLD homologues, LinB is the only HLD known to be involved in the utilization of γ-hexachlorocyclohexane (γ-HCH), a man-made chlorinated pesticide. In this study, to gain insight into the functional evolution of HLDs toward γ-HCH utilization, the linB gene in γ-HCH-degrading Sphingobium japonicum strain UT26 was replaced by 7 other HLD or HLD homologous genes, including 3 putative ancestral enzymes. Interestingly, strains carrying genes for DmmA from a marine metagenome or Rluc_anc, an ancestor of Renilla-luciferin 2-monooxygenase and LinB, produced 2,5-dichlorophenol and 2,5-dichlorohydroquinone from γ-HCH, which are indicators of LinB activity, and grew in minimal medium supplied with γ-HCH as a sole carbon source. These results indicated that other HLDs, in addition to LinB, can play roles in the γ-HCH utilization.

Reinforcing Tunnel Network Exploration in Proteins Using Gaussian Accelerated Molecular Dynamics

Tunnels are structural conduits in biomolecules responsible for transporting chemical compounds and solvent molecules from the active site. They have been shown to be present in a wide variety of enzymes across all functional and structural classes. However, the study of such pathways is experimentally challenging, because they are typically transient. Computational methods, such as molecular dynamics (MD) simulations, have been successfully proposed to explore tunnels. Conventional MD (cMD) provides structural details to characterize tunnels but suffers from sampling limitations to capture rare tunnel openings on longer time scales. Therefore, in this study, we explored the potential of Gaussian accelerated MD (GaMD) simulations to improve the exploration of complex tunnel networks in enzymes. We used the haloalkane dehalogenase LinB and its two variants with engineered transport pathways, which are not only well-known for their application potential but have also been extensively studied experimentally and computationally regarding their tunnel networks and their importance in multistep catalytic reactions. Our study demonstrates that GaMD efficiently improves tunnel sampling and allows the identification of all known tunnels for LinB and its two mutants. Furthermore, the improved sampling provided insight into a previously unknown transient side tunnel (ST). The extensive conformational landscape explored by GaMD simulations allowed us to investigate in detail the mechanism of ST opening. We determined variant-specific dynamic properties of ST opening, which were previously inaccessible due to limited sampling of cMD. Our comprehensive analysis supports multiple indicators of the functional relevance of the ST, emphasizing its potential significance beyond structural considerations. In conclusion, our research proves that the GaMD method can overcome the sampling limitations of cMD for the effective study of tunnels in enzymes, providing further means for identifying rare tunnels in enzymes with the potential for drug development, precision medicine, and rational protein engineering.

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.

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.

Lessons from the genomes of lindane-degrading sphingomonads

Bacterial strains capable of degrading man-made xenobiotic compounds are good materials to study bacterial evolution towards new metabolic functions. Lindane (γ-hexachlorocyclohexane, γ-HCH, or γ-BHC) is an especially good target compound for the purpose, because it is relatively recalcitrant but can be degraded by a limited range of bacterial strains. A comparison of the complete genome sequences of lindane-degrading sphingomonad strains clearly demonstrated that (i) lindane-degrading strains emerged from a number of different ancestral hosts that have recruited lin genes encoding enzymes that are able to channel lindane to central metabolites, (ii) in sphingomonads lin genes have been acquired by horizontal gene transfer mediated by different plasmids and in which IS6100 plays a role in recruitment and distribution of genes, and (iii) IS6100 plays a role in dynamic genome rearrangements providing genetic diversity to different strains and ability to evolve to other states. Lindane-degrading bacteria whose genomes change so easily and quickly are also fascinating starting materials for tracing the bacterial evolution process experimentally in a relatively short time period. As the origin of the specific lin genes remains a mystery, such genes will be useful probes for exploring the cryptic 'gene pool' available to bacteria.

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.

Structure-Based Engineering of Amidase from Pantoea sp. for Efficient 2-Chloronicotinic Acid Biosynthesis

2-Chloronicotinic acid is a key intermediate of pharmaceuticals and pesticides. Amidase-catalyzed hydrolysis provides a promising enzymatic method for 2-chloronicotinic acid production from 2-chloronicotinamide. However, biocatalytic hydrolysis of 2-chloronicotinamide is difficult due to the strong steric and electronic effect caused by 2-position chlorine substituent of the pyridine ring. In this study, an amidase from a Pantoea sp. (Pa-Ami) was designed and engineered to have improved catalytic properties. Single mutant G175A and double mutant G175A/A305T strains exhibited 3.2- and 3.7-fold improvements in their specific activity for 2-chloronicotinamide, and the catalytic efficiency was significantly increased, with k _cat/K_m values 3.1 and 10.0 times higher than that of the wild type, respectively. Structure-function analysis revealed that the distance between Oγ of Ser177 (involved in the catalytic triad) and the carbonyl carbon of 2-chloronicotinamide was shortened in the G175A mutant, making the nucleophilic attack on the Oγ of Ser177 easier by virtue of proper orientation. In addition, the A305T mutation contributed to a suitable tunnel formation to facilitate the substrate entry and product release, resulting in improved catalytic efficiency. With the G175A/A305T double mutant as a biocatalyst, a maximum of 1,220 mM 2-chloronicotinic acid was produced with a 94% conversion, and the space-time yield reached as high as 575 g_product liter^-1 day^-1 These results provide not only a novel robust biocatalyst for the production of 2-chloronicotinic acid but also new insights into amidase structure-function relationships.IMPORTANCE In recent years, the demand for 2-chloronicotinic acid has been greatly increased. To date, several chemical methods have been used for the synthesis of 2-chloronicotinic acid, but all include tedious steps and/or drastic reaction conditions, resulting in both economic and environmental issues. It is requisite to develop an efficient and green synthesis route. We recently screened Pa-Ami and demonstrated its potential for synthesis of 2-chloronicotinic acid from 2-chloronicotinamide. However, chlorine substitution on the pyridine ring of nicotinamide significantly affected the activity of Pa-Ami. Especially for 2-chloronicotinamide, the enzyme activity and catalytic efficiency were relatively low. In this study, based on structure-function analysis, we succeeded in engineering the amidase by structure-guided saturation mutagenesis. The engineered Pa-Ami exhibited quite high catalytic activity toward 2-chloronicotinamide and could serve as a promising biocatalyst for the biosynthesis of 2-chloronicotinic acid.

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.

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.

Single-walled carbon nanotube release affects the microbial enzyme-catalyzed oxidation processes of organic pollutants and lignin model compounds in nature

The question how microbial enzyme-catalyzed oxidation processes of organic pollutants and lignin model compounds (LMCs) are affected by the release of single-walled carbon nanotube (SWCNT) into the environment remains to be addressed at the molecular level. We have, therefore concentrated the effects of SWCNT on some important properties associated with enzyme activity and function during microbial oxidation of polycyclic aromatic hydrocarbons (benzo(a)pyrene, acenaphthene and anthracene), LMCs (2,6-dimethoxyphenol, guaiacol and veratryl alcohol) and β-hexachlorocyclohexane, including the behaviour of water molecules, hydrogen bonds (HBs) and hydrophobic interactions (HYs) between ligand and the enzyme, and conformational dynamics in N- and C-terminus. Our study revealed that SWCNT significantly affected the behaviour of water molecules within 5 Å of both these substrates and their respective enzymes during oxidation (p < 0.01), by increasing or decreasing the water number near them. SWCNT tended to significantly enhance or reduce the stability of atom pairs that formed the HBs and HYs (p < 0.01). N- and C-terminus conformations underwent transitions between positive and negative states or between positive state or between negative state in all analyzed complexes. Significant conformational transitions were found for all C-terminus, but only for a part of N-terminus after the inclusion of the SWCNT. These results showed that SWCNT release would significantly affect the microbial enzyme-catalyzed processes of organic pollutants and LMCs in nature.

Catalytic Promiscuity of Ancestral Esterases and Hydroxynitrile Lyases

Catalytic promiscuity is a useful, but accidental, enzyme property, so finding catalytically promiscuous enzymes in nature is inefficient. Some ancestral enzymes were branch points in the evolution of new enzymes and are hypothesized to have been promiscuous. To test the hypothesis that ancestral enzymes were more promiscuous than their modern descendants, we reconstructed ancestral enzymes at four branch points in the divergence hydroxynitrile lyases (HNL's) from esterases ∼ 100 million years ago. Both enzyme types are α/β-hydrolase-fold enzymes and have the same catalytic triad, but differ in reaction type and mechanism. Esterases catalyze hydrolysis via an acyl enzyme intermediate, while lyases catalyze an elimination without an intermediate. Screening ancestral enzymes and their modern descendants with six esterase substrates and six lyase substrates found higher catalytic promiscuity among the ancestral enzymes (P < 0.01). Ancestral esterases were more likely to catalyze a lyase reaction than modern esterases, and the ancestral HNL was more likely to catalyze ester hydrolysis than modern HNL's. One ancestral enzyme (HNL1) along the path from esterase to hydroxynitrile lyases was especially promiscuous and catalyzed both hydrolysis and lyase reactions with many substrates. A broader screen tested mechanistically related reactions that were not selected for by evolution: decarboxylation, Michael addition, γ-lactam hydrolysis and 1,5-diketone hydrolysis. The ancestral enzymes were more promiscuous than their modern descendants (P = 0.04). Thus, these reconstructed ancestral enzymes are catalytically promiscuous, but HNL1 is especially so.

Dehydrochlorination mechanism of γ-hexachlorocyclohexane degraded by dehydrochlorinase LinA from Sphingomonas paucimobilis UT26

The biotransformation pathway from γ-HCH to 1,3,4,6-TCDN catabolized by dehydrochlorinase LinA contains two discontinuous dehydrochlorination reactions and a conformational transition for the product of the first dehydrochlorination reaction.

Evolution of a Catalytic Mechanism

The means by which superfamilies of specialized enzymes arise by gene duplication and functional divergence are poorly understood. The escape from adaptive conflict hypothesis, which posits multiple copies of a gene encoding a primitive inefficient and highly promiscuous generalist ancestor, receives support from experiments showing that resurrected ancestral enzymes are indeed more substrate-promiscuous than their modern descendants. Here, we provide evidence in support of an alternative model, the innovation-amplification-divergence hypothesis, which posits a single-copied ancestor as efficient and specific as any modern enzyme. We argue that the catalytic mechanisms of plant esterases and descendent acetone cyanohydrin lyases are incompatible with each other (e.g., the reactive substrate carbonyl must bind in opposite orientations in the active site). We then show that resurrected ancestral plant esterases are as catalytically specific as modern esterases, that the ancestor of modern acetone cyanohydrin lyases was itself only very weakly promiscuous, and that improvements in lyase activity came at the expense of esterase activity. These observations support the innovation-amplification-divergence hypothesis, in which an ancestor gains a weak promiscuous activity that is improved by selection at the expense of the ancestral activity, and not the escape from adaptive conflict in which an inefficient generalist ancestral enzyme steadily loses promiscuity throughout the transition to a highly active specialized modern enzyme.

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.

Kinetic and sequence-structure-function analysis of LinB enzyme variants with β- and δ-hexachlorocyclohexane

Organochlorine insecticide hexachlorocyclohexane (HCH) has recently been classified as a ‘Persistent Organic pollutant’ by the Stockholm Convention. The LinB haloalkane dehalogenase is a key upstream enzyme in the recently evolved Lin pathway for the catabolism of HCH in bacteria. Here we report a sequence-structure-function analysis of ten naturally occurring and thirteen synthetic mutants of LinB. One of the synthetic mutants was found to have ∼80 fold more activity for β- and δ-hexachlorocyclohexane. Based on detailed biophysical calculations, molecular dynamics and ensemble docking calculations, we propose that the latter variant is more active because of alterations to the shape of its active site and increased conformational plasticity.

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.