Shattering of Peptide ions on self-assembled monolayer surfaces

Time- and collision energy-resolved surface-induced dissociation (SID) of des-Arg(1)- and des-Arg(9)-bradykinin on a fluorinated self-assembled monolayer (SAM) surface was studied by use of a novel Fourier transform ion cyclotron resonance mass spectrometer (FT-ICR MS) specially equipped to perform SID experiments. Time-resolved fragmentation efficiency curves (TFECs) were modeled by an RRKM-based approach developed in our laboratory that utilizes a very flexible analytical expression for the internal energy deposition function capable of reproducing both single- and multiple-collision activation in the gas phase and excitation by collisions with a surface. Both experimental observations and modeling establish a very sharp transition in the dynamics of ion-surface interaction: the shattering transition. The experimental signature for this transition is the appearance of prompt (time-independent) fragmentation, which becomes dominant at high collision energies. Shattering opens a variety of dissociation pathways that are not accessible to slow collisional and thermal ion activation. This results in much better sequence coverage for the singly protonated peptides than dissociation patterns obtained with any of the slow activation methods. Modeling demonstrated that, for short reaction delays, dissociation of these peptides is solely determined by shattering. Internal energies required for shattering transition are approximately the same for des-Arg(1) and des-Arg(9)-bradykinin, resulting in the overlap of fragmentation efficiency curves obtained at short reaction delays. At longer delay times, parent ions depletion is mainly determined by a slow decay rate and fragmentation efficiency curves for des-Arg(1) and des-Arg(9)-bradykinin diverge. Dissociation thresholds of 1.17 and 1.09 eV and activation entropies of -22.2 and -23.3 cal/(mol K) were obtained for des-Arg(1) and des-Arg(9)-bradykinin from RRKM modeling of time-resolved data. Dissociation parameters for des-Arg(1)-bradykinin are in good agreement with parameters derived from thermal experiments. However, there is a significant discrepancy between the thermal data and dissociation parameters for des-Arg(9)-bradykinin obtained in this study. The difference is attributed to the differences in conformations that undergo thermal activation and activation by ion-surface collisions prior to dissociation.     

Size matters! Fragmentation chemistry of [Cu(L)n]2+ complexes of diacylglycerophosphocholines as a function of coordination number (n = 2-7)

[Cu(L)(n)](2+) complexes of 1,2-dihexanoyl-sn-glycero-3-phosphocholine (L = D6PC) are formed upon electrospray ionization mass spectrometry (ESI-MS) of an 8 mM solution of D6PC with 4 mM CuCl(2) in 10 mM ammonium acetate buffer, pH 6.1. The collision-induced dissociation (CID) reactions of the [Cu(L)(n)](2+) complexes were examined in a linear ion trap mass spectrometer. A rich fragmentation chemistry was observed, including: loss of a neutral ligand; intermolecular ligand-ligand S(N)2 methylation; metal ion induced ligand fragmentation via carboxylate abstraction; and phosphate abstraction. The dominant reaction channel depends on the size (n) of the complex. Thus loss of neutral ligand(s) is the sole reaction channel for n = 5-7. At n = 4, S(N)2 methylation and carboxylate abstraction start to compete with neutral ligand loss. At n = 2 the carboxylate abstraction and phosphate abstraction reactions dominate the CID spectrum. The carboxylate abstraction and phosphate abstraction reactions are likely to be driven via neighboring group pathways. PM3 calculations, carried out to compare competing neighboring pathways based on the relative stabilities of the product ions, suggest a preference for five-membered ring formation for ligand fragmentation involving both carboxylate and phosphate abstraction.     

A quantitative study of intrinsic non-covalent interactions within complexes of α-cyclodextrin and benzoate derivatives

A novel deconvolution method for energy-resolved reaction cross sections is applied to determine intrinsic gas-phase dissociation energies for non-covalent α-cyclodextrin host-guest complexes. M06-2X//M06-L/6-31+G(d,p) calculations reproduce the experimental results and enable us to quantify the contribution of intermolecular hydrogen bonding.     

Gas-phase binding of non-covalent protein complexes between bovine pancreatic trypsin inhibitor and its target enzymes studied by electrospray ionization tandem mass spectrometry

