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.     

How optimal are the binding energetics of barnase and barstar?

The extracellular ribonuclease barnase and its intracellular inhibitor barstar bind fast and with high affinity. Although extensive experimental and theoretical studies have been carried out on this system, it is unclear what the relative importance of different contributions to the high affinity is and whether binding can be improved through point mutations. In this work, we first applied Poisson-Boltzmann electrostatic calculations to 65 barnase-barstar complexes with mutations in both barnase and barstar. The continuum electrostatic calculations with a van der Waals surface dielectric boundary definition result in the electrostatic interaction free energy providing the dominant contribution favoring barnase-barstar binding. The results show that the computed electrostatic binding free energy can be improved through mutations at W44/barstar and E73/barnase. Furthermore, the determinants of binding affinity were quantified by applying COMparative BINding Energy (COMBINE) analysis to derive quantitative structure-activity relationships (QSARs) for the 65 complexes. The COMBINE QSAR model highlights approximately 20 interfacial residue pairs as responsible for most of the differences in binding affinity between the mutant complexes, mainly due to electrostatic interactions. Based on the COMBINE model, together with Brownian dynamics simulations to compute diffusional association rate constants, several mutants were designed to have higher binding affinities than the wild-type proteins.     

The nature of intermolecular interactions between aromatic amino acid residues

The nature of intermolecular interactions between aromatic amino acid residues has been investigated by a combination of molecular dynamics and ab initio methods. The potential energy surface of various interacting pairs, including tryptophan, phenilalanine, and tyrosine, was scanned for determining all the relevant local minima by a combined molecular dynamics and conjugate gradient methodology with the AMBER force field. For each of these minima, single-point correlated ab initio calculations of the binding energy were performed. The agreement between empirical force field and ab initio binding energies of the minimum energy structures is excellent. Aromatic-aromatic interactions can be rationalized on the basis of electrostatic and van der Waals interactions, whereas charge transfer or polarization phenomena are small for all intermolecular complexes and, particularly, for stacked structures. Proteins 2002;48:117-125.     

Heterologous Expression inEscherichia coliof Soluble Active-Site Random Mutants of Haloalkane Dehalogenase fromXanthobacter autotrophicusGJ10 by Coexpression of Molecular Chaperonins GroEL/ES

A system for heterologous expression inEscherichia coliof dehaloalkane dehalogenase DhlA fromXanthobacter autotrophicusstrain GJ10 is presented. The strategy involved overexpression ofE. colichaperonins GroEL/ES which facilitated the production of soluble DhlA. When active-site mutant forms were constructed they could not to any detectable degree be expressed in a soluble state in the absence of overproduced GroEL/ES. However, with the described expression system, wild-type DhlA as well as variant forms randomly mutated in the active-site residues Phe172 and Trp175 were reliably produced. An introduced C-terminal (His)₅-tag provided an immunological handle as well as a site for metal ion coordination utilized in affinity chromatography for the purification of recombinant DhlA. The purified His-tagged enzyme, DhlA-5His, was confirmed to be catalytically fully active when measuring the dehalogenase activity with dichloroethane as substrate.     

Exploring the binding of BACE-1 inhibitors using comparative binding energy analysis (COMBINE)

ABSTRACT: BACKGROUND: The inhibition of the activity of beta-secretase (BACE-1) is a potentially important approach for the treatment of Alzheimer disease. To explore the mechanism of inhibition, we describe the use of 46 X-ray crystallographic BACE-1/inhibitor complexes to derive quantitative structure-activity relationship (QSAR) models. The inhibitors were aligned by superimposing 46 X-ray crystallographic BACE-1/inhibitor complexes, and gCOMBINE software was used to perform COMparative BINding Energy (COMBINE) analysis on these 46 minimized BACE-1/inhibitor complexes. The major advantage of the COMBINE analysis is that it can quantitatively extract key residues involved in binding the ligand and identify the nature of the interactions between the ligand and receptor. RESULTS: By considering the contributions of the protein residues to the electrostatic and van der Waals intermolecular interaction energies, two predictive and robust COMBINE models were developed: (i) the 3-PC distance-dependent dielectric constant model (built from a single X-ray crystal structure) with a q2 value of 0.74 and an SDEC value of 0.521; and (ii) the 5-PC sigmoidal electrostatic model (built from the actual complexes present in the Brookhaven Protein Data Bank) with a q2 value of 0.79 and an SDEC value of 0.41. CONCLUSIONS: These QSAR models and the information describing the inhibition provide useful insights into the design of novel inhibitors via the optimization of the interactions between ligands and those key residues of BACE-1.     

