Functionally relevant motions of haloalkane dehalogenases occur in the specificity-modulating cap domains

One-nanosecond molecular dynamics trajectories of three haloalkane dehalogenases (DhlA, LinB, and DhaA) are compared. The main domain was rigid in all three dehalogenases, whereas the substrate specificity-modulating cap domains showed considerably higher mobility. The functionally relevant motions were spread over the entire cap domain in DhlA, whereas they were more localized in LinB and DhaA. The highest amplitude of essential motions of DhlA was noted in the alpha4'-helix-loop-alpha4-helix region, formerly proposed to participate in the large conformation change needed for product release. The highest amplitude of essential motions of LinB and DhaA was observed in the random coil before helix 4, linking two domains of these proteins. This flexibility is the consequence of the modular composition of haloalkane dehalogenases. Two members of the catalytic triad, that is, the nucleophile and the base, showed a very high level of rigidity in all three dehalogenases. This rigidity is essential for their function. One of the halide-stabilizing residues, important for the catalysis, shows significantly higher flexibility in DhlA compared with LinB and DhaA. Enhanced flexibility may be required for destabilization of the electrostatic interactions during the release of the halide ion from the deeply buried active site of DhlA. The exchange of water molecules between the enzyme active site and bulk solvent was very different among the three dehalogenases. The differences could be related to the flexibility of the cap domains and to the number of entrance tunnels.     

Catalytic mechanism of the maloalkane dehalogenase LinB from Sphingomonas paucimobilis UT26

Haloalkane dehalogenases are bacterial enzymes capable of carbon-halogen bond cleavage in halogenated compounds. To obtain insights into the mechanism of the haloalkane dehalogenase from Sphingomonas paucimobilis UT26 (LinB), we studied the steady-state and presteady-state kinetics of the conversion of the substrates 1-chlorohexane, chlorocyclohexane, and bromocyclohexane. The results lead to a proposal of a minimal kinetic mechanism consisting of three main steps: (i) substrate binding, (ii) cleavage of the carbon-halogen bond with simultaneous formation of an alkyl-enzyme intermediate, and (iii) hydrolysis of the alkyl-enzyme intermediate. Release of both products, halide and alcohol, is a fast process that was not included in the reaction mechanism as a distinct step. Comparison of the kinetic mechanism of LinB with that of haloalkane dehalogenase DhlA from Xantobacter autotrophicus GJ10 and the haloalkane dehalogenase DhaA from Rhodococcus rhodochrous NCIMB 13064 shows that the overall mechanisms are similar. The main difference is in the rate-limiting step, which is hydrolysis of the alkylenzyme intermediate in LinB, halide release in DhlA, and liberation of an alcohol in DhaA. The occurrence of different rate-limiting steps for three enzymes that belong to the same protein family indicates that extrapolation of this important catalytic property from one enzyme to another can be misleading even for evolutionary closely related proteins. The differences in the rate-limiting step were related to: (i) number and size of the entrance tunnels, (ii) protein flexibility, and (iii) composition of the halide-stabilizing active site residues based on comparison of protein structures.     

Halide-stabilizing residues of haloalkane dehalogenases studied by quantum mechanic calculations and site-directed mutagenesis

