Generative artificial intelligence (AI) models trained on natural protein sequences have been used to design functional enzymes. However, their ability to predict individual reaction steps in enzyme catalysis remains unclear, limiting the potential use of sequence information for enzyme engineering. In this study, we demonstrated that sequence information can predict the rate of the SN2 step of a haloalkane dehalogenase using a generative maximum-entropy (MaxEnt) model. We then designed lower-order protein variants of haloalkane dehalogenase using the model. Kinetic measurements confirmed the successful design of protein variants that enhance catalytic activity, above that of the wild type, in the overall reaction and in particular in the SN2 step. On the simulation side, we provided molecular insights into these designs for the SN2 step using the empirical valence bond (EVB) and metadynamics simulations. The EVB calculations showed activation barriers consistent with experimental reaction rates, while examining the effect of amino acid replacements on the electrostatic effect on the activation barrier and the consequence of water penetration, as well as the extent of ground state destabilization/stabilization. Metadynamics simulations emphasize the importance of the substrate positioning in enzyme catalysis. Overall, our AI-guided approach successfully enabled the design of a variant with a faster rate for the SN2 step than the wild-type enzyme, despite haloalkane dehalogenase being extensively optimized through natural evolution.

• MeSH
• Hydrolases * chemistry   metabolism   genetics   MeSH
• Kinetics MeSH
• Protein Engineering MeSH
• Molecular Dynamics Simulation MeSH
• Artificial Intelligence * MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• haloalkane dehalogenase MeSH Browser
• Hydrolases * MeSH

Haloalkane dehalogenase (HLD) enzymes employ an SN 2 nucleophilic substitution mechanism to erase halogen substituents in diverse organohalogen compounds. Subfamily I and II HLDs are well-characterized enzymes, but the mode and purpose of multimerization of subfamily III HLDs are unknown. Here we probe the structural organization of DhmeA, a subfamily III HLD-like enzyme from the archaeon Haloferax mediterranei, by combining cryo-electron microscopy (cryo-EM) and x-ray crystallography. We show that full-length wild-type DhmeA forms diverse quaternary structures, ranging from small oligomers to large supramolecular ring-like assemblies of various sizes and symmetries. We optimized sample preparation steps, enabling three-dimensional reconstructions of an oligomeric species by single-particle cryo-EM. Moreover, we engineered a crystallizable mutant (DhmeAΔGG ) that provided diffraction-quality crystals. The 3.3 Å crystal structure reveals that DhmeAΔGG forms a ring-like 20-mer structure with outer and inner diameter of ~200 and ~80 Å, respectively. An enzyme homodimer represents a basic repeating building unit of the crystallographic ring. Three assembly interfaces (dimerization, tetramerization, and multimerization) were identified to form the supramolecular ring that displays a negatively charged exterior, while its interior part harboring catalytic sites is positively charged. Localization and exposure of catalytic machineries suggest a possible processing of large negatively charged macromolecular substrates.

• Keywords
• DhmeA, Haloferax mediterranei, catalysis, cryo-EM, haloalkane dehalogenase, multimerization, x-ray crystallography,
• MeSH
• Cryoelectron Microscopy methods   MeSH
• Hydrolases * chemistry   MeSH
• Crystallography, X-Ray MeSH
• Substrate Specificity MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• haloalkane dehalogenase MeSH Browser
• Hydrolases * MeSH

Haloalkane dehalogenases are enzymes that catalyze the cleavage of carbon-halogen bonds in halogenated compounds. They serve as model enzymes for studying structure-function relationships of >100.000 members of the α/β-hydrolase superfamily. Detailed kinetic analysis of their reaction is crucial for understanding the reaction mechanism and developing novel concepts in protein engineering. Fluorescent substrates, which change their fluorescence properties during a catalytic cycle, may serve as attractive molecular probes for studying the mechanism of enzyme catalysis. In this work, we present the development of the first fluorescent substrates for this enzyme family based on coumarin and BODIPY chromophores. Steady-state and pre-steady-state kinetics with two of the most active haloalkane dehalogenases, DmmA and LinB, revealed that both fluorescent substrates provided specificity constant two orders of magnitude higher (0.14-12.6 μM-1 s-1) than previously reported representative substrates for the haloalkane dehalogenase family (0.00005-0.014 μM-1 s-1). Stopped-flow fluorescence/FRET analysis enabled for the first time monitoring of all individual reaction steps within a single experiment: (i) substrate binding, (ii-iii) two subsequent chemical steps and (iv) product release. The newly introduced fluorescent molecules are potent probes for fast steady-state kinetic profiling. In combination with rapid mixing techniques, they provide highly valuable information about individual kinetic steps and mechanism of haloalkane dehalogenases. Additionally, these molecules offer high specificity and efficiency for protein labeling and can serve as probes for studying protein hydration and dynamics as well as potential markers for cell imaging.

