The catalytic mechanism of fluoroacetate dehalogenase: a computational exploration of biological dehalogenation

Defluorination of Organofluorine Compounds Using Dehalogenase Enzymes from Delftia acidovorans (D4B)

Organofluorine compounds have been widely used as pharmaceuticals, agricultural pesticides, and water-resistant coatings for decades; however, these compounds are recognized as environmental pollutants. The capability of microorganisms and enzymes to defluorinate organofluorine compounds is both rare and highly desirable to facilitate environmental remediation efforts. Recently, a strain of Delftia acidovorans (D4B) was identified with potential biodegradation activity toward perfluoroalkyl substances (PFAS) and other organofluorine compounds. Genomic analysis found haloacid and fluoroacetate dehalogenases as enzymes associated with Delftia acidovorans. Here, defluorination activity of these enzymes toward different fluorinated substrates was investigated after their recombinant expression and purification from E. coli. Using an electrochemical fluoride probe, 19F NMR, and mass spectrometry to monitor defluorination, we identified two dehalogenases, DeHa2 (a haloacid dehalogenase) and DeHa4 (a fluoroacetate dehalogenase), with activity toward mono- and difluoroacetate. Of the two dehalogenases, DeHa4 demonstrated a low pH optimum compared to DeHa2, which lost catalytic activity under acidic conditions. DeHa2 and DeHa4 are relatively small proteins, operate under aerobic conditions, and remain active for days in the presence of substrates. Significantly, while there have been many reports on dehalogenation of monofluoroacetate by dehalogenases, this study adds to the relatively small list of enzymes reported to carry out enzymatic defluorination of the more recalcitrant disubstituted carbon in an organofluorine compound. Thus, DeHa2 and DeHa4 represent organofluorine dehalogenases that may be used in the future to design and engineer robust defluorination agents for environmental remediation efforts.

Toward the development of a molecular toolkit for the microbial remediation of per-and polyfluoroalkyl substances

Per- and polyfluoroalkyl substances (PFAS) are highly fluorinated synthetic organic compounds that have been used extensively in various industries owing to their unique properties. The PFAS family encompasses diverse classes, with only a fraction being commercially relevant. These substances are found in the environment, including in water sources, soil, and wildlife, leading to human exposure and fueling concerns about potential human health impacts. Although PFAS degradation is challenging, biodegradation offers a promising, eco-friendly solution. Biodegradation has been effective for a variety of organic contaminants but is yet to be successful for PFAS due to a paucity of identified microbial species capable of transforming these compounds. Recent studies have investigated PFAS biotransformation and fluoride release; however, the number of specific microorganisms and enzymes with demonstrable activity with PFAS remains limited. This review discusses enzymes that could be used in PFAS metabolism, including haloacid dehalogenases, reductive dehalogenases, cytochromes P450, alkane and butane monooxygenases, peroxidases, laccases, desulfonases, and the mechanisms of microbial resistance to intracellular fluoride. Finally, we emphasize the potential of enzyme and microbial engineering to advance PFAS degradation strategies and provide insights for future research in this field.

Mechanisms and Opportunities for Rational In Silico Design of Enzymes to Degrade Per- and Polyfluoroalkyl Substances (PFAS)

Per and polyfluoroalkyl substances (PFAS) present a unique challenge to remediation techniques because their strong carbon-fluorine bonds make them difficult to degrade. This review explores the use of in silico enzymatic design as a potential PFAS degradation technique. The scope of the enzymes included is based on currently known PFAS degradation techniques, including chemical redox systems that have been studied for perfluorooctanesulfonic acid (PFOS) and perfluorooctanoic acid (PFOA) defluorination, such as those that incorporate hydrated electrons, sulfate, peroxide, and metal catalysts. Bioremediation techniques are also discussed, namely the laccase and horseradish peroxidase systems. The redox potential of known reactants and enzymatic radicals/metal-complexes are then considered and compared to potential enzymes for degrading PFAS. The molecular structure and reaction cycle of prospective enzymes are explored. Current knowledge and techniques of enzyme design, particularly radical-generating enzymes, and application are also discussed. Finally, potential routes for bioengineering enzymes to enable or enhance PFAS remediation are considered as well as the future outlook for computational exploration of enzymatic in situ bioremediation routes for these highly persistent and globally distributed contaminants.

