The Catalytic Mechanism of Fluoroacetate Dehalogenase: A Computational Exploration of Biological Dehalogenation

The biological dehalogenation of fluoroacetate carried out by fluoroacetate dehalogenase is discussed by using quantum mechanical/molecular mechanical (QM/MM) calculations for a whole‐enzyme model of 10 800 atoms. Substrate fluoroacetate is anchored by a hydrogen‐bonding network with water molecules and the surrounding amino acid residues of Arg105, Arg108, His149, Trp150, and Tyr212 in the active site in a similar way to haloalkane dehalogenase. Asp104 is likely to act as a nucleophile to attack the α‐carbon of fluoroacetate, resulting in the formation of an ester intermediate, which is subsequently hydrolyzed by the nucleophilic attack of a water molecule to the carbonyl carbon atom. The cleavage of the strong C? F bond is greatly facilitated by the hydrogen‐bonding interactions between the leaving fluorine atom and the three amino acid residues of His149, Trp150, and Tyr212. The hydrolysis of the ester intermediate is initiated by a proton transfer from the water molecule to His271 and by the simultaneous nucleophilic attack of the water molecule. The transition state and produced tetrahedral intermediate are stabilized by Asp128 and the oxyanion hole composed of Phe34 and Arg105.     

Bacterial hydrolytic dehalogenases and related enzymes: occurrences, reaction mechanisms, and applications

Dehalogenases catalyze the cleavage of the carbon-halogen bond of organohalogen compounds. They have been attracting a great deal of attention partly because of their potential applications in the chemical industry and bioremediation. In this personal account, we describe occurrences, reaction mechanisms, and applications of bacterial hydrolytic dehalogenases and related enzymes, particularly L-2-haloacid dehalogenase, DL-2-haloacid dehalogenase, fluoroacetate dehalogenase, and 2-haloacrylate reductase. L-2-Haloacid dehalogenase is a representative enzyme of the haloacid dehalogenase (HAD) superfamily, which includes the P-type ATPases and other hydrolases. Structural and mechanistic analyses of this enzyme have yielded important insights into the mode of action of the HAD superfamily proteins. Fluoroacetate dehalogenase is unique in that it catalyzes the cleavage of the highly stable C--F bond of a fluorinated aliphatic compound. In the reactions of L-2-haloacid dehalogenase and fluoroacetate dehalogenase, the carboxylate group of Asp performs a nucleophilic attack on the alpha-carbon atom of the substrate, displacing the halogen atom. This mechanism is common to haloalkane dehalogenase and 4-chlorobenzoyl-CoA dehalogenase. DL-2-Haloacid dehalogenase is unique in that a water molecule directly attacks the substrate, displacing the halogen atom. The occurrence of 2-haloacrylate reductase was recently reported, revealing a new pathway for the degradation of unsaturated aliphatic organohalogen compounds.     

Tunneling motions of argon on chlorofluoromethane

The rotational supersonic jet Fourier transform microwave spectra of the 35Cl and 37Cl species of the molecular complex chlorofluoromethane-argon show that, in its equilibrium conformation, the argon atom is located out of the ClCF plane, interacting with the F and Cl atoms. All rotational transitions are split into several quadrupole components, each of them further split into two lines, due to the tunneling motion of the Ar atom between two equivalent positions, below and above the ClCF plane. The feasible low energy pathway between the structurally degenerate conformations is described, in a first approximation, by a circular motion around the C-Cl bond, with barriers estimated to be about 61 and 100 cm(-1).     

Reaction mechanism of soluble epoxide hydrolase: insights from molecular dynamics simulations