• Keywords
• Enzyme kinetics, Fluorescent substrate, Haloalkane dehalogenase, Mechanism, Protein labeling,
• Publication type
• Journal Article MeSH

Rational enzyme design presents a major challenge that has not been overcome by computational approaches. One of the key challenges is the difficulty in assessing the magnitude of the maximum possible catalytic activity. In an attempt to overcome this challenge, we introduce a strategy that takes an active enzyme (assuming that its activity is close to the maximum possible activity), design mutations that reduce the catalytic activity, and then try to restore that catalysis by mutating other residues. Here we take as a test case the enzyme haloalkane dehalogenase (DhlA), with a 1,2-dichloroethane substrate. We start by demonstrating our ability to reproduce the results of single mutations. Next, we design mutations that reduce the enzyme activity and finally design double mutations that are aimed at restoring the activity. Using the computational predictions as a guide, we conduct an experimental study that confirms our prediction in one case and leads to inconclusive results in another case with 1,2-dichloroethane as substrate. Interestingly, one of our predicted double mutants catalyzes dehalogenation of 1,2-dibromoethane more efficiently than the wild-type enzyme.

• Keywords
• EVB, dehalogenase, enzyme design, nucleophilic substitution, transient kinetics,
• MeSH
• Models, Chemical * MeSH
• Ethylene Dichlorides chemistry   MeSH
• Hydrolases chemistry   MeSH
• Catalytic Domain MeSH
• Models, Molecular * MeSH
• Computer Simulation * MeSH
• Substrate Specificity MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Research Support, N.I.H., Extramural MeSH
• Names of Substances
• Ethylene Dichlorides MeSH
• ethylene dichloride MeSH Browser
• haloalkane dehalogenase MeSH Browser
• Hydrolases MeSH

Haloalkane dehalogenases catalyze the hydrolysis of halogen-carbon bonds in organic halogenated compounds and as such are of great utility as biocatalysts. The crystal structures of the haloalkane dehalogenase DhlA from the bacterium from Xanthobacter autotrophicus GJ10, specifically adapted for the conversion of the small 1,2-dichloroethane (DCE) molecule, display the smallest catalytic site (110 Å3) within this enzyme family. However, during a substrate-specificity screening, we noted that DhlA can catalyze the conversion of far bulkier substrates, such as the 4-(bromomethyl)-6,7-dimethoxy-coumarin (220 Å3). This large substrate cannot bind to DhlA without conformational alterations. These conformational changes have been previously inferred from kinetic analysis, but their structural basis has not been understood. Using molecular dynamic simulations, we demonstrate here the intrinsic flexibility of part of the cap domain that allows DhlA to accommodate bulky substrates. The simulations displayed two routes for transport of substrates to the active site, one of which requires the conformational change and is likely the route for bulky substrates. These results provide insights into the structure-dynamics function relationships in enzymes with deeply buried active sites. Moreover, understanding the structural basis for the molecular adaptation of DhlA to 1,2-dichloroethane introduced into the biosphere during the industrial revolution provides a valuable lesson in enzyme design by nature.

• Keywords
• active site, conformational change, dichloroethane degradation, enzyme catalysis, enzyme kinetics, enzyme mechanism, ethylene dichloride, haloalkane dehalogenase, molecular dynamics, molecular evolution, organic halogen, organohalogen, protein conformation,
• MeSH
• Ethylene Dichlorides metabolism   MeSH
• Halogenation MeSH
• Hydrolases chemistry   metabolism   MeSH
• Catalytic Domain MeSH
• Kinetics MeSH
• Protein Conformation MeSH
• Crystallography, X-Ray MeSH
• Coumarins chemistry   metabolism   MeSH
• Methylation MeSH
• Molecular Dynamics Simulation MeSH
• Molecular Docking Simulation MeSH
• Substrate Specificity MeSH
• Xanthobacter chemistry   enzymology   metabolism   MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• Ethylene Dichlorides MeSH
• ethylene dichloride MeSH Browser
• haloalkane dehalogenase MeSH Browser
• Hydrolases MeSH
• Coumarins MeSH