Structural insights into hydrolytic defluorination of difluoroacetate by microbial fluoroacetate dehalogenases

Fluorine forms the strongest single bond to carbon with the highest bond dissociation energy among natural products. However, fluoroacetate dehalogenases (FADs) have been shown to hydrolyze this bond in fluoroacetate under mild reaction conditions. Furthermore, two recent studies demonstrated that the FAD RPA1163 from Rhodopseudomonas palustris can also accept bulkier substrates. In this study, we explored the substrate promiscuity of microbial FADs and their ability to defluorinate polyfluorinated organic acids. Enzymatic screening of eight purified dehalogenases with reported fluoroacetate defluorination activity revealed significant hydrolytic activity against difluoroacetate in three proteins. Product analysis using liquid chromatography-mass spectrometry identified glyoxylic acid as the final product of enzymatic DFA defluorination. The crystal structures of DAR3835 from Dechloromonas aromatica and NOS0089 from Nostoc sp. were determined in the apo-state along with the DAR3835 H274N glycolyl intermediate. Structure-based site-directed mutagenesis of DAR3835 demonstrated a key role for the catalytic triad and other active site residues in the defluorination of both fluoroacetate and difluoroacetate. Computational analysis of the dimer structures of DAR3835, NOS0089, and RPA1163 indicated the presence of one substrate access tunnel in each protomer. Moreover, protein-ligand docking simulations suggested similar catalytic mechanisms for the defluorination of both fluoroacetate and difluoroacetate, with difluoroacetate being defluorinated via two consecutive defluorination reactions producing glyoxylate as the final product. Thus, our findings provide molecular insights into substrate promiscuity and catalytic mechanism of FADs, which are promising biocatalysts for applications in synthetic chemistry and bioremediation of fluorochemicals.

Strategies for the Biodegradation of Polyfluorinated Compounds

Many cite the strength of C–F bonds for the poor microbial biodegradability of polyfluorinated organic compounds (PFCs). However, commercial PFCs almost invariably contain more functionality than fluorine. The additional functionality provides a weak entry point for reactions that activate C–F bonds and lead to their eventual cleavage. This metabolic activation strategy is common in microbial biodegradation pathways and is observed with aromatic hydrocarbons, chlorinated compounds, phosphonates and many other compounds. Initial metabolic activation precedes critical bond breakage and assimilation of nutrients. A similar strategy with commercial PFCs proceeds via initial attack at the non-fluorinated functionalities: sulfonates, carboxylates, chlorines, phenyl rings, or phosphonates. Metabolic transformation of these non-fluorinated groups can activate the C–F bonds, allowing more facile cleavage than a direct attack on the C–F bonds. Given that virtually all compounds denoted as “PFAS” are not perfluorinated and are not alkanes, it is posited here that considering their individual chemical classes is more useful for both chemical and microbiological considerations of their fate.

Defluorination Capability of l-2-Haloacid Dehalogenases in the HAD-Like Hydrolase Superfamily Correlates with Active Site Compactness

l-2-Haloacid dehalogenases, industrially and environmentally important enzymes that catalyse cleavage of the carbon-halogen bond in S-2-halocarboxylic acids, were known to hydrolyse chlorinated, brominated and iodinated substrates but no activity towards fluorinated compounds had been reported. A screen for novel dehalogenase activities revealed four l-2-haloacid dehalogenases capable of defluorination. We now report crystal structures for two of these enzymes, Bpro0530 and Rha0230, as well as for the related proteins PA0810 and RSc1362, which hydrolyse chloroacetate but not fluoroacetate, all at ∼2.2 Å resolution. Overall structure and active sites of these enzymes are highly similar. In molecular dynamics (MD) calculations, only the defluorinating enzymes sample more compact conformations, which in turn allow more effective interactions with the small fluorine atom. Structural constraints, based on X-ray structures and MD calculations, correctly predict the defluorination activity of the homologous enzyme ST2570.