Energetics of complex formation of the DNA hairpin structure d(GCGAAGC) with aromatic ligands

The energy contributions of various physical interactions to the total Gibbs energy of complex formation of the biologically important DNA hairpin d(GCGAAGC) with aromatic antitumor antibiotics daunomycin and novantron and the mutagens ethidium and proflavine have been calculated. It has been shown that the relatively small value of the total energy of binding of the ligands to the hairpin is the sum of components great in absolute value and different in sign. The contributions of van der Waals interactions and both intra- and intermolecular hydrogen bonds and bonds with aqueous environment have been studied. According to the calculations, the hydrophobic and van der Waals components are energetically favorable in complex formation of the ligands with the DNA pairpin d(GCGAAGC), whereas the electrostatic (with consideration of hydrogen bonds) and entropic components are unfavorable.     

Energetics of complex formation of the dna hairpin structure d(GCGAAGC) with aromatic ligands

The energy contributions of various physical interactions to the total Gibbs energy of complex formation of the biologically important DNA hairpin d(GCGAAGC) with aromatic antitumor antibiotics daunomycin and novantrone and mutagens ethidium bromide and proflavin have been calculated. It has been shown that the relatively small value of the total energy of binding of the ligands to the hairpin is the sum of components great in absolute value and different in sign. The contributions of van der Waals interactions and both intra- and intermolecular hydrogen bonds and bonds with aqueous environment have been studied. According to the calculations, the hydrophobic and van der Waals components are energetically favorable in complex formation of the ligands with the hairpin d(GCGAAGC), whereas the electrostatic (with consideration of hydrogen bonds) and entropic components are unfavorable.     

阳离子-π相互作用在两种典型蛋白质折叠型中的偏好性

阳离子-π相互作用是一种在阳离子和芳香性体系之间形成的一种作用力。在蛋白质中,带正电荷的氨基酸(Lys、Arg)和芳香族氨基酸(Phe、Tyr、Trp)之间可以形成阳离子-π相互作用。对α/β类蛋白中两种典型折叠类型――单绕和双绕的研究表明:(1)单绕结构中阳离子-π相互作用的分布密度是双绕结构的2.3倍。(2)Arg-Phe组合偏好在双绕中出现,Arg-Tyr组合偏好在单绕中出现。(3)在单绕中除Lys-Phe组合外,其余5种组合的阳离子-π相互作用能量高于双绕的对应组合,其中以Arg-Trp组合的能量最高。(4)在单绕结构中,样本所含氨基酸残基数量和样本中阳离子-π的数量有明显的相关性,在双绕结构中没有发现类似的相关性。(5)在单绕和双绕结构当中,把阳离子-π相互作用能量分解为静电能和范德华力能揭示出静电能与范德华力能之比接近2∶1,静电作用在阳离子-π相互作用中起主要作用。     

Probing the electronic structures and properties of neutral and anionic ScSi n (0,−1) (n = 1–6) clusters using ccCA-TM and G4 theory

We apply molecular docking, molecular dynamics (MD) simulation, and binding free energy calculation to investigate and reveal the binding mechanism between five xanthine inhibitors and DPP-4. The electrostatic and van der Waals interactions of the five inhibitors with DPP-4 are analyzed and discussed. The computed binding free energies using MM-PBSA method are in qualitatively agreement with experimental inhibitory potency of five inhibitors. The hydrogen bonds of inhibitors with Ser630 and Asp663 can stabilize the inhibitors in binding sites. The van der Waals interactions, especially the key contacts with His740, Asn710, Trp629, and Tyr666 have larger contributions to the binding free energy and play important roles in distinguishing the variant bioactivity of five inhibitors.     

Protein docking using continuum electrostatics and geometric fit

