Search

The first steps of hydration of the protonated aromatic amino acids phenylalanine, tryptophan, and tyrosine were studied experimentally employing a mass spectrometer equipped with a drift cell to examine the sequential addition of individual water molecules in equilibrium experiments and theoretically by a combination of molecular mechanics and electronic structure calculations (B3LYP/6-311++G**) on the three amino acid systems including up to five water molecules. It is found that both the ammonium and carboxyl groups offer good water binding sites with binding energies of the order of 13 kcal/mol for the first water molecule. Subsequent water molecules bind less strongly, in the range of 7-11 kcal/mol for the second through fifth water molecules. The ammonium group is able to host up to three water molecules and the carboxyl group one water molecule before additional water molecules bind either to the amino acid side chain as in tyrosine or to already-bound water in a second solvation shell around the ammonium group. Reasons for the surprisingly high water affinity of the neutral carboxyl group, comparable to that of the charge-carrying ammonium group, are found to be high intrinsic hydrophilicity, favorable charge-dipole alignment, and--for the case of multiply hydrated species--favorable dipole-dipole interaction among water molecules and the lack of alternative fully exposed hydration sites.

The sequential addition of water molecules to protonated and deprotonated forms of the four mononucleotides dAMP, dCMP, dGMP, and dTMP was studied experimentally by equilibrium measurements using an electrospray mass spectrometer equipped with a drift cell and theoretically by computational methods including molecular modeling and density functional theory calculations. Experiments were carried out in positive and negative ion mode, and calculations included the protonated and deprotonated forms of the four nucleotides. For deprotonated anionic nucleotides the experimental enthalpies of hydration (DeltaH degrees n) were found to be similar for all four systems and varied between -10.1 and -11.5 kcal mol-1 for the first water molecule (n = 1) and -8.3 and -9.6 kcal mol-1 for additional water molecules (n = 2-4). Theory indicated that the first water molecule binds to the charge-carrying phosphate group. Simulations of deprotonated mononucleotides with four water molecules yielded a large number of structures with similar energies. In some of the structures all four water molecules cluster around the phosphate group, and in other structures the four water molecules each hydrate a different functional group of the nucleotide. These include the phosphate group, the deoxyribose hydroxyl group, and various functional groups on the nucleobases. Experimental DeltaH degrees 1 values for the protonated cationic mononucleotides ranged from -10.5 to -13.5 kcal mol-1 with more negative values (< or =-12 kcal mol-1) for dCMP, dGMP, and dTMP and the least negative value for dAMP. For n = 2-4 DeltaH degrees n values varied from -6.9 to -9.7 kcal/mol and were similar in value to the deprotonated nucleotides except for dAMP. Theory on the protonated nucleotides indicated that the first water molecule binds to the charge-carrying group for dCMP, dGMP, and dTMP. For protonated dAMP, on the other hand, the charge-carrying N3 group is well self-solvated by the phosphate group and not readily available for a hydrogen bond with the water molecule. The insight gained on nucleotide stabilization by individual water molecules is used to discuss the competition between hydration of individual nucleotides and Watson-Crick base pairing.

Noncovalent peptide-peptide and peptide-water interactions in small model systems were examined using an electrospray mass spectrometer equipped with a high-pressure drift cell. The results of these aggregation and hydration experiments were interpreted with the aid of molecular mechanics (MM) and density functional theory (DFT) calculations. The systems investigated include bare deprotonated monomers and dimers [P(1,2)-H](-) and hydrated deprotonated monomers and dimers [P(1,2)-H](-).(H(2)O)(n)() for the peptides dialanine (P = AA) and diglycine (P = GG). Mass spectra indicated that both peptides AA and GG form exclusively dimer ions in the electrospray process. Monomeric ions were generated by high-energy injection of the dimers into the drift cell. Temperature-dependent hydration equilibrium experiments carried out in the drift cell yielded water binding energies ranging from 11.7 (first water molecule) to 7.1 kcal/mol (fourth water) for [AA-H](-) and 11.0 to 7.4 kcal/mol for [GG-H](-). The first water molecule adding to the dimer ions [AA-H](-).(AA) and [GG-H](-).(GG) is bound by 8.4 and 7.5 kcal/mol, respectively. The hydration mass spectra for the monomers and dimers provide a means to compare the ability of water and a neutral peptide to solvate a deprotonated peptide [P-H](-). The data indicate that a similar degree of solvation is achieved by four water molecules, [P-H](-).(H(2)O)(4), or one neutral peptide, [P-H](-).(P). Temperature-dependent kinetics experiments yielded activation energies for dissociation of the dimers [AA-H](-).(AA) and [GG-H](-).(GG) of 34.9 and 32.2 kcal/mol, respectively. MM and DFT calculations carried out for the dialanine system indicated that the dimer binding energy is 24.3 kcal/mol, when the [AA-H](-) and AA products are relaxed to their global minimum structures. However, a value of 38.9 kcal/mol is obtained if [AA-H](-) and AA dissociate but retain the structures of the moieties in the dimer, suggesting the transition state occurs early in the dissociation process. Similar results were found for the diglycine dimer.

