MPW1K Performs Much Better than B3LYP in DFT Calculations on Reactions that Proceed by Proton-Coupled Electron Transfer (PCET)

DFT calculations have been performed with the B3LYP and MPW1K functional on the hydrogen atom abstraction reactions of ethenoxyl with ethenol and of phenoxyl with both phenol and alpha-naphthol. Comparison with the results of G3 calculations shows that B3LYP seriously underestimates the barrier heights for the reaction of ethenoxyl with ethenol by both proton-coupled electron transfer (PCET) and hydrogen atom transfer (HAT) mechanisms. The MPW1K functional also underestimates the barrier heights, but by much less than B3LYP. Similarly, comparison with the results of experiments on the reaction of phenoxyl radical with alpha-naphthol indicates that the barrier height for the preferred PCET mechanism is calculated more accurately by MPW1K than by B3LYP. These findings indicate that the MPW1K functional is much better suited than B3LYP for calculations on hydrogen abstraction reactions by both HAT and PCET mechanisms.     

Are 1,5-disubstituted semibullvalenes that have C2v equilibrium geometries necessarily bishomoaromatic?

Unrestricted density functional theory (UB3LYP), CASSCF, and CASPT2 calculations have been employed to compute the relative energies of the C(s) and C(2v) geometries of several 1,5-disubstituted semibullvalenes. Substitution at these positions with R = F, -CH(2)-, or -O- affords semibullvalenes that are predicted to have C(2v) equilibrium geometries. Calculated singlet-triplet energy splittings and the energies of isodesmic reactions are used to assess the amount of bishomoaromatic character at these geometries. The results of these calculations show that employing strain to destabilize the C(s) geometries of semibullvalenes can lead to a significant decrease in the amount of bishomoaromatic stabilization of the C(2v) geometries, due to reduced through-space interaction between the two allyl groups. However, the C(2v) equilibrium geometries of the 1,5-disubstituted semibullvalenes with R = F and -RR- = -O- do benefit from stabilizing through-bond interactions between the two allyl groups. These interactions involve mixing of the bisallyl HOMO with the low-lying C-F or C-O sigma orbital combinations of the same symmetry. In contrast, for -RR- = -CH(2)-, through-bond interactions destabilize the bisallyl HOMO and are predicted to make the ground state of this semibullvalene a triplet.     

Ligand-assisted reduction of osmium tetroxide with molecular hydrogen via a [3+2] mechanism

Osmium tetroxide is reduced by molecular hydrogen in the presence of ligands in both polar and nonpolar solvents. In CHCl3 containing pyridine (py) or 1,10-phenanthroline (phen), OsO4 is reduced by H2 to the known Os(VI) dimers L2Os(O)2(mu-O)2Os(O)2L2 (L2 = py2, phen). However, in the absence of ligands in CHCl3 and other nonpolar solvents, OsO4 is unreactive toward H2 over a week at ambient temperatures. In basic aqueous media, H2 reduces OsO4(OH)n(n-) (n = 0, 1, 2) to the isolable Os(VI) complex, OsO2(OH)4(2-), at rates close to that found in py/CHCl3. Depending on the pH, the aqueous reactions are exergonic by deltaG = -20 to -27 kcal mol(-1), based on electrochemical data. The second-order rate constants for the aqueous reactions are larger as the number of coordinated hydroxide ligands increases, k(OsO4) = 1.6(2) x 10(-2) M(-1) s(-1) < k(OsO4(OH)-) = 3.8(4) x 10(-2) M(-1) s(-1) < k(OsO4(OH)2(2-)) = 3.8(4) x 10(-1) M(-1) s(-1). The observation of primary deuterium kinetic isotope effects, k(H2)/k(D2) = 3.1(3) for OsO4 and 3.6(4) for OsO4(OH)-, indicates that the rate-determining step in each case involves H-H bond cleavage. Density functional calculations and thermochemical arguments favor a concerted [3+2] addition of H2 across two oxo groups of OsO4(L)n and argue against H* or H- abstraction from H2 or [2+2] addition of H2 across one Os=O bond. The [3+2] mechanism is analogous to that of alkene addition to OsO4(L)n to form diolates, for which acceleration by added ligands has been extensively documented. The observation that ligands also accelerate H2 addition to OsO4(L)n highlights the analogy between these two reactions.     

Structure and torsional flexibility of the linkage between guanine and fluorene residues in the deoxyguanosine–aminofluorene and deoxyguanosine–acetylaminofluorene carcinogenic adducts

