Comparison of the Accurate Kohn−Sham Solution with the Generalized Gradient Approximations (GGAs) for the SN2 Reaction F- + CH3F → FCH3 + F-: A Qualitative Rule To Predict Success or Failure of GGAs

Structure and fragmentation chemistry of the peptide radical cations of glycylphenylalanylglycine (GFG)

Herein, we explore the generation and characterization of the radical cations of glycylphenylalanylglycine, or [GFG]•+, formed via dissociative electron-transfer reaction from the tripeptide to copper(II) within a ternary complex. A comprehensive investigation employing isotopic labeling, infrared multiple-photon dissociation (IRMPD) spectroscopy, and density functional theory (DFT) calculations elucidated the details and energetics in formation of the peptide radical cations as well as their dissociation products. Unlike conventional aromatic-containing peptide radical cations that primarily form canonical π-radicals, our findings reveal that 75% of the population of the experimentally produced [GFG]•+ precursors are [GFα•G]+, where the radical resides on the middle α-carbon of the phenylalanyl residue. This unexpected isomeric ion has an enthalpy of 6.8 kcal/mol above the global minimum, which has an N-terminal captodative structure, [Gα•FG]+, comprising 25% of the population. The [b₂-H]•+ product ions are also present in a ratio of 75/25 from [GFα•G]+/ [Gα•FG]+, the results of which are obtained from matches between the IRMPD action spectrum and predicted IR absorption spectra of the [b₂-H]•+ candidate structures, as well as from IRMPD isomer population analyses.

Phosphate Migration versus the Loss of Phosphoric Acid in Protonated Phosphopeptides: A Computational Study

Residue-specific phosphorylation is a protein post-translational modification that regulates cellular functions. Experimental determination of the exact sites of protein phosphorylation provides an understanding of the signaling and processes at work for a given cellular state. Any experimental artifact that involves migration of the phosphate group during measurement is a concern, as the outcome can lead to erroneous conclusions that may confound studies on cellular signal transduction. Herein, we examine computationally the mechanism by which a phosphate group migrates from one serine residue to another serine in monoprotonated pentapeptides [BA-pSer-Gly-Ser-BB + H]+ → [BA-Ser-Gly-pSer-BB + H]+ (where BA and BB are different combinations of the three basic amino acids, histidine, lysine, and arginine). In addition to moving the phosphate group, the overall mechanism involves transferring a proton from the N-terminal amino acid, BA, to the C-terminal amino acid, BB. This is not a synchronous process, and there is a key high-energy intermediate, structure C, that is zwitterionic with both the basic amino acids protonated and the phosphate group attached to both serine residues and carrying a negative charge. The barriers to moving the phosphate group are calculated to be in the range of 219-274 kJ mol-1 at the B3LYP/6-31G(d) level. These barriers are systematically slightly lower and in good agreement with single-point energy calculations at both M06-2X/6-311++G(d,p) and MP2/6-31++G(d,p) levels. The competitive reaction, loss of phosphoric acid from the protonated pentapeptides, has a barrier in the range of 176-202 kJ mol-1 at the B3LYP/6-31G(d) level. Extension of the theory to M06-2X/6-311++G(d,p)//B3LYP/6-31G(d) and MP2/6-31++G(d,p)// B3LYP/6-31G(d) gives higher values for the loss of phosphoric acid, falling in the range of 196-226 kJ mol-1; these are comparable to the barriers against phosphate migration at the same levels of theory. For larger peptides His-pSer-(Gly)n-Ser-His, where n has values from 2 to 5, the barriers against the loss of phosphoric acid are higher than those against the phosphate group migration. This difference is most pronounced and significant when n = 4 and 5 (the differences are approximately 80 kJ mol-1 under the single-point energy calculations at the M06-2X and MP2 levels). Energy differences using two more recent functionals, M08-HX and MN15, on His-pSer-(Gly)n-Ser-His, where n = 1 and 5, are in good agreement with the M06-2X and MP2 calculations. These results provide the mechanistic rationale for phosphate migration versus other competing reactions in the gas phase under tandem mass spectrometry conditions.

Protonated α-N-Acetyl Galactose Glycopeptide Dissociation Chemistry

