Electron attachment step in electron capture dissociation (ECD) and electron transfer dissociation (ETD)

We have made use of classical dynamics trajectory simultions and ab initio electronic structure calculations to estimate the cross sections with which electrons are attached (in electron capture dissociation (ECD)) or transferred (in electron transfer dissociation (ETD)) to a model system that contained both an S-S bond that is cleaved and a -NH(3)(+) positively charged site. We used a Landau-Zener-Stueckelberg curve-crossing approximation to estimate the ETD rates for electron transfer from a CH(3)(-) anion to the -NH(3)(+) Rydberg orbital or the S-S sigma* orbital. We draw conclusions about ECD from our ETD results and from known experimental electron-attachment cross sections for cations and sigma-bonds. We predict the cross section for ETD at the positive site of our model compound to be an order of magnitude larger than that for transfer to the Coulomb-stabilized S-S bond site. We also predict that, in ECD, the cross section for electron capture at the positive site will be up to 3 orders of magnitude larger than that for capture at the S-S bond site. These results seem to suggest that attachment to such positive sites should dominate in producing S-S bond cleavage in our compound. However, we also note that cleavage induced by capture at the positive site will be diminished by an amount that is related to the distance from the positive site to the S-S bond. This dimunition can render cleavage through Coulomb-assisted S-S sigma* attachment competitive for our model compound. Implications for ECD and ETD of peptides and proteins in which SS or N-C(alpha) bonds are cleaved are also discussed, and we explain that such events are most likely susceptible to Coulomb-assisted attachment, because the S-S sigma* and C=O pi* orbitals are the lowest-lying antibonding orbitals in most peptides and proteins.     

The role of excited Rydberg States in electron transfer dissociation

Ab initio electronic structure methods are used to estimate the cross sections for electron transfer from donor anions having electron binding energies ranging from 0.001 to 0.6 eV to each of three sites in a model disulfide-linked molecular cation. The three sites are (1) the S-S sigma(*) orbital to which electron attachment is rendered exothermic by Coulomb stabilization from the nearby positive site, (2) the ground Rydberg orbital of the -NH(3)(+) site, and (3) excited Rydberg orbitals of the same -NH(3)(+) site. It is found that attachment to the ground Rydberg orbital has a somewhat higher cross section than attachment to either the sigma orbital or the excited Rydberg orbital. However, it is through attachment either to the sigma(*) orbital or to certain excited Rydberg orbitals that cleavage of the S-S bond is most likely to occur. Attachment to the sigma(*) orbital causes prompt cleavage because the sigma energy surface is repulsive (except at very long range). Attachment to the ground or excited Rydberg state causes the S-S bond to rupture only once a through-bond electron transfer from the Rydberg orbital to the S-S sigma(*) orbital takes place. For the ground Rydberg state, this transfer requires surmounting an approximately 0.4 eV barrier that renders the S-S bond cleavage rate slow. However, for the excited Rydberg state, the intramolecular electron transfer has a much smaller barrier and is prompt.     

Side-chain fragmentation of alkylated cysteine residues in electron capture dissociation mass spectrometry

The recent development of novel fragmentation processes based on either electron capture directly or transfer from an anion show great potential for solving problems in proteomics that are intractable by the more widely employed thermal-based fragmentation processes such as collision induced dissociation. The dominant fragmentation occurring upon electron capture dissociation of peptides is cleavage of N-C alpha bonds in the peptide backbone to form c and z* ions. In the case of disulfide-linked peptides, it has also been shown that electron capture on one of the cystine sulfur atoms is favored, resulting in cleavage of the disulfide bond. In this study, we report that electron capture on the sulfur of alkylated cysteine residues is also a dominant process, causing cysteine side-chain loss from z* ions.     

Comments on some chemical properties of diphenyl disulfide and its derivatives: II.Correlating the studies of reductive cleavage of disulfide bonds with the studies of anti-human immunodeficiency virus activity with lowest unoccupied molecular orbital properties

The higher anti-human immunodeficiency virus activity of a symmetrical 2,2′-disubstitued derivative of diphenyl disulfide (DPDS) has been explained by the lower energy of the lowest unoccupied molecular orbital (LUMO), resulted from a better hydrogen bond stabilization of the σ*_SS bond orbital (BO). This conclusion entails the participation of σ*_SS BO in constructing the LUMO. The higher content of σ*_SS BO, compared to π*_CC BOs of phenyl groups, in LUMO of DPDS has been found through analysis of the LUMO of DPDS expanded in the BO space. The high content of σ*_SS BO (%σ*_SS) in the LUMO of DPDS has laid the foundation for the formation of σ-type radical anion intermediate in the stepwise reductive cleavage of disulfide bond in the symmetrical 4,4′-disubstitued DPDS derivatives. For the nine 4,4′-disubstituted DPDS-derivatives under reductive cleavage studies, the increasing %σ*_SS in the LUMOs is parallel to the increasing value of inner reorganization energy.     