The potential of electrospray ionization (ESI) mass spectrometry (MS) to detect non-covalent protein complexes has been demonstrated repeatedly. However, questions about correlation of the solution and gas-phase structures of these complexes still produce vigorous scientific discussion. Here, we demonstrate the evaluation of the gas-phase binding of non-covalent protein complexes formed between bovine pancreatic trypsin inhibitor (BPTI) and its target enzymes over a wide range of dissociation constants. Non-covalent protein complexes were detected by ESI-MS. The abundance of the complex ions in the mass spectra is less than expected from the values of the dissociation constants of the complexes in solution. Collisionally activated dissociation (CAD) tandem mass spectrometry (MS/MS) and a collision model for ion activation were used to evaluate the binding of non-covalent complexes in the gas phase. The internal energy required to induce dissociation was calculated for three collision gases (Ne, Ar, Kr) over a wide range of collision gas pressures and energies using an electrospray ionization source. The order of binding energies of the gas-phase ions for non-covalent protein complexes formed by the ESI source and assessed using CAD-MS/MS appears to differ from that of the solution complexes. The implication is that solution structure of these complexes was not preserved in the gas phase.     

Modeling metal cation-phosphate interactions in nucleic acids in the gas phase via alkali metal cation-triethyl phosphate complexes

Threshold collision-induced dissociation techniques are employed to determine the bond dissociation energies (BDEs) of complexes of alkali metal cations, Na+, K+, Rb+, and Cs+, to triethyl phosphate (TEP). The primary and lowest energy dissociation pathway in all cases is the endothermic loss of the neutral TEP ligand. Theoretical electronic structure calculations at the B3LYP/6-311+G(2d,2p)//B3LYP/6-31G* level of theory are used to determine the structures, molecular parameters, and theoretical estimates for the BDEs of these complexes. For the complexes to Rb+ and Cs+, theoretical calculations were performed using hybrid basis sets in which the effective core potentials and valence basis sets of Hay and Wadt were used to describe the alkali metal cation, while the standard basis sets were used for all other atoms. The agreement between theory and experiment is excellent for the complexes to Na+ and K+ and is somewhat less satisfactory for the complexes to the heavier alkali metal cations, Rb+ and Cs+, where effective core potentials were used to describe the cation. The trends in the binding energies are examined. The binding of alkali metal cations to triethyl phosphate is compared with that to trimethylphosphate.     

Dissociation Energy Computations for Saline Bonds Implied in Interactions Mediated by Peptidoglicans

Summary: Dissociation energy and hydration energy calculations, in water solution, are presented for saline bonds mediated by Ca²⁺ and Mg²⁺ ions with Brőnstedt type bases ( COO ,  OSO₃ ,  OH). A computationally intensive method, Polarisable Continuum Model (PCM) using 6-31G^* basis set, was applied. Hydration energies were computed by various methods, as well as dissociation energies of some L₂M complexes. L₂Ca complexes result as more stable against dissociation than L₂Mg complexes. Hydration energy calculation results, for some of the methods, here used, seem rather reliable as compared to experimental results.     

A simple model for metal cation-phosphate interactions in nucleic acids in the gas phase: Alkali metal cations and trimethyl phosphate

Threshold collision-induced dissociation techniques are employed to determine the bond dissociation energies (BDEs) of complexes of alkali metal cations to trimethyl phosphate, TMP. Endothermic loss of the intact TMP ligand is the only dissociation pathway observed for all complexes. Theoretical calculations at the B3LYP/6-31G* level of theory are used to determine the structures, vibrational frequencies, and rotational constants of neutral TMP and the M+(TMP) complexes. Theoretical BDEs are determined from single point energy calculations at the B3LYP/6-311+G(2d,2p) level using the B3LYP/6-31G* optimized geometries. The agreement between theory and experiment is reasonably good for all complexes except Li+(TMP). The absolute M+-(TMP) BDEs are found to decrease monotonically as the size of the alkali metal cation increases. No activated dissociation was observed for alkali metal cation binding to TMP. The binding of alkali metal cations to TMP is compared with that to acetone and methanol.     

