Carnosine and Carcinine Derivatives Rapidly React with Hypochlorous Acid to Form Chloramines and Dichloramines.

Hypochlorous acid (HOCl) is a highly reactive, toxic species generated by neutrophils via the action of myeloperoxidase in order to destroy invading pathogens. However, when HOCl is produced inappropriately, it can damage host tissue and proteins and plays a role in the initiation and progression of disease. Carnosine, a peptide of ß-alanine and histidine, has been shown to react rapidly with HOCl yielding monochloramines and can undergo intramolecular transchlorination. The current study examines the kinetics and pH dependence of the reactions of carnosine and novel structural derivatives with HOCl and the occurrence of intra- and intermolecular transchlorination processes. We demonstrate that the transchlorination reactions of carnosine are pH dependent, with intramolecular transfer favored at higher pH. Carcinine, having a structure identical to carnosine though lacking the carboxylic acid group of the histidine residue, reacts with HOCl and forms monochloramines though intramolecular transfer reactions are not observed, and this is supported by computational modeling. Novel analogues with one (carnosine+1) and two (carnosine+2) methylene groups in the alkyl chain of the ß-alanine react with HOCl to yield monochloramines that undergo transchlorinations to yield a mixture of mono- and dichloramines. The latter are stable over 24 h. The ability of carnosine and derivatives to react rapidly with HOCl to give long-lived, poorly reactive, species may prevent damage to proteins and other targets at sites of inflammation.

An ONIOM investigation of the effect of conformation on bond dissociation energies in peptides.

In the present study, we use the ONIOM strategy of Morokuma and coworkers to examine the various CH bond dissociation energies (BDEs) of a small peptide (2ONW) and compare these with values obtained for its component individual amino acid residues. To evaluate suitable methods for ONIOM-based geometry optimizations, we test an "internal consistency" approach against full B3-LYP//B3-LYP results, and find B3-LYP/6-31G(d):AM1 to be appropriate. We find that the BDEs at the α-carbon in 2ONW are generally larger than the corresponding values for the individual residues on their own. This is attributed to the constraints of the peptide backbone leading to conformations that are not ideal for captodative stabilization of the resulting α-radicals. At the more flexible ß- and γ-positions, the differences between the BDEs in 2ONW and the individual residues are smaller. Overall, the α-BDEs are smaller than the ß- and γ-BDEs in most cases. Thus, to rationalize the inertness of peptide backbones with respect to α-hydrogen abstraction that is frequently found experimentally, it is necessary to consider alternative protection mechanisms such as the polar effect. © 2018 Wiley Periodicals, Inc.

Solvation of the Glycyl Radical.

The effect of adding explicit water molecules to the neutral (N) and zwitterionic (Z) forms of the glycyl radical has been examined. The results show that a minimum of three water molecules is required to stabilize the Z radical as a local minimum, with an energy gap of 123 kJ mol-1 between the N and Z forms at this point, in favor of the N form. Increasing the number of water molecules to ∼20 leads to a converged Z-N energy difference of ∼50 kJ mol-1 still in favor of the N form, even though the radical is not considered fully solvated from a structural point of view. Thus, energetic convergence is determined mainly by solvation of the polar functional groups, and a complete coverage of the entire molecule is not necessary. Because aqueous closed-shell glycine exists as a zwitterion while aqueous glycyl radical prefers the neutral form, the conversion between the two necessitates a change along the hydrogen-abstraction reaction pathway. In this regard, the transition structure for α-hydrogen abstraction by the ·OH radical has greater resemblance to glycine than to the glycyl radical. Overall, the barrier for hydrogen abstraction from Z glycine is larger than that from the N isomer, and this might act to provide some protection against radical damage to the free amino acid in the (aqueous) biological environment.

Computational Tale of Two Enzymes: Glycerol Dehydration With or Without B12.