Comments on chemical properties reported for diphenyl disulfide and its derivatives: The merit of the phenyl groups

The symmetrical 2,2′-disubstitued derivatives of diphenyl disulfide showing widely spanning rates of electrophilic attack of the HIV-1 nucleocapsid protein p7 zinc fingers have been rationalized, based on the lowest unoccupied molecular orbital (LUMO)-lowering approach, by the substituents' π-effects and the hydrogen bond stabilization effects. In the 2,2′-amide- and 4,4′-N-amide-substituted derivatives, the extent of LUMO lowering has been reduced by the destabilization of lone-pair bond orbital, lp(N), present on the nitrogen atom of N-amide. From the natural bond orbital viewpoint, hydrogen bond stabilization of LUMO is mainly governed by stabilization of the σ*_SS bond orbital.     

The Pi‐Nature of Methyl Obtained by the Natural Bond Orbital Method: Orbital‐Based Rationalizations of Site‐Dependent Substitution Effects on Fine Color‐Tuning of Luminescence

Both C‐H bonding and antibonding (σ_CH and σ*_CH) of a methyl group would contribute to the highest occupied or lowest unoccupied molecular orbitals (HOMO or LUMO) in methylated derivatives of Ir(ppz)₂ 3 iq (ppz = 1‐phenylpyrazole and 3iq = isoquinoline‐3‐carboxylate). This is found by analysis of HOMO (or LUMO) formed by linear combination of bond orbitals using the natural bond orbital (NBO) method. The elevated level of HOMO (or LUMO) uniformly found for each methylated derivative, indicating the σ_CH‐destabilization outweighs the σ*_CH‐stabilization. To broaden the HOMO‐LUMO gap, methylation at a carbon having smaller contribution to HOMO and/or larger contribution to LUMO is suggested.     

Ab initio through‐space/bond interaction analysis of the long C–C bonds in Bi(anthracene‐9,10‐dimethylene) photoisomers

The bi(anthracene‐9,10‐dimethylene) photoisomer has remarkably long C–C single bonds. To examine the lengthening of the C–C bond, we propose a novel procedure for quantitatively analyzing orbital interactions in a molecule at the level of the ab initio molecular orbital method. In this procedure, we can cut off the specific through‐space/bond interactions in a molecule by artificially increasing the absolute magnitude of the exponents in a Gaussian function. Then, the spatial orbital interactions were perfectly cut off, and, each term that makes up the total energy, that is, the nuclear–electron attractions, the electron–electron repulsions, and the nuclear–nuclear repulsions cancel each other. Several model molecules of the photoisomer were analyzed by this procedure. It was found that the orbital interaction between the p orbital on the benzene ring and the σ* orbital on the C–C bond in question, σ→σ* electron transfer through π orbital, weakens the C–C bond efficiently when these orbitals were located in the “periplanar” conformation. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Coulomb-assisted dissociative electron attachment: application to a model peptide

The fragmentation of positively charged gas-phase samples of peptides is used to infer the primary structure of such molecules. In electron capture dissociation (ECD) experiments, very low-energy electrons attach to the sample and rupture bonds to effect the fragmentation. It turns out that ECD fragmentation tends to produce cleavage of very specific types of bonds. In earlier works by this group, it has been suggested that the presence of positive charges produces stabilizing Coulomb potentials that allow low-energy electrons to exothermically attach to sigma orbitals of certain bonds and thus to cleave those bonds. In the present effort, the stabilizing effects of Coulomb potentials due to proximal positive charges are examined for a small model peptide molecule that contains a wide range of bond types. Direct attachment of an electron to the sigma orbitals of eight different bonds as well as indirect sigma bond cleavage, in which an electron first binds to a carbonyl C=O pi orbital, are examined using ab initio methods. It is found that direct attachment to and subsequent cleavage of any of the eight sigma bonds is not likely except for highly positively charged samples. It is also found that attachment to a C=O pi orbital followed by cleavage of the nitrogen-to-alpha-carbon bond is the most likely outcome. Interestingly, this bond cleavage is the one that is seen most commonly in ECD experiments. So, the results presented here seem to offer good insight into one aspect of the ECD process, and they provide a means by which one can estimate (on the basis of a simple Coulomb energy formula) which bonds may be susceptible to cleavage by low-energy electron attachment.     

