Investigation of Noncovalent Interactions in Deprotonated Peptides: Structural and Energetic Competition between Aggregation and Hydration

Cold Denaturation of Proteins in the Absence of Solvent: Implications for Protein Storage

The effect of temperature on the stability of proteins is well explored above 298 K, but harder to track experimentally below 273 K. Variable-temperature ion mobility mass spectrometry (VT IM-MS) allows us to measure the structure of molecules at sub-ambient temperatures. Here we monitor conformational changes that occur to two isotypes of monoclonal antibodies (mAbs) on cooling by measuring their collision cross sections (CCS) at discrete drift gas temperatures from 295 to 160 K. The CCS at 250 K is larger than predicted from collisional theory and experimental data at 295 K. This restructure is attributed to change in the strength of stabilizing intermolecular interactions. Below 250 K the CCS of the mAbs increases in line with prediction implying no rearrangement. Comparing data from isotypes suggest disulfide bridging influences thermal structural rearrangement. These findings indicate that in vacuo deep-freezing minimizes denaturation and maintains the native fold and VT IM-MS measurements at sub ambient temperatures provide new insights to the phenomenon of cold denaturation.

Cold Denaturation of Proteins in the Absence of Solvent: Implications for Protein Storage

The effect of temperature on the stability of proteins is well explored above 298 K, but harder to track experimentally below 273 K. Variable-temperature ion mobility mass spectrometry (VT IM-MS) allows us to measure the structure of molecules at sub-ambient temperatures. Here we monitor conformational changes that occur to two isotypes of monoclonal antibodies (mAbs) on cooling by measuring their collision cross sections (CCS) at discrete drift gas temperatures from 295 to 160 K. The CCS at 250 K is larger than predicted from collisional theory and experimental data at 295 K. This restructure is attributed to change in the strength of stabilizing intermolecular interactions. Below 250 K the CCS of the mAbs increases in line with prediction implying no rearrangement. Comparing data from isotypes suggest disulfide bridging influences thermal structural rearrangement. These findings indicate that in vacuo deep-freezing minimizes denaturation and maintains the native fold and VT IM-MS measurements at sub ambient temperatures provide new insights to the phenomenon of cold denaturation.

Deamidation reaction network mapping of pharmacologic and related proteins: impact of solvation dielectric on the degradation energetics of asparagine dipeptides

Asn-X deamidation pathways in the FV region of the monoclonal antibody (mAb).

Assisted dipeptide bond formation: glycine as a case study

Peptide bond formation is a crucial chemical process that dominates most biological mechanisms and is claimed to be a governing factor in the origin of life. Dipeptides made from glycine are studied computationally via Density Functional Theory (DFT) using two different basis sets. This reaction was investigated from both a thermodynamic and kinetic point of view. The effect of explicit assistance via the introduction of discrete solvent molecules was investigated. Water, methanol, and cyclohexane were all employed as solvent media in addition to gas to investigate their effects on the mechanism of peptide bond formation. This computational investigation revealed that methanol is slightly better than water to leverage peptide bond formation both kinetically and thermodynamically, while cyclohexane, a non-polar and non-protic solvent, is the least effective after gas as a medium of solvation. Energetic results in the gas environment are very close to those obtained in polar and protic solvents, suggesting that peptide bonds can be formed under interstellar conditions.

Probing Gas-Phase-Clustering Thermodynamics with Ion Mobility-Mass Spectrometry: Association Energies of Phenylalanine Ions with Gas-Phase Alcohols