We recently provided mass spectrometric, H/D labeling, and computational evidence of pyranose to furanose N-acetylated ion isomerization reactions that occurred prior to glycosidic bond cleavage in both O- and N-linked glycosylated amino acid model systems (Guan et al. Phys. Chem. Chem. Phys., 2021, 23, 23256-23266). These reactions occurred irrespective of the glycosidic linkage stereochemistry (α or β) and the N-acetylated hexose structure (GlcNAc or GalNAc). In the present article, we test the generality of the preceding findings by examining threonyl α-GalNAc-glycosylated peptides. We utilize computational chemistry to compare the various dissociation and isomerization pathways accessible with collisional activation. We then interrogate the structure(s) of the resulting charged glycan and peptide fragments with infrared "action" spectroscopy. Isomerization of the original pyranose, the protonated glycopeptide [AT(GalNAc)A+H]⁺, is predicted to be facile compared to direct dissociation, as is the glycosidic bond cleavage of the newly formed furanose form, i.e., furanose oxazolinium ion structures are predicted to predominate. IR action spectra for the m/z 204, C₈H₁₄N₁O₅⁺, glycan fragment population support this prediction. The IR action spectra of the complementary m/z 262 peptide fragment were assigned as a mixture of the lowest-energy structures of [ATA+H]⁺ consistent with the literature. If general, the change to a furanose m/z 204 product ion structure fundamentally alters the ion population available for MS³ dissociation and glycopeptide sequence identification.

Spin-component-scaled and dispersion-corrected second-order Møller-Plesset perturbation theory: A path toward chemical accuracy

Second-order Moller-Plesset perturbation theory (MP2) provides a valuable alternative to density functional theory for modeing problems in organic and biological chemistry. However, MP2 suffers from known limitations in the description...

Generation, Characterization, and Dissociation of Radical Cations Derived from Prolyl-glycyl-glycine

Radical cations of an aliphatic tripeptide prolyl-glycyl-glycine (PGG^•+) and its sequence ions [a₃ + H]^•+ and [b₂ - H]^•+ have been generated by collision-induced dissociation of the [Cu(Phen)(PGG)]^•2+ complex, where Phen = 1,10-phenanthroline. Infrared multiple photon dissociation spectroscopy, ion-molecule reaction experiments, and theoretical calculations have been used to investigate the structures of these ions. The unpaired electron in these three radical cations is located at different α-carbons. The PGG^•+ radical cation has a captodative structure with the radical at the α-carbon of the proline residue and the proton on the oxygen of the first amide group. This structure is at the global minimum on the potential energy surface (PES). By contrast, the [a₃ + H]^•+ and [b₂ - H]^•+ ions are not the lowest-energy structures on their respective PESs, and their radicals are formally located at the C-terminal and second α-carbons, respectively. Density functional theory calculations on the structures of the ternary copper(II) complex ion suggest that the charge-solvated isomer of the metal complex is the precursor ion that dissociates to give the PGG^•+ radical cation. The isomer of the complex in which PGG is bound as a zwitterion dissociates to give the [a₃ + H]^•+ and [b₂ - H]^•+ ions.

Size Dependent Fragmentation Chemistry of Short Doubly Protonated Tryptic Peptides

Tandem mass spectrometry of electrospray ionized multiply charged peptide ions is commonly used to identify the sequence of peptide(s) and infer the identity of source protein(s). Doubly protonated peptide ions are consistently the most efficiently sequenced ions following collision-induced dissociation of peptides generated by tryptic digestion. While the broad characteristics of longer (N ≥ 8 residue) doubly protonated peptides have been investigated, there is comparatively little data on shorter systems where charge repulsion should exhibit the greatest influence on the dissociation chemistry. To address this gap and further understand the chemistry underlying collisional-dissociation of doubly charged tryptic peptides, two series of analytes ([G_xR+2H]²⁺ and [A_xR+2H]²⁺, x = 2-5) were investigated experimentally and with theory. We find distinct differences in the preference of bond cleavage sites for these peptides as a function of size and to a lesser extent composition. Density functional calculations at two levels of theory predict that the threshold relative energies required for bond cleavages at the same site for peptides of different size are quite similar (for example, b₂-y_N-2). In isolation, this finding is inconsistent with experiment. However, the predicted extent of entropy change of these reactions is size dependent. Subsequent RRKM rate constant calculations provide a far clearer picture of the kinetics of the competing bond cleavage reactions enabling rationalization of experimental findings. The M06-2X data were substantially more consistent with experiment than were the B3LYP data.

Are Anions of Cyclobutane Beryllium Derivatives Stabilized through Four-Center One-Electron Bonds?

High-level G4 ab initio calculations allowed us to show that C₄H₄(BeX)₄ (X = H, Cl) derivatives behave as rather efficient electron capturers due to their ability to trap the extra electron through the formation of a four-membered beryllium ring. This finding is in agreement with previous work showing the ability of highly electron-deficient atoms, such as beryllium, to lead to multicenter one-electron bonds. In our particular case, the formation of the four-center bond is characterized, in very good harmony, by different topological methods such as quantum theory of atoms in molecules (QTAIM), the electron localization function (ELF), and the noncovalent interactions (NCI) approach and is accompanied by large electron affinity values, around 300 kJ·mol^-1, in the gas phase. Preliminary results may anticipate that the ability of groups of beryllium atoms to trap electrons decays on going to bigger systems.

