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.

Phylogenetic analysis of haloalkane dehalogenases

Haloalkane dehalogenases (HLDs) are enzymes that catalyze the cleavage of carbon-halogen bonds by a hydrolytic mechanism. Although comparative biochemical analyses have been published, no classification system has been proposed for HLDs, to date, that reconciles their phylogenetic and functional relationships. In the study presented here, we have analyzed all sequences and structures of genuine HLDs and their homologs detectable by database searches. Phylogenetic analyses revealed that the HLD family can be divided into three subfamilies denoted HLD-I, HLD-II, and HLD-III, of which HLD-I and HLD-III are predicted to be sister-groups. A mismatch between the HLD protein tree and the tree of species, as well as the presence of more than one HLD gene in a few genomes, suggest that horizontal gene transfers, and perhaps also multiple gene duplications and losses have been involved in the evolution of this family. Most of the biochemically characterized HLDs are found in the HLD-II subfamily. The dehalogenating activity of two members of the newly identified HLD-III subfamily has only recently been confirmed, in a study motivated by this phylogenetic analysis. A novel type of the catalytic pentad (Asp-His-Asp+Asn-Trp) was predicted for members of the HLD-III subfamily. Calculation of the evolutionary rates and lineage-specific innovations revealed a common conserved core as well as a set of residues that characterizes each HLD subfamily. The N-terminal part of the cap domain is one of the most variable regions within the whole family as well as within individual subfamilies, and serves as a preferential site for the location of relatively long insertions. The highest variability of discrete sites was observed among residues that are structural components of the access channels. Mutations at these sites modify the anatomy of the channels, which are important for the exchange of ligands between the buried active site and the bulk solvent, thus creating a structural basis for the molecular evolution of new substrate specificities. Our analysis sheds light on the evolutionary history of HLDs and provides a structural framework for designing enzymes with new specificities.