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

Comprehensive review on Haloalkane dehalogenase (LinB): a β-hexachlorocyclohexane (HCH) degrading enzyme

Haloalkane dehalogenase, LinB, is a member of the α/β hydrolase family of enzymes. It has a wide range of halogenated substrates, but, has been mostly studied in context of degradation of hexachlorocyclohexane (HCH) isomers, especially β-HCH (5-12% of total HCH isomers), which is the most recalcitrant and persistent among all the HCH isomers. LinB was identified to directly act on β-HCH in a one or two step transformation which decreases its toxicity manifold. Thereafter, many studies focused on LinB including its structure determination using X-ray crystallographic studies, structure comparison with other haloalkane dehalogenases, substrate specificity and kinetic studies, protein engineering and site-directed mutagenesis studies in search of better catalytic activity of the enzyme. LinB was mainly identified and characterized in bacteria belonging to sphingomonads. Detailed sequence comparison of LinB from different sphingomonads further revealed the residues critical for its activity and ability to catalyze either one or two step transformation of β-HCH. Association of LinB with IS6100 elements is also being discussed in detail in sphingomonads. In this review, we summarized vigorous efforts done by different research groups on LinB for developing better bioremediation strategies against HCH contamination. Also, kinetic studies, protein engineering and site directed mutagenesis studies discussed here forms the basis of further exploration of LinB's role as an efficient enzyme in bioremediation projects.

Epistasis facilitates functional evolution in an ancient transcription factor

A protein's genetic architecture - the set of causal rules by which its sequence produces its functions - also determines its possible evolutionary trajectories. Prior research has proposed that the genetic architecture of proteins is very complex, with pervasive epistatic interactions that constrain evolution and make function difficult to predict from sequence. Most of this work has analyzed only the direct paths between two proteins of interest - excluding the vast majority of possible genotypes and evolutionary trajectories - and has considered only a single protein function, leaving unaddressed the genetic architecture of functional specificity and its impact on the evolution of new functions. Here, we develop a new method based on ordinal logistic regression to directly characterize the global genetic determinants of multiple protein functions from 20-state combinatorial deep mutational scanning (DMS) experiments. We use it to dissect the genetic architecture and evolution of a transcription factor's specificity for DNA, using data from a combinatorial DMS of an ancient steroid hormone receptor's capacity to activate transcription from two biologically relevant DNA elements. We show that the genetic architecture of DNA recognition consists of a dense set of main and pairwise effects that involve virtually every possible amino acid state in the protein-DNA interface, but higher-order epistasis plays only a tiny role. Pairwise interactions enlarge the set of functional sequences and are the primary determinants of specificity for different DNA elements. They also massively expand the number of opportunities for single-residue mutations to switch specificity from one DNA target to another. By bringing variants with different functions close together in sequence space, pairwise epistasis therefore facilitates rather than constrains the evolution of new functions.

Genome evolution related to γ-hexachlorocyclohexane metabolic function in the soil microbial population

γ-Hexachlorocyclohexane (γ-HCH)-degrading strain, Sphingobium sp. TA15, was newly isolated from an experimental field soil from which the archetypal γ-HCH-degrading strain, S. japonicum UT26, was isolated previously. Comparison of the complete genome sequences of these 2 strains revealed that TA15 shares the same basic genome backbone with UT26, but also has the variable regions that are presumed to have changed either from UT26 or from a putative common ancestor. Organization and localization of lin genes of TA15 were different from those of UT26. It was inferred that transposition of IS6100 had played a crucial role in these genome rearrangements. The accumulation of toxic dead-end products in TA15 was lower than in UT26, suggesting that TA15 utilizes γ-HCH more effectively than UT26. These results suggested that genome evolution related to the γ-HCH metabolic function in the soil microbial population is ongoing.

A mechanistic view of enzyme evolution

New enzyme functions often evolve through the recruitment and optimization of latent promiscuous activities. How do mutations alter the molecular architecture of enzymes to enhance their activities? Can we infer general mechanisms that are common to most enzymes, or does each enzyme require a unique optimization process? The ability to predict the location and type of mutations necessary to enhance an enzyme's activity is critical to protein engineering and rational design. In this review, via the detailed examination of recent studies that have shed new light on the molecular changes underlying the optimization of enzyme function, we provide a mechanistic perspective of enzyme evolution. We first present a global survey of the prevalence of activity-enhancing mutations and their distribution within protein structures. We then delve into the molecular solutions that mediate functional optimization, specifically highlighting several common mechanisms that have been observed across multiple examples. As distinct protein sequences encounter different evolutionary bottlenecks, different mechanisms are likely to emerge along evolutionary trajectories toward improved function. Identifying the specific mechanism(s) that need to be improved upon, and tailoring our engineering efforts to each sequence, may considerably improve our chances to succeed in generating highly efficient catalysts in the future.

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.

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.

Comparison of the complete genome sequences of four γ-hexachlorocyclohexane-degrading bacterial strains: insights into the evolution of bacteria able to degrade a recalcitrant man-made pesticide

γ-Hexachlorocyclohexane (γ-HCH) is a recalcitrant man-made chlorinated pesticide. Here, the complete genome sequences of four γ-HCH-degrading sphingomonad strains, which are most unlikely to have been derived from one ancestral γ-HCH degrader, were compared. Together with several experimental data, we showed that (i) all the four strains carry almost identical linA to linE genes for the conversion of γ-HCH to maleylacetate (designated “specific” lin genes), (ii) considerably different genes are used for the metabolism of maleylacetate in one of the four strains, and (iii) the linKLMN genes for the putative ABC transporter necessary for γ-HCH utilization exhibit structural divergence, which reflects the phylogenetic relationship of their hosts. Replicon organization and location of the lin genes in the four genomes are significantly different with one another, and that most of the specific lin genes are located on multiple sphingomonad-unique plasmids. Copies of IS6100, the most abundant insertion sequence in the four strains, are often located in close proximity to the specific lin genes. Analysis of the footprints of target duplication upon IS6100 transposition and the experimental detection of IS6100 transposition strongly suggested that the IS6100 transposition has caused dynamic genome rearrangements and the diversification of lin-flanking regions in the four strains.

How mutational epistasis impairs predictability in protein evolution and design

There has been much debate about the extent to which mutational epistasis, that is, the dependence of the outcome of a mutation on the genetic background, constrains evolutionary trajectories. The degree of unpredictability introduced by epistasis, due to the non-additivity of functional effects, strongly hinders the strategies developed in protein design and engineering. While many studies have addressed this issue through systematic characterization of evolutionary trajectories within individual enzymes, the field lacks a consensus view on this matter. In this work, we performed a comprehensive analysis of epistasis by analyzing the mutational effects from nine adaptive trajectories toward new enzymatic functions. We quantified epistasis by comparing the effect of mutations occurring between two genetic backgrounds: the starting enzyme (for example, wild type) and the intermediate variant on which the mutation occurred during the trajectory. We found that most trajectories exhibit positive epistasis, in which the mutational effect is more beneficial when it occurs later in the evolutionary trajectory. Approximately half (49%) of functional mutations were neutral or negative on the wild-type background, but became beneficial at a later stage in the trajectory, indicating that these functional mutations were not predictable from the initial starting point. While some cases of strong epistasis were associated with direct interaction between residues, many others were caused by long-range indirect interactions between mutations. Our work highlights the prevalence of epistasis in enzyme adaptive evolution, in particular positive epistasis, and suggests the necessity of incorporating mutational epistasis in protein engineering and design to create highly efficient catalysts.