We present a series of QM/MM calculations aimed at understanding the mechanism of the biological dehydration of glycerol. Strikingly and unusually, this process is catalyzed by two different radical enzymes, one of which is a coenzyme-B12-dependent enzyme and the other which is a coenzyme-B12-independent enzyme. We show that glycerol dehydration in the presence of the coenzyme-B12-dependent enzyme proceeds via a 1,2-OH shift, which benefits from a significant catalytic reduction in the barrier. In contrast, the same reaction in the presence of the coenzyme-B12-independent enzyme is unlikely to involve the 1,2-OH shift; instead, a strong preference for direct loss of water from a radical intermediate is indicated. We show that this preference, and ultimately the evolution of such enzymes, is strongly linked with the reactivities of the species responsible for abstracting a hydrogen atom from the substrate. It appears that the hydrogen-reabstraction step involving the product-related radical is fundamental to the mechanistic preference. The unconventional 1,2-OH shift seems to be required to generate a product-related radical of sufficient reactivity to cleave the relatively inactive C-H bond arising from the B12 cofactor. In the absence of B12, it is the relatively weak S-H bond of a cysteine residue that must be homolyzed. Such a transformation is much less demanding, and its inclusion apparently enables a simpler overall dehydration mechanism.

Effect of Hydrogen Bonding and Partial Deprotonation on the Oxidation of Peptides.

In a recent computational study, we found that hydrogen bonding/partial deprotonation facilitates subsequent electron transfer from amides to HO•. We have now analyzed these and related reactions with a glycine derivative as a model peptide, investigating not only reaction energies but also barriers for the individual steps. We find that partial deprotonation not only assists subsequent electron transfer (a sequential proton-loss electron-transfer (SPLET)-type reaction pathway) but also promotes sequential hydrogen-atom transfer (HAT, in a sequential proton-loss hydrogen-atom-transfer (SPLHAT)-type process), both being potential alternatives to direct HAT as routes for peptide oxidation. Each of these alternative pathways is calculated to have energy requirements that make them accessible and competitive. These oxidative processes may produce α-carbon-centered peptide radicals that, through deprotonation, are readily oxidized to the corresponding imines. We have also examined the possibility of competing reactions of amino acid side chains by comparing reactions of the glycine model with those of an analogous valine derivative. We find that, while the side chains of amino acids are more reactive toward direct HAT, a preceding partial deprotonation instead continues to favor the SPLET- and SPLHAT-type reactions, leading to the production of α-carbon-centered peptide radicals. Taken together, these processes have broad implications that impact many aspects of the science and utility of peptides.

Watson-Crick Base Pair Radical Cation as a Model for Oxidative Damage in DNA.

The deleterious cellular effects of ionizing radiation are well-known, but the mechanisms causing DNA damage are poorly understood. The accepted molecular events involve initial oxidation and deprotonation at guanine sites, triggering hydrogen atom abstraction reactions from the sugar moieties, causing DNA strand breaks. Probing the chemistry of the initially formed radical cation has been challenging. Here, we generate, spectroscopically characterize, and examine the reactivity of the Watson-Crick nucleobase pair radical cation in the gas phase. We observe rich chemistry, including proton transfer between the bases and propagation of the radical site in deoxyguanosine from the base to the sugar, thus rupturing the sugar. This first example of a gas-phase model system providing molecular-level details on the chemistry of an ionized DNA base pair paves the way toward a more complete understanding of molecular processes induced by radiation. It also highlights the role of radical propagation in chemistry, biology, and nanotechnology.

Impact of Hydrogen Bonding on the Susceptibility of Peptides to Oxidation.

The tendency of peptides to be oxidized is intimately connected with their function and even their ability to exist in an oxidative environment. Here we report high-level theoretical studies that show that hydrogen bonding can alter the susceptibility of peptides to oxidation, with complexation to a hydrogen-bond acceptor facilitating oxidation, and vice versa, impacting the feasibility of a diverse range of biological processes. It can even provide an energetically viable mechanistic alternative to direct hydrogen-atom abstraction. We find that hydrogen bonding to representative reactive groups leads to a broad (≈400 kJ mol-1 ) spectrum of ionization energies in the case of model amide, thiol and phenol systems. While some of the oxidative processes at the extreme ends of the spectrum are energetically prohibitive, subtle environmental and solvent effects could potentially mitigate the situation, leading to a balance between hydrogen bonding and oxidative susceptibility.

Hydrogen-adduction to open-shell graphene fragments: spectroscopy, thermochemistry and astrochemistry.