Comprehensive Understanding of Fluoroacetate Dehalogenase-Catalyzed Degradation of Fluorocarboxylic Acids: A QM/MM Approach

Fluorochemicals are persistent, bioaccumulative, and toxic compounds that are widely tributed in the environment. Developing efficient biodegradation strategies to decompose the fluorochemicals via breaking the inert C-F bonds presents a holistic challenge. As a promising biodegradation enzyme candidate, fluoroacetate dehalogenase (FAcD) has been reported as the only non-metallic enzyme to catalyze the cleavage of the strong C-F bond. Here, we systematically investigated the catalytic actions of FAcD toward its natural substrate fluoroacetate using molecular dynamics simulations and quantum mechanism/molecular mechanism calculations. We propose that the enzymatic transformation involves four elementary steps, (I) C-F bond activation, (II) nucleophilic attack, (III) C-O bond cleavage, and (IV) proton transfer. Our results show that nucleophilic attack is the rate-determining step. However, for difluoroacetate and trifluoroacetate, C-F bond activation, instead of nucleophilic attack, becomes the rate-determining step. We show that FAcD, originally recognized as α-fluorocarboxylic acid degradation enzyme, can catalyze the defluorination of difluoroacetate to glyoxylate, which is captured by our high-resolution mass spectrometry experiments. In addition, we employed amino acid electrostatic analysis method to screen potential mutation hotspots for tuning FAcD's electrostatic environment to favor substrate conversion. The comprehensive understanding of catalytic mechanism will inform a rational enzyme engineering strategy to degrade fluorochemicals for benefits of environmental sustainability.

Fluoroacetate dehalogenase catalyzed dehalogenation of halogenated carboxylic acids: A QM/MM approach

Dehalogenation is one of the most important reactions in environmental pollution control, for instance, the degradation of persistent organic pollutants (POPs). Recently, fluoroacetate dehalogenase (FAcD) has been reported to catalyze the dehalogenation reactions, which shows great potential in treating halogenated pollutants. Here the dehalogenation mechanism catalyzed by FAcD was fully deciphered with the aid of quantum mechanics/molecular mechanics method. The results show that FAcD catalyzed dehalogenation efficiency follows the order of defluorination > dechlorination > debromination. The corresponding Boltzmann-weighted average barriers are 10.1, 19.7, and 20.9 kcal mol^-1. Positive/negative correlations between activation barriers and structural parameters (e.g. distance and angle) for FAcD catalyzed dechlorination and debromination were established. Based on the structure-energy relationship, we propose that mutation of the binding pocket amino acids (e.g. His155, Trp156, Tyr219) to smaller proton donor amino acids (e.g. Serine, Threonine, Cysteine, Asparagine) may increase the efficiency for dechlorination and debromination. The results may of practical value for the efficient degradation of chlorined and bromined pollutants by harnessing FAcD.

Metabolic and Pharmaceutical Aspects of Fluorinated Compounds

The applications of fluorine in drug design continue to expand, facilitated by an improved understanding of its effects on physicochemical properties and the development of synthetic methodologies that are providing access to new fluorinated motifs. In turn, studies of fluorinated molecules are providing deeper insights into the effects of fluorine on metabolic pathways, distribution, and disposition. Despite the high strength of the C-F bond, the departure of fluoride from metabolic intermediates can be facile. This reactivity has been leveraged in the design of mechanism-based enzyme inhibitors and has influenced the metabolic fate of fluorinated compounds. In this Perspective, we summarize the literature associated with the metabolism of fluorinated molecules, focusing on examples where the presence of fluorine influences the metabolic profile. These studies have revealed potentially problematic outcomes with some fluorinated motifs and are enhancing our understanding of how fluorine should be deployed.