Bonding between electron-deficient atoms: strong Lewis-acid character preserved in X–Y–X (X = B, Al; Y = Be, Mg) bridges

Beryllium bis(diazaborolyl) derivatives and their Mg and Al-containing analogues are stable compounds stabilized through covalent bonds between electron-deficient atoms, and behave as good Lewis acids.

Degradation and mechanism of microcystin-LR by PbCrO4 nanorods driven by visible light

This work focuses on the photocatalytic removal of recalcitrant organic pollutants in water treatment. Based on facile precipitation reaction, we fabricated a photocatalyst (PbCrO₄) in single crystals that present evident response to visible light and employed the catalyst in the photocatalytic decomposition of microcystin-LR (MC-LR). In the degradation test using the nanorods with prepared PbCrO₄ photocatalyst, a 100% removal efficiency (27 min reaction) and a kinetics constant of 0.1356 min^-1 were achieved. Such a high performance of PbCrO₄ in photocatalytic conversion of MC-LR was ascribed to its high carrier separation efficiency, positive valence band (VB) position, and good delocalization of VB and conduction band (CB). The test of electron spin-resonance resonance (ESR) demonstrated that excessive free OH radicals were produced during the PbCrO₄ photocatalysis of MC-LR. The density functional theory (DFT) and LC/MS/MS technology were employed to ascertain the intermediates during the MC-LR photocatalytic degradation. The major intermediates were resulted from the attack of hydroxyl radicals to the ADDA side chains of MC-LR structure. This study provides a proof-of-concept strategy to develop effective photocatalysts to efficiently produce OH radicals for the visible-light induced photocatalytic degradation of MC-LR in water.

Improved agreement between experimental and computational results for collision-induced dissociation mass spectrometry of cation-tagged hexoses

Low exact exchange DFT methods underestimate the dehydration barriers as the charge localization in such transition states is challenging to describe.

Sequence Ion Structures and Dissociation Chemistry of Deprotonated Sucrose Anions

We investigate the tandem mass spectrometry of regiospecifically labeled, deprotonated sucrose analytes. We utilize density functional theory calculations to model the pertinent gas-phase fragmentation chemistry of the prevalent glycosidic bond cleavages (B₁-Y₁ and C₁-Z₁ reactions) and compare these predictions to infrared spectroscopy experiments on the resulting B₁ and C₁ product anions. For the C₁ anions, barriers to interconversion of the pyranose [α-glucose-H]^-, C₁ anions to entropically favorable ring-open aldehyde-terminated forms were modest (41 kJ mol^-1) consistent with the observation of a band assigned to a carbonyl stretch at ~ 1680-1720 cm^-1. For the B₁ anions, our transition structure calculations predict the presence of both deprotonated 1,6-anhydroglucose and carbon 2-ketone ((4S,5S,6R)-4,5-dihydroxy-6-(hydroxymethyl)dihydro-2H-pyran-3(4H)-one) anion structures, with the latter predominating. This hypothesis is supported by our spectroscopic data which show diagnostic bands at 1600, 1674, and 1699 cm^-1 (deprotonated carbon 2-ketone structures), and at ~ 1541 cm^-1 (both types of structure) and RRKM rate calculations. The deprotonated carbon 2-ketone structures are also the lowest energy product B₁ anions. Graphical Abstract ᅟ.

Multiscale Simulation of the Interaction and Adsorption of Ions on a Hydrophobic Graphene Surface

The adsorption of ions on a graphene surface is very important to control relevant graphene-based processes. In this work, a multiscale simulation was carried out to study the adsorption of Na⁺ /Cl^- ions on graphene by combining quantum mechanics calculations and molecular dynamics (MD) simulations. The interaction energies of the ions with graphene were computed using density functional theory (DFT). It was found that the ions show strong interaction with a graphene cluster and the overwhelming portion of the interaction energy is the ion-π orbital interaction. The large orbital interaction can be ascribed to the two contributions arising from the ion-induced polarization of graphene and the charge transfer between ion and graphene. Their different contribution degrees reveal that the polarization effect plays a main role on the orbital interaction for ion adsorption. Comparatively, for Na/Cl atom adsorption, the charge transfer shows large part to the orbital interaction with weak atom-induced polarization. The obtained interaction energies were applied to develop new interaction potentials between ion and graphene, and then MD simulations were used to study the interfacial adsorption behavior of Na⁺ /Cl^- aqueous solution onto the graphene surface. Due to enhanced ion-π interactions, Na⁺ /Cl^- cooperatively demonstrates a strong ion adsorption layer through direct contact with the hydrophobic graphene surface. Our simulation result presents a new understanding of ion-graphene interactions.