The computer program DOT quickly finds low-energy docked structures for two proteins by performing a systematic search over six degrees of freedom. A novel feature of DOT is its energy function, which is the sum of both a Poisson-Boltzmann electrostatic energy and a van der Waals energy, each represented as a grid-based correlation function. DOT evaluates the energy of interaction for many orientations of the moving molecule and maintains separate lists scored by either the electrostatic energy, the van der Waals energy or the composite sum of both. The free energy is obtained by summing the Boltzmann factor over all rotations at each grid point. Three important findings are presented. First, for a wide variety of protein-protein interactions, the composite-energy function is shown to produce larger clusters of correct answers than found by scoring with either van der Waals energy (geometric fit) or electrostatic energy alone. Second, free-energy clusters are demonstrated to be indicators of binding sites. Third, the contributions of electrostatic and attractive van der Waals energies to the total energy term appropriately reflect the nature of the various types of protein-protein interactions studied.     

Cloning,functional expression,biochemical characterization,and structural analysis of a haloalkane dehalogenase from <Emphasis Type="Italic">Plesiocystis pacifica</Emphasis> SIR-1

A haloalkane dehalogenase (DppA) from Plesiocystis pacifica SIR-1 was identified by sequence comparison in the NCBI database, cloned, functionally expressed in Escherichia coli, purified, and biochemically characterized. The three-dimensional (3D) structure was determined by X-ray crystallography and has been refined at 1.95 Å resolution to an R-factor of 21.93%. The enzyme is composed of an α/β-hydrolase fold and a cap domain and the overall fold is similar to other known haloalkane dehalogenases. Active site residues were identified as Asp123, His278, and Asp249 and Trp124 and Trp163 as halide-stabilizing residues. DppA, like DhlA from Xanthobacter autotrophicus GJ10, is a member of the haloalkane dehalogenase subfamily HLD-I. As a consequence, these enzymes have in common the relative position of their catalytic residues within the structure and also show some similarities in the substrate specificity. The enzyme shows high preference for 1-bromobutane and does not accept chlorinated alkanes, halo acids, or halo alcohols. It is a monomeric protein with a molecular mass of 32.6 kDa and exhibits maximum activity between 33 and 37°C with a pH optimum between pH 8 and 9. The K_m and k_cat values for 1-bromobutane were 24.0 mM and 8.08 s^?1. Furthermore, from the 3D-structure of DppA, it was found that the enzyme possesses a large and open active site pocket. Docking experiments were performed to explain the experimentally determined substrate preferences.     

Critical contribution of aromatic rings to specific recognition of polyether rings. The case of ciguatoxin CTX3C-ABC and its specific antibody 1C49

To address how proteins recognize polyether toxin compounds, we focused on the interaction between the ABC ring compound of ciguatoxin 3C and its specific antibody, 1C49. Surface plasmon resonance analyses indicated that Escherichia coli-expressed variable domain fragments (Fv) of 1C49 had the high affinity constants and slow dissociation constants typical of antigen-antibody interactions. Linear van't Hoff analyses suggested that the interaction is enthalpy-driven. We resolved the crystal structure of 1C49 Fv bound to ABC ring compound of ciguatoxin 3C at a resolution of 1.7A. The binding pocket of the antibody had many aromatic rings and bound the antigen by shape complementarity typical of hapten-antibody interactions. Three hydrogen bonds and many van der Waals interactions were present. We mutated several residues of the antibody to Ala, and we used surface plasmon resonance to analyze the interactions between the mutated antibodies and the antigen. This analysis identified Tyr-91 and Trp-96 in the light chain as hot spots for the interaction, and other residues made incremental contributions by conferring enthalpic advantages and reducing the dissociation rate constant. Systematic mutation of Tyr-91 indicated that CH-pi and pi-pi interactions between the aromatic ring at this site and the antigen made substantial contributions to the association, and van der Waals interactions inhibited dissociation, suggesting that aromaticity and bulkiness are critical for the specific recognition of polyether compounds by proteins.     

Utilization of trihalogenated propanes by Agrobacterium radiobacter AD1 through heterologous expression of the haloalkane dehalogenase from Rhodococcus sp. strain M15-3.