Harnessing fluoroacetate dehalogenase for defluorination of fluorocarboxylic acids: in silico and in vitro approach

Widely distributed fluorocarboxylic acids have aroused worldwide environmental concerns due to its toxicity, persistence, and bioaccumulation. Enzyme-based eco-friendly biodegradation techniques have become increasingly important in treating fluorocarboxylic acids. Here we utilized in silico and in vitro approaches to investigate the defluorination mechanism of fluoroacetate dehalogenase (FAcD) toward monofluoropropionic acids at atomic-level. The experimentally determined k_cat and k_M for defluorination of 2-fluoropropionic acid are 330 ± 60 min^-1 and 6.12 ± 0.13 mM. The in silico results demonstrated positive/negative correlations between activation barriers and structural parameters (e.g. distance and angle) under different enzymatic conformations. We also screened computationally and tested in vitro (enzyme assay and kinetic study) the catalytic proficiency of FAcD toward polyfluoropropionic acids and perfluoropropionic acids which are known to be challenging for enzymatic degradation. The results revealed potential degradation activity of FAcD enzyme toward 2,3,3,3-tetrafluoropropionic acids. Our work will initiate the development of a new "integrated approach" for enzyme engineering to degrade environmentally persistent fluorocarboxylic acids.

Time-resolved crystallography reveals allosteric communication aligned with molecular breathing

A comprehensive understanding of protein function demands correlating structure and dynamic changes. Using time-resolved serial synchrotron crystallography, we visualized half-of-the-sites reactivity and correlated molecular-breathing motions in the enzyme fluoroacetate dehalogenase. Eighteen time points from 30 milliseconds to 30 seconds cover four turnover cycles of the irreversible reaction. They reveal sequential substrate binding, covalent-intermediate formation, setup of a hydrolytic water molecule, and product release. Small structural changes of the protein mold and variations in the number and placement of water molecules accompany the various chemical steps of catalysis. Triggered by enzyme-ligand interactions, these repetitive changes in the protein framework’s dynamics and entropy constitute crucial components of the catalytic machinery.

In silico design of potentially functional artificial metallo-haloalkane dehalogenase containing catalytic zinc

Artificial metalloenzymes are unique as they combine the good features of homogeneous and enzymatic catalysts, and they can potentially improve some difficult catalytic assays. This study reports a method that can be used to create an artificial metal-binding site prior to proving it to be functional in a wet lab. Haloalkane dehalogenase was grafted into a metal-binding site to form an artificial metallo-haloalkane dehalogenase and was studied for its potential functionalities in silico. Computational protocols regarding dynamic metal docking were studied using native metalloenzymes and functional artificial metalloenzymes. Using YASARA Structure, a simulation box covering template structure was created to be filled with water molecules followed by one mutated water molecule closest to the metal-binding site to metal ion. A simple energy minimization step was subsequently run using an AMBER force field to allow the metal ion to interact with the metal-binding residues. Long molecular dynamic simulation using YASARA Structure was performed to analyze the stability of the metal-binding site and the distance between metal-binding residues. Metal ions fluctuating around 2.0 Å across a 20 ns simulation indicated a stable metal-binding site. Metal-binding energies were predicted using FoldX, with a native metalloenzyme (carbonic anhydrase) scoring 18.0 kcal/mol and the best mutant model (C1a) scoring 16.4 kcal/mol. Analysis of the metal-binding site geometry was performed using CheckMyMetal, and all scores for the metalloenzymes and mutant models were in an acceptable range. Like native metalloenzymes, the metal-binding site of C1a was supported by residues in the second coordination shell to maintain a more coordinated metal-binding site. Short-chain multihalogenated alkanes (1,2-dibromoethane and 1,2,3-trichloropropane) were able to dock in the active site of C1a. The halides of the substrate were in contact with both the metal and halide-stabilizing residues, thus indicating a better stabilization of the substrate. The simple catalytic mechanism proposed is that the metal ion interacted with halogen and polarized the carbon-halogen bond, thus making the alpha carbon susceptible to attack by nucleophilic hydroxide. The interaction between halogen in the metal ion and halide-stabilizing residues may help to improve the stabilization of the substrate-enzyme complex and reduce the activation energy. This study reports a modified dynamic metal-docking protocol and validation tests to verify the metal-binding site. These approaches can be applied to design different kinds of artificial metalloenzymes or metal-binding sites.