Noncovalent interactions of Cu+ with N-donor ligands (pyridine, 4,4-dipyridyl, 2,2-dipyridyl, and 1,10-phenanthroline): collision-induced dissociation and theoretical studies

Collision-induced dissociation of complexes of Cu+ bound to a variety of N-donor ligands (N-L) with Xe is studied using guided ion beam tandem mass spectrometry. The N-L ligands examined include pyridine, 4,4-dipyridyl, 2,2-dipyridyl, and 1,10-phenanthroline. In all cases, the primary and lowest-energy dissociation channel observed corresponds to the endothermic loss of a single intact N-L ligand. Sequential dissociation of additional N-L ligands is observed at elevated energies for the pyridine and 4,4-dipyridyl complexes containing more than one ligand. Ligand exchange processes to produce Cu+Xe are also observed as minor reaction pathways in several systems. The primary cross section thresholds are interpreted to yield 0 and 298 K bond dissociation energies (BDEs) after accounting for the effects of multiple ion-neutral collisions, the kinetic and internal energy distributions of the reactants, and dissociation lifetimes. Density functional theory calculations at the B3LYP/6-31G* level are performed to obtain model structures, vibrational frequencies, and rotational constants for the neutral N-L ligands and the Cu+(N-L)x complexes. The relative stabilities of the various conformations of these N-L ligands and Cu+(N-L)x complexes as well as theoretical BDEs are determined from single point energy calculations at the B3LYP/6-311+G(2d,2p) level of theory using B3LYP/6-31G* optimized geometries. Excellent agreement between theory and experiment is observed for all complexes containing one or two N-L ligands, while theory systematically underestimates the strength of binding for complexes containing more than two N-L ligands. The ground-state structures of these complexes and the trends in the sequential BDEs are explained in terms of stabilization gained from sd-hybridization and repulsive ligand-ligand interactions. The nature of the binding interactions in the Cu+(N-L)x complexes are examined via natural bond orbital analyses.     

Studies on the non-covalent complexes between oleanolic acid and cyclodextrins using electrospray ionization tandem mass spectrometry

Non-covalent inclusion complexes formed between an anti-inflammatory drug, oleanolic acid (OA), and alpha-, beta- and gamma-cyclodextrins (CDs) were investigated by means of solubility studies and electrospray ionization tandem mass spectrometry (ESI-MS(n)). The order of calculated association constants (K(1 : 1)) of complexes between OA and different CDs in solution is in good agreement with the order of their relative peak intensities and the relative CID energies of the complexes under the same ESI-MS(n) conditions. These results indicate a direct correlation between the behaviors of solution- and gas-phase complexes. ESI-MS can thus be used to evaluate solution-phase non-covalent complexes successfully. The experimental results show that the most stable 1 : 1 inclusion complexes between three CDs and OA can be formed, but 2 : 1 CD-OA complexes can be formed with beta- and gamma-CDs. Multi-component complexes of alpha-CD-OA-beta-CD (1 : 1 : 1), alpha-CD-OA-gamma-CD (1 : 1 : 1) and beta-CD-OA-gamma-CD (1 : 1 : 1) were found in equimolar CD mixtures with excess OA. The formation of 2 : 1 and multi-component 1 : 1 : 1 non-covalent CD-OA complexes indicates that beta- and gamma-CD are able to form sandwich-type inclusion non-covalent complexes with OA. The above results can be partly supported by the relative sizes of OA and CD cavities by molecular modeling calculations. All the complexes allow the detection of gaseous deprotonated CD-OA complexes in the negative ion mode at high abundances. The relative stabilities of the CDs-OA inclusion complexes in the gas phase can be evaluated from the relative CID energies in the ion trap (alpha-CD-OA < beta-CD-OA < gamma-CD-OA) in the negative ion mode.     

Binding Energies of Proton-Bound Dimers of Imidazole and n-Acetylalanine Methyl Ester Obtained by Blackbody Infrared Radiative Dissociation