Trihalogenated propanes are toxic and recalcitrant organic compounds. Attempts to obtain pure bacterial cultures able to use these compounds as sole carbon and energy sources were unsuccessful. Both the haloalkane dehalogenase from Xanthobacter autotrophicus GJ10 (DhlA) and that from Rhodococcus sp. strain m15-3 (DhaA) were found to dehalogenate trihalopropanes to 2,3-dihalogenated propanols, but the kinetic properties of the latter enzyme are much better. Broad-host-range dehalogenase expression plasmids, based on RSF1010 derivatives, were constructed with the haloalkane dehalogenase from Rhodococcus sp. strain m15-3 under the control of the heterologous promoters P(lac), P(dhlA), and P(trc). The resulting plasmids yielded functional expression in several gram-negative bacteria. A catabolic pathway for trihalopropanes was designed by introducing these broad-host-range dehalogenase expression plasmids into Agrobacterium radiobacter AD1, which has the ability to utilize dihalogenated propanols for growth. The recombinant strain AD1(pTB3), expressing the haloalkane dehalogenase gene under the control of the dhlA promoter, was able to utilize both 1,2,3-tribromopropane and 1,2-dibromo-3-chloropropane as sole carbon sources. Moreover, increased expression of the haloalkane dehalogenase resulted in elevated resistance to trihalopropanes.     

Importance of van der Waals contact between Glu 35 and Trp 109 to the catalytic action of human lysozyme.

The importance of van der Waals contact between Glu 35 and Trp 109 to the active-site structure and the catalytic properties of human lysozyme (HL) has been investigated by site-directed mutagenesis. The X-ray analysis of mutant HLs revealed that both the replacement of Glu 35 by Asp or Ala, and the replacement of Trp 109 by Phe or Ala resulted in a significant but localized change in the active-site cleft geometry. A prominent movement of the backbone structure was detected in the region of residues 110 to 120 and in the region of residues 100 to 115 for the mutations concerning Glu 35 and Trp 109, respectively. Accompanied by the displacement of the main-chain atoms with a maximal deviation of C alpha atom position ranging from 0.7 A to 1.0 A, the mutant HLs showed a remarkable change in the catalytic properties against Micrococcus luteus cell substrate as compared with native HL. Although the replacement of Glu 35 by Ala completely abolished the lytic activity, HL-Asp 35 mutant retained a weak but a certain lytic activity, showing the possible involvement of the side-chain carboxylate group of Asp 35 in the catalytic action. The kinetic consequence derived from the replacement of Trp 109 by Phe or Ala together with the result of the structural change suggested that the structural detail of the cleft lobe composed of the residues 100 to 115 centered at Ala 108 was responsible for the turnover in the reaction of HL against the bacterial cell wall substrate. The results revealed that the van der Waals contact between Glu 35 and Trp 109 was an essential determinant in the catalytic action of HL.     

Molecular dynamics simulation and steered molecular dynamics simulation on irisin dimers

Irisin is found closely associated with promoting the browning of beige fat cells in white adipose tissue. The crystal structure reveals that irisin forms a continuous inter-subunit β-sheet dimer. Here, molecular dynamics (MD) simulation and steered molecular dynamics (SMD) simulation were performed to investigate the dissociation process and the intricate interactions between the two irisin monomers. In the process of MD, the interactions between the monomers were roughly analyzed through the average numbers of both hydrophobic contacts and H-bonds. Then, SMD was performed to investigate the accurate interaction energy between the monomers. By the analysis of dissociation energy, the van der Waals (vdW) force was identified as the major energy to maintain the dimer structure, which also verified the results of MD simulation. Meanwhile, 11 essential residues were discovered by the magnitude of rupture force during dissociation. Among them, residues Arg75, Glu79, Ile77, Ala88, and Trp90 were reported in a previous study using the method of mutagenesis and size exclusion chromatography, and several new important residues (Arg72, Leu74, Phe76, Gln78, Val80, and Asp91) were also identified. Interestingly, the new important residues that we discovered and the important residues that were reported are located in the opposite side of the β-sheet of the dimer.     

Steered molecular dynamics simulation of the binding of the β2 and β3 regions in domain-swapped human cystatin C dimer