Dehalogenases: From Improved Performance to Potential Microbial Dehalogenation Applications

The variety of halogenated substances and their derivatives widely used as pesticides, herbicides and other industrial products is of great concern due to the hazardous nature of these compounds owing to their toxicity, and persistent environmental pollution. Therefore, from the viewpoint of environmental technology, the need for environmentally relevant enzymes involved in biodegradation of these pollutants has received a great boost. One result of this great deal of attention has been the identification of environmentally relevant bacteria that produce hydrolytic dehalogenases—key enzymes which are considered cost-effective and eco-friendly in the removal and detoxification of these pollutants. These group of enzymes catalyzing the cleavage of the carbon-halogen bond of organohalogen compounds have potential applications in the chemical industry and bioremediation. The dehalogenases make use of fundamentally different strategies with a common mechanism to cleave carbon-halogen bonds whereby, an active-site carboxylate group attacks the substrate C atom bound to the halogen atom to form an ester intermediate and a halide ion with subsequent hydrolysis of the intermediate. Structurally, these dehalogenases have been characterized and shown to use substitution mechanisms that proceed via a covalent aspartyl intermediate. More so, the widest dehalogenation spectrum of electron acceptors tested with bacterial strains which could dehalogenate recalcitrant organohalides has further proven the versatility of bacterial dehalogenators to be considered when determining the fate of halogenated organics at contaminated sites. In this review, the general features of most widely studied bacterial dehalogenases, their structural properties, basis of the degradation of organohalides and their derivatives and how they have been improved for various applications is discussed.

Fluoroacetate in plants - a review of its distribution, toxicity to livestock and microbial detoxification

Fluoroacetate producing plants grow worldwide and it is believed they produce this toxic compound as a defence mechanism against grazing by herbivores. Ingestion by livestock often results in fatal poisonings, which causes significant economic problems to commercial farmers in many countries such as Australia, Brazil and South Africa. Several approaches have been adopted to protect livestock from the toxicity with limited success including fencing, toxic plant eradication and agents that bind the toxin. Genetically modified bacteria capable of degrading fluoroacetate have been able to protect ruminants from fluoroacetate toxicity under experimental conditions but concerns over the release of these microbes into the environment have prevented the application of this technology. Recently, a native bacterium from an Australian bovine rumen was isolated which can degrade fluoroacetate. This bacterium, strain MFA1, which belongs to the Synergistetes phylum degrades fluoroacetate to fluoride ions and acetate. The discovery and isolation of this bacterium provides a new opportunity to detoxify fluoroacetate in the rumen. This review focuses on fluoroacetate toxicity in ruminant livestock, the mechanism of fluoroacetate toxicity, tolerance of some animals to fluoroaceate, previous attempts to mitigate toxicity, aerobic and anaerobic microbial degradation of fluoroacetate, and future directions to overcome fluoroacetate toxicity.

Catalytic mechanism of C–F bond cleavage: insights from QM/MM analysis of fluoroacetate dehalogenase

The catalytic mechanisms of fluoroacetate dehalogenase (FAcD) toward substrates fluoroacetate and chloroacetate were studied by a combined quantum mechanics/molecular mechanics (QM/MM) method.

Mechanistic insights into the dehalogenation reaction of fluoroacetate/fluoroacetic acid