Solvent Effect on the Potential Energy Surfaces of the F- + CH3CH2Br Reaction

Although substantial work has been undertaken on reaction pathways involved in base-promoted elimination reactions and bimolecular nucleophilic substitution reaction of F^- on CH₃CH₂X (X = Cl, Br, I), the effect of solvents with varying dielectric constants on the stereochemistry of each of the reaction species involved across the reaction profile have not yet been clearly understood. The present investigation reports the effect of solvents on the potential energy surfaces (PES) and structures of the species appearing in the reaction pathway of F^- with bromoethane. The PESs in the gas phase have been computed at MP2 level and CCSD(T) level. The performance of several hybrid density functional, such as B3LYP, M06, M06L, BHandH, X3LYP, M05, M05-2X, and M06-2X have also been investigated toward describing the elimination and nucleophilic substitution reactions. With respect to MAE values and to make the computation cost-effective, we have explored the implicit continuum solvent model, CPCM in solvents like cyclohexane, methanol, acetonitrile, dimethyl sulfoxide and water. The reactant complexes proceed through the subsequent steps to produce fluoroethane as the substitution product and ethylene as one of the elimination products. For elimination reaction both syn and anti elimination have been explored. The calculated relatives energies values, which are negative in the gas phase, are found to be positive in polar solvents since the point charge in the separated reactants are more stabilized than the dispersed charge in the transient complex, which has also been analyzed through NBO analysis.

Electronic structure of BN-aromatics: Choice of reliable computational tools

The importance of having reliable calculation tools to interpret and predict the electronic properties of BN-aromatics is directly linked to the growing interest for these very promising new systems in the field of materials science, biomedical research, or energy sustainability. Ionization energy (IE) is one of the most important parameters to approach the electronic structure of molecules. It can be theoretically estimated, but in order to evaluate their persistence and propose the most reliable tools for the evaluation of different electronic properties of existent or only imagined BN-containing compounds, we took as reference experimental values of ionization energies provided by ultra-violet photoelectron spectroscopy (UV-PES) in gas phase-the only technique giving access to the energy levels of filled molecular orbitals. Thus, a set of 21 aromatic molecules containing B-N bonds and B-N-B patterns has been merged for a comparison between experimental IEs obtained by UV-PES and various theoretical approaches for their estimation. Time-Dependent Density Functional Theory (TD-DFT) methods using B3LYP and long-range corrected CAM-B3LYP functionals are used, combined with the ΔSCF approach, and compared with electron propagator theory such as outer valence Green's function (OVGF, P3) and symmetry adapted cluster-configuration interaction ab initio methods. Direct Kohn-Sham estimation and "corrected" Kohn-Sham estimation are also given. The deviation between experimental and theoretical values is computed for each molecule, and a statistical study is performed over the average and the root mean square for the whole set and sub-sets of molecules. It is shown that (i) ΔSCF+TDDFT(CAM-B3LYP), OVGF, and P3 are the most efficient way for a good agreement with UV-PES values, (ii) a CAM-B3LYP range-separated hybrid functional is significantly better than B3LYP for the purpose, especially for extended conjugated systems, and (iii) the "corrected" Kohn-Sham result is a fast and simple way to predict IEs.

Constrained-Orbital Density Functional Theory. Computational Method and Applications to Surface Chemical Processes

We present a method for performing density-functional theory (DFT) calculations in which one or more Kohn-Sham orbitals are constrained to be localized on individual atoms. This constrained-orbital DFT (CO-DFT) approach can be used to tackle two prevalent shortcomings of DFT: the lack of transparency with regard to the governing electronic structure in large (planewave based) DFT calculations and the limitations of semilocal DFT in describing systems with localized electrons or a large degree of static correlation. CO-DFT helps to address the first of these issues by decomposing complex orbital transformations occurring during elementary chemical processes into simpler and more intuitive transformations. The second issue is addressed by using the CO-DFT method to generate configuration states for multiconfiguration Kohn-Sham calculations. We demonstrate both of these applications for elementary reaction steps involved in the oxygen evolution reaction.

Isomerization versus dissociation of phenylalanylglycyltryptophan radical cations

