Non-enzymatic and enzymatic hydrolysis of alkyl halides: A haloalkane dehalogenation enzyme evolved to stabilize the gas-phase transition state of an SN2 displacement reaction

The semiempirical PM3 method, calibrated against ab initio HF/6–31+G(d) theory, has been used to elucidate the reaction of 1,2-dichloroethane (DCE) with the carboxylate of Asp-124 at the active site of haloalkane dehalogenase of Xanthobacter autothropicus. Asp-124 and 13 other amino acid side chains that make up the active site cavity (Glu-56, Trp-125, Phe-128, Phe-172, Trp-175, Leu-179, Val-219, Phe-222, Pro-223, Val-226, Leu-262, Leu-263, and His-289) were included in the calculations. The three most significant observations of the present study are that: (i) the DCE substrate and Asp-124 carboxylate, in the reactive ES complex, are present as an ion-molecule complex with a structure similar to that seen in the gas-phase reaction of AcO− with DCE; (ii) the structures of the transition states in the gas-phase and enzymatic reaction are much the same where the structure formed at the active site is somewhat exploded; and (iii) the enthalpies in going from ground states to transition states in the enzymatic and gas-phase reactions differ by only a couple kcal/mol. The dehalogenase derives its catalytic power from: (i) bringing the electrophile and nucleophile together in a low-dielectric environment in an orientation that allows the reaction to occur without much structural reorganization; (ii) desolvation; and (iii) stabilizing the leaving chloride anion by Trp-125 and Trp-175 through hydrogen bonding.

A family of phase-variable restriction enzymes with differing specificities generated by high-frequency gene rearrangements

The hsd genes of Mycoplasma pulmonis encode restriction and modification enzymes exhibiting a high degree of sequence similarity to the type I enzymes of enteric bacteria. The S subunits of type I systems dictate the DNA sequence specificity of the holoenzyme and are required for both the restriction and the modification reactions. The M. pulmonis chromosome has two hsd loci, both of which contain two hsdS genes each and are complex, site-specific DNA inversion systems. Embedded within the coding region of each hsdS gene are a minimum of three sites at which DNA inversions occur to generate extensive amino acid sequence variations in the predicted S subunits. We show that the polymorphic hsdS genes produced by gene rearrangement encode a family of functional S subunits with differing DNA sequence specificities. In addition to creating polymorphisms in hsdS sequences, DNA inversions regulate the phase-variable production of restriction activity because the other genes required for restriction activity (hsdR and hsdM) are expressed only from loci that are oriented appropriately in the chromosome relative to the hsd promoter. These data cast doubt on the prevailing paradigms that restriction systems are either selfish or function to confer protection from invasion by foreign DNA.

RNA analysis by ion-pair reversed-phase high performance liquid chromatography

Ion-pair reversed-phase high performance liquid chromatography (IP RP HPLC) is presented as a new, superior method for the analysis of RNA. IP RP HPLC provides a fast and reliable alternative to classical methods of RNA analysis, including separation of different RNA species, quantification and purification. RNA is stable under the analysis conditions used; degradation of RNA during the analyses was not observed. The versatility of IP RP HPLC for RNA analysis is demonstrated. Components of an RNA ladder, ranging in size from 155 to 1770 nt, were resolved. RNA transcripts of up to 5219 nt were analyzed, their integrity determined and they were quantified and purified. Purification of mRNA from total RNA is described, separating mouse rRNA from poly(A)+ mRNA. IP RP HPLC is also suitable for the separation and purification of DIG-labeled from unlabeled RNA. RNA purified by IP RP HPLC exhibits improved stability.

Structure of the ripple phase in lecithin bilayers.

The phases of the x-ray form factors are derived for the ripple (Pbeta') thermodynamic phase in the lecithin bilayer system. By combining these phases with experimental intensity data, the electron density map of the ripple phase of dimyristoyl-phosphatidylcholine is constructed. The phases are derived by fitting the intensity data to two-dimensional electron density models, which are created by convolving an asymmetric triangular ripple profile with a transbilayer electron density profile. The robustness of the model method is indicated by the result that many different models of the transbilayer profile yield essentially the same phases, except for the weaker, purely ripple (0,k) peaks. Even with this residual ambiguity, the ripple profile is well determined, resulting in 19 angstroms for the ripple amplitude and 10 degrees and 26 degrees for the slopes of the major and the minor sides, respectively. Estimates for the bilayer head-head spacings show that the major side of the ripple is consistent with gel-like structure, and the minor side appears to be thinner with lower electron density.

Progress toward chemical accuracy in the computer simulation of condensed phase reactions.

We describe a procedure for the generation of chemically accurate computer-simulation models to study chemical reactions in the condensed phase. The process involves (i) the use of a coupled semiempirical quantum and classical molecular mechanics method to represent solutes and solvent, respectively; (ii) the optimization of semiempirical quantum mechanics (QM) parameters to produce a computationally efficient and chemically accurate QM model; (iii) the calibration of a quantum/classical microsolvation model using ab initio quantum theory; and (iv) the use of statistical mechanical principles and methods to simulate, on massively parallel computers, the thermodynamic properties of chemical reactions in aqueous solution. The utility of this process is demonstrated by the calculation of the enthalpy of reaction in vacuum and free energy change in aqueous solution for a proton transfer involving methanol, methoxide, imidazole, and imidazolium, which are functional groups involved with proton transfers in many biochemical systems. An optimized semiempirical QM model is produced, which results in the calculation of heats of formation of the above chemical species to within 1.0 kcal/mol (1 kcal = 4.18 kJ) of experimental values. The use of the calibrated QM and microsolvation QM/MM (molecular mechanics) models for the simulation of a proton transfer in aqueous solution gives a calculated free energy that is within 1.0 kcal/mol (12.2 calculated vs. 12.8 experimental) of a value estimated from experimental pKa values of the reacting species.

Gas-phase nitronium ion affinities.

Evaluation of nitronium ion-transfer equilibria, L1NO2+ + L2 = L2NO2+ + L1 (where L1 and L2 are ligands 1 and 2, respectively) by Fourier-transform ion cyclotron resonance mass spectrometry and application of the kinetic method, based on the metastable fragmentation of L1(NO2+)L2 nitronium ion-bound dimers led to a scale of relative gas-phase nitronium ion affinities. This scale, calibrated to a recent literature value for the NO2+ affinity of water, led for 18 ligands, including methanol, ammonia, representative ketones, nitriles, and nitroalkanes, to absolute NO2+ affinities, that fit a reasonably linear general correlation when plotted vs. the corresponding proton affinities (PAs). The slope of the plot depends to a certain extent on the specific nature of the ligands and, hence, the correlations between the NO2+ affinities, and the PAs of a given class of compounds display a better linearity than the general correlation and may afford a useful tool for predicting the NO2+ affinity of a molecule based on its PA. The NO2+ binding energies are considerably lower than the corresponding PAs and well below the binding energies of related polyatomic cations, such as NO+, a trend consistent with the available theoretical results on the structure and the stability of simple NO2+ complexes. The present study reports an example of extension of the kinetic method to dimers, such as L1(NO2+)L2, bound by polyatomic ions, which may considerably widen its scope. Finally, measurement of the NO2+ affinity of ammonia allowed evaluation of the otherwise inaccessible PA of the amino group of nitramide and, hence, direct experimental verification of previous theoretical estimates.