We apply a combination of state-of-the-art experimental and quantum-chemical methods to elucidate the electronic and chemical energetics of hydrogen adduction to a model open-shell graphene fragment. The lowest-energy adduct, 1H-phenalene, is determined to have a bond dissociation energy of 258.1 kJ mol-1, while other isomers exhibit reduced or in some cases negative bond dissociation energies, the metastable species being bound by the emergence of a conical intersection along the high-symmetry dissociation coordinate. The gas-phase excitation spectrum of 1H-phenalene and its radical cation are recorded using laser spectroscopy coupled to mass-spectrometry. Several electronically excited states of both species are observed, allowing the determination of the excited-state bond dissociation energy. The ionization energy of 1H-phenalene is determined to be 7.449(17) eV, consistent with high-level W1X-2 calculations.

Restricted-Open-Shell G4(MP2)-Type Procedures.

In the present study, we have reformulated the G4(MP2) and G4(MP2)-6X procedures for use with a restricted-open-shell (RO) formalism. We find that the resulting ROG4(MP2) and ROG4(MP2)-6X procedures generally perform comparably to the original unrestricted (U) variants, including their performance on radicals. Our analysis suggests that this is due mainly to the inclusion of empirical parameters that overcome the slightly less good performance of the U variants. However, a major practical advantage of ROG4(MP2) and ROG4(MP2)-6X is that they can be used in a wider range of computational chemistry software packages than the U analogs. We have demonstrated the importance of this aspect with a large-scale ROG4(MP2)-6X computation for the dissociation of the dodecahedryl dimer (C20H19)2.

Reactivity of disulfide bonds is markedly affected by structure and environment: implications for protein modification and stability.

Disulfide bonds play a key role in stabilizing protein structures, with disruption strongly associated with loss of protein function and activity. Previous data have suggested that disulfides show only modest reactivity with oxidants. In the current study, we report kinetic data indicating that selected disulfides react extremely rapidly, with a variation of 104 in rate constants. Five-membered ring disulfides are particularly reactive compared with acyclic (linear) disulfides or six-membered rings. Particular disulfides in proteins also show enhanced reactivity. This variation occurs with multiple oxidants and is shown to arise from favorable electrostatic stabilization of the incipient positive charge on the sulfur reaction center by remote groups, or by the neighboring sulfur for conformations in which the orbitals are suitably aligned. Controlling these factors should allow the design of efficient scavengers and high-stability proteins. These data are consistent with selective oxidative damage to particular disulfides, including those in some proteins.

Role of Hydrogen Bonding on the Reactivity of Thiyl Radicals: A Mass Spectrometric and Computational Study Using the Distonic Radical Ion Approach.

Experimental and computational quantum chemistry investigations of the gas-phase ion-molecule reactions between the distonic ions +H3N(CH2)nS• (n = 2-4) and the reagents dimethyl disulfide, allyl bromide, and allyl iodide demonstrate that intramolecular hydrogen bonding can modulate the reactivity of thiyl radicals. Thus, the 3-ammonium-1-propanethiyl radical (n = 3) exhibits the lowest reactivity of these distonic ions toward all substrates. Theoretical calculations on this distonic ion highlight that its most stable conformation involves a six-membered ring configuration, and that it has the strongest intramolecular hydrogen bond. In addition, the calculations indicate that the barrier heights for radical abstraction by this hydrogen-bond-stabilized 3-ammonium-1-propanethiyl radical are the highest among the systems examined, consistent with the experimental observations.

Frequency Scale Factors for Some Double-Hybrid Density Functional Theory Procedures: Accurate Thermochemical Components for High-Level Composite Protocols.

In the present study, we have obtained geometries and frequency scale factors for a number of double-hybrid density functional theory (DH-DFT) procedures. We have evaluated their performance for obtaining thermochemical quantities [zero-point vibrational energies (ZPVE) and thermal corrections for 298 K enthalpies (ΔH298) and 298 K entropies (S298)] to be used within high-level composite protocols (using the W2X procedure as a probe). We find that, in comparison with the previously prescribed protocol for optimization and frequency calculations (B3-LYP/cc-pVTZ+d), the use of contemporary DH-DFT methods such as DuT-D3 and DSD-type procedures leads to a slight overall improved performance compared with B3-LYP. A major strength of this approach, however, lies in the better robustness of the DH-DFT methods in that the largest deviations are notably smaller than those for B3-LYP. In general, the specific choices of the DH-DFT procedure and the associated basis set do not drastically change the performance. Nonetheless, we find that the DSD-PBE-P86/aug'-cc-pVTZ+d combination has a very slight edge over the others that we have examined, and we recommend its general use for geometry optimization and vibrational frequency calculations, in particular within high-level composite methods such as the higher-level members of the WnX series of protocols. The scale factors determined for DSD-PBE-P86/aug'-cc-pVTZ+d are 0.9830 (ZPVE), 0.9876 (ΔH298), and 0.9923 (S298).