Two haloalkane dehalogenases, LinBUT and LinBMI, each with 296 amino acid residues, exhibit only seven amino acid residue differences between them, but LinBMI's catalytic performance towards β-hexachlorocyclohexane (β-HCH) is considerably higher than LinBUT's. To elucidate the molecular basis governing this difference, intermediate mutants between LinBUT and LinBMI were constructed and kinetically characterized. The activities of LinBUT-based mutants gradually increased by cumulative mutations into LinBUT, and the effects of the individual amino acid substitutions depended on combination with other mutations. These results indicated that LinBUT's β-HCH degradation activity can be enhanced in a stepwise manner by the accumulation of point mutations.

• Keywords
• Biodegradation, Haloalkane dehalogenase, Protein evolution, Xenobiotics, β-Hexachlorocyclohexane,
• Publication type
• Journal Article MeSH

Haloalkane dehalogenases are enzymes that catalyze the cleavage of the carbon-halogen bond by a hydrolytic mechanism. Genomes of Mycobacterium tuberculosis and M. bovis contain at least two open reading frames coding for the polypeptides showing a high sequence similarity with biochemically characterized haloalkane dehalogenases. We describe here the cloning of the haloalkane dehalogenase genes dmbA and dmbB from M. bovis 5033/66 and demonstrate the dehalogenase activity of their translation products. Both of these genes are widely distributed among species of the M. tuberculosis complex, including M. bovis, M. bovis BCG, M. africanum, M. caprae, M. microti, and M. pinnipedii, as shown by the PCR screening of 48 isolates from various hosts. DmbA and DmbB proteins were heterologously expressed in Escherichia coli and purified to homogeneity. The DmbB protein had to be expressed in a fusion with thioredoxin to obtain a soluble protein sample. The temperature optimum of DmbA and DmbB proteins determined with 1,2-dibromoethane is 45 degrees C. The melting temperature assessed by circular dichroism spectroscopy of DmbA is 47 degrees C and DmbB is 57 degrees C. The pH optimum of DmbA depends on composition of a buffer with maximal activity at 9.0. DmbB had a single pH optimum at pH 6.5. Mycobacteria are currently the only genus known to carry more than one haloalkane dehalogenase gene, although putative haloalkane dehalogenases can be inferred in more then 20 different bacterial species by comparative genomics. The evolution and distribution of haloalkane dehalogenases among mycobacteria is discussed.

• MeSH
• Bacterial Proteins chemistry   genetics   isolation & purification   metabolism   MeSH
• Hydrolases chemistry   genetics   isolation & purification   metabolism   MeSH
• Cloning, Molecular * MeSH
• Hydrogen-Ion Concentration MeSH
• Humans MeSH
• Molecular Sequence Data MeSH
• Mycobacterium bovis enzymology   genetics   MeSH
• Mycobacterium classification   enzymology   genetics   MeSH
• Sequence Analysis, DNA MeSH
• Cattle MeSH
• Enzyme Stability MeSH
• Temperature MeSH
• Animals MeSH
• Check Tag
• Humans MeSH
• Cattle MeSH
• Animals MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• Bacterial Proteins MeSH
• haloalkane dehalogenase MeSH Browser
• Hydrolases MeSH

We evaluate the applicability of automated molecular docking techniques and quantum mechanical calculations to the construction of a set of structures of enzyme-substrate complexes for use in Comparative binding energy (COMBINE) analysis to obtain 3D structure-activity relationships. The data set studied consists of the complexes of eighteen substrates docked within the active site of haloalkane dehalogenase (DhlA) from Xanthobacter autotrophicus GJ10. The results of the COMBINE analysis are compared with previously reported data obtained for the same dataset from modelled complexes that were based on an experimentally determined structure of the DhlA-dichloroethane complex. The quality of fit and the internal predictive power of the two COMBINE models are comparable, but better external predictions are obtained with the new approach. Both models show a similar composition of the principal components. Small differences in the relative contributions that are assigned to important residues for explaining binding affinity differences can be directly linked to structural differences in the modelled enzyme-substrate complexes: (i) rotation of all substrates in the active site about their longitudinal axis, (ii) repositioning of the ring of epihalohydrines and the halogen substituents of 1,2-dihalopropanes, and (iii) altered conformation of the long-chain molecules (halobutanes and halohexanes). For external validation, both a novel substrate not included in the training series and two different mutant proteins were used. The results obtained can be useful in the future to guide the rational engineering of substrate specificity in DhlA and other related enzymes.