Haloalkane dehalogenases catalyze cleavage of the carbon-halogen bond in halogenated aliphatic compounds, resulting in the formation of an alcohol, a halide, and a proton as the reaction products. Three structural features of haloalkane dehalogenases are essential for their catalytic performance: (i) a catalytic triad, (ii) an oxyanion hole, and (iii) the halide-stabilizing residues. Halide-stabilizing residues are not structurally conserved among different haloalkane dehalogenases. The level of stabilization of the transition state structure of S(N)2 reaction and halide ion provided by each of the active site residues in the enzymes DhlA, LinB, and DhaA was quantified by quantum mechanic calculations. The residues that significantly stabilize the halide ion were assigned as the primary (essential) or the secondary (less important) halide-stabilizing residues. Site-directed mutagenesis was conducted with LinB enzyme to confirm location of its primary halide-stabilizing residues. Asn38Asp, Asn38Glu, Asn38Phe, Asn38Gln, Trp109Leu, Phe151Leu, Phe151Trp, Phe151Tyr, and Phe169Leu mutants of LinB were constructed, purified, and kinetically characterized. The following active site residues were classified as the primary halide-stabilizing residues: Trp125 and Trp175 of DhlA; Asn38 and Trp109 of LinB; and Asn41 and Trp107 of DhaA. All these residues make a hydrogen bond with the halide ion released from the substrate molecule, and their substitution results in enzymes with significantly modified catalytic properties. The following active site residues were classified as the secondary halide-stabilizing residues: Phe172, Pro223, and Val226 of DhlA; Trp207, Pro208, and Ile211 of LinB; and Phe205, Pro206, and Ile209 of DhaA. The differences in the halide stabilizing residues of three haloalkane dehalogenases are discussed in the light of molecular adaptation of these enzymes to their substrates.     

Exploring the binding site structure of the PPAR gamma ligand-binding domain by computational solvent mapping

Solvent mapping moves molecular probes, small organic molecules containing various functional groups, around the protein surface, finds favorable positions, clusters the conformations, and ranks the clusters based on the average free energy. Using at least six different solvents as probes, the probes cluster in major pockets of the functional site, providing detailed and reliable information on the amino acid residues that are important for ligand binding. Solvent mapping was applied to 12 structures of the peroxisome proliferator activated receptor gamma (PPARgamma) ligand-binding domain (LBD), including 2 structures without a ligand, 2 structures with a partial agonist, and 8 structures with a PPAR agonist bound. The analysis revealed 10 binding "hot spots", 4 in the ligand-binding pocket, 2 in the coactivator-binding region, 1 in the dimerization domain, 2 around the ligand entrance site, and 1 minor site without a known function. Mapping is a major source of information on the role and cooperativity of these sites. It shows that large portions of the ligand-binding site are already formed in the PPARgamma apostructure, but an important pocket near the AF-2 transactivation domain becomes accessible only in structures that are cocrystallized with strong agonists. Conformational changes were seen in several other sites, including one involved in the stabilization of the LBD and two others at the region of the coactivator binding. The number of probe clusters retained by these sites depends on the properties of the bound agonist, providing information on the origin of correlations between ligand and coactivator binding.     

Comparative surface geometry of the protein kinase family

Identifying conserved pockets on the surfaces of a family of proteins can provide insight into conserved geometric features and sites of protein–protein interaction. Here we describe mapping and comparison of the surfaces of aligned crystallographic structures, using the protein kinase family as a model. Pockets are rapidly computed using two computer programs, FADE and Crevasse. FADE uses gradients of atomic density to locate grooves and pockets on the molecular surface. Crevasse, a new piece of software, splits the FADE output into distinct pockets. The computation was run on 10 kinase catalytic cores aligned on the αF‐helix, and the resulting pockets spatially clustered. The active site cleft appears as a large, contiguous site that can be subdivided into nucleotide and substrate docking sites. Substrate specificity determinants in the active site cleft between serine/threonine and tyrosine kinases are visible and distinct. The active site clefts cluster tightly, showing a conserved spatial relationship between the active site and αF‐helix in the C‐lobe. When the αC‐helix is examined, there are multiple mechanisms for anchoring the helix using spatially conserved docking sites. A novel site at the top of the N‐lobe is present in all the kinases, and there is a large conserved pocket over the hinge and the αC‐β4 loop. Other pockets on the kinase core are strongly conserved but have not yet been mapped to a protein–protein interaction. Sites identified by this algorithm have revealed structural and spatially conserved features of the kinase family and potential conserved intermolecular and intramolecular binding sites.     

Crystallographic and kinetic evidence of a collision complex formed during halide import in haloalkane dehalogenase.