The crystal structure of the human cystatin C (hCC) dimer revealed that a stable twofold-symmetric dimer was formed via 3D domain swapping. Domain swapping with the need for near-complete unfolding has been proposed as a possible route for amyloid fibril initiation. Thus, the interesting interactions that occur between the two molecules may be important for the further aggregation of the protein. In this work, we performed steered molecular dynamics (SMD) simulations to investigate the dissociation of the β2 and β3 strands in the hCC dimer. The energy changes observed during the SMD simulations showed that electrostatic interactions were the dominant interactions involved in stabilizing the two parts of the dimer during the early stages of SMD simulation, whereas van der Waals (VDW) interactions and electrostatic interactions were equally matched during the latter stages. Furthermore, our data indicated that the two parts of the dimer are stabilized by intermolecular hydrogen bonds among the residues Arg51 (β2), Gln48 (β2), Asp65 (β3), and Glu67 (β3), salt bridges among the residues Arg53 (β2), Arg51 (β2), and Asp65 (β3), and VDW interactions among the residues Gln48 (β2), Arg51 (β2), Glu67 (β3), Asp65 (β3), Phe63 (β3), and Asn61 (β3). The residues Gln48 (β2), Arg51 (β2), Asp65 (β3) and Glu67 (β3) appear to be crucial, as they play important roles in both electrostatic and VDW interactions. Thus, the present study determined the key residues involved in the stabilization of the domain-swapped dimer structure, and also provided molecular-level insights into the dissociation process of the hCC dimer.     

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.     

Application of molecular dynamics and free energy perturbation methods to metalloporphyrin-ligand systems II: CO and dioxygen binding to myoglobin.

The protein contribution to the relative binding affinity of the ligands CO and O2 toward myoglobin (Mb) has been simulated using free energy perturbation calculations. The tautomers of the His E7 residue are different for the oxymyoglobin (MbO2) and carboxymyoglobin (MbCO) systems. This was modeled by performing two-step calculations that mutate the ligand and mutate the His E7 tautomers in separate steps. Differences in hydrogen bonding to the O2 and CO ligands were incorporated into the model. The O2 complex was calculated to be 2-3 kcal/mol more stable than the corresponding CO complex when compared to the same difference in an isolated heme control. This value agrees well with the experimental value of 2.0 kcal/mol. In qualitative agreement with experiments, the Fe-C-O bond is found to be bent (theta = 159.8 degrees) with a small tilt (theta = 6.2 degrees). The contributions made by each of the 29 residues--within the 9.0-A radius of the iron atom--to the free energy difference are separated into van der Waals and electrostatic contributions; the latter contributions are dominant. Aside from the proximal histidine and the heme group, the residues having the largest difference in free energy in mutating MbO2-->MbCO are His E7, Phe CD1, Phe CD4, Val E11, and Thr E10.     

Utilization of Trihalogenated Propanes by Agrobacterium radiobacter AD1 through Heterologous Expression of the Haloalkane Dehalogenase from Rhodococcus sp. Strain m15-3

Trihalogenated propanes are toxic and recalcitrant organic compounds. Attempts to obtain pure bacterial cultures able to use these compounds as sole carbon and energy sources were unsuccessful. Both the haloalkane dehalogenase from Xanthobacter autotrophicus GJ10 (DhlA) and that from Rhodococcus sp. strain m15-3 (DhaA) were found to dehalogenate trihalopropanes to 2,3-dihalogenated propanols, but the kinetic properties of the latter enzyme are much better. Broad-host-range dehalogenase expression plasmids, based on RSF1010 derivatives, were constructed with the haloalkane dehalogenase from Rhodococcus sp. strain m15-3 under the control of the heterologous promoters P_lac, P_dhlA, and P_trc. The resulting plasmids yielded functional expression in several gram-negative bacteria. A catabolic pathway for trihalopropanes was designed by introducing these broad-host-range dehalogenase expression plasmids into Agrobacterium radiobacter AD1, which has the ability to utilize dihalogenated propanols for growth. The recombinant strain AD1(pTB3), expressing the haloalkane dehalogenase gene under the control of the dhlA promoter, was able to utilize both 1,2,3-tribromopropane and 1,2-dibromo-3-chloropropane as sole carbon sources. Moreover, increased expression of the haloalkane dehalogenase resulted in elevated resistance to trihalopropanes.