• MeSH
• Principal Component Analysis MeSH
• Models, Chemical MeSH
• Databases, Protein MeSH
• Hydrocarbons, Halogenated chemistry   metabolism   MeSH
• Hydrolases chemistry   metabolism   MeSH
• Protein Conformation MeSH
• Quantitative Structure-Activity Relationship MeSH
• Least-Squares Analysis MeSH
• Molecular Conformation MeSH
• Models, Molecular MeSH
• Molecular Structure MeSH
• Computer Simulation MeSH
• Static Electricity MeSH
• Substrate Specificity MeSH
• Thermodynamics MeSH
• Protein Binding MeSH
• Binding Sites MeSH
• Xanthobacter enzymology   MeSH
• Imaging, Three-Dimensional MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• haloalkane dehalogenase MeSH Browser
• Hydrocarbons, Halogenated MeSH
• Hydrolases MeSH

Haloalkane dehalogenases are microbial enzymes that catalyze cleavage of the carbon-halogen bond by a hydrolytic mechanism. Until recently, these enzymes have been isolated only from bacteria living in contaminated environments. In this report we describe cloning of the dehalogenase gene dhmA from Mycobacterium avium subsp. avium N85 isolated from swine mesenteric lymph nodes. The dhmA gene has a G+C content of 68.21% and codes for a polypeptide that is 301 amino acids long and has a calculated molecular mass of 34.7 kDa. The molecular masses of DhmA determined by sodium dodecyl sulfate-polyacrylamide gel electrophoresis and by gel permeation chromatography are 34.0 and 35.4 kDa, respectively. Many residues essential for the dehalogenation reaction are conserved in DhmA; the putative catalytic triad consists of Asp123, His279, and Asp250, and the putative oxyanion hole consists of Glu55 and Trp124. Trp124 should be involved in substrate binding and product (halide) stabilization, while the second halide-stabilizing residue cannot be identified from a comparison of the DhmA sequence with the sequences of three dehalogenases with known tertiary structures. The haloalkane dehalogenase DhmA shows broad substrate specificity and good activity with the priority pollutant 1,2-dichloroethane. DhmA is significantly less stable than other currently known haloalkane dehalogenases. This study confirms that a hydrolytic dehalogenase is present in the facultative pathogen M. avium. The presence of dehalogenase-like genes in the genomes of other mycobacteria, including the obligate pathogens Mycobacterium tuberculosis and Mycobacterium bovis, as well as in other bacterial species, including Mesorhizobium loti, Xylella fastidiosa, Photobacterium profundum, and Caulobacter crescentus, led us to speculate that haloalkane dehalogenases have some other function besides catalysis of hydrolytic dehalogenation of halogenated substances.

• MeSH
• Alkanes metabolism   MeSH
• Phylogeny MeSH
• Transcription, Genetic MeSH
• Hydrolases chemistry   genetics   isolation & purification   metabolism   MeSH
• Cloning, Molecular methods   MeSH
• Molecular Sequence Data MeSH
• Mycobacterium avium enzymology   genetics   growth & development   MeSH
• Swine MeSH
• Amino Acid Sequence MeSH
• Sequence Analysis, DNA MeSH
• Sequence Alignment MeSH
• Animals MeSH
• Check Tag
• Animals MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• Alkanes MeSH
• haloalkane dehalogenase MeSH Browser
• Hydrolases MeSH

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.

• MeSH
• Hydrolases chemistry   metabolism   MeSH
• Solvents MeSH
• Protein Structure, Secondary MeSH
• Substrate Specificity MeSH
• Temperature MeSH
• Water chemistry   MeSH
• Xanthobacter enzymology   MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• haloalkane dehalogenase MeSH Browser
• Hydrolases MeSH
• Solvents MeSH
• Water MeSH