The dissociation kinetics of protonated n-acetyl-L-alanine methyl ester dimer (AcAlaME(d)), imidazole dimer, and their cross dimer were measured using blackbody infrared radiative dissociation (BIRD). Master equation modeling of these data was used to extract threshold dissociation energies (E(o)) for the dimers. Values of 1.18 +/- 0.06, 1.11 +/- 0.04, and 1.12 +/- 0.08 eV were obtained for AcAlaME(d), imidazole dimer, and the cross dimer, respectively. Assuming that the reverse activation barrier for dissociation of the ion-molecule complex is negligible, the value of E(o) can be compared to the dissociation enthalpy (DeltaH(d) degrees ) from HPMS data. The E(o) values obtained for the imidazole dimer and the cross dimer are in agreement with HPMS values; the value for AcAlaME(d) is somewhat lower. Radiative rate constants used in the master equation modeling were determined using transition dipole moments calculated at the semiempirical (AM1) level for all dimers and compared to ab initio (RHF/3-21G*) calculations where possible. To reproduce the experimentally measured dissociation rates using master equation modeling, it was necessary to multiply semiempirical transition dipole moments by a factor between 2 and 3. Values for transition dipole moments from the ab initio calculations could be used for two of the dimers but appear to be too low for AcAlaME(d). These results demonstrate that BIRD, in combination with master equation modeling, can be used to determine threshold dissociation energies for intermediate size ions that are in neither the truncated Boltzmann nor the rapid energy exchange limit.     

Retention of bioactive ligand conformation in a gaseous protein-trisaccharide complex

Arrhenius parameters, obtained with the blackbody infrared radiative dissociation technique, are reported for the dissociation of gaseous protonated complexes of a single-chain variable fragment (scFv) of the monoclonal antibody Se155-4 with structurally constrained trisaccharide ligands that resemble the bioactive conformer. The similarity in the dissociation activation energies measured for the +10 charge-state complexes of the constrained ligands and the native trisaccharide is evidence that the bioactive conformation of the native ligand is retained in the gas phase.     

The effect of the secondary structure on dissociation of peptide radical cations: fragmentation of angiotensin III and its analogues

Fragmentation of protonated RVYIHPF and RVYIHPF-OMe and the corresponding radical cations was studied using time- and collision energy-resolved surface-induced dissociation (SID) in a Fourier transform ion cyclotron resonance mass spectrometer (FT-ICR MS) specially equipped to perform SID experiments. Peptide radical cations were produced by gas-phase fragmentation of Co (III)(salen)-peptide complexes. Both the energetics and the mechanisms of dissociation of even-electron and odd-electron angiotensin III ions are quite different. Protonated molecules are much more stable toward fragmentation than the corresponding radical cations. RRKM modeling of the experimental data suggests that this stability is largely attributed to differences in threshold energies for dissociation, while activation entropies are very similar. Detailed analysis of the experimental data obtained for radical cations demonstrated the presence of two distinct structures separated by a high free-energy barrier. The two families of structures were ascribed to the canonical and zwitterionic forms of the radical cations produced in our experiments.     

Noncovalent interactions of Zn+ with N-donor ligands (pyridine, 4,4'-dipyridyl, 2,2'-dipyridyl, and 1,10-phenanthroline): collision-induced dissociation and theoretical studies

The binding interactions in complexes of Zn(+) with nitrogen donor ligands, (N-L) = pyridine (x = 1-4), 4,4'-dipyridyl (x = 1-3), 2,2'-dipyridyl (x = 1-2), and 1,10-phenanthroline (x = 1-2), are examined in detail. The bond dissociation energies (BDEs) for loss of an intact ligand from the Zn(+)(N-L)(x) complexes are reported. Experimental BDEs are obtained from thermochemical analyses of the threshold regions of the collision-induced dissociation cross sections of Zn(+)(N-L)(x) complexes. Density functional theory calculations at the B3LYP/6-31G* level of theory are performed to determine stable structures of these species and to provide molecular parameters needed for the thermochemical analysis of experimental data. Relative stabilities of the various conformations of these N-donor ligands and their complexes to Zn(+) as well as theoretical BDEs are determined from single point energy calculations at the B3LYP/6-311+G(2d,2p) and M06/6-311+G(2d,2p) levels of theory using the B3LYP/6-31G* optimized geometries. The experimental BDEs for the Zn(+)(N-L)(x) complexes are in reasonably good agreement with values derived from density functional theory calculations. BDEs derived from M06 calculations provide better agreement with the measured values than those based on B3LYP calculations. Trends in the sequential BDEs are explained in terms of sp polarization of Zn(+) and repulsive ligand-ligand interactions. Comparisons are made to the analogous Cu(+)(N-L)(x) and Ni(+)(N-L)(x) complexes previously studied.     