Molecular dynamics simulations have been performed to gain insights into the catalytic mechanism of the hydrolysis of epoxides to vicinal diols by soluble epoxide hydrolase (sEH). The binding of a substrate, 1S,2S-trans-methylstyrene oxide, was studied in two conformations in the active site of the enzyme. It was found that only one is likely to be found in the active enzyme. In the preferred conformation the phenyl group of the substrate is pi-sandwiched between two aromatic residues, Tyr381 and His523, whereas the other conformation is pi-stacked with only one aromatic residue, Trp334. Two simulations were carried out to 1 ns for each conformation to evaluate the protonation state of active site residue His523. It was found that a protonated histidine is essential for keeping the active site from being disrupted. Long time scale, 4 ns, molecular dynamics simulation was done for the structure with the most likely combination of binding conformation and protonation state of His523. Near Attack Conformers (NACs) are present 5.3% of the time and nucleophilic attack on either epoxide carbon atom, approximately 75% on C(1) and approximately 25% on C(2), is found. A maximum of one hydrogen bond between the epoxide oxygen and either of the active site tyrosines, Tyr465 and Tyr381, is present, in agreement with experimental mutagenesis results that reveal a slight loss in activity if one tyrosine is mutated and essential loss of all activity upon double mutation of the two tyrosines in question. It was found that a hydrogen bond from Tyr465 to the substrate oxygen is essential for controlling the regioselectivity of the reaction. Furthermore, a relationship between the presence of this hydrogen bond and the separation of reactants was found. Two groups of amino acid segments were identified each as moving collectively. Furthermore, an overall anti-correlation was found between the movements of these two individually collectively moving groups, made up by parts of the cap-region, including the two tyrosines, and the site of the catalytic triad, respectively. This overall anti-correlated collective domain motion is, perhaps, involved in the conversion of E.NAC to E.TS.     

Reaction of N-chloromethyldiethylamine,N-chloromethylbenzamide,N-chloromethylphthalimide with amides of trivalent phosphorus acids

Conclusions In the series of N-chloromethylphthalimide, N-chloromethylbenzamide, N-chloromethyldiethylamine, in reactions with amides of trivalent phosphorus acids, the contribution of the reaction at the phosphorus atom increases with increase in the polarity of the C-Cl bond. N-Chloromethylphthalimide, whose C-Cl bond is least polar in this series, reacts with hexaethyltriamidophosphite and diethylamidodiphenylphosphinite reacts at the N and P atoms to form the corresponding quasiphosphonium salts and acid chlorides of trivalent phosphorus acids.Translated from Izvestiya Akademii Nauk SSSR, Seriya Khimicheskaya, No. 1, pp. 173–178, January, 1985.     

Improving the Thermostability of a Methyl Parathion Hydrolase by Adding the Ionic Bond on Protein Surface

The thermostability of the methyl parathion hydrolase (MPH_OCH) from Ochrobactrum sp. M231 was improved using site-directed mutagenesis. Two prolines (Pro76 and Pro78) located on the protein surface were selected for mutations after inspection of the sequence alignment of MPH_OCH and OPHC2, a thermostable organophosphorus hydrolase from Pseudomonas pseudoalcaligenes C2-1. The temperature of the double-point mutant (P76D/P78K) at which the mutant lost 50% of its activity (T50) was approximately 68 °C, which is higher than that of WT enzyme (64 °C), P76D (67 °C), and P78K (59 °C). Structural analysis of P76D/P78K indicated that the substituted residues (Asp76 and Lys78) could generate an ionic bond and increase the structural electrostatic energy, which could then increase the stability of the protein. These results also suggest that the thermal stability of proteins could be improved by adding the ionic bond on protein surface.     

Structures,vibrational frequencies,topologies, and energies of hydrogen bonds in cysteine‐formaldehyde complexes