Vapor assisted mobility shift measurements were made with atmospheric pressure drift-tube ion mobility-mass spectrometry (IM-MS) to determine the thermodynamic properties of weakly bound ion-molecule clusters formed from protonated phenylalanine and neutral vapor molecules with hydroxyl functional groups. Relative binding energies and gas-phase association energies of amino acid ions clustered with small organic molecules have been established previously using high-pressure mass spectrometry. However, the issue of volatility largely prohibits the use of high-pressure mass spectrometry for the determination of gas-phase associations of amino acid ions clustered with neutral vapor molecules in many instances. In contrast, ion mobility measurements can be made at atmospheric pressure with volatile vapor additives near and above their boiling points, providing access to clustering equilibria not possible using high-vacuum techniques. In this study, we report the gas-phase association energies, enthalpies, and entropies for a protonated phenylalanine ion clustered with three neutral vapor molecules: 2-propanol, 1-butanol, and 2-pentanol based upon measurements at temperatures ranging from 120 to 180 °C. The gas-phase enthalpy and entropy changes ranged between -4 to -7 kcal/mol and -3 to 6 cal/(mol K), respectively. We found enthalpically favored ion-neutral cluster reactions for phenylalanine with entropic barriers for the formation of phenylalanine-1-butanol and phenylalanine-2-pentanol cluster ions, while phenylalanine-2-propanol cluster ion formation is both enthalpically and (weakly) entropically favorable. Under the measurement conditions examined, phenylalanine-vapor modifier cluster ion formation is clearly observed via shifts in the drift time for the three test vapor molecules. In comparison, negligible shifts in mobility are observed for protonated arginine exposed to the same vapor modifiers.

Water proton NMR detection of amide hydrolysis and diglycine dimerization

The transverse relaxation rate of water protons R2(1H2O) is found to be sensitive to amide hydrolysis and diglycine dimerization. The results demonstrate the feasibility of using R2(1H2O) as a diagnostic tool to detect chemical changes in aqueous solutions. Potential applications include drug product formulation and inspection.

Sequential water molecule binding enthalpies for aqueous nanodrops containing a mono-, di- or trivalent ion and between 20 and 500 water molecules

Sequential water molecule binding enthalpies, ΔH_n,n-1, are important for a detailed understanding of competitive interactions between ions, water and solute molecules, and how these interactions affect physical properties of ion-containing nanodrops that are important in aerosol chemistry. Water molecule binding enthalpies have been measured for small clusters of many different ions, but these values for ion-containing nanodrops containing more than 20 water molecules are scarce. Here, ΔH_n,n-1 values are deduced from high-precision ultraviolet photodissociation (UVPD) measurements as a function of ion identity, charge state and cluster size between 20-500 water molecules and for ions with +1, +2 and +3 charges. The ΔH_n,n-1 values are obtained from the number of water molecules lost upon photoexcitation at a known wavelength, and modeling of the release of energy into the translational, rotational and vibrational motions of the products. The ΔH_n,n-1 values range from 36.82 to 50.21 kJ mol^-1. For clusters containing more than ∼250 water molecules, the binding enthalpies are between the bulk heat of vaporization (44.8 kJ mol^-1) and the sublimation enthalpy of bulk ice (51.0 kJ mol^-1). These values depend on ion charge state for clusters with fewer than 150 water molecules, but there is a negligible dependence at larger size. There is a minimum in the ΔH_n,n-1 values that depends on the cluster size and ion charge state, which can be attributed to the competing effects of ion solvation and surface energy. The experimental ΔH_n,n-1 values can be fit to the Thomson liquid drop model (TLDM) using bulk ice parameters. By optimizing the surface tension and temperature change of the logarithmic partial pressure for the TLDM, the experimental sequential water molecule binding enthalpies can be fit with an accuracy of ±3.3 kJ mol^-1 over the entire range of cluster sizes.

Microsolvation of 2-azetidinone: a model for the peptide group–water interactions

The geometries of the 2-azetidinone–(H₂O)_n clusters, determined by rotational spectroscopy, show the preference of water to interact with the CO group and the effects of cooperative hydrogen bonding.

Host-guest chemistry in the gas phase: complex formation of cucurbit[6]uril with proton-bound water dimer

The hydration of cucurbit[6]uril (CB[6]) in the gas phase is investigated using electrospray ionization traveling wave ion mobility mass spectrometry (ESI-TWIM-MS). Highly abundant dihydrated and tetrahydrated species of diprotonated CB[6] are found in the ESI-TWIM-MS spectrum. The hydration patterns of the CB[6] ion and the dissociation patterns of the hydrated CB[6] ion indicate that two water molecules are bound to each other, forming a water dimer in the CB[6] complex. Ion mobility studies combined with the structures calculated by density functional theory suggest that the proton-bound water dimer is present as a Zundel-like structure in the CB[6] portal, forming a hydrogen bond network with carbonyl groups of the CB[6]. When a large guest molecule is bound to a CB[6] portal, water molecules cannot bind to the portal. In addition, the strong binding energy of the water dimer blocks the portal, hindering the insertion of the long alkyl chain of the guest molecule into the CB[6] cavity. With small alkali metal cations, such as Li(+) and Na(+), a single water molecule interacts with the CB[6] portal, forming hydrogen bonds with the carbonyl groups of CB[6]. A highly stable Zundel-like structure of the proton-bound water dimer or a metal-bound water molecule at the CB[6] portal is suggested as an initial hydration process for CB[6], which is only dissolved in aqueous solution with acid or alkali metal ions.