Four isomers of the radical cation of tripeptide phenylalanylglycyltryptophan, in which the initial location of the radical center is well defined, have been isolated and their collision-induced dissociation (CID) spectra examined. These ions, the π-centered [FGW_π˙]⁺, α-carbon- [FG_α˙W]⁺, N-centered [FGW_N˙]⁺ and ζ-carbon- [F_ζ˙GW]⁺ radical cations, were generated via collision-induced dissociation (CID) of transition metal-ligand-peptide complexes, side chain fragmentation of a π-centered radical cation, homolytic cleavage of a labile nitrogen-nitrogen single bond, and laser induced dissociation of an iodinated peptide, respectively. The π-centered and tryptophan N-centered peptide radical cations produced almost identical CID spectra, despite the different locations of their initial radical sites, which indicated that interconversion between the π-centered and tryptophan N-centered radical cations is facile. By contrast, the α-carbon-glycyl radical [FG_α˙W]⁺, and ζ-phenyl radical [F_ζ˙GW]⁺, featured different dissociation product ions, suggesting that the interconversions among α-carbon, π-centered (or tryptophan N-centered) and ζ-carbon-radical cations have higher barriers than those to dissociation. Density functional theory calculations have been used to perform systematic mechanistic investigations on the interconversions between these isomers and to study selected fragmentation pathways for these isomeric peptide radical cations. The results showed that the energy barrier for interconversion between [FGW_π˙]⁺ and [FGW_N˙]⁺ is only 31.1 kcal mol^-1, much lower than the barriers to their dissociation (40.3 kcal mol^-1). For the [FGW_π˙]⁺, [FG_α˙W]⁺, and [F_ζ˙GW]⁺, the barriers to interconversion are higher than those to dissociation, suggesting that interconversions among these isomers are not competitive with dissociations. The [z₃ - H]˙⁺ ions isolated from [FG_α˙W]⁺ and [F_ζ˙GW]⁺ show distinctly different fragmentation patterns, indicating that the structures of these ions are different and this result is supported by the DFT calculations.

Cationized Carbohydrate Gas-Phase Fragmentation Chemistry

We investigate the fragmentation chemistry of cationized carbohydrates using a combination of tandem mass spectrometry, regioselective labeling, and computational methods. Our model system is D-lactose. Barriers to the fundamental glyosidic bond cleavage reactions, neutral loss pathways, and structurally informative cross-ring cleavages are investigated. The most energetically favorable conformations of cationized D-lactose were found to be similar. In agreement with the literature, larger group I cations result in structures with increased cation coordination number which require greater collision energy to dissociate. In contrast with earlier proposals, the B _n -Y _m fragmentation pathways of both protonated and sodium-cationized analytes proceed via protonation of the glycosidic oxygen with concerted glycosidic bond cleavage. Additionally, for the sodiated congeners our calculations support sodiated 1,6-anhydrogalactose B _n ion structures, unlike the preceding literature. This affects the subsequent propensity of formation and prediction of B _n /Y _m branching ratio. The nature of the anomeric center (α/β) affects the relative energies of these processes, but not the overall ranking. Low-energy cross-ring cleavages are observed for the metal-cationized analytes with a retro-aldol mechanism producing the ^0,2 A ₂ ion from the sodiated forms. Theory and experiment support the importance of consecutive fragmentation processes, particularly for the protonated congeners at higher collision energies. Graphical Abstract ᅟ.

On the Influence of the Protonation States of Active Site Residues on AChE Reactivation: A QM/MM Approach

Acetylcholinesterase (AChE), an enzyme of the serine hydrolase superfamily, is a mediator of signal transmission at cholinergic synapses by catalyzing acetylcholine cleavage into acetate and choline. This enzyme is vulnerable to covalent inhibition by organophosphate compounds (like VX). Covalent inhibition of AChE does not revert spontaneously. Known reactivator compounds have limited action in restoring catalytic activity. QM/MM simulations of VX-inhibited AChE reactivation by pralidoxime (2-PAM), a classical reactivator, were performed. These afforded a broad view of the effect of protonation states of active-site residues, and provide evidence for the role of Glu202, which needs to be protonated for reactivation to occur. In situ deprotonation of 2-PAM for both protonation states of Glu202 showed that His447 is able to deprotonate 2-PAM with the assistance of Glu202. Because the active site of serine hydrolases is highly conserved, this work provides new insights on the interplay between the catalytic triad residues and this glutamate, newly identified as protonatable.

E2 and SN2 Reactions of X(-) + CH3CH2X (X = F, Cl); an ab Initio and DFT Benchmark Study