The hydrogen bonding interactions between cysteine (Cys) and formaldehyde (FA) were studied with density functional theory regarding their geometries, energies, vibrational frequencies, and topological features of the electron density. The quantum theory of atoms in molecules and natural bond orbital analyses were employed to elucidate the interaction characteristics in the Cys‐FA complexes. The intramolecular hydrogen bonds (H‐bonds) formed between the hydroxyl and the N atom of cysteine moiety in some Cys‐FA complexes were strengthened because of the cooperativity. Most of intermolecular H‐bonds involve the O atom of cysteine/FA moiety as proton acceptors, while the strongest H‐bond involves the O atom of FA moiety as proton acceptor, which indicates that FA would rather accept proton than providing one. The H‐bonds formed between the CH group of FA and the S atom of cysteine in some complexes are so weak that no hydrogen bonding interactions exist among them. In most of complexes, the orbital interaction of H‐bond is predominant during the formation of complex. The electron density (ρ_b) and its Laplace (?²ρ_b) at the bond critical point significantly correlate with the H‐bond parameter δR, while a linearly relationship between the second‐perturbation energy E(2) and ρ_b has been found as well. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Why substituting the asparagine at position 35 in Bacillus circulans xylanase with an aspartic acid remarkably improves the enzymatic catalytic activity? A quantum chemistry-based calculation study

Xylanases from Bacillus circulans (BCX) are known as configuration-retaining glycoside hydrolases, which hydrolyze xylans with two glutamic acid residues (Glu78 and Glu172) serving as catalytic active residues according to a double displacement mechanism. Existing experimental researches show that mutating the asparagines (Asn) to aspartic acid (Asp) at position 35 next to Glu172 can obviously improve the catalytic activity of BCX. To better understand the inherent mechanism for the experimental finding, we performed quantum chemistry calculations on two model systems to mimic the catalyses of wild-type and mutant BCXs. Geometrical structures and relative energies of intermediates and transition states involved in the hydrolysis reactions are given in detail. It is found that in the wild-type model system Asn35 interacts with Glu172 via a loose hydrogen bond, while in the mutant model system Asp35 forms a very tight hydrogen bond with Glu172. The glycosidic bond cleavage is proposed to be the rate-determining step for the hydrolysis reaction, whose barrier varies from 98 to 65 kJ mol⁻¹ when Asn35 is replaced by Asp35, showing the presence of Asp35 remarkably reduces the energy demand for the hydrolysis reaction. The present result provides a theoretical elucidation for why a single amino acid substitution can importantly influences catalytic activity of BCX.     

Reaction of 1,1-dichloroethane with tellurium in the system hydrazine hydrate-potassium hydroxide

The reaction of 1,1-dichloroethane with tellurium in the system hydrazine hydrate-potassium hydroxide gave ethyltellanyl derivatives as a result of replacement of one chlorine atom in the substrate by tellurium, and of the other, by hydrogen. Probable mechanisms of reduction of the C-Cl bond and mass spectra of the products were considered. The mass spectra of ditelluroacetals revealed unusual rearrangement of the molecular ion, leading to the formation of Te-Te and new C-C bond.     

Reaction dynamics following electron capture of chlorofluorocarbon adsorbed on water cluster: a direct density functional theory molecular dynamics study

The electron capture dynamics of halocarbon and its water complex have been investigated by means of the full dimensional direct density functional theory molecular dynamics method in order to shed light on the mechanism of electron capture of a halocarbon adsorbed on the ice surface. The CF(2)Cl(2) molecule and a cyclic water trimer (H(2)O)(3) were used as halocarbon and water cluster, respectively. The dynamics calculation of CF(2)Cl(2) showed that both C-Cl bonds are largely elongated after the electron capture, while one of the Cl atoms is dissociated from CF(2)Cl(2) (-) as a Cl(-) ion. Almost all total available energy was transferred into the internal modes of the parent CF(2)Cl radical on the product state, while the relative translational energy of Cl(-) was significantly low due to the elongation of two C-Cl bonds. In the case of a halocarbon-water cluster system, the geometry optimization of neutral complex CF(2)Cl(2)(H(2)O)(3) showed that one of the Cl atoms interacts with n orbital of water molecules of trimer and the other Cl atom existed as a dangling Cl atom. After the electron capture, only one C-Cl bond (dangling Cl atom) was rapidly elongated, whereas the other C-Cl bond is silent during the reaction. The dangling Cl atom was directly dissociated from CF(2)Cl(2) (-)(H(2)O)(3) as Cl(-). The fast Cl(-) ion was generated from CF(2)Cl(2) (-)(H(2)O)(3) on the water cluster. The mechanism of the electron capture of halocarbon on water ice was discussed on the basis of the theoretical results.     