Ions in size-selected aqueous nanodrops: sequential water molecule binding energies and effects of water on ion fluorescence

The effects of water on ion fluorescence were investigated, and average sequential water molecule binding energies to hydrated ions, M(z)(H(2)O)(n), at large cluster size were measured using ion nanocalorimetry. Upon 248-nm excitation, nanodrops with ~25 or more water molecules that contain either rhodamine 590(+), rhodamine 640(+), or Ce(3+) emit a photon with average energies of approximately 548, 590, and 348 nm, respectively. These values are very close to the emission maxima of the corresponding ions in solution, indicating that the photophysical properties of these ions in the nanodrops approach those of the fully hydrated ions at relatively small cluster size. As occurs in solution, these ions in nanodrops with 8 or more water molecules fluoresce with a quantum yield of ~1. Ce(3+) containing nanodrops that also contain OH(-) fluoresce, whereas those with NO(3)(-) do not. This indirect fluorescence detection method has the advantages of high sensitivity, and both the size of the nanodrops as well as their constituents can be carefully controlled. For ions that do not fluoresce in solution, such as protonated tryptophan, full internal conversion of the absorbed 248-nm photon occurs, and the average sequential water molecule binding energies to the hydrated ions can be accurately obtained at large cluster sizes. The average sequential water molecule binding energies for TrpH(+)(H(2)O)(n) and a doubly protonated tripeptide, [KYK + 2H](2+)(H(2)O)(n), approach asymptotic values of ~9.3 (n ≥ 11) and ~10.0 kcal/mol (n ≥ 25), respectively, consistent with a liquidlike structure of water in these nanodrops.

Effect of cation complexation on the structure of a conformationally flexible multiply charged anion: stabilization of excess charge in the Na+ x adenosine 5'-triphosphate dianion ion-pair complex

We report a computational study of the conformationally and tautomerically flexible cation-dianion complex of Na(+) with doubly deprotonated adenosine 5'-triphosphate (ATP) using a hierarchical selection method. The method uses molecular dynamics to generate initial conformeric structures, followed by a classification process that groups conformers into five "families" to ensure that a representative sample of structures is retained for further analysis, while very similar conformational structures are eliminated. Hierarchical ab initio calculations (DFT and MP2) of typical conformers of the families are then performed to identify the lowest-energy conformeric structures. The procedure described should provide a useful methodology for conducting higher-level ab initio calculations of medium-sized gas-phase biological molecules for interpreting contemporary laser spectroscopy measurements. For Na(+) x [ATP-2H](2) (considering tautomers where the phosphate chain of ATP is doubly deprotonated), the calculations reveal that the sodium cation interacts directly with the negatively charged phosphates (maximum distance = 2.54 A) in all of the low-energy conformers, while a number of the structures also display close cation-adenine interactions producing compact ball-like structures. These compact structures generally correspond to the lowest-energy conformers. The structural variation between the bare [ATP-2H](2-) molecular ion (Burke et al. J. Phys. Chem. A 2005 , 109 , 9775-9785) and the Na(+) x [ATP-2H](2-) cluster is discussed in detail, including the effect of sodiation on the intramolecular hydrogen-bonding network within ATP in a gas-phase environment.

Ab initio investigation of the hydration of deprotonated amino acids

