101
|
Wood GPF, Radom L, Petersson GA, Barnes EC, Frisch MJ, Montgomery JA. A restricted-open-shell complete-basis-set model chemistry. J Chem Phys 2006; 125:094106. [PMID: 16965071 DOI: 10.1063/1.2335438] [Citation(s) in RCA: 164] [Impact Index Per Article: 9.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A restricted-open-shell model chemistry based on the complete basis set-quadratic Becke3 (CBS-QB3) model is formulated and denoted ROCBS-QB3. As the name implies, this method uses spin-restricted wave functions, both for the direct calculations of the various components of the electronic energy and for extrapolating the correlation energy to the complete-basis-set limit. These modifications eliminate the need for empirical corrections that are incorporated in standard CBS-QB3 to compensate for spin contamination when spin-unrestricted wave functions are used. We employ an initial test set of 19 severely spin-contaminated species including doublet radicals and both singlet and triplet biradicals. The mean absolute deviation (MAD) from experiment for the new ROCBS-QB3 model (3.6+/-1.5 kJ mol(-1)) is slightly smaller than that of the standard unrestricted CBS-QB3 version (4.8+/-1.5 kJ mol(-1)) and substantially smaller than the MAD for the unrestricted CBS-QB3 before inclusion of the spin correction (16.1+/-1.5 kJ mol(-1)). However, when applied to calculate the heats of formation at 298 K for the moderately spin-contaminated radicals in the G2/97 test set, ROCBS-QB3 does not perform quite as well as the standard unrestricted CBS-QB3, with a MAD from experiment of 3.8+/-1.6 kJ mol(-1) (compared with 2.9+/-1.6 kJ mol(-1) for standard CBS-QB3). ROCBS-QB3 performs marginally better than standard CBS-QB3 for the G2/97 set of ionization energies with a MAD of 4.1+/-0.1 kJ mol(-1) (compared with 4.4+/-0.1 kJ mol(-1)) and electron affinities with a MAD of 3.9+/-0.2 kJ mol(-1) (compared with 4.3+/-0.2 kJ mol(-1)), but the differences in MAD values are comparable to the experimental uncertainties. Our overall conclusion is that ROCBS-QB3 eliminates the spin correction in standard CBS-QB3 with no loss in accuracy.
Collapse
|
102
|
Corral I, Mó O, Yáñez M, Salpin JY, Tortajada J, Moran D, Radom L. An Experimental and Theoretical Investigation of Gas-Phase Reactions of Ca2+ with Glycine. Chemistry 2006; 12:6787-96. [PMID: 16807970 DOI: 10.1002/chem.200600127] [Citation(s) in RCA: 54] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
The gas-phase reactions between Ca(2+) and glycine ([Ca(gly)](2+)) have been investigated through the use of mass spectrometry techniques and B3-LYP/cc-pWCVTZ density functional theory computations. The major peaks observed in the electrospray MS/MS spectrum of [Ca(gly)](2+) correspond to the formation of the [Ca,C,O(2),H](+), NH(2)CH(2) (+), CaOH(+), and NH(2)CH(2)CO(+) fragment ions, which are produced in Coulomb explosion processes. The computed potential energy surface (PES) shows that not only are these species the most stable product ions from a thermodynamic point of view, but they may be produced with barriers lower than for competing processes. Carbon monoxide is a secondary product, derived from the unimolecular decomposition of some of the primary ions formed in the Coulomb explosions. In contrast to what is found for the reactions of Ca(2+) with urea ([Ca(urea)](2+)), minimal unimolecular losses of neutral fragments are observed for the gas-phase fragmentation processes of [Ca(gly)](2+), which is readily explained in terms of the topological differences between their respective PESs.
Collapse
|
103
|
Wood GPF, Easton CJ, Rauk A, Davies MJ, Radom L. Effect of Side Chains on Competing Pathways for β-Scission Reactions of Peptide-Backbone Alkoxyl Radicals. J Phys Chem A 2006; 110:10316-23. [PMID: 16928124 DOI: 10.1021/jp062916j] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
High-level quantum chemistry calculations have been carried out to investigate beta-scission reactions of alkoxyl radicals located at the alpha-carbon of a peptide backbone. This type of alkoxyl radical may undergo three possible beta-scission reactions, namely C-C beta-scission of the backbone, C-N beta-scission of the backbone, and C-R beta-scission of the side chain. We find that the rates for the C-C beta-scission reactions are all very fast, with rate constants of the order 10(12) s(-1) that are essentially independent of the side chain. The C-N beta-scission reactions are all slow, with rate constants that range from 10(-0.7) to 10(-4.5) s(-1). The rates of the C-R beta-scission reactions depend on the side chain and range from moderately fast (10(7) s(-1)) to very fast (10(12) s(-1)). The rates of the C-R beta-scission reactions correlate well with the relative stabilities of the resultant side-chain product radicals (*R), as reflected in calculated radical stabilization energies (RSEs). The order of stabilities for the side-chain fragment radicals for the natural amino acids is found to be Ala < Glu < Gln approximately Leu approximately Met approximately Lys approximately Arg < Asp approximately Ile approximately Asn approximately Val < Ser approximately Thr approximately Cys < Phe approximately Tyr approximately His approximately Trp. We predict that for side-chain C-R beta-scission reactions to effectively compete with the backbone C-C beta-scission reactions, the side-chain fragment radicals would generally need an RSE greater than approximately 30 kJ mol(-1). Thus, the residues that may lead to competitive side-chain beta-scission reactions are Ser, Thr, Cys, Phe, Tyr, His, and Trp.
Collapse
|
104
|
Hemelsoet K, Moran D, Van Speybroeck V, Waroquier M, Radom L. An Assessment of Theoretical Procedures for Predicting the Thermochemistry and Kinetics of Hydrogen Abstraction by Methyl Radical from Benzene. J Phys Chem A 2006; 110:8942-51. [PMID: 16836458 DOI: 10.1021/jp061823z] [Citation(s) in RCA: 57] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The reaction enthalpy (298 K), barrier (0 K), and activation energy and preexponential factor (600-800 K) have been examined computationally for the abstraction of hydrogen from benzene by the methyl radical, to assess their sensitivity to the applied level of theory. The computational methods considered include high-level composite procedures, including W1, G3-RAD, G3(MP2)-RAD, and CBS-QB3, as well as conventional ab initio and density functional theory (DFT) methods, with the latter two classes employing the 6-31G(d), 6-31+G(d,p) and/or 6-311+G(3df,2p) basis sets, and including ZPVE/thermal corrections obtained from 6-31G(d) or 6-31+G(d,p) calculations. Virtually all the theoretical procedures except UMP2 are found to give geometries that are suitable for subsequent calculation of the reaction enthalpy and barrier. For the reaction enthalpy, W1, G3-RAD, and URCCSD(T) give best agreement with experiment, while the large-basis-set DFT procedures slightly underestimate the endothermicity. The reaction barrier is slightly more sensitive to the choice of basis set and/or correlation level, with URCCSD(T) and the low-cost BMK method providing values in close agreement with the benchmark G3-RAD value. Inspection of the theoretically calculated rate parameters reveals a minor dependence on the level of theory for the preexponential factor. There is more sensitivity for the activation energy, with a reasonable agreement with experiment being obtained for the G3 methods and the hybrid functionals BMK, BB1K, and MPW1K, especially in combination with the 6-311+G(3df,2p) basis set. Overall, the high-level G3-RAD composite procedure, URCCSD(T), and the cost-effective DFT methods BMK, BB1K, and MPW1K give the best results among the methods assessed for calculating the thermochemistry and kinetics of hydrogen abstraction by the methyl radical from benzene.
Collapse
|
105
|
Sandala GM, Smith DM, Coote ML, Golding BT, Radom L. Insights into the hydrogen-abstraction reactions of diol dehydratase: relevance to the catalytic mechanism and suicide inactivation. J Am Chem Soc 2006; 128:3433-44. [PMID: 16522124 DOI: 10.1021/ja057902q] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
High-level quantum chemistry calculations have been used to examine the hydrogen-abstraction reactions of diol dehydratase (DDH) in the context of both the catalytic mechanism and the enzyme dysfunction phenomenon termed suicide inactivation. The barriers for the catalytic hydrogen-abstraction reactions of ethane-1,2-diol and propane-1,2-diol are examined in isolation, as well as in the presence of various Brønsted acids and bases. Modest changes in the magnitudes of the initial and final abstraction barriers are seen, depending on the strength of the acid or base, and on whether these effects are considered individually or together. The most significant changes (ca. 20 kJ mol(-1)) are found for the initial abstraction barrier when the spectator OH group is partially deprotonated. Kinetic isotope effects including Eckart tunneling corrections (KIEs) have also been calculated for these model systems. We find that contributions from tunneling are of a magnitude similar to that of the contributions from semiclassical theory alone, meaning that quantum effects serve to significantly accelerate the rate of hydrogen transfer. The calculated KIEs for the partially deprotonated system are in qualitative agreement with experimentally determined values. In complementary investigations, the ability of DDH to become deactivated by certain substrate analogues is examined. In all cases, the formation of a stable radical intermediate causes the hydrogen re-abstraction step to become an extremely endothermic process. The consequent inability of 5'-deoxyadenosyl radical to be regenerated breaks the catalytic cycle, resulting in the suicide inactivation of DDH.
Collapse
|
106
|
Barlow CK, Moran D, Radom L, McFadyen WD, O'Hair RAJ. Metal-Mediated Formation of Gas-Phase Amino Acid Radical Cations,. J Phys Chem A 2006; 110:8304-15. [PMID: 16821814 DOI: 10.1021/jp056471v] [Citation(s) in RCA: 62] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The results from an investigation of the collision-induced dissociation (CID) of the ternary complexes [Cu(II)(terpy)(AA)](2+) are presented (terpy = 2,2':6',2' '-terpyridine; AA = one of the twenty common amino acids). These complexes show a rich gas-phase chemistry, which depends on the identity of the amino acid. For the histidine-, lysine- and tryptophan-containing complexes, oxidative dissociation of the amino acid is observed, yielding the amino acid radical cation. The results of further mass selection and CID of these amino acid radical cations are presented. The CID of the series [Fe(III)(salen)(AA)](+) (where salen = N,N'-ethylenebis(salicylideneaminato)) is also examined. These complexes undergo loss of the neutral amino acid in all cases, although the radical cation of arginine is also produced and its subsequent fragmentation examined. B3-LYP/6-31G(d) computations were carried out to test aspects of the proposed fragmentation mechanism of the histidine and arginine radical cations.
Collapse
|
107
|
Chan B, Radom L. Design of Effective Zeolite Catalysts for the Complete Hydrogenation of CO2. J Am Chem Soc 2006; 128:5322-3. [PMID: 16620086 DOI: 10.1021/ja0602492] [Citation(s) in RCA: 52] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Ab initio molecular orbital calculations have been applied to the study of the three-stage zeolite-catalyzed hydrogenation of CO2 to methanol. The results present strong evidence that appropriate chemical modifications to ZSM-5 can lead to significantly lower energy barriers for the three component reactions, that is, hydrogenation of CO2, HCO2H, and CH2O. Zeolites incorporating either Na+ or Ge are more effective catalysts than conventional acidic zeolites for the hydrogenation of CO2 to give HCO2H, but amine-based zeolites do not lead to significantly lower barriers for any of the three hydrogenation reactions. However, we predict that when all three features, namely, Na+, N, and Ge, are incorporated in the zeolite, there is a dramatic improvement in catalytic activity for all three reactions.
Collapse
|
108
|
Wee S, Mortimer A, Moran D, Wright A, Barlow CK, O'Hair RAJ, Radom L, Easton CJ. Gas-phase regiocontrolled generation of charged amino acid and peptide radicals. Chem Commun (Camb) 2006:4233-5. [PMID: 17031442 DOI: 10.1039/b608724h] [Citation(s) in RCA: 51] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The combined use of advanced mass spectrometry experiments, condensed-phase synthesis of serine and homoserine nitrate ester radical precursors, and high-level ab initio calculations provides a powerful way of examining the fundamental reactivity of radicals derived from peptides.
Collapse
|
109
|
Haworth NL, Sullivan MB, Wilson AK, Martin JML, Radom L. Structures and Thermochemistry of Calcium-Containing Molecules. J Phys Chem A 2005; 109:9156-68. [PMID: 16332025 DOI: 10.1021/jp052889h] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A variety of theoretical procedures, including the high-level ab initio methods G3, G3[CC](dir,full), and W2C//ACQ, have been used to predict the structures and heats of formation of several small calcium-containing molecules (CaH, CaH2, CaO, CaOH, Ca(OH)2, CaF, CaF2, CaS, CaCl, and CaCl2). B3-LYP and CCSD(T) with both the (aug-)cc-pWCVQZ and (aug-)cc-pWCVQ+dZ basis sets are found to give molecular geometries that agree well with the experimental results. The CCSD(T)(riv)/(aug-)cc-pWCVQ+dZ results are found to be the most accurate, with a mean absolute deviation from experiment of just 0.008 angstroms. Zero-point vibrational energies (ZPVEs) and thermochemical corrections are found to be relatively insensitive to the level of theory, except in the case of molecules with highly anharmonic calcium-centered bending modes (CaH2, Ca(OH)2, CaF2, CaCl2), where special procedures need to be employed in order to obtain satisfactory results. Several potential improvements to the W2C method were investigated, most of which do not produce significant changes in the heats of formation. It was observed, however, that for CaO and CaS the scalar relativistic corrections are unexpectedly large and highly basis set dependent. In these cases, Douglas-Kroll CCSD(T)/(aug-)cc-pWCV5Z calculations appear to give a converged result. The G3[CC](dir,full) and best W2C-type heats of formation are both found generally to agree well with experimental values recommended in recent critical compendia. However, in some cases (CaO, Ca(OH)2, and CaF2), they differ from one another by more than their predicted error margins. The available experimental data are not sufficiently precise to distinguish definitively between the two sets of results although, in general, when discrepancies exist the W2C heats of formation are lower in energy and tend to be in better agreement with experiment. In the case of CaO, the W2C heat of formation (20.7 kJ mol(-1)) is approximately 20 kJ mol(-1) lower than the G3[CC](dir, full) result and most of the experimental data. Extensive investigation of possible refinements of the W2C method has failed to reveal any weaknesses that could account for this discrepancy. We therefore believe that the heat of formation of CaO is likely to lie closer to the more recent direct experimental determination of 27 kJ mol(-1) than to the value of approximately 40 kJ mol(-1) recommended in recent thermochemical reviews.
Collapse
|
110
|
Sandala GM, Smith DM, Radom L. Divergent mechanisms of suicide inactivation for ethanolamine ammonia-lyase. J Am Chem Soc 2005; 127:8856-64. [PMID: 15954794 DOI: 10.1021/ja051527k] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Ab initio molecular orbital calculations have been used to study the mechanism of suicide inactivation of ethanolamine ammonia-lyase induced by three different substrate analogues. Analysis of the normal catalytic mechanism with 2-aminoethanol (ethanolamine) as substrate predicts that both the hydrogen-abstraction and hydrogen-reabstraction steps involving the B(12)-cofactor are likely to be exothermic. On the other hand, the proposed inactivation mechanism for the first substrate analogue, glycolaldehyde, leads to a highly stabilized radical that results in a very endothermic (by ca. 90 kJ mol(-)(1)) hydrogen-reabstraction step, which is thought to halt the normal function of the enzyme. Curiously, the energy requirements for a catalytically imposed mechanism in the case of the second substrate analogue, 2-hydroxyethylhydrazine (HEH), parallel those for the catalytic substrate, despite the fact that HEH is found to inactivate EAL experimentally. However, further analysis reveals the presence of a lower energy pathway for HEH that leads to the formation of the highly stabilized hydrazinium radical cation. In a manner similar to when glycolaldehyde is the substrate analogue, this results in an endothermicity for the hydrogen-reabstraction step that is prohibitively large. In contrast to these related inactivation mechanisms, the third substrate analogue, 2-aminoacetaldehyde, apparently accomplishes the inactivation of EAL in an entirely different manner. A pathway for the experimentally observed formation of acetic acid and ammonium cation has been identified and appears catalytic in the sense that 5'-deoxyadenosyl radical is regenerated. However, mechanisms to account for the subsequent formation of 4',5'-anhydroadenosine and degradation of the corrinoid ring of the cofactor have not been elucidated.
Collapse
|
111
|
Izgorodina EI, Coote ML, Radom L. Trends in R−X Bond Dissociation Energies (R = Me, Et, i-Pr, t-Bu; X = H, CH3, OCH3, OH, F): A Surprising Shortcoming of Density Functional Theory. J Phys Chem A 2005; 109:7558-66. [PMID: 16834125 DOI: 10.1021/jp052021r] [Citation(s) in RCA: 181] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The performance of a variety of high-level composite procedures, as well as lower-cost density functional theory (DFT)- and second-order perturbation theory (MP2)-based methods, for the prediction of absolute and relative R-X bond dissociation energies (BDEs) was examined for R = Me, Et, i-Pr and t-Bu, and X = H, CH(3), OCH(3), OH and F. The methods considered include the high-level G3(MP2)-RAD and G3-RAD procedures, a variety of pure and hybrid DFT methods (B-LYP, B3-LYP, B3-P86, KMLYP, B1B95, MPW1PW91, MPW1B95, BB1K, MPW1K, MPWB1K and BMK), standard restricted (open-shell) MP2 (RMP2), and two recently introduced variants of MP2, namely spin-component-scaled MP2 (SCS-MP2) and scaled-opposite-spin MP2 (SOS-MP2). The high-level composite procedures show very good agreement with experiment and are used to evaluate the performance of the lower-level DFT- and MP2-based procedures. The best DFT methods (KMLYP and particularly BMK) provide very reasonable predictions for the absolute heats of formation and R-X BDEs for the systems studied. However, all of the DFT methods overestimate the stabilizing effect on BDEs in going from R = Me to R = t-Bu, leading in some cases to incorrect qualitative behavior. In contrast, the MP2-based methods generally show larger errors (than the best DFT methods) in the absolute heats of formation and BDEs, but better behavior for the relative BDEs, although they do tend to underestimate the stabilizing effect on BDEs in going from R = Me to R = t-Bu. The potentially less computationally expensive SOS-MP2 method offers particular promise as a reliable method that might be applicable to larger systems.
Collapse
|
112
|
Corral I, Mó O, Yáñez M, Radom L. Why Are the Ca2+ and K+ Binding Energies of Formaldehyde and Ammonia Reversed with Respect to Their Proton Affinities? J Phys Chem A 2005; 109:6735-42. [PMID: 16834027 DOI: 10.1021/jp051052j] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The binding energies (BEs) of alkali metal monocations and alkaline-earth metal dications to a series of small oxygen and nitrogen bases have been evaluated by means of CCSD(T) calculations on B3-LYP optimized structures. These calculations were carried out both using all-electron basis sets, and additionally using an effective-core potential (ECP) to describe the inner electrons of the metal. A theoretical model aiming at analyzing the effects on the binding energy trends of electrostatic, polarization, and covalent contributions, as well as geometry distortion, was employed. From this analysis, we conclude that although the neutral-ion interaction energy for alkali and alkaline-earth metal cations is dominated by electrostatic contributions, in many cases the correct basicity trends are only attained once polarization effects are also included in the model. This is indeed the case when Ca2+ and K+ are bound to ammonia and formaldehyde. Geometry distortions triggered by polarization are also necessary, in some cases, to obtain the correct basicity trends. In addition, in particular for alkaline-earth dications, the energy associated with covalent interactions sometimes dictates the basicity trend. Our observations imply that simple models based on ion-dipole interactions, that are frequently used in the literature to explain affinity trends in ion-molecule reactions, are generally not likely to be reliable.
Collapse
|
113
|
Wood GPF, Moran D, Jacob R, Radom L. Bond Dissociation Energies and Radical Stabilization Energies Associated with Model Peptide-Backbone Radicals. J Phys Chem A 2005; 109:6318-25. [PMID: 16833974 DOI: 10.1021/jp051860a] [Citation(s) in RCA: 65] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Bond dissociation energies (BDEs) and radical stabilization energies (RSEs) have been calculated for a series of models that represent a glycine-containing peptide-backbone. High-level methods that have been used include W1, CBS-QB3, U-CBS-QB3, and G3X(MP2)-RAD. Simpler methods used include MP2, B3-LYP, BMK, and MPWB1K in association with the 6-311+G(3df,2p) basis set. We find that the high-level methods produce BDEs and RSEs that are in good agreement with one another. Of the simpler methods, RBMK and RMPWB1K achieve good accuracy for BDEs and RSEs for all the species that were examined. For monosubstituted carbon-centered radicals, we find that the stabilizing effect (as measured by RSEs) of carbonyl substituents (CX=O) ranges from 24.7 to 36.9 kJ mol(-1), with the largest stabilization occurring for the CH=O group. Amino groups (NHY) also stabilize a monosubstituted alpha-carbon radical, with the calculated RSEs ranging from 44.5 to 49.5 kJ mol(-1), the largest stabilization occurring for the NH2 group. In combination, NHY and CX=O substituents on a disubstituted carbon-centered radical produce a large stabilizing effect ranging from 82.0 to 125.8 kJ mol(-1). This translates to a captodative (synergistic) stabilization of 12.8 to 39.4 kJ mol(-1). For monosubstituted nitrogen-centered radicals, we find that the stabilizing effect of methyl and related (CH2Z) substituents ranges from 25.9 to 31.7 kJ mol(-1), the largest stabilization occurring for the CH3 group. Carbonyl substituents (CX=O) destabilize a nitrogen-centered radical relative to the corresponding closed-shell molecule, with the calculated RSEs ranging from -30.8 to -22.3 kJ mol(-1), the largest destabilization occurring for the CH=O group. In combination, CH2Z and CX=O substituents at a nitrogen radical center produce a destabilizing effect ranging from -19.0 to -0.2 kJ mol(-1). This translates to an additional destabilization associated with disubstitution of -18.6 to -7.8 kJ mol(-1).
Collapse
|
114
|
Chan B, Del Bene JE, Elguero J, Radom L. On the Relationship between the Preferred Site of Hydrogen Bonding and Protonation. J Phys Chem A 2005; 109:5509-17. [PMID: 16839079 DOI: 10.1021/jp0516994] [Citation(s) in RCA: 49] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
Ab initio molecular orbital calculations have been employed to investigate the interactions between a set of basic substrates (B) with H+ and HF, and the interaction between acids of varying strength (AH+) with two bases, vinylamine and furan. The preferred site for protonation of the substrates appears to be determined primarily by the ability of the protonated species (BH+) to delocalize the acquired positive charge. On the other hand, localization of a pair of electrons at a proton-acceptor site of B tends to be more important in determining the preferred site for hydrogen bonding with HF. The behavior of acids stronger than HF lies between these extremes. Consistent with a previously proposed Hammond postulate for complexes, when a substrate (B) interacts with a range of acids (AH+), proton transfer is generally found to occur when the proton affinity of A is significantly less than that of B. When the proton affinity of A is greater than that of B, a hydrogen-bonded complex is generally formed without proton transfer. Strongest binding (relative to the lowest energy components) occurs when the proton affinities of A and B are comparable. Proton transfer from AH+ is found to take place in some cases when this would not be predicted on the basis of protonation energies alone, because of specific interactions in the resulting complexes.
Collapse
|
115
|
Chan B, Radom L. Base-catalyzed hydrogenation: rationalizing the effects of catalyst and substrate structures and solvation. J Am Chem Soc 2005; 127:2443-54. [PMID: 15724999 DOI: 10.1021/ja0450253] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Ab initio molecular orbital calculations have been used to study the base-catalyzed hydrogenation of carbonyl compounds. It is found that these hydrogenation reactions share many common features with S(N)2 reactions. Both types of reactions are described by double-well energy profiles, with deep wells and a low or negative overall energy barrier in the gas phase, while the solution-phase profiles show very shallow wells and much higher barriers. For the hydrogenation reactions, the assembly of the highly ordered transition structure is found to be a major limiting factor to the rate of reaction. In the gas phase, the overall barriers for reactions catalyzed by Group I methoxides increase steadily down the group, due to the decreasing charge density on the metal. On the other hand, for Group II and Group III metals, the overall barriers decrease down the group, which is attributed to the increasing ionic character of the metal-oxygen bond. The reaction with B(OCH(3))(3) has an exceptionally high barrier, which is attributed to pi-electron donation from the oxygen lone pairs of the methoxy groups to the formally vacant p orbital on B, as well as to the high covalent character of the B-O bonds. In solution, these reactivity trends are generally the opposite of the corresponding gas-phase trends. While similar barriers are obtained for reactions catalyzed by methoxides and by tert-butoxides, reactions with benzyloxides have somewhat higher barriers. Aromatic ketones are found to be more reactive than purely aliphatic ketones. Moreover, comparison between catalytic hydrogenation of 2,2,5,5-tetramethylcyclopentanone and pivalophenone shows that factors such as steric effects may also be important in differentiating their reactivity. Solvation studies with a wide range of solvents indicate a steady decrease in barrier with decreasing solvent dielectric constant, with nonpolar solvents generally leading to considerably lower barriers than polar solvents. In practice, a good balance between polarity and catalyst solubility is required in selecting the most suitable solvent for the base-catalyzed hydrogenation reaction.
Collapse
|
116
|
Sandala GM, Smith DM, Coote ML, Radom L. Suicide inactivation of dioldehydratase by glycolaldehyde and chloroacetaldehyde: an examination of the reaction mechanism. J Am Chem Soc 2004; 126:12206-7. [PMID: 15453720 DOI: 10.1021/ja047377f] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
High-level ab initio calculations have been used to study the mechanism for the inactivation of diol dehydratase (DDH) by glycolaldehyde or 2-chloroacetaldehyde. As in the case of the catalytic substrates of DDH, e.g., ethane-1,2-diol, the 5'-deoxyadenosyl radical (Ado*) is able to abstract a hydrogen atom from both substrate analogues in the initial step on the reaction pathway, as evidenced by comparable energy barriers. However, in subsequent step(s), each substrate analogue produces the highly stable glycolaldehyde radical. The barrier for hydrogen atom reabstraction by the glycolaldehyde radical is calculated to be too high ( approximately 110 kJ mol-1) to allow Ado* to be regenerated and recombine with the cob(II)alamin radical, the latter therefore remaining tightly bound to DDH. As a consequence, the catalytic pathway is disrupted, and DDH becomes an impotent enzyme. Interconversion of equivalent structures of the glycolaldehyde radical via the symmetrical cis-ethanesemidione radical is calculated to require 38 kJ mol-1. EPR indications of a symmetrical cis-ethanesemidione structure are likely to be the result of formation of an equilibrium mixture of glycolaldehyde radical structures, this equilibration being facilitated by partial deprotonation of the glycolaldehyde radical by the carboxylate of an amino acid residue within the active site of DDH.
Collapse
|
117
|
Corral I, Mó O, Yáñez M, Salpin JY, Tortajada J, Radom L. Gas-Phase Reactions between Urea and Ca2+: The Importance of Coulomb Explosions. J Phys Chem A 2004. [DOI: 10.1021/jp046624z] [Citation(s) in RCA: 44] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
118
|
Barratt BJW, Easton CJ, Henry DJ, Li IHW, Radom L, Simpson JS. Inhibition of Peptidylglycine α-Amidating Monooxygenase by Exploitation of Factors Affecting the Stability and Ease of Formation of Glycyl Radicals. J Am Chem Soc 2004; 126:13306-11. [PMID: 15479085 DOI: 10.1021/ja046204n] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Peptidylglycine alpha-amidating monooxygenase catalyzes the biosynthesis of peptide hormones through radical cleavage of the C-terminal glycine residues of the corresponding prohormones. We have correlated ab initio calculations of radical stabilization energies and studies of free radical brominations with the extent of catalysis displayed by peptidylglycine alpha-amidating monooxygenase, to identify classes of inhibitors of the enzyme. In particular we find that, in closely related systems, the substitution of glycolate for glycine reduces the calculated radical stabilization energy by 34.7 kJ mol(-1), decreases the rate of bromination with N-bromosuccinimide at reflux in carbon tetrachloride by a factor of at least 2000, and stops catalysis by the monooxygenase, while maintaining binding to the enzyme.
Collapse
|
119
|
Topf M, Sandala GM, Smith DM, Schofield CJ, Easton CJ, Radom L. The Unusual Bifunctional Catalysis of Epimerization and Desaturation by Carbapenem Synthase. J Am Chem Soc 2004; 126:9932-3. [PMID: 15303862 DOI: 10.1021/ja047899v] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
High-level ab initio calculations have been used to study the mechanism for the conversion of (3S,5S)-carbapenam to the biologically active beta-lactam antibiotic, (5R)-carbapenem, catalyzed by carbapenem synthase. This process involves epimerization at C5 and desaturation at C2/C3. Our calculations suggest that the reaction proceeds via initial abstraction of the C5 hydrogen atom, followed by epimerization. In addition, we have identified an attractive mechanism for coupling the epimerization and desaturation in thermodynamically favorable steps with the aid of an external reductant. Other mechanisms that have been examined have significantly higher energy requirements or do not appear to be consistent with available experimental evidence.
Collapse
|
120
|
Gill PMW, Gordon MS, Head-Gordon M, Radom L. Remembrance: John A. Pople (1925–2004). J Chem Phys 2004. [DOI: 10.1063/1.1757682] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
121
|
|
122
|
Henry DJ, Coote ML, Gómez-Balderas R, Radom L. Comparison of the Kinetics and Thermodynamics for Methyl Radical Addition to CC, CO, and CS Double Bonds. J Am Chem Soc 2004; 126:1732-40. [PMID: 14871104 DOI: 10.1021/ja039139a] [Citation(s) in RCA: 55] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The barriers, enthalpies, and rate constants for the addition of methyl radical to the double bonds of a selection of alkene, carbonyl, and thiocarbonyl species (CH(2)=Z, CH(3)CH=Z, and (CH(3))(2)C=Z, where Z = CH(2), O, or S) and for the reverse beta-scission reactions have been investigated using high-level ab inito calculations. The results are rationalized with the aid of the curve-crossing model. The addition reactions proceed via early transition structures in all cases. The barriers for addition of methyl radical to C=C bonds are largely determined by the reaction exothermicities. Addition to the unsubstituted carbon center of C=C double bonds is favored over addition to the substituted carbon center, both kinetically (lower barriers) and thermodynamically (greater exothermicities). The barriers for addition to C=O bonds are influenced by both the reaction exothermicity and the singlet-triplet gap of the substrate. Addition to the carbon center is favored over addition to the oxygen, also both thermodynamically and kinetically. For the thiocarbonyl systems, addition to the carbon center is thermodynamically favored over addition to sulfur. However, in this case, the reaction is contrathermodynamic, addition to the sulfur center having a lower barrier due to spin density considerations. Entropic differences among corresponding addition and beta-scission reactions are relatively minor, and the differences in reaction rates are thus dominated by differences in the respective reaction barriers.
Collapse
|
123
|
Chan B, Radom L. Understanding Metal-Free Catalytic Hydrogenation: A Systematic Theoretical Study of the Hydrogenation of Ethene. Aust J Chem 2004. [DOI: 10.1071/ch04031] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
Abstract
Metal-free catalytic hydrogenation of ethene has been examined using high-level [G3(MP2)-RAD] ab initio molecular orbital theory. The dependence of the catalytic activity on the nature of the catalyst Z–X–H has been explored. We find that the catalytic activity is generally greater as Z–X–H becomes more acidic, both for first- and second-row atoms X. Molecules in which X is a second-row atom generally lead to more effective catalysis than the corresponding first-row analogues. The proton affinity at X of Z–X–H also contributes significantly to the catalysis in some cases (e.g. amines).
Collapse
|
124
|
Coote ML, Radom L. Substituent Effects in Xanthate-Mediated Polymerization of Vinyl Acetate: Ab Initio Evidence for an Alternative Fragmentation Pathway. Macromolecules 2003. [DOI: 10.1021/ma035477k] [Citation(s) in RCA: 92] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
125
|
Coote ML, Pross A, Radom L. Variable trends in R-X bond dissociation energies (R = Me, Et, i-Pr, t-Bu). Org Lett 2003. [PMID: 14627416 DOI: 10.1021/ol035860] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
[structure: see text] High level ab initio molecular orbital calculations confirm experimental indications that the effect of alkyl substituents (R = Me, Et, i-Pr, t-Bu) on R-X bond dissociation energies varies considerably according to the nature of X. A simple qualitative explanation in terms of valence-bond theory is presented, highlighting the increasing importance of the stabilization of R-X by the ionic R(+)X(-) configuration for electronegative X substituents (such as F, OH, and OCH(3)).
Collapse
|