Haloalkane dehalogenase (DhlA) converts haloalkanes to their corresponding alcohols and halide ions. The rate-limiting step in the reaction of DhlA is the release of the halide ion. The kinetics of halide release have been analyzed by measuring halide binding with stopped-flow fluorescence experiments. At high halide concentrations, halide import occurs predominantly via the rapid formation of a weak initial collision complex, followed by transport of the ion to the active site. To obtain more insight in this collision complex, we determined the X-ray structure of DhlA in the presence of bromide and investigated the kinetics of mutants that were constructed on the basis of this structure. The X-ray structure revealed one bromide ion firmly bound in the active site and two bromide ions weakly bound on the surface of the enzyme. One of the weakly bound ions is close to Thr197 and Phe294, near the entrance of the earlier proposed tunnel for substrate import. Kinetic analysis of bromide import by the Thr197Ala and Phe294Ala mutants of DhlA at high halide concentration showed that the rate constants for halide binding no longer displayed a wild-type-like parabolic increase with increasing bromide concentrations. This is in agreement with an elimination or a decrease in affinity of the surface-located halide-binding site. Likewise, chloride binding kinetics of the mutants indicated significant differences with wild-type enzyme. The results indicate that Thr197 and Phe294 are involved in the formation of an initial collision complex for halide import in DhlA and provide experimental evidence for the role of the tunnel in substrate and product transport.     

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.

This is a PLOS Computational Biology Software Article

     

Thermodynamic analysis of halide binding to haloalkane dehalogenase suggests the occurrence of large conformational changes

Haloalkane dehalogenase (DhlA) hydrolyzes short-chain haloalkanes to produce the corresponding alcohols and halide ions. Release of the halide ion from the active-site cavity can proceed via a two-step and a three-step route, which both contain slow enzyme isomerization steps. Thermodynamic analysis of bromide binding and release showed that the slow unimolecular isomerization steps in the three-step bromide export route have considerably larger transition state enthalpies and entropies than those in the other route. This suggests that the three-step route involves different and perhaps larger conformational changes than the two-step export route. We propose that the three-step halide export route starts with conformational changes that result in a more open configuration of the active site from which the halide ion can readily escape. In addition, we suggest that the two-step route for halide release involves the transfer of the halide ion from the halide-binding site in the cavity to a binding site somewhere at the protein surface, where a so-called collision complex is formed in which the halide ion is only weakly bound. No large structural rearrangements are necessary for this latter process.     

Identification of substrate binding sites in enzymes by computational solvent mapping

Enzyme structures determined in organic solvents show that most organic molecules cluster in the active site, delineating the binding pocket. We have developed algorithms to perform solvent mapping computationally, rather than experimentally, by placing molecular probes (small molecules or functional groups) on a protein surface, and finding the regions with the most favorable binding free energy. The method then finds the consensus site that binds the highest number of different probes. The probe-protein interactions at this site are compared to the intermolecular interactions seen in the known complexes of the enzyme with various ligands (substrate analogs, products, and inhibitors). We have mapped thermolysin, for which experimental mapping results are also available, and six further enzymes that have no experimental mapping data, but whose binding sites are well characterized. With the exception of haloalkane dehalogenase, which binds very small substrates in a narrow channel, the consensus site found by the mapping is always a major subsite of the substrate-binding site. Furthermore, the probes at this location form hydrogen bonds and non-bonded interactions with the same residues that interact with the specific ligands of the enzyme. Thus, once the structure of an enzyme is known, computational solvent mapping can provide detailed and reliable information on its substrate-binding site. Calculations on ligand-bound and apo structures of enzymes show that the mapping results are not very sensitive to moderate variations in the protein coordinates.     

Oligosaccharide binding to barley alpha-amylase 1