Ab initio assessment of the electronic structure of 5α-reduced progestins

Progesterone (P) yields to 5α-reduced progestins, namely 5α-pregnanedione (DHP), tetrahydroprogesterone (THP), and allopregnanolone (ALLO-P). The geometries and electronic structure of these steroids were assessed by ab initio calculations using the 6-31G* basis set. The parameters measured were bond distances, valence angles, and dihedral angles. Likewise, the following were calculated: total energy; frontier orbitals, i.e., highest occupied molecular orbital (HOMO); lowest unoccupied molecular orbital (LUMO); dipole moment; atomic charges; and electrostatic potentials. The frontier orbitals of P were located at the π-double bond. However, the HOMO of the 5α-progestins was extended into the molecule, while the LUMO was confined at the C20 carbonyl group. The atomic charges, electronic density surfaces and electrostatic potentials showed patterns according to the stereochemical arrangement of the C3 and C20 carbonyl and hydroxyl functional groups. Interestingly, P and THP showed the larger dipole moment and high electronic density at the A-ring because the double bond and the 3α-hydroxy group, respectively. The present results might explain to some extent the metabolism of the studied progestins. Similarly, some physicochemical properties, such as dipole moments and electrostatic potentials, seem related with important biological actions such as uterine contractility and control of gonadotropin secretion. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 67: 329–338, 1998     

Theory of electron localization and its application to blue-shifting hydrogen bonds

We propose a theory of electron localization or stabilization by electron localization through the interactions between occupied (i) and vacant (j*) orbitals under certain conditions, which have been believed so far to cause only electron delocalization. Electrons localize when the electrons redistributed by the interaction are more stable in the i-th occupied orbital than in the overlap region: h(ij*) > s(ij*)h(ii) for s(ij*) > 0. Electron delocalization occurs when h(ij*) < s(ij*)h(ii) for s(ij*) > 0. The h(ij*) and s(ij*)h(ii) terms represent the energy of the electrons in the overlap region and the energy of the redistributed electrons in the occupied orbital, respectively. The theory of electron localization is substantiated by the correlation of the C-H bond lengths of fluorinated methanes H(4-n)CF(n) (n = 1, 2, 3) to the electron population of the σ(CH) bonding orbital, and successfully applied to understanding blue-shifting hydrogen bonds in F(3)CH···X (X = CO, N(2), OC, Ne, OC(CH(3))(2)) and designing some proton donors, HCO(2)CH(3) and hypervalent molecules HPF(4) and HSF(5), for blue-shifting hydrogen bonds.     

The effects of C-S and C-Se bonds on torquoselectivity: stereoselective olefination of α-thio and α-selenoketones with ynolates

Highly Z-selective olefination of acyclic α-thio and α-selenoketones with ynolates has been achieved, and the theoretical calculations of the transition states in the ring-opening of the intermediates, the β-lactone enolates, revealed that the torquoselectivity was controlled by the secondary orbital interactions between the σ orbital of the C-S bond or a lone pair orbital on the S and σ^∗ orbitals of the breaking C-O bond, and the σ orbital of the breaking C-O bond or a lone pair orbital on the O on the ring and the σ^∗ orbitals of the C-S bond. The synthetic applications of the resulting olefins are also shown.     

The Bond Order of C2 from a Strictly N‐Representable Natural Orbital Energy Functional Perspective

The bond order of the ground electronic state of the carbon dimer has been analyzed in the light of natural orbital functional theory calculations carried out with an approximate, albeit strictly N‐representable, energy functional. Three distinct solutions have been found from the Euler equations of the minimization of the energy functional with respect to the natural orbitals and their occupation numbers, which expand upon increasing values of the internuclear coordinate. In the close vicinity of the minimum energy region, two of the solutions compete around a discontinuity point. The former, corresponding to the absolute minimum energy, features two valence natural orbitals of each of the following symmetries, σ, σ*, π and π*, and has three bonding interactions and one antibonding interaction, which is very suggestive of a bond order large than two but smaller than three. The latter, features one σ–σ* linked pair of natural orbitals and three degenerate pseudo‐bonding like orbitals, paired each with one triply degenerate pseudo‐antibonding orbital, which points to a bond order larger than three. When correlation effects, other than Hartree–Fock for example, between the paired natural orbitals are accounted for, this second solution vanishes yielding a smooth continuous dissociation curve. Comparison of the vibrational energies and electron ionization energies, calculated on this curve, with their corresponding experimental marks, lend further support to a bond order for C ₂ intermediate between acetylene and ethylene.     