Hydrogen-atom attack on phenol and toluene is ortho-directed.

The reaction of H + phenol and H/D + toluene has been studied in a supersonic expansion after electric discharge. The (1 + 1') resonance-enhanced multiphoton ionization (REMPI) spectra of the reaction products, at m/z = parent + 1, or parent + 2 amu, were measured by scanning the first (resonance) laser. The resulting spectra are highly structured. Ionization energies were measured by scanning the second (ionization) laser, while the first laser was tuned to a specific transition. Theoretical calculations, benchmarked to the well-studied H + benzene → cyclohexadienyl radical reaction, were performed. The spectrum arising from the reaction of H + phenol is attributed solely to the ortho-hydroxy-cyclohexadienyl radical, which was found in two conformers (syn and anti). Similarly, the reaction of H/D + toluene formed solely the ortho isomer. The preference for the ortho isomer at 100-200 K in the molecular beam is attributed to kinetic, not thermodynamic effects, caused by an entrance channel barrier that is ∼5 kJ mol(-1) lower for ortho than for other isomers. Based on these results, we predict that the reaction of H + phenol and H + toluene should still favour the ortho isomer under elevated temperature conditions in the early stages of combustion (200-400 °C).

α-Hydrogen Abstraction by •OH and •SH Radicals from Amino Acids and Their Peptide Derivatives.

We have used computational quantum chemistry to investigate the thermochemistry of α-hydrogen abstraction from the full set of amino acids normally found in proteins, as well as their peptide forms, by •OH and •SH radicals. These reactions, with their reasonable complexity in the electronic structure (at the α-carbon), are chosen as a consistent set of models for conducting a fairly robust assessment of theoretical procedures. Our benchmarking investigation shows that, in general, the performance for the various classes of theoretical methods improves in the order nonhybrid DFT → hybrid DFT → double-hybrid DFT → composite procedures. More specifically, we find that the DSD-PBE-P86 double-hybrid DFT procedure yields the best agreement with our high-level W1X-2 vibrationless barriers and reaction energies for this particular set of systems. A significant observation is that, when one considers relative instead of absolute values for the vibrationless barriers and reaction energies, even nonhybrid DFT procedures perform fairly well. To exploit this feature in a cost-effective manner, we have examined a number of multilayer schemes for the calculation of reaction energies and barriers for the abstraction reactions. We find that accurate values can be obtained when a "core" of glycine plus the abstracting radical is treated by DSD-PBE-P86, and the substituent effects are evaluated with M06-2X. Inspection of the set of calculated thermochemical data shows that the correlation between the free energy barriers and reaction free energies is strongest when the reactions are either endergonic or nearly thermoneutral.

Factors influencing the formation of polybromide monoanions in solutions of ionic liquid bromide salts.

Six different bromide salts - tetraethylammonium bromide ([N2,2,2,2]Br, Br), 1-ethyl-1-methylpiperidinium bromide ([C2MPip]Br, Br), 1-ethyl-1-methylpyrrolidinium bromide ([C2MPyrr]Br, Br), 1-ethyl-3-methylimidazolium bromide ([C2MIm]Br, Br), 1-ethylpyridinium bromide ([C2Py]Br, Br), and 1-(2-hydroxyethyl)pyridinium bromide ([C2OHPy]Br, Br) - were studied in regards to their capacity to form polybromide monoanion products on addition of molecular bromine in acetonitrile solutions. Using complementary spectroscopic and computational methods for the examination of tribromide and pentabromide anion formation, key factors influencing polybromide sequestration were identified. Here, we present criteria for the targeted synthesis of highly efficient bromine sequestration agents.

Beyond the Halogen Bond: Examining the Limits of Extended Polybromide Networks through Quantum-Chemical Investigations.

The bonding environments of some polybromide monoanions and networks were examined by quantum-chemical methods to investigate electronic interactions between dibromine-dibromine contacts. Examination of thermodynamic parameters and a bond critical point analysis give strong evidence for such bonding modes, which have been previously treated disparately in the literature. The thermodynamic stability of large polybromides up to [Br37 ](-) was also predicted by these methods.

On the inclusion of post-MP2 contributions to double-hybrid density functionals.

In this study, we explore the effect of supplementing the DuT spin-component-scaled double-hybrid density functional method with post-second-order Møller-Plesset-type theory (MP2) correlation terms. We find that the inclusion of additional MP3 correlation energies has almost no effect on the performance. Further addition of correlation effects from MP4 generally leads to a small improvement in the performance. However, we find that the inclusion of the higher-order perturbative correlation effects does not rectify some major shortcomings of DuT for more challenging systems, and the use of MP4, in fact, leads to a significant deterioration in the performance in some cases. We also find that the use of correlation energies from CCSD(T) instead of those from MP3 and MP4 does not lead to a substantial improvement over the MP4-based method, both in general and in some difficult cases that we have examined. An additional observation is that, for large systems that are dominated by noncovalent interactions, DuT and the two MPn-based post-MP2 double-hybrid density functional theory (DFT) procedures all benefit from the inclusion of dispersion corrections. Overall, our investigation suggests that the current generation of MP2-based double-hybrid DFT methods may already be providing close to the optimal performance that can be achieved with the double-hybrid methodology paired with spin-component-scaling. Development of even better double hybrids is an active research field, and we hope that our study provides valuable insights. We recommend the continuing use of existing MP2-based double-hybrid methods as a bridging level between hybrid density functional procedures and high-level wave-function-based procedures.

Preparation of an ion with the highest calculated proton affinity: ortho-diethynylbenzene dianion.

Owing to the increased proton affinity that results from additional negative charges, multiply-charged anions have been proposed as one route to prepare and access a range of new and powerful "superbases". Paradoxically, while the additional electrons in polyanions increase basicity they serve to diminish the electron binding energy and thus, it had been thought, hinder experimental synthesis. We report the synthesis and isolation of the ortho-diethynylbenzene dianion (ortho-DEB2-) and present observations of this novel species undergoing gas-phase proton-abstraction reactions. Using a theoretical model based on Marcus-Hush theory, we attribute the stability of ortho-DEB2- to the presence of a barrier that prevents spontaneous electron detachment. The proton affinity of 1843 kJ mol-1 calculated for this dianion superbase using high-level quantum chemistry calculations significantly exceeds that of the lithium monoxide anion, the most basic system previously prepared. The ortho-diethynylbenzene dianion is therefore the strongest base that has been experimentally observed to date.

W2X and W3X-L: Cost-Effective Approximations to W2 and W4 with kJ mol(-1) Accuracy.

We have formulated the W2X and W3X-L protocols as cost-effective alternatives to W2 and W3/W4, respectively, and to supplement our previously developed set of W1X and W3X procedures. The W2X procedure provides an accurate approximation to the all-electron scalar-relativistic CCSD(T)/CBS energy, with a mean absolute deviation (MAD) of 0.6 kJ mol(-1) from benchmark energies provided by the CCSD(T) component in the W4 protocol. Such a performance is comparable to that of W2w (0.5 kJ mol(-1)) but comes at a significantly lower cost. Comparison of computational requirements shows that W2X should be applicable to systems that can be treated by the W1w method. Thus, W2X provides an accurate means for the treatment of medium-sized systems such as naphthalene. For the calculation of post-CCSD(T) effects, we propose a slight modification to the method used in our previously devised W3X procedure. Our new W3-type protocol (W3X-L) combines this new post-CCSD(T) treatment with our new W2X procedure. It has an MAD from benchmark values of 0.8 kJ mol(-1) for the W4-11 set, which is comparable to that for the computationally more demanding W3.2 method (0.6 kJ mol(-1)). However, the use of the even relatively modest post-CCSD(T) calculations in W3X-L still represents a computational bottleneck, and this currently restricts its application to systems up to the size of benzene with our current computing resources.