The complexation of five deprotonated anionic amino acids (glycine, L-alanine, L-valine, L-Aspartic acid, and L-glutamine) with one water molecule, has been investigated using a MP2/63-11++G(d,p) approach fully accounting for the basis set superposition errors. For each amino acid, several energetic minima have been identified, and we provide spectroscopic information allowing to discriminate them. Our results strongly suggest that two complexes should coexist under the experimental conditions for [Ala - H](-), [Val - H](-), and [Asp - H](-). Comparisons with the experimental enthalpies, entropies, and Gibbs free energies recently obtained by Wincel [J. Am. Soc. Mass Spectrom. 2008, 19, 1091-1097] show that our simulation reproduces the most significant structure/energy experimental trends, though the entropic changes induced by hydration are slightly overestimated.

Electron capture by a hydrated gaseous peptide: effects of water on fragmentation and molecular survival

The effects of water on electron capture dissociation products, molecular survival, and recombination energy are investigated for diprotonated Lys-Tyr-Lys solvated by between zero and 25 water molecules. For peptide ions with between 12 and 25 water molecules attached, electron capture results in a narrow distribution of product ions corresponding to primarily the loss of 10-12 water molecules from the reduced precursor. From these data, the recombination energy (RE) is determined to be equal to the energy that is lost by evaporating on average 10.7 water molecules, or 4.3 eV. Because water stabilizes ions, this value is a lower limit to the RE of the unsolvated ion, but it indicates that the majority of the available RE is deposited into internal modes of the peptide ion. Plotting the fragment ion abundances for ions formed from precursors with fewer than 11 water molecules as a function of hydration extent results in an energy resolved breakdown curve from which the appearance energies of the b 2 (+), y 2 (+), z 2 (+*), c 2 (+), and (KYK + H) (+) fragment ions formed from this peptide ion can be obtained; these values are 78, 88, 42, 11, and 9 kcal/mol, respectively. The propensity for H atom loss and ammonia loss from the precursor changes dramatically with the extent of hydration, and this change in reactivity can be directly attributed to a "caging" effect by the water molecules. These are the first experimental measurements of the RE and appearance energies of fragment ions due to electron capture dissociation of a multiply charged peptide. This novel ion nanocalorimetry technique can be applied more generally to other exothermic reactions that are not readily accessible to investigation by more conventional thermochemical methods.

Hydration energies of deprotonated amino acids from gas phase equilibria measurements

Singly hydrated clusters of deprotonated amino acids were studied using an electrospray high-pressure mass spectrometer equipped with a pulsed ion-beam reaction chamber. Thermochemical data, DeltaH(o), DeltaS(o), and DeltaG(o), for the hydration reaction [AA - H](-) + H(2)O = [AA - H](-).(H(2)O) were obtained from gas-phase equilibria determinations for AA = Gly, Ala, Val, Pro, Phe, Lys, Met, Trp, Gln, Arg, and Asp. The hydration free-energy changes are found to depend significantly on the side-chain substituents. The water binding energy in [AA - H](-).(H(2)O) increases with the gas-phase acidity of AA. The anionic hydrogen bond strengths in [AA - H](-).(H(2)O) are compared with those of the cationic bonds in the corresponding AAH(+).(H(2)O) systems.

Stabilization of excess charge in isolated adenosine 5'-triphosphate and adenosine 5'-diphosphate multiply and singly charged anions

Multiply charged anions (MCAs) represent highly energetic species in the gas phase but can be stabilized through formation of molecular clusters with solvent molecules or counterions. We explore the intramolecular stabilization of excess negative charge in gas-phase MCAs by probing the intrinsic stability of the [adenosine 5'-triphosphate-2H](2-) ([ATP-2H](2-)), [adenosine 5'-diphosphate-2H](2-) ([ADP-2H](2-)), and H(3)P(3)O(10)(2-) dianions and their protonated monoanionic analogues. The relative activation barriers for decay of the dianions via electron detachment or ionic fragmentation are investigated using resonance excitation of ions isolated within a quadrupole trap. All of the dianions decayed via ionic fragmentation demonstrating that the repulsive Coulomb barriers (RCB) for ionic fragmentation lie below the RCBs for electron detachment. Both the electrospray ionization mass spectra (ESI-MS) and total fragmentation energies for [ATP-2H](2-), [ADP-2H](2-), and H(3)P(3)O(10)(2-) indicate that the multiply charged H(3)P(3)O(10)(2-) phosphate moiety is stabilized by the presence of the adenosine group and the stability of the dianions increases in the order H(3)P(3)O(10)(2-) < [ADP-2H](2-) < [ATP-2H](2-). Fully optimized, B3LYP/6-31+G* minimum energy structures illustrate that the excess charges in all of the phosphate anions are stabilized by intramolecular hydrogen bonding either within the phosphate chain or between the phosphate and the adenosine. We develop a model to illustrate that the relative magnitudes of the RCBs and hence the stability of these ions is dominated by the extent of intramolecular hydrogen bonding.