Binding energies of water to doubly hydrated cationized glutamine and structural analogues in the gas phase

The modes of metal-ion and water binding in doubly hydrated complexes of lithiated and sodiated glutamine (Gln) are probed using blackbody infrared radiative dissociation experiments and density functional theory calculations. Threshold dissociation energies, E0, for loss of a water molecule from these complexes are obtained from master-equation modeling of these data. The values of E0 are 36 +/- 1 and 38 +/- 2 kJ/mol for the lithiated and sodiated glutamine complexes, respectively, and are consistent with calculated water binding energies for the nonzwitterionic form of these complexes. Calculated water binding energies for the zwitterionic forms of these complexes are significantly higher. In contrast, calculations indicate that the zwitterionic form of Gln in these complexes is more stable than the nonzwitterionic form by 8 and 15 kJ/mol when lithiated and sodiated, respectively. Doubly hydrated lithiated and sodiated complexes of asparagine methyl ester (AsnOMe), asparagine ethyl ester (AsnOEt), and glutamine methyl ester (GlnOMe) were also studied for comparison to Gln. Although these clusters lack the acidic group of Gln and therefore have different water coordination behavior, these results further support the conclusion that Gln is nonzwitterionic in these clusters. Surprisingly, the complexes containing sodium are more stable than those containing lithium, a result that is attributed to subtle differences in how these two metal ions bind to the amino acid esters in these complexes.     

The collaborative role of molecular conformation and energetics in the binding of gas-phase non-covalent polymer/amine complexes

A large series of similar non-covalent complexes were probed using ion mobility spectrometry, molecular mechanics/molecular dynamics (MM/MD), electrospray-tandem mass spectrometry (ESI-MS/MS) and RRKM theory in order to determine the effects of charge state and charge location upon the conformation, the 0 K activation energy (E(0)) and the entropy of activation (ΔS(?)) of the dissociation of these complexes. The non-covalent complexes consisted of poly(methylmethacrylate) oligomers and singly and doubly charged diaminoalkanes of varying length. This allowed for control of the charge separation within the complexes, as well as the size of the complex. A destabilizing effect was observed in complexes containing protons in close proximity, and/or short oligomers. Interestingly, a multiple charge stabilizing effect was observed when charge sites were sufficiently separated and/or when the polymer moiety of the complex was large. ΔS(?) values of doubly charged complexes showed a greater increase with increasing polymer size in comparison to singly charged complexes. This entropic observation is explained by structure, where IMS and MM/MD determined that the charge location was the determining factor of the overall conformation of these complexes and multiple charging resulted in more rigid arrangements. Dissociation of a tightly bound complex is more entropically favorable than a loosely bound complex. Also presented is a MM/MD refinement regime derived from IMS measurements.     

The role of phosphorylated residues in peptide-peptide noncovalent complexes formation

Electrospray mass spectrometry (ESI-MS) has become the tool of choice for the study of noncovalent complexes. Our previous work has highlighted the role of phosphorylated amino acid residues in the formation of noncovalent complexes through electrostatic interaction with arginine residues’ guanidinium groups. In this study, we employ tandem mass spectrometry to investigate the gas-phase stability and dissociation pathways of these noncovalent complexes. The only difference in the three phosphopeptides tested is the nature of the phosphorylated amino acid residue. In addition the absence of acidic residues and an amidated carboxyl terminus insured that the only negative charge came from the phosphate, which allowed for the comparison of the noncovalent bond between arginine residues and each of the different phosphorylated residues. Dissociation curves were generated by plotting noncovalent complex ion intensities as a function of the nominal energy given to the noncovalent complex ion before entering the collision cell. These results showed that noncovalent complexes formed with phosphorylated tyrosine were the most stable, followed by serine and threonine, which had similar stability.     

Sequential bond energies of Fe+ (CO2)n, n = 1-5, determined by threshold collision-induced dissociation and ab initio theory