Enzymatic subsite mapping earlier predicted 10 binding subsites in the active site substrate binding cleft of barley alpha-amylase isozymes. The three-dimensional structures of the oligosaccharide complexes with barley alpha-amylase isozyme 1 (AMY1) described here give for the first time a thorough insight into the substrate binding by describing residues defining 9 subsites, namely -7 through +2. These structures support that the pseudotetrasaccharide inhibitor acarbose is hydrolyzed by the active enzymes. Moreover, sugar binding was observed to the starch granule-binding site previously determined in barley alpha-amylase isozyme 2 (AMY2), and the sugar binding modes are compared between the two isozymes. The "sugar tongs" surface binding site discovered in the AMY1-thio-DP4 complex is confirmed in the present work. A site that putatively serves as an entrance for the substrate to the active site was proposed at the glycone part of the binding cleft, and the crystal structures of the catalytic nucleophile mutant (AMY1D180A) complexed with acarbose and maltoheptaose, respectively, suggest an additional role for the nucleophile in the stabilization of the Michaelis complex. Furthermore, probable roles are outlined for the surface binding sites. Our data support a model in which the two surface sites in AMY1 can interact with amylose chains in their naturally folded form. Because of the specificities of these two sites, they may locate/orient the enzyme in order to facilitate access to the active site for polysaccharide chains. Moreover, the sugar tongs surface site could also perform the unraveling of amylose chains, with the aid of Tyr-380 acting as "molecular tweezers."     

Reconstruction of mycobacterial dehalogenase Rv2579 by cumulative mutagenesis of haloalkane dehalogenase LinB

The homology model of protein Rv2579 from Mycobacterium tuberculosis H37Rv was compared with the crystal structure of haloalkane dehalogenase LinB from Sphingomonas paucimobilis UT26, and this analysis revealed that 6 of 19 amino acid residues which form an active site and entrance tunnel are different in LinB and Rv2579. To characterize the effect of replacement of these six amino acid residues, mutations were introduced cumulatively into the six amino acid residues of LinB. The sixfold mutant, which was supposed to have the active site of Rv2579, exhibited haloalkane dehalogenase activity with the haloalkanes tested, confirming that Rv2579 is a member of the haloalkane dehalogenase protein family.     

Computational solvent mapping reveals the importance of local conformational changes for broad substrate specificity in mammalian cytochromes P450

Computational solvent mapping moves small organic molecules as probes around a protein surface, finds favorable binding positions, clusters the conformations, and ranks the clusters on the basis of their average free energy. Prior mapping studies of enzymes, crystallized in either substrate-free or substrate-bound form, have shown that the largest number of solvent probe clusters invariably overlaps in the active site. We have applied this method to five cytochromes P450. As expected, the mapping of two bacterial P450s, P450 cam (CYP101) and P450 BM-3 (CYP102), identified the substrate-binding sites in both ligand-bound and ligand-free P450 structures. However, the mapping finds the active site only in the ligand-bound structures of the three mammalian P450s, 2C5, 2C9, and 2B4. Thus, despite the large cavities seen in the unbound structures of these enzymes, the features required for binding small molecules are formed only in the process of substrate binding. The ability of adjusting their binding sites to substrates that differ in size, shape, and polarity is likely to be responsible for the broad substrate specificity of these mammalian P450s. Similar behavior was seen at "hot spots" of protein-protein interfaces that can also bind small molecules in grooves created by induced fit. In addition, the binding of S-warfarin to P450 2C9 creates a high-affinity site for a second ligand, which may help to explain the prevalence of drug-drug interactions involving this and other mammalian P450s.     

Reconstruction of Mycobacterial Dehalogenase Rv2579 by Cumulative Mutagenesis of Haloalkane Dehalogenase LinB

The homology model of protein Rv2579 from Mycobacterium tuberculosis H37Rv was compared with the crystal structure of haloalkane dehalogenase LinB from Sphingomonas paucimobilis UT26, and this analysis revealed that 6 of 19 amino acid residues which form an active site and entrance tunnel are different in LinB and Rv2579. To characterize the effect of replacement of these six amino acid residues, mutations were introduced cumulatively into the six amino acid residues of LinB. The sixfold mutant, which was supposed to have the active site of Rv2579, exhibited haloalkane dehalogenase activity with the haloalkanes tested, confirming that Rv2579 is a member of the haloalkane dehalogenase protein family.     