Computational study of ketosteroid isomerase: insights from molecular dynamics simulation of enzyme bound substrate and intermediate

Delta(5)-3-Ketosteroid Isomerase (KSI) catalyzes the isomerization of 5,6-unsaturated ketosteroids to their 4,5-unsaturated isomers at a rate approaching the diffusion limit. The isomerization reaction follows a two-step general acid-base mechanism starting with Asp38-CO(2)(-) mediated proton abstraction from a sp(3)-hybridized carbon atom, alpha to carbonyl group, providing a dienolate intermediate. In the second step, Asp38-CO(2)H protonates the C6 of the intermediate providing a 4,5-unsaturated ketosteroid. The details of the mechanism have been highly controversial despite several experimental and computational studies of this enzyme. The general acid-base catalysis has been proposed to involve either a catalytic diad or a cooperative hydrogen bond mechanism. In this paper, we report our results from the 1.5 nanosecond molecular dynamics (MD) simulation of enzyme bound natural substrate (E.S) and enzyme bound intermediate (E.In) solvated in a TIP3P water box. The final coordinates from our MD simulation strongly support the cooperative hydrogen bond mechanism. The MD simulation of E.S and E.In shows that both Tyr14 and Asp99 are hydrogen bonded to the O3 of the substrate or intermediate. The average hydrogen bonding distance between Tyr14-OH and O3 becomes shorter and exhibits less fluctuation on E.S --> E.In. We also observe dynamic motions of water moving in and out of the active site in the E.S structures. This free movement of water disappears in the E.In structures. The active site is shielded by hydrophobic residues, which come together and squeeze out the waters from the active site in the E.In complex.     

Characterizing the protonation states of the catalytic residues in apo and substrate-bound human T-cell leukemia virus type 1 protease

Human T-cell leukemia virus type 1 (HTLV-1) protease is an attractive target when developing inhibitors to treat HTLV-1 associated diseases. To study the catalytic mechanism and design novel HTLV-1 protease inhibitors, the protonation states of the two catalytic aspartic acid residues must be determined. Free energy simulations have been conducted to study the proton transfer reaction between the catalytic residues of HTLV-1 protease using a combined quantum mechanical and molecular mechanical (QM/MM) molecular dynamics simulation. The free energy profiles for the reaction in the apo-enzyme and in an enzyme – substrate complex have been obtained. In the apo-enzyme, the two catalytic residues are chemically equivalent and are expected to be both unprotonated. Upon substrate binding, the catalytic residues of HTLV-1 protease evolve to a singly protonated state, in which the OD1 of Asp32 is protonated and forms a hydrogen bond with the OD1 of Asp32′, which is unprotonated. The HTLV-1 protease–substrate complex structure obtained from this simulation can serve as the Michaelis complex structure for further mechanistic studies of HTLV-1 protease while providing a receptor structure with the correct protonation states for the active site residues toward the design of novel HTLV-1 protease inhibitors through virtual screening.     

Reaction of N-chloromethyldiethylamine,N-chloromethylbenzamide,and N-chloromethylphthalimide with esters of trivalent phosphorus thioacids

Conclusions N-Chlororaethylbenzamide and N-chloromethylphthalimide react with S-butyldiphenylthiophosphinite at the P atom. The reaction of S,S-dibutylphenyldithiophosphonite with N-chloromethylphthalimide proceeds preferentially at the P atom, and with N-chloromethylbenzamide at the S atom. With increase in the polarity of the C-Cl bond in the N-chloromethyldiethylamine > N-chloromethylbenzamide > N-chloromethylphthalimide series, the tendency to react at the S atom increases.Translated from Izvestiya Akademli Nauk SSSR, Seriya Khimicheskaya, No. 1, pp. 178–182, January, 1985.     

How dihydrofolate reductase facilitates protonation of dihydrofolate

