Eva Chovancová

CAVER 3.0: A Tool for the Analysis of Transport Pathways in Dynamic Protein Structures

Tunnels and channels facilitate the transport of small molecules, ions and water solvent in a large variety of proteins. Characteristics of individual transport pathways, including their geometry, physico-chemical properties and dynamics are instrumental for understanding of structure-function relationships of these proteins, for the design of new inhibitors and construction of improved biocatalysts. CAVER is a software tool widely used for the identification and characterization of transport pathways in static macromolecular structures. Herein we present a new version of CAVER enabling automatic analysis of tunnels and channels in large ensembles of protein conformations. CAVER 3.0 implements new algorithms for the calculation and clustering of pathways. A trajectory from a molecular dynamics simulation serves as the typical input, while detailed characteristics and summary statistics of the time evolution of individual pathways are provided in the outputs. To illustrate the capabilities of CAVER 3.0, the tool was applied for the analysis of molecular dynamics simulation of the microbial enzyme haloalkane dehalogenase DhaA. CAVER 3.0 safely identified and reliably estimated the importance of all previously published DhaA tunnels, including the tunnels closed in DhaA crystal structures. Obtained results clearly demonstrate that analysis of molecular dynamics simulation is essential for the estimation of pathway characteristics and elucidation of the structural basis of the tunnel gating. CAVER 3.0 paves the way for the study of important biochemical phenomena in the area of molecular transport, molecular recognition and enzymatic catalysis. The software is freely available as a multiplatform command-line application at http://www.caver.cz.

HotSpot Wizard: a web server for identification of hot spots in protein engineering

HotSpot Wizard is a web server for automatic identification of ‘hot spots’ for engineering of substrate specificity, activity or enantioselectivity of enzymes and for annotation of protein structures. The web server implements the protein engineering protocol, which targets evolutionarily variable amino acid positions located in the active site or lining the access tunnels. The ‘hot spots’ for mutagenesis are selected through the integration of structural, functional and evolutionary information obtained from: (i) the databases RCSB PDB, UniProt, PDBSWS, Catalytic Site Atlas and nr NCBI and (ii) the tools CASTp, CAVER, BLAST, CD-HIT, MUSCLE and Rate4Site. The protein structure and e-mail address are the only obligatory inputs for the calculation. In the output, HotSpot Wizard lists annotated residues ordered by estimated mutability. The results of the analysis are mapped on the enzyme structure and visualized in the web browser using Jmol. The HotSpot Wizard server should be useful for protein engineers interested in exploring the structure of their favourite protein and for the design of mutations in site-directed mutagenesis and focused directed evolution experiments. HotSpot Wizard is available at http://loschmidt.chemi.muni.cz/hotspotwizard/.

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. Proteins 2007.

Substrate specificity of haloalkane dehalogenases

An enzymes substrate specificity is one of its most important characteristics. The quantitative comparison of broad-specificity enzymes requires the selection of a homogenous set of substrates for experimental testing, determination of substrate-specificity data and analysis using multivariate statistics. We describe a systematic analysis of the substrate specificities of nine wild-type and four engineered haloalkane dehalogenases. The enzymes were characterized experimentally using a set of 30 substrates selected using statistical experimental design from a set of nearly 200 halogenated compounds. Analysis of the activity data showed that the most universally useful substrates in the assessment of haloalkane dehalogenase activity are 1-bromobutane, 1-iodopropane, 1-iodobutane, 1,2-dibromoethane and 4-bromobutanenitrile. Functional relationships among the enzymes were explored using principal component analysis. Analysis of the untransformed specific activity data revealed that the overall activity of wild-type haloalkane dehalogenases decreases in the following order: LinB~DbjA>DhlA~DhaA~DbeA~DmbA>DatA~DmbC~DrbA. After transforming the data, we were able to classify haloalkane dehalogenases into four SSGs (substrate-specificity groups). These functional groups are clearly distinct from the evolutionary subfamilies, suggesting that phylogenetic analysis cannot be used to predict the substrate specificity of individual haloalkane dehalogenases. Structural and functional comparisons of wild-type and mutant enzymes revealed that the architecture of the active site and the main access tunnel significantly influences the substrate specificity of these enzymes, but is not its only determinant. The identification of other structural determinants of the substrate specificity remains a challenge for further research on haloalkane dehalogenases.

Software tools for identification, visualization and analysis of protein tunnels and channels.

Protein structures contain highly complex systems of voids, making up specific features such as surface clefts or grooves, pockets, protrusions, cavities, pores or channels, and tunnels. Many of them are essential for the migration of solvents, ions and small molecules through proteins, and their binding to the functional sites. Analysis of these structural features is very important for understanding of structure-function relationships, for the design of potential inhibitors or proteins with improved functional properties. Here we critically review existing software tools specialized in rapid identification, visualization, analysis and design of protein tunnels and channels. The strengths and weaknesses of individual tools are reported together with examples of their applications for the analysis and engineering of various biological systems. This review can assist users with selecting a proper software tool for study of their biological problem as well as highlighting possible avenues for further development of existing tools. Development of novel descriptors representing not only geometry, but also electrostatics, hydrophobicity or dynamics, is needed for reliable identification of biologically relevant tunnels and channels.

Characterization of protein glycosylation in Francisella tularensis subsp. holarctica; identification of a novel glycosylated lipoprotein required for virulence

FTH_0069 is a previously uncharacterized strongly immunoreactive protein that has been proposed to be a novel virulence factor in Francisella tularensis. Here, the glycan structure modifying two C-terminal peptides of FTH_0069 was identified utilizing high resolution, high mass accuracy mass spectrometry, combined with in-source CID tandem MS experiments. The glycan observed at m/z 1156 was determined to be a hexasaccharide, consisting of two hexoses, three N-acetylhexosamines, and an unknown monosaccharide containing a phosphate group. The monosaccharide sequence of the glycan is tentatively proposed as X-P-HexNAc-HexNAc-Hex-Hex-HexNAc, where X denotes the unknown monosaccharide. The glycan is identical to that of DsbA glycoprotein, as well as to one of the multiple glycan structures modifying the type IV pilin PilA, suggesting a common biosynthetic pathway for the protein modification. Here, we demonstrate that the glycosylation of FTH_0069, DsbA, and PilA was affected in an isogenic mutant with a disrupted wbtDEF gene cluster encoding O-antigen synthesis and in a mutant with a deleted pglA gene encoding pilin oligosaccharyltransferase PglA. Based on our findings, we propose that PglA is involved in both pilin and general F. tularensis protein glycosylation, and we further suggest an inter-relationship between the O-antigen and the glycan synthesis in the early steps in their biosynthetic pathways.

A Single Mutation in a Tunnel to the Active Site Changes the Mechanism and Kinetics of Product Release in Haloalkane Dehalogenase LinB

Background: Tunnel properties affect ligand passage in enzymes with buried active sites. Results: A tunnel mutation from leucine to tryptophan changes the mechanism of bromide ion release from haloalkane dehalogenase LinB. Conclusion: Interactions of the bromide ion with the tryptophan increase free energy barrier for its passage, causing the reaction mechanism change. Significance: The results provide guidelines for enzyme engineering. 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.

Biochemical characterization of a novel haloalkane dehalogenase from a cold-adapted bacterium.

ABSTRACT A haloalkane dehalogenase, DpcA, from Psychrobacter cryohalolentis K5, representing a novel psychrophilic member of the haloalkane dehalogenase family, was identified and biochemically characterized. DpcA exhibited a unique temperature profile with exceptionally high activities at low temperatures. The psychrophilic properties of DpcA make this enzyme promising for various environmental applications.