Computational selection, identification and structural analysis of ω-aminotransferases with various substrate specificities from the genome sequence of Mesorhizobium loti MAFF303099

ω-Aminotransferase (ω-AT) is an important class of enzymes for the synthesis of chiral amines or β-amino acids. Family profile analysis was applied to screen putative ω-ATs from Mesorhizobium loti MAFF303099, a nitrogen fixation bacterium that has a larger number of ATs than other microorganisms. By family profile analysis, we selected 10 putative ω-ATs according to E-value. The functions of the putative ω-ATs were investigated by examining activities towards amines and/or β-amino acids. 10 putative proteins were found to have ω-AT activity with narrow or broad substrate specificity. Structure analysis using crystal structure of mll7127 and homology models of mll1632 and mll3663 indicated that the structures of active sites of the enzymes were very similar and highly conserved, but their substrate specificities appeared to be determined by residues positioned at the entrance region of the active site binding pockets.     

Anatomy of protein pockets and cavities: measurement of binding site geometry and implications for ligand design.

Identification and size characterization of surface pockets and occluded cavities are initial steps in protein structure-based ligand design. A new program, CAST, for automatically locating and measuring protein pockets and cavities, is based on precise computational geometry methods, including alpha shape and discrete flow theory. CAST identifies and measures pockets and pocket mouth openings, as well as cavities. The program specifies the atoms lining pockets, pocket openings, and buried cavities; the volume and area of pockets and cavities; and the area and circumference of mouth openings. CAST analysis of over 100 proteins has been carried out; proteins examined include a set of 51 monomeric enzyme-ligand structures, several elastase-inhibitor complexes, the FK506 binding protein, 30 HIV-1 protease-inhibitor complexes, and a number of small and large protein inhibitors. Medium-sized globular proteins typically have 10-20 pockets/cavities. Most often, binding sites are pockets with 1-2 mouth openings; much less frequently they are cavities. Ligand binding pockets vary widely in size, most within the range 10(2)-10(3)A3. Statistical analysis reveals that the number of pockets and cavities is correlated with protein size, but there is no correlation between the size of the protein and the size of binding sites. Most frequently, the largest pocket/cavity is the active site, but there are a number of instructive exceptions. Ligand volume and binding site volume are somewhat correlated when binding site volume is < or =700 A3, but the ligand seldom occupies the entire site. Auxiliary pockets near the active site have been suggested as additional binding surface for designed ligands (Mattos C et al., 1994, Nat Struct Biol 1:55-58). Analysis of elastase-inhibitor complexes suggests that CAST can identify ancillary pockets suitable for recruitment in ligand design strategies. Analysis of the FK506 binding protein, and of compounds developed in SAR by NMR (Shuker SB et al., 1996, Science 274:1531-1534), indicates that CAST pocket computation may provide a priori identification of target proteins for linked-fragment design. CAST analysis of 30 HIV-1 protease-inhibitor complexes shows that the flexible active site pocket can vary over a range of 853-1,566 A3, and that there are two pockets near or adjoining the active site that may be recruited for ligand design.     

Crystal structure of the haloalkane dehalogenase from Sphingomonas paucimobilis UT26

The haloalkane dehalogenase from Sphingomonas paucimobilis UT26 (LinB) is the enzyme involved in the degradation of the important environmental pollutant gamma-hexachlorocyclohexane. The enzyme hydrolyzes a broad range of halogenated cyclic and aliphatic compounds. Here, we present the 1.58 A crystal structure of LinB and the 2.0 A structure of LinB with 1,3-propanediol, a product of debromination of 1,3-dibromopropane, in the active site of the enzyme. The enzyme belongs to the alpha/beta hydrolase family and contains a catalytic triad (Asp108, His272, and Glu132) in the lipase-like topological arrangement previously proposed from mutagenesis experiments. The LinB structure was compared with the structures of haloalkane dehalogenase from Xanthobacter autotrophicus GJ10 and from Rhodococcus sp. and the structural features involved in the adaptation toward xenobiotic substrates were identified. The arrangement and composition of the alpha-helices in the cap domain results in the differences in the size and shape of the active-site cavity and the entrance tunnel. This is the major determinant of the substrate specificity of this haloalkane dehalogenase.     