We have computed consistent benchmark potential energy surfaces (PESs) for the anti-E2, syn-E2, and SN2 pathways of X(-) + CH3CH2X with X = F and Cl. This benchmark has been used to evaluate the performance of 31 popular density functionals, covering local-density approximation, generalized gradient approximation (GGA), meta-GGA, and hybrid density-functional theory (DFT). The ab initio benchmark has been obtained by exploring the PESs using a hierarchical series of ab initio methods [up to CCSD(T)] in combination with a hierarchical series of Gaussian-type basis sets (up to aug-cc-pVQZ). Our best CCSD(T) estimates show that the overall barriers for the various pathways increase in the order anti-E2 (X = F) < SN2 (X = F) < SN2 (X = Cl) ∼ syn-E2 (X = F) < anti-E2 (X = Cl) < syn-E2 (X = Cl). Thus, anti-E2 dominates for F(-) + CH3CH2F, and SN2 dominates for Cl(-) + CH3CH2Cl, while syn-E2 is in all cases the least favorable pathway. Best overall agreement with our ab initio benchmark is obtained by representatives from each of the three categories of functionals, GGA, meta-GGA, and hybrid DFT, with mean absolute errors in, for example, central barriers of 4.3 (OPBE), 2.2 (M06-L), and 2.0 kcal/mol (M06), respectively. Importantly, the hybrid functional BHandH and the meta-GGA M06-L yield incorrect trends and qualitative features of the PESs (in particular, an erroneous preference for SN2 over the anti-E2 in the case of F(-) + CH3CH2F) even though they are among the best functionals as measured by their small mean absolute errors of 3.3 and 2.2 kcal/mol in reaction barriers. OLYP and B3LYP have somewhat higher mean absolute errors in central barriers (5.6 and 4.8 kcal/mol, respectively), but the error distribution is somewhat more uniform, and as a consequence, the correct trends are reproduced.

Formation of a(1) ions directly from oxazolone b(2) ions: an energy-resolved and computational study

It is well-known that oxazolone b2 ions fragment extensively by elimination of CO to form a2 ions, which often fragment further to form a1 ions. Less well-known is that some oxazolone b2 ions may fragment directly to form a1 ions. The present study uses energy-resolved collision-induced dissociation experiments to explore the occurrence of the direct b2→a1 fragmentation reaction. The experimental results show that the direct b2→a1 reaction is generally observed when Gly is the C-terminal residue of the oxazolone. When the C-terminal residue is more complex, it is able to provide increased stability of the a2 product in the b2→a2 fragmentation pathway. Our computational studies of the relative critical reaction energies for the b2→a2 reaction compared with those for the b2→a1 reaction provide support that the critical reaction energies are similar for the two pathways when the C-terminal residue of the oxazolone is Gly. By contrast, when the nitrogen of the oxazolone ring in the b2 ion does not bear a hydrogen, as in the Ala-Sar and Tyr-Sar (Sar = N-methylglycine) oxazolone b2 ions, a1 ions are not formed but rather neutral imine elimination from the N-terminus of the b2 ion becomes a dominant fragmentation reaction. The M06-2X/6-31+G(d,p) density functional theory calculations are in general agreement with the experimental data for both types of reaction. In contrast, the B3LYP/6-31+G(d,p) model systematically underestimates the barriers of these SN2-like b2→a1 reaction. The difference between the two methods of barrier calculation are highly significant (P < 0.001) for the b2→a1 reaction, but only marginally significant (P = 0.05) for the b2→a2 reaction. The computations provide further evidence of the limitations of the B3LYP functional when describing SN2-like reactions.

Assessment of DFT methods for studying acid gas capture by ionic liquids

For the first time, this work reports an analysis of the performance of Density Functional methods for studying acid gas capture (CO₂ and SO₂) by ionic liquids (ILs).

A simplified Tamm-Dancoff density functional approach for the electronic excitation spectra of very large molecules

Two approximations in the Tamm-Dancoff density functional theory approach (TDA-DFT) to electronically excited states are proposed which allow routine computations for electronic ultraviolet (UV)- or circular dichroism (CD) spectra of molecules with 500-1000 atoms. Speed-ups compared to conventional time-dependent DFT (TD-DFT) treatments of about two to three orders of magnitude in the excited state part at only minor loss of accuracy are obtained. The method termed sTDA ("s" for simplified) employs atom-centered Löwdin-monopole based two-electron repulsion integrals with the asymptotically correct 1∕R behavior and perturbative single excitation configuration selection. It is formulated generally for any standard global hybrid density functional with given Fock-exchange mixing parameter ax. The method performs well for two standard benchmark sets of vertical singlet-singlet excitations for values of ax in the range 0.2-0.6. The mean absolute deviations from reference data are only 0.2-0.3 eV and similar to those from standard TD-DFT. In three cases (two dyes and one polypeptide), good mutual agreement between the electronic spectra (up to 10-11 eV excitation energy) from the sTDA method and those from TD(A)-DFT is obtained. The computed UV- and CD-spectra of a few typical systems (e.g., C60, two transition metal complexes, [7]helicene, polyalanine, a supramolecular aggregate with 483 atoms and about 7000 basis functions) compare well with corresponding experimental data. The method is proposed together with medium-sized double- or triple-zeta type atomic-orbital basis sets as a quantum chemical tool to investigate the spectra of huge molecular systems at a reliable DFT level.