Dihydrofolate Reductase (DHFR) catalyzes the reduction of dihydrofolate (H2F) to tetrahydrofolate. On the basis of 10-12.5 ns molecular dynamics simulations of two conformations (closed and occluded) of the ternary DHFR/NADPH/H2F complex from Escherichia coli and a free energy perturbation approach, we have calculated the pKa value for the N5 atom in H2F. Our results suggest that the N5 atom in H2F is responsible for the pH dependency of the catalyzed reaction, meaning that DHFR facilitates protonation of H2F by approximately 4 pKa units. The mechanism behind this increase is due to favorable electrostatic interactions between the Asp27 residue and a proton at the N5 atom. The electrostatic interactions are enhanced by a hydrophobic active site, which to a large extent is made hydrophobic by the M20 loop in DHFR. Moreover, we find that the conformation imposed on H2F by DHFR to some extent also favors protonation of the N5 atom. Our results add support to previous findings and suggestions by Callender and co-workers [e.g., Deng, J.; Callender, R. J. Am. Chem. Soc. 1998, 120, 7730-7737] and explain why mutation of Asp27 may lead to severely reduced activity at neutral pH.     

Low energy (0-10 eV) electron driven reactions in the halogenated organic acids CCl3COOH, CClF2COOH, and CF3CHNH2COOH (trifluoroalanine)

Negative ion formation following resonant electron attachment to the three title molecules is studied by means of a beam experiment with mass spectrometric detection of the anions. All three molecules exhibit a pronounced resonance in the energy range around 1 eV which decomposes by the loss of a neutral hydrogen atom thereby generating the closed shell anion (M-H)(-) (or RCOO(-)), a reaction which is also a common feature in the non-substituted organic acids. The two chlorine containing molecules CCl(3)COOH and CClF(2)COOH exhibit an additional strong and narrow resonance at very low energy (close to 0 eV) which decomposes by the cleavage of the C-Cl bond with the excess charge finally localised on either of the two fragments Cl(-) and (M-Cl)(-). This reaction is by two to three orders of magnitude more effective than hydrogen loss. Apart from these direct bond cleavages (C-Cl, O-H) resonant attachment of subexcitation electrons trigger additional remarkably complex unimolecular decompositions leading, e.g., to the formation of the bihalide ions ClHCl(-) and ClHF(-) from CCl(3)COOH and CClF(2)COOH, respectively, or the loss of a neutral CF(2) unit from trifluoroalanine thereby generating the fluoroglycine radical anion. These reactions require substantial rearrangement in the transitory negative ion, i.e., the cleavage of different bonds and formation of new bonds. F(-) from both chlorodifluoroacetic acid and trifluoroalanine is formed at comparatively low intensity (more than three orders of magnitude less than Cl(-) from the chlorine containing molecules) and predominantly within a broad resonant feature around 7-8 eV characterised as core excited resonance.     

Novel Bicyclic Lactams as XaaPro Type VI beta Turn Mimics: Design, Synthesis, and Evaluation