Fluoroacetate is a toxic compound whose environmental accumulation may represent an important contamination problem, its elimination is therefore a challenging issue. Fluoroacetate dehalogenase catalyzes its degradation through a two step process initiated by an S(N)2 reaction in which the aspartate residue performs a nucleophilic attack on the carbon bonded to the fluorine; the second step is hydrolysis that releases the product as glycolate. In this paper, we present a study based on density functional theory calculations of the S(N)2 initiation reaction modeled through the interaction between the substrate and the propionate anion as the nucleophile. Results are analyzed within the framework of the reaction force and using the reaction electronic flux to identify and characterize the electronic activity that drives the reaction. Our results reveal that the selective protonation of the substrate catalyzes the reaction by decreasing the resistance of the structural and electronic reorganization needed to reach the transition state. Finally, the reaction energy is modulated by the degree of stabilization of the fluoride anion formed after the S(N)2 reaction. In this way, a site-induced partial protonation acts as a chemical switch in a key process that determines the output of the reaction.

Insights into enzymatic halogenation from computational studies

The halogenases are a group of enzymes that have only come to the fore over the last 10 years thanks to the discovery and characterization of several novel representatives. They have revealed the fascinating variety of distinct chemical mechanisms that nature utilizes to activate halogens and introduce them into organic substrates. Computational studies using a range of approaches have already elucidated many details of the mechanisms of these enzymes, often in synergistic combination with experiment. This Review summarizes the main insights gained from these studies. It also seeks to identify open questions that are amenable to computational investigations. The studies discussed herein serve to illustrate some of the limitations of the current computational approaches and the challenges encountered in computational mechanistic enzymology.

Theoretical study of the catalytic mechanism of E1 subunit of pyruvate dehydrogenase multienzyme complex from Bacillus stearothermophilus

Pyruvate dehydrogenase multienzyme complex (PDHc) is a member of a family of 2-oxo acid dehydrogenase (OADH) multienzyme complexes involved in several central points of oxidative metabolism, and the E1 subunit is the most important component in the entire PDHc catalytic system, which catalyzes the reversible transfer of an acetyl group from a pyruvate to the lipoyl group of E2 subunit lipoly domain. In this article, the catalytic mechanism of the E1 subunit has been systematically studied using density functional theory (DFT). Four possible pathways with different general acid/base catalysts in decarboxylation and reductive acylation processes were explored. Our calculation results indicate that the 4'-amino pyrimidine of ThDP and residue His128 are the most likely proton donors in the decarboxylation and reductive acylation processes, respectively. During the reaction, each C-C and C-S bond formation or cleavage process, except for the liberation of CO2, is always accompanied by a proton transfer between the substrates and proton donors. The liberation of CO2 is calculated to be the rate-limiting step for the overall reaction, with an energy barrier of 13.57 kcal/mol. The decarboxylation process is endothermic by 5.32 kcal/mol, whereas the reductive acylation process is exothermic with a value of 5.74 kcal/mol. The assignment of protonation states of the surrounding residues can greatly influence the reaction. Residues His128 and His271 play roles in positioning the first substrate pyruvate and second substrate lipoyl group, respectively.

Substrate specificity of fluoroacetate dehalogenase: an insight from crystallographic analysis, fluorescence spectroscopy, and theoretical computations