A series of gas-phase experiments and extensive theoretical modeling was done on the family of singly protonated peptides AARAA, Ac-AARAA, and AARAA-OMe. (AARAA)H(+) underwent extensive H/D exchange with D(2)O, whereas the other two peptides with blocked termini did not, implying that a salt bridge was involved in the H/D exchange process. Ion mobility measurements and complementary molecular modeling unambiguously identified the 300 K structures of all three protonated peptides as charge solvation structures, not salt bridges. High-level density functional theory calculations indicated the global minimum of (AARAA)H(+) was a charge solvation structure with the lowest-energy salt bridge structure 4.8 kcal/mol higher in energy. Uptake of the first five water molecules of hydration at 260 K showed near identical propensities for all three peptides consistent with a common structural motif. Quantitative measurements of Delta H degrees and Delta S degrees for the first two waters of hydration were very similar for all three peptides, again suggestive of a common structure. A detailed search of the potential energy surface for the singly hydrated (AARAA)H(+) using molecular mechanics and density functional theory approaches indicated a charge solvation structure was the global minimum, but now the lowest-energy salt bridge structure was only 1.8 kcal/mol higher in energy. Importantly, a low-energy transition state connecting the charge solvation and the salt bridge structures was found where the D(2)O molecule facilitated H/D exchange via the relay mechanism. This "relay" transition state was 7 kcal/mol below the (AARAA)H(+) + D(2)O asymptotic energy, suggesting that facile H/D exchange could occur in this system. There was no equivalent low-lying relay mechanism transition state for the (Ac-AARAA)H(+) and (AARAA-OMe)H(+) peptides, consistent with the fact that H/D exchange was not observed. Hence, the combined experimental and theoretical methods confirmed that a salt bridge was involved in the H/D exchange by D(2)O of (AARAA)H(+), but it existed only as a kinetic intermediate, not as a global minimum structure. These findings suggest that caution must be observed in drawing structural conclusions from H/D exchange only. A prescription is given here for understanding both the structural and H/D exchange mechanistic aspects of bare and singly hydrated peptides.

The sequential addition of water molecules to a series of small protonated peptides was studied by equilibrium experiments using electrospray ionization combined with drift cell techniques. The experimental data were compared to theoretical structures of selected hydrated species obtained by molecular mechanics simulations. The sequential water binding energies were measured to be of the order of 7-15 kcal/mol, with the largest values for the first water molecule adding to either a small nonarginine containing peptide (e.g., protonated dialanine) or to a larger peptide in a high charge state (e.g., triply protonated neurotensin). General trends are (a) that the first water molecules are more strongly bound than the following water molecules, (b) that very small peptides (2-3 residues) bind the first few water molecules more strongly than larger peptides, (c) that the first few water molecules bind more strongly to higher charge states than to lower charge states, and (d) that water binds less strongly to a protonated guanidino group (arginine containing peptides) than to a protonated amino group. Experimental differential entropies of hydration were found to be of the order of -20 cal/mol/K although values vary from system to system. At constant experimental conditions the number of water molecules adding to any peptide ion is strongly dependent on the peptide charge state (with higher charge states adding proportionally more water molecules) and only weakly dependent on the choice of peptide. For small peptides molecular mechanics calculations indicate that the first few water molecules add preferentially to the site of protonation until a complete solvation shell is formed around the charge. Subsequent water molecules add either to water molecules of the first solvation shell or add to charge remote functional groups of the peptide. In larger peptides, charge remote sites generally compete more effectively with charge proximate sites even for the first few water molecules.