Simulation studies of the Cl- + CH3I SN2 nucleophilic substitution reaction: comparison with ion imaging experiments

In the previous work of Mikosch et al. [Science 319, 183 (2008)], ion imaging experiments were used to study the Cl(-) + CH3I → ClCH3 + I(-) reaction at collision energies E(rel) of 0.39, 0.76, 1.07, and 1.9 eV. For the work reported here MP2(fc)/ECP/d direct dynamics simulations were performed to obtain an atomistic understanding of the experiments. There is good agreement with the experimental product energy and scattering angle distributions for the highest three E(rel), and at these energies 80% or more of the reaction is direct, primarily occurring by a rebound mechanism with backward scattering. At 0.76 eV there is a small indirect component, with isotropic scattering, involving formation of the pre- and post-reaction complexes. All of the reaction is direct at 1.07 eV. Increasing E(rel) to 1.9 eV opens up a new indirect pathway, the roundabout mechanism. The product energy is primarily partitioned into relative translation for the direct reactions, but to CH3Cl internal energy for the indirect reactions. The roundabout mechanism transfers substantial energy to CH3Cl rotation. At E(rel) = 0.39 eV both the experimental product energy partitioning and scattering are statistical, suggesting the reaction is primarily indirect with formation of the pre- and post-reaction complexes. However, neither MP2 nor BhandH/ECP/d simulations agree with experiment and, instead, give reaction dominated by direct processes as found for the higher collision energies. Decreasing the simulation E(rel) to 0.20 eV results in product energy partitioning and scattering which agree with the 0.39 eV experiment. The sharp transition from a dominant direct to indirect reaction as E(rel) is lowered from 0.39 to 0.20 eV is striking. The lack of agreement between the simulations and experiment for E(rel) = 0.39 eV may result from a distribution of collision energies in the experiment and/or a shortcoming in both the MP2 and BhandH simulations. Increasing the reactant rotational temperature from 75 to 300 K for the 1.9 eV collisions, results in more rotational energy in the CH3Cl product and a larger fraction of roundabout trajectories. Even though a ClCH3-I(-) post-reaction complex is not formed and the mechanistic dynamics are not statistical, the roundabout mechanism gives product energy partitioning in approximate agreement with phase space theory.

Corrected small basis set Hartree-Fock method for large systems

A quantum chemical method based on a Hartree-Fock calculation with a small Gaussian AO basis set is presented. Its main area of application is the computation of structures, vibrational frequencies, and noncovalent interaction energies in huge molecular systems. The method is suggested as a partial replacement of semiempirical approaches or density functional theory (DFT) in particular when self-interaction errors are acute. In order to get accurate results three physically plausible atom pair-wise correction terms are applied for London dispersion interactions (D3 scheme), basis set superposition error (gCP scheme), and short-ranged basis set incompleteness effects. In total nine global empirical parameters are used. This so-called Hartee-Fock-3c (HF-3c) method is tested for geometries of small organic molecules, interaction energies and geometries of noncovalently bound complexes, for supramolecular systems, and protein structures. In the majority of realistic test cases good results approaching large basis set DFT quality are obtained at a tiny fraction of computational cost.

Dispersion corrected hartree-fock and density functional theory for organic crystal structure prediction

We present and evaluate dispersion corrected Hartree-Fock (HF) and Density Functional Theory (DFT) based quantum chemical methods for organic crystal structure prediction. The necessity of correcting for missing long-range electron correlation, also known as van der Waals (vdW) interaction, is pointed out and some methodological issues such as inclusion of three-body dispersion terms are discussed. One of the most efficient and widely used methods is the semi-classical dispersion correction D3. Its applicability for the calculation of sublimation energies is investigated for the benchmark set X23 consisting of 23 small organic crystals. For PBE-D3 the mean absolute deviation (MAD) is below the estimated experimental uncertainty of 1.3 kcal/mol. For two larger π-systems, the equilibrium crystal geometry is investigated and very good agreement with experimental data is found. Since these calculations are carried out with huge plane-wave basis sets they are rather time consuming and routinely applicable only to systems with less than about 200 atoms in the unit cell. Aiming at crystal structure prediction, which involves screening of many structures, a pre-sorting with faster methods is mandatory. Small, atom-centered basis sets can speed up the computation significantly but they suffer greatly from basis set errors. We present the recently developed geometrical counterpoise correction gCP. It is a fast semi-empirical method which corrects for most of the inter- and intramolecular basis set superposition error. For HF calculations with nearly minimal basis sets, we additionally correct for short-range basis incompleteness. We combine all three terms in the HF-3c denoted scheme which performs very well for the X23 sublimation energies with an MAD of only 1.5 kcal/mol, which is close to the huge basis set DFT-D3 result.