FTMove: A Web Server for Detection and Analysis of Cryptic and Allosteric Binding Sites by Mapping Multiple Protein Structures

Protein mapping distributes many copies of different molecular probes on the surface of a target protein in order to determine binding hot spots, regions that are highly preferable for ligand binding. While mapping of X-ray structures by the FTMap server is inherently static, this limitation can be overcome by the simultaneous analysis of multiple structures of the protein. FTMove is an automated web server that implements this approach. From the input of a target protein, by PDB code, the server identifies all structures of the protein available in the PDB, runs mapping on them, and combines the results to form binding hot spots and binding sites. The user may also upload their own protein structures, bypassing the PDB search for similar structures. Output of the server consists of the consensus binding sites and the individual mapping results for each structure - including the number of probes located in each binding site, for each structure. This level of detail allows the users to investigate how the strength of a binding site relates to the protein conformation, other binding sites, and the presence of ligands or mutations. In addition, the structures are clustered on the basis of their binding properties. The use of FTMove is demonstrated by application to 22 proteins with known allosteric binding sites; the orthosteric and allosteric binding sites were identified in all but one case, and the sites were typically ranked among the top five. The FTMove server is publicly available at https://ftmove.bu.edu.     

Computational Selection,Identification and Structural Analysis of ω-Aminotransferases with Various Substrate Specificities from the Genome Sequence of Mesorhizobium loti MAFF303099

ω-Aminotransferase (ω-AT) is an important class of enzymes for the synthesis of chiral amines or β-amino acids. Family profile analysis was applied to screen putative ω-ATs from Mesorhizobium loti MAFF303099, a nitrogen fixation bacterium that has a larger number of ATs than other microorganisms. By family profile analysis, we selected 10 putative ω-ATs according to E-value. The functions of the putative ω-ATs were investigated by examining activities towards amines and/or β-amino acids. 10 putative proteins were found to have ω-AT activity with narrow or broad substrate specificity. Structure analysis using crystal structure of mll7127 and homology models of mll1632 and mll3663 indicated that the structures of active sites of the enzymes were very similar and highly conserved, but their substrate specificities appreared to be determined by residues positioned at the entrance region of the active site binding pockets.     

Unveiling Prolyl Oligopeptidase Ligand Migration by Comprehensive Computational Techniques

Prolyl oligopeptidase (POP) is a large 80 kDa protease, which cleaves oligopeptides at the C-terminal side of proline residues and constitutes an important pharmaceutical target. Despite the existence of several crystallographic structures, there is an open debate about migration (entrance and exit) pathways for ligands, and their coupling with protein dynamics. Recent studies have shown the capabilities of molecular dynamics and classical force fields in describing spontaneous binding events and nonbiased ligand migration pathways. Due to POP’s size and to the buried nature of its active site, an exhaustive sampling by means of conventional long enough molecular dynamics trajectories is still a nearly impossible task. Such a level of sampling, however, is possible with the breakthrough protein energy landscape exploration technique. Here, we present an exhaustive sampling of POP with a known inhibitor, Z-pro-prolinal. In >3000 trajectories Z-pro-prolinal explores all the accessible surface area, showing multiple entrance events into the large internal cavity through the pore in the β-propeller domain. Moreover, we modeled a natural substrate binding and product release by predicting the entrance of an undecapeptide substrate, followed by manual active site cleavage and nonbiased exit of one of the products (a dipeptide). The product exit shows preference from a flexible 18-amino acid residues loop, pointing to an overall mechanism where entrance and exit occur in different sites.