Collision-induced dissociation of the Fe+ (CO2)n complexes for n = 1-5 is studied using kinetic energy dependent guided ion beam mass spectrometry. In all cases, the primary products are endothermic loss of an intact neutral ligand from the complex. The cross section thresholds are interpreted to yield 0 K bond energies after accounting for the effects of multiple ion-molecule collisions, internal energy of the complexes, and unimolecular decay rates. These values are compared with density functional theoretical values for all five complexes. Theory provides bond energies in reasonable agreement with experiment for n = 1-4 and predictions for the infrared spectroscopy of these complexes that agree nicely with experimental results of Gregoire and Duncan (J. Chem. Phys. 2002, 117, 2120). Our thermochemical results are also compared with the Fe+ (CO)n and Fe+ (N2)n complexes, previously studied.     

Hydrogen Bonding versus Halogen Bonding: Spectroscopic Investigation of Gas-Phase Complexes Involving Bromide and Chloromethanes

Hydrogen bonding and halogen bonding are important non-covalent interactions that are known to occur in large molecular systems, such as in proteins and crystal structures. Although these interactions are important on a large scale, studying hydrogen and halogen bonding in small, gas-phase chemical species allows for the binding strengths to be determined and compared at a fundamental level. In this study, anion photoelectron spectra are presented for the gas-phase complexes involving bromide and the four chloromethanes, CH₃Cl, CH₂Cl₂, CHCl₃, and CCl₄. The stabilisation energy and electron binding energy associated with each complex are determined experimentally, and the spectra are rationalised by high-level CCSD(T) calculations to determine the non-covalent interactions binding the complexes. These calculations involve nucleophilic bromide and electrophilic bromine interactions with chloromethanes, where the binding motifs, dissociation energies and vertical detachment energies are compared in terms of hydrogen bonding and halogen bonding.     

Dipole effects on cation-pi interactions: absolute bond dissociation energies of complexes of alkali metal cations to N-methylaniline and N,N-dimethylaniline

Threshold collision-induced dissociation of M (+)( nMA) x with Xe is studied using guided ion beam mass spectrometry, where nMA = N-methylaniline and N, N-dimethylaniline and x = 1 and 2. M (+) includes the following alkali metal cations: Li (+), Na (+), K (+), Rb (+), and Cs (+). In all cases, the primary dissociation pathway corresponds to the endothermic loss of an intact nMA ligand. The primary cross section thresholds are interpreted to yield 0 and 298 K bond dissociation energies (BDEs) for ( nMA) x-1 M (+)-( nMA) after accounting for the effects of multiple ion-neutral collisions, the internal and kinetic energy distributions of the reactants, and the dissociation lifetimes. Density functional theory calculations at the B3LYP/6-31G* level of theory are used to determine the structures of these complexes, which are also used in single-point calculations at the MP2(full)/6-311+G(2d,2p) level to determine theoretical BDEs. The results of these studies are compared to previous studies of the analogous M (+)(aniline) x complexes to examine the effects of methylation of the amino group on the binding interactions. Comparisons are also made to a wide variety of cation-pi complexes previously studied to elucidate the contributions that ion-dipole, ion-induced-dipole, and ion-quadrupole interactions make to the overall binding.     

CASPT2 investigation of ethane dissociation and methyl recombination using canonical variational transition state theory

Ethane dissociation and the reverse recombination reactions were investigated based on CASPT2 and CASSCF calculations. The CASPT2 (partial) geometry optimization calculations and the CASSCF frequency calculations provided geometrical parameters, potential energies, and vibrational frequencies along the reaction pathway. For determining dissociation and recombination rate constants at a temperature range from 200 to 2000 K, two models (models 1 and 2) were used on the basis of the canonical variational transition state theory. The different methods for accounting for the five transitional modes were proposed in the two models. Dissociation activation parameters evaluated using the two models are in good agreement with the data available in the literature. Model 1 predicts reasonable rate constant values for methyl recombination at high temperature, and model 2 predicts reasonable values in both high‐ and low‐temperature ranges. © 2008 Wiley Periodicals, Inc. Int J Chem Kinet 40: 161–173, 2008