Unpaired electrons at the second‐order reduced density matrix level: Covalent bonding,and coulomb and fermi correlations in closed shell systems

The usual one‐electron populations in atomic orbitals of closed shell systems are split into unpaired and paired at the (spin‐dependent) second‐order reduced density matrix level. The unpaired electron in an orbital is defined as the “simultaneous occurrence of an electron and an electron hole of opposite spins in the same spatial orbital,” which for simplicity is called “electropon.” The electropon population in a given orbital reveals whether and to what degree the Coulomb correlations, and hence, the chemical bonding between this orbital and the remaining orbitals of the system are globally favorable or unfavorable. The interaction of two electropons in two target orbitals reveals the quality (favorable or unfavorable) and the strength of the covalent bonding between these orbitals; this establish a bridge between the notion of “unpaired electrons” and the traditional covalent structure of valence‐bond (VB) theory. Favorable/unfavorable bonding between two orbitals is characterized by the positive/negative (Coulomb) correlation of two electropons of opposite spins, or alternatively, by the negative/positive (Fermi) correlation of two parallel spin electropons. A spin‐free index is defined, and the relationship between the electropon viewpoint for chemical bonding and the well‐known two‐electron Coulomb and Fermi correlations is established. Benchmark calculations are achieved for ethylene, hexatriene, benzene, pyrrole, methylamine, and ammonia molecules on the basis of physically meaningful natural orbitals. The results, obtained in the framework of both orthogonal and nonorthogonal population analysis methods, provide the same conceptual pictures, which are in very good agreement with elementary chemical knowledge and VB theory. © 2013 Wiley Periodicals, Inc.     

S K-edge XAS and DFT calculations on SAM dependent pyruvate formate-lyase activating enzyme: nature of interaction between the Fe4S4 cluster and SAM and its role in reactivity

S K-edge X-ray absorption spectroscopy on the resting oxidized and the S-adenosyl-l-methionine (SAM) bound forms of pyruvate formate-lyase activating enzyme are reported. The data show an increase in pre-edge intensity, which is due to additional contributions from sulfide and thiolate of the Fe(4)S(4) cluster into the C-S σ* orbital. This experimentally demonstrates that there is a backbonding interaction between the Fe(4)S(4) cluster and C-S σ* orbitals of SAM in this inner sphere complex. DFT calculations that reproduce the data indicate that this backbonding is enhanced in the reduced form and that this configurational interaction between the donor and acceptor orbitals facilitates the electron transfer from the cluster to the SAM, which otherwise has a large outer sphere electron transfer barrier. The energy of the reductive cleavage of the C-S bond is sensitive to the dielectric of the protein in the immediate vicinity of the site as a high dielectric stabilizes the more charge separated reactant increasing the reaction barrier. This may provide a mechanism for generation of the 5'-deoxyadenosyl radical upon substrate binding.     

Orbital interactions in large molecular systems: Adsorption of H2 on Ni surfaces

The concept of orbital interactions is applied to the adsorption of H₂ on to the Ni (110) and (111) surfaces. We calculate first two orbitals of a Ni cluster one of which forms an orbital pair with the σ MO and the other with the σ*MO of a H₂ molecule. Contributions of these paired orbitals of fragments to the density of states of the surface-adsorbate extended system are then examined. It is shown that the orbital of the surface that participates in electron delocalization to σ* of the H₂ molecule is located significantly below the Fermi level both in the (110) and in the (111) adsorption models. The σ MO of H₂ and its counterpart of the surface represent mainly overlap repulsion which is shown to be stronger on the (111) surface than on the (110) surface. It is feasible to understand chemical interactions of large systems by using the paired interacting orbitals.     

Theoretical Study and NBO Analysis of the Decomposition Kinetics of Oxetane, 2-Methyloxetane and 2,2-Dimethyloxetane