The design, enantioselective synthesis, and structural characterization of novel bicyclic lactams as peptide mimics of the type VI beta turn is described. The mimics duplicate the conformation of the backbone and disposition of the side-chain atoms of the central two residues of the turn. The Gly L-Pro mimic, lactam 6, was prepared in good overall yield starting from (S)-2-(2'-propenyl)proline. (1)H NMR spectroscopy defined the relative stereochemistry of the substituents and conformational characteristics of the six-membered ring of the lactam; X-ray crystallographic analysis confirmed the conformational and stereochemical assignment. Examination of the crystal structure of lactam 6 revealed that the central amide bond was twisted appreciably out of planarity. The twisting of the amide bond was attributed to angle strain resulting from the presence of the sp(2)-hybridized nitrogen atom at the junction of the two rings. Alkylation of the enolate of the N,N-dimethylformamidine derivative of lactam 6 with benzyl bromide afforded stereoselectively the formamidine 11, a mimic of an L-Phe L-Pro dipeptide in the type VI turn conformation. The efficient synthetic route to highly functionalized peptidomimetics such as 11 will prove highly useful in peptide structure-function studies.     

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.     

C-Cl bond fission dynamics and angular momentum recoupling in the 235 nm photodissociation of allyl chloride

The photodissociation dynamics of allyl chloride at 235 nm producing atomic Cl((2)P(J);J=1/2,3/2) fragments is investigated using a two-dimensional photofragment velocity ion imaging technique. Detection of the Cl((2)P(1/2)) and Cl((2)P(3/2)) products by [2+1] resonance enhanced multiphoton ionization shows that primary C-Cl bond fission of allyl chloride generates 66.8% Cl((2)P(3/2)) and 33.2% Cl((2)P(1/2)). The Cl((2)P(3/2)) fragments evidenced a bimodal translational energy distribution with a relative weight of low kinetic energy Cl((2)P(3/2))/high kinetic energy Cl((2)P(3/2)) of 0.097/0.903. The minor dissociation channel for C-Cl bond fission, producing low kinetic energy chlorine atoms, formed only chlorine atoms in the Cl((2)P(3/2)) spin-orbit state. The dominant C-Cl bond fission channel, attributed to an electronic predissociation that results in high kinetic energy Cl atoms, produced both Cl((2)P(1/2)) and Cl((2)P(3/2)) atomic fragments. The relative branching for this dissociation channel is Cl((2)P(1/2))/[Cl((2)P(1/2))+Cl((2)P(3/2))]=35.5%. The average fraction of available energy imparted into product recoil for the high kinetic energy products was found to be 59%, in qualitative agreement with that predicted by a rigid radical impulsive model. Both the spin-orbit ground and excited chlorine atom angular distributions were close to isotropic. We compare the observed Cl((2)P(1/2))/[Cl((2)P(1/2))+Cl((2)P(3/2))] ratio produced in the electronic predissociation channel of allyl chloride with a prior study of the chlorine atom spin-orbit states produced from HCl photodissociation, concluding that angular momentum recoupling in the exit channel at long interatomic distance determines the chlorine atom spin-orbit branching.     

Enzymatic coupling of specific peptides at nonspecific ligation sites: effect of Asp189Glu mutation in trypsin on substrate mimetic-mediated reactions

Two main drawbacks seriously restrict the synthetic value of proteases as reagents in peptide fragment coupling: (i) native proteolytic activity and, thus, risk of undesired peptide cleavage; (ii) limited enzyme specificities restricting the amino acid residues between which a peptide bond can be formed. While the latter can be overcome by the use of substrate mimetics achieving peptide bond formation at nonspecific ligation sites, the risk of proteolytic cleavage still remains and hinders the wide acceptance of this powerful strategy for peptide coupling. This paper reports on the effect of the trypsin point mutant Asp189Glu on substrate mimetic-mediated reactions. The effect of this mutation on the steady-state hydrolysis of substrate mimetics of the 4-guanidinophenyl ester type and on trypsin-specific Lys- and Arg-containing peptides was investigated. The results were confirmed by enzymatic coupling reactions using substrate mimetics as the acyl donor and specific amino acid-containing peptides as the acyl acceptor. The competition assay verifies the predicted shift in substrate preference from Lys and Arg to the substrate mimetics and, thus, from cleavage to synthesis of peptide bonds. The combination of results obtained qualifies the trypsin mutant D189E as the first substrate mimetic-specific peptide ligase.     

喹诺酮药物与血清蛋白相互作用的三维荧光光谱研究

应用三维荧光光谱和三维荧光偏振光谱研究了数种喹诺酮药物与牛血清白蛋白(BSA)分子间的相互作用。由三维荧光(偏振)光谱得到的指纹信息说明了喹诺酮药物与BSA结合反应对BSA分子构象的影响。通过研究喹诺酮药物发生相互作用前后BSA荧光偏振度及各向异性的变化,定量说明了喹诺酮药物-BSA所发生的结合反应。