Selectivity of labeled bromoethylamine for protein alkylation

Alkylation of cysteine residues has been used extensively for characterization of proteins and their mode of action in biological systems, research endeavors that are at the core of proteomics. Treatment with a simple alkylating agent such as [2-(13)C] bromoethylamine would result in labeled thialysine at the ε-position. This chemical modification of proteins would allow investigations via both (13)C NMR spectroscopy and mass spectrometry. However [2-(13)C] labeled bromoethylamine is not available commercially. We investigated its synthesis at acid pH with the goal of obtaining singly labeled bromoethylamine and understanding the mechanistic details of the reaction. Based on our experimental and theoretical results, bromination of [2-(13)C] labeled ethanolamine in acidic conditions takes place via exclusive attack of the nucleophile (HBr) at the hydroxyl bearing C. Moreover, hydrogen bonding guides the nucleophilic attack, resulting in no label scrambling of the bromoethylamine product. Protein alkylation at cysteine residue with the synthesized Br(13)CH(2)CH(2)NH(2)-HBr is successful. Ab initio calculations in which CH(3)SH serves as a model for the cysteine residue suggest that in gas phase intermolecular attack by the sulfur bearing nucleophile is favored over the intramolecular substitution by the amino group by 15.4 kJ mol(-1). Solution modeling shows that the trend is preserved at basic pH, which is the experimental one, but is reversed at neutral pH.

On the PES for the interaction of an H atom with an H chemisorbate on a graphenic platelet

Motivated by the problem of H(2) formation in diffuse clouds of the interstellar medium (ISM), we study the effect of including van der Waals-type corrections in DFT calculations on the entrance PES of the Eley-Rideal reaction H(b) + H(a)-GR → H(b)-H(a) + GR for a graphenic surface GR. The present calculations make use of the PBE-D3 dispersion corrected functional of Grimme et al. (2010) and are carried out on cluster models of graphenic surfaces: C(24)H(12) and C(54)H(18). To assess the soundness of the chosen functional we start by revisiting the H-GR adsorption potential. We find a satisfactory on top physisorption well (43-48 meV) correctly located at an H-GR distance of 3 Å. We then revisit the H(b)-H(a)-GR system using both the PW91 and PBE functionals. Our calculations do not reproduce the tiny potential barrier reported earlier for large H(b)distances from the surface. The barrier in the calculations of Sidis et al. (2000) and Morisset et al. (2003, 2004) has been traced to their previous use of an LSDA + POSTSCF PW91 procedure rather than the genuine PW91 one. The new PBE-D3 PES for the H(b)-H(a)-GR system is reported as a function of the H(b) distance to the surface and its impact parameter relative to the H(a) chemisorbate for the so-called "fixed puckered" ("diabatic" or "sudden") approach. The results are discussed in relation to recent experimental and theoretical work.

Gas-phase structure and fragmentation pathways of singly protonated peptides with N-terminal arginine

The gas-phase structures and fragmentation pathways of the singly protonated peptide arginylglycylaspartic acid (RGD) are investigated by means of collision-induced-dissociation (CID) and detailed molecular mechanics and density functional theory (DFT) calculations. It is demonstrated that despite the ionizing proton being strongly sequestered at the guanidine group, protonated RGD can easily be fragmented on charge directed fragmentation pathways. This is due to facile mobilization of the C-terminal or aspartic acid COOH protons thereby generating salt-bridge (SB) stabilized structures. These SB intermediates can directly fragment to generate b(2) ions or facilely rearrange to form anhydrides from which both b(2) and b(2)+H(2)O fragments can be formed. The salt-bridge stabilized and anhydride transition structures (TSs) necessary to form b(2) and b(2)+H(2)O are much lower in energy than their traditional charge solvated counterparts. These mechanisms provide compelling evidence of the role of SB and anhydride structures in protonated peptide fragmentation which complements and supports our recent findings for tryptic systems (Bythell, B. J.; Suhai, S.; Somogyi, A.; Paizs, B. J. Am. Chem. Soc. 2009, 131, 14057-14065.). In addition to these findings we also report on the mechanisms for the formation of the b(1) ion, neutral loss (H(2)O, NH(3), guanidine) fragment ions, and the d(3) ion.