A theoretical study of the thermal decomposition kinetics of oxetane (1), 2-methyloxetane (2), and 2,2-dimethyloxetane (3) has been carried out at the B3LYP/6-311+G**, B3PW91/6-311+G**, and MPW1PW91/6-311+G** levels of theory. The MPW1PW91/6-311+G** method was found to give a reasonable good agreement with the experimental kinetics and thermodynamic parameters. The decomposition reaction of compounds 1～3 yields formaldehyde and the corresponding substituted olefin. Based on the optimized ground state geometries using MPW1PW91/6-311+G** method, the natural bond orbital (NBO) analysis of donor-acceptor (bond-antibond) interactions revealed that the stabilization energies associated with the electronic delocalization from σC3-C4 bonding to σ*O1-C2 antibonding orbitals decrease from compounds 1 to 3. The σC3-C4→σO1-C2 resonance energies for compounds 1～3 are 2.63, 2.59 and 2.45 kcal mol-1, respectively. Further, the results showed that the energy gaps between σC3-C4 bonding and σ*O1-C2 antibonding orbitals decrease from compounds 1 to 3. Also, the decomposition process in these compounds are controlled by σ→σ* resonance energies. Moreover, the obtained order of energy barriers could be explained by the number of electron-releasing methyl groups substituted to the Csp3 atom (which is attached to oxygen atom). NBO analysis shows that the occupancies of σCsp3-O bonds decrease for compounds 1～3 as 3<2<1, and those of σCsp3-O bonds increase in the opposite order (3 > 2 > 1). This fact illustrates a comparatively easier thermal decomposition of the sCsp3-O bond in compound 3 compared to compound 2, and in compound 2 compared to compound 1. NBO results indicate that these reactions are occurring through a concerted and asynchronous four-membered cyclic transition state type of mechanism.     

Density functional study of platinum polyyne monomer,oligomer, and polymer: Ground state geometrical and electronic structures

Geometries of monomers and oligomers of a platinum polyyne and its free ligands were optimized using density functional theory with B3LYP hybrid functional. The LANL2DZ basis set was used for Pt and the 6‐31G* for other atoms in geometry optimizations. The electronic structures of these compounds were analyzed using Stuttgart/Dresden ECPs (SDD) basis set for metal atoms and 6‐311G* for others. The polymerization has very little effect on the bond lengths and by introducing the metal, the acetylide bond length increases slightly. The strong overlap between metal sp_x orbitals and σp_x orbitals of acetylides results in localized σ bonding. The hybridization between the ligand pπ orbitals and the platinum dπ orbital resulted in the π‐conjugation enhancement. This conjugation enhancement causes some effects such as the highest‐occupied molecular orbital–lowest‐unoccupied molecular orbital gap reduction and charge transfer characteristic of low‐energy vertical transitions. © 2013 Wiley Periodicals, Inc.     

Strong electron correlation in the decomposition reaction of dioxetanone with implications for firefly bioluminescence

Dioxetanone, a key component of the bioluminescence of firefly luciferin, is itself a chemiluminescent molecule due to two conical intersections on its decomposition reaction surface. While recent calculations of firefly luciferin have employed four electrons in four active orbitals [(4,4)] for the dioxetanone moiety, a study of dioxetanone [F. Liu et al., J. Am. Chem. Soc. 131, 6181 (2009)] indicates that a much larger active space is required. Using a variational calculation of the two-electron reduced-density-matrix (2-RDM) [D. A. Mazziotti, Acc. Chem. Res. 39, 207 (2006)], we present the ground-state potential energy surface as a function of active spaces from (4,4) to (20,17) to determine the number of molecular orbitals required for a correct treatment of the strong electron correlation near the conical intersections. Because the 2-RDM method replaces exponentially scaling diagonalizations with polynomially scaling semidefinite optimizations, we readily computed large (18,15) and (20,17) active spaces that are inaccessible to traditional wave function methods. Convergence of the electron correlation with active-space size was measured with complementary RDM-based metrics, the von Neumann entropy of the one-electron RDM as well as the Frobenius and infinity norms of the cumulant 2-RDM. Results show that the electron correlation is not correctly described until the (14,12) active space with small variations present through the (20,17) space. Specifically, for active spaces smaller than (14,12), we demonstrate that at the first conical intersection, the electron in the σ(?) orbital of the oxygen-oxygen bond is substantially undercorrelated with the electron of the σ orbital and overcorrelated with the electron of the carbonyl oxygen's p orbital. Based on these results, we estimate that in contrast to previous treatments, an accurate calculation of the strong electron correlation in firefly luciferin requires an active space of 28 electrons in 25 orbitals, beyond the capacity of traditional multireference wave function methods.