Intermolecular Interactions in Biomolecular Systems Examined by Mass Spectrometry

With the development of electrospray and matrix-assisted laser desorption ionization, mass spectrometry (MS) evolved into a powerful tool in the field of biochemistry. Whereas MS is primarily analytical in nature, an increasing number of MS research groups employ the method to address fundamental biochemical questions. Probing the interaction of noncovalently bound molecules in the mass spectrometer is one of the most interesting MS-based experiments possible today, with the potential of making a significant contribution to the basic understanding of the structure and function of biochemical complexes. Here we review a number of current research efforts employing primarily MS techniques to investigate intermolecular interactions in biochemical systems. Examples chosen include the interaction of biomolecules with solvent molecules; interactions between nucleic-acid molecules, in particular, interactions in duplex and quadruplex structures; and interactions between proteins involved in neurodegenerative diseases. Finally we conclude by presenting a few examples of very large biomolecular assemblies in the mega-Dalton range analyzed by MS.

C alpha-C backbone fragmentation dominates in electron detachment dissociation of gas-phase polypeptide polyanions

Fragmentation of peptide polyanions by electron detachment dissociation (EDD) has been induced by electron irradiation of deprotonated polypeptides [M-nH](n-) with >10 eV electrons. EDD has been found to lead preferentially to a* and x fragment ions (C(alpha)-C backbone cleavage) arising from the dissociation of oxidized radical anions [M-nH]((n-1)-*. We demonstrate that C(alpha)-C cleavages, which are otherwise rarely observed in tandem mass spectrometry, can account for most of the backbone fragmentation, with even-electron x fragments dominating over radical a* ions. Ab initio calculations at the B3 LYP level of theory with the 6-311+G(2 p,2 d)//6-31+G(d,p) basis set suggested a unidirectional mechanism for EDD (cleavage always N-terminal to the radical site), with a*, x formation being favored over a, x* fragmentation by 74.2 kJ mol(-1). Thus, backbone C(alpha)-C bonds N-terminal to proline residues should be immune to EDD, in agreement with the observations. EDD may find application in mass spectrometry for such tasks as peptide sequencing and localization of labile post-translational modifications, for example, those introduced by sulfation and phosphorylation. EDD can now be performed not only in Fourier transform mass spectrometry, but also in far more widely used quadrupole (Paul) ion traps.

Multidimensional separations of ubiquitin conformers in the gas phase: relating ion cross sections to H/D exchange measurements

Investigating gas-phase structures of protein ions can lead to an improved understanding of intramolecular forces that play an important role in protein folding. Both hydrogen/deuterium (H/D) exchange and ion mobility spectrometry provide insight into the structures and stabilities of different gas-phase conformers, but how best to relate the results from these two methods has been hotly debated. Here, high-field asymmetric waveform ion mobility spectrometry (FAIMS) is combined with Fourier-transform ion cyclotron resonance mass spectrometry (FT/ICR MS) and is used to directly relate ubiquitin ion cross sections and H/D exchange extents. Multiple conformers can be identified using both methods. For the 9+ charge state of ubiquitin, two conformers (or unresolved populations of conformers) that have cross sections differing by 10% are resolved by FAIMS, but only one conformer is apparent using H/D exchange at short times. For the 12+ charge state, two conformers (or conformer populations) have cross sections differing by <1%, yet H/D exchange of these conformers differ significantly (6 versus 25 exchanges). These and other results show that ubiquitin ion collisional cross sections and H/D exchange distributions are not strongly correlated and that factors other than surface accessibility appear to play a significant role in determining rates and extents of H/D exchange. Conformers that are not resolved by one method could be resolved by the other, indicating that these two methods are highly complementary and that more conformations can be resolved with this combination of methods than by either method alone.