The high substrate specificity of fluoroacetate dehalogenase was explored by using crystallographic analysis, fluorescence spectroscopy, and theoretical computations. A crystal structure for the Asp104Ala mutant of the enzyme from Burkholderia sp. FA1 complexed with fluoroacetate was determined at 1.2 Å resolution. The orientation and conformation of bound fluoroacetate is different from those in the crystal structure of the corresponding Asp110Asn mutant of the enzyme from Rhodopseudomonas palustris CGA009 reported recently (J. Am. Chem. Soc. 2011, 133, 7461). The fluorescence of the tryptophan residues of the wild-type and Trp150Phe mutant enzymes from Burkholderia sp. FA1 incubated with fluoroacetate and chloroacetate was measured to gain information on the environment of the tryptophan residues. The environments of the tryptophan residues were found to be different between the fluoroacetate- and chloroacetate-bound enzymes; this would come from different binding modes of these two substrates in the active site. Docking simulations and QM/MM optimizations were performed to predict favorable conformations and orientations of the substrates. The F atom of the substrate is oriented toward Arg108 in the most stable enzyme-fluoroacetate complex. This is a stable but unreactive conformation, in which the small O-C-F angle is not suitable for the S(N)2 displacement of the F(-) ion. The cleavage of the C-F bond is initiated by the conformational change of the substrate to a near attack conformation (NAC) in the active site. The second lowest energy conformation is appropriate for NAC; the C-O distance and the O-C-F angle are reasonable for the S(N) 2 reaction. The activation energy is greatly reduced in this conformation because of three hydrogen bonds between the leaving F atom and surrounding amino acid residues. Chloroacetate cannot reach the reactive conformation, due to the longer C-Cl bond; this results in an increase of the activation energy despite the weaker C-Cl bond.

Mapping the reaction coordinates of enzymatic defluorination

The carbon-fluorine bond is the strongest covalent bond in organic chemistry, yet fluoroacetate dehalogenases can readily hydrolyze this bond under mild physiological conditions. Elucidating the molecular basis of this rare biocatalytic activity will provide the fundamental chemical insights into how this formidable feat is achieved. Here, we present a series of high-resolution (1.15-1.80 Å) crystal structures of a fluoroacetate dehalogenase, capturing snapshots along the defluorination reaction: the free enzyme, enzyme-fluoroacetate Michaelis complex, glycolyl-enzyme covalent intermediate, and enzyme-product complex. We demonstrate that enzymatic defluorination requires a halide pocket that not only supplies three hydrogen bonds to stabilize the fluoride ion but also is finely tailored for the smaller fluorine halogen atom to establish selectivity toward fluorinated substrates. We have further uncovered dynamics near the active site which may play pivotal roles in enzymatic defluorination. These findings may ultimately lead to the development of novel defluorinases that will enable the biotransformation of more complex fluorinated organic compounds, which in turn will assist the synthesis, detoxification, biodegradation, disposal, recycling, and regulatory strategies for the growing markets of organofluorines across major industrial sectors.

The assessment of density functionals for DNA–protein stacked and T-shaped complexes

The present work uses 129 nucleobase – amino acid CCSD(T)/CBS stacking and T-shaped interaction energies as reference data to test the ability of various density functionals with double-zeta quality basis sets, as well as some semi-empirical and molecular mechanics methods, to accurately describe noncovalent DNA–protein π–π and π+–π interactions. The goal of this work is to identify methods that can be used in hybrid approaches (QM/MM, ONIOM) for large-scale modeling of enzymatic systems involving active-site (substrate) π–π contacts. Our results indicate that AMBER is a more appropriate choice for the lower-level method in hybrid techniques than popular semi-empirical methods (AM1, PM3), and suggest that AMBER accurately describes the π–π interactions found throughout DNA–protein complexes. The M06–2X and PBE-D density functionals were found to provide very promising descriptions of the 129 nucleobase – amino acid interaction energies, which suggests that these may be the most suitable methods for describing high-level regions. Therefore, M06–2X and PBE-D with both the 6–31G(d) and 6–31+G(d,p) basis sets were further examined through potential-energy surface scans to better understand how these techniques describe DNA–protein π–π interactions in both minimum and nonminimum regions of the potential-energy surfaces, which is critical information when modeling enzymatic reaction pathways. Our results suggest that studies of stacked nucleobase – amino acid systems should implement the PBE-D/6–31+G(d,p) method. However, if T-shaped contacts are involved and (or) smaller basis sets must be considered due to limitations in computational resources, then M06–2X/6–31G(d) provides an overall excellent description of both nucleobase – amino acid stacking and T-shaped interactions for a range of DNA–protein π–π and π+–π interactions.