Potential energy surfaces for rotations around two central CN bonds in N-(deoxyguanosin-8-yl)-2-acetylaminofluorene (AAF–dG) and its deacetylated derivative (AF–dG) were studied using Amber 95 molecular mechanics. Both of these adducts are known to be strong mutagens and carcinogens. New Amber 95 force field parameters were derived for the linkage connecting guanine and fluorene moieties in AAF–dG and AF–dG. For this purpose, we determined ab initio MP2/cc-pVDZ//B3-LYP/6-31G* and polarized continuum model Hartree–Fock/6-31G* potential energy surfaces of smaller model systems that included the N-methylimidazole–acetylaniline and N-methylimidazole–aniline adducts. The molecular mechanics parameters were adjusted to minimize differences between the gas-phase ab initio and molecular mechanics surfaces of these model systems. The resulting parameters were transferred to AF–dG and AAF–dG. The barrier for the rotation of the fluorene residue in AF–dG was found to be less than 2 kcal/mol. Such a small barrier renders the fluorene moiety freely rotatable at room temperature. In contrast, the fluorene rotation in AAF–dG is hindered by a significantly larger barrier of 10 kcal/mol. This barrier corresponds to conformations in which the fluorene and acetyl groups lie in the same plane, and is largely due to steric repulsion. Similarly, the coplanar arrangement of guanine and the bridging amino or acetyl groups is disfavored by 5–10 kcal/mol, with AAF–dG again being the more rigid of the two molecules. Energy minima for a rotation around a bond between guanine and the bridging nitrogen are found at ±80° in AAF–dG, and at 120° and –90° for AF–dG. Overall, the fluorene–dG linkages in AF–dG and AAF–dG adducts have significantly different equilibrium structures and torsional flexibilities. These differences may be contributing factors for the observed disparity in mutagenic effects of these adducts.Electronic Supplementary Material: Supplementary material is available in the online version of this article at Acknowledgements. This work was supported by the NSF REU grant no. CHE-0243825 to Loyola University Chicago. We thank to Tom Ellenberger and Shuchismita Dutta for providing us with their results prior to publication.     

Surface phase behavior and microstructure of lipid/PEG-emulsifier monolayer-coated microbubbles

Langmuir trough methods and fluorescence microscopy were combined to investigate the phase behavior and microstructure of monolayer shells coating micron-scale bubbles (microbubbles) typically used in biomedical applications. The monolayer shell consisted of a homologous series of saturated acyl chain phospholipids and an emulsifier containing a single hydrophobic stearate chain and polyethylene glycol (PEG) head group. PEG-emulsifier was fully miscible with expanded phase lipids and phase separated from condensed phase lipids. Phase coexistence was observed in the form of dark condensed phase lipid domains surrounded by a sea of bright, emulsifier-rich expanded phase. A rich assortment of condensed phase area fractions and domain morphologies, including networks and other novel structures, were observed in each batch of microbubbles. Network domains were reproduced in Langmuir monolayers under conditions of heating–cooling followed by compression–expansion, as well as in microbubble shells that underwent surface flow with slight compression. Domain size decreased with increased cooling rate through the phase transition temperature, and domain branching increased with lipid acyl chain length at high cooling rates. Squeeze-out of the emulsifier at a surface pressure near 35 mN/m was indicated by a plateau in Langmuir isotherms and directly visualized with fluorescence microscopy, although collapse of the solid lipid domains occurred at much higher surface pressures. Compression of the monolayer past the PEG-emulsifier squeeze-out surface pressure resulted in a dark shell composed entirely of lipid. Under certain conditions, the PEG-emulsifier was reincorporated upon subsequent expansion. Factors that affect shell formation and evolution, as well as implications for the rational design of microbubbles in medical applications, are discussed.     

Do deviations from bond enthalpy additivity define the thermodynamic stabilities of diradicals?

Deviations from bond enthalpy additivity (DeltaBEA) are frequently used to assess the thermodyamic stabilities of diradicals. (U)B3LYP/6-31G calculations have been performed in order to determine how well DeltaBEA values actually do reflect the thermodynamic stabilities of the triplet states of diradicals in which one or both nonbonding electrons occupy a delocalized pi orbital. The calculations find that different pathways for forming sigma,pi-diradicals, such as alpha,2- and alpha,4-dehydrotoluene (4 and 6), give DeltaBEA values that differ by ca. 1 kcal/mol. The path dependency of the DeltaBEA values is computed to be one order of magnitude larger for non-Kekulé hydrocarbon diradicals, such as m-benzoquinodimethane (12) and 1,3-dimethylenecyclobutane-2,4-diyl (15), than for sigma,pi-diradicals. Since the DeltaBEA values for forming 4, 6, 12, and 15 are all path dependent, we conclude that DeltaBEA values for diradicals with one or two delocalized, nonbonding pi electrons do not, in general, uniquely define the thermodynamic stabilities of the diradicals. Hence, DeltaBEA values should not be used for this purpose, especially for non-Kekulé hydrocarbon diradicals.