Conformational Structure of Tyrosine, Tyrosyl-glycine, and Tyrosyl-glycyl-glycine by Double Resonance Spectroscopy

Enantioselective Complexation of Protonated Tyrosine by a Chiral Crown-ether: The Nature of the Hydrogen Bonds Makes the Difference

The complexes formed between (18-crown-6)-tetracarboxylic acid, denoted as 18C6TA, and the two enantiomers of protonated tyrosine, L- and D-Tyr, are studied in a cryogenic ion trap by combining mass spectrometry, laser spectroscopy, and DFT calculations. Both UV and IR photodissociation spectra indicate the formation of multiple isomers for each complex, some of them interconverting upon IR irradiation. Conformer-selective vibrational spectroscopy reveals that all structures involve an internally hydrogen-bonded folded structure of the crown ether. The complexes formed with L-Tyr involve two NH…O interactions with the ether oxygen atoms, accompanied by two hydrogen bonds to the crown ether carboxylic function, one from the NH group and the other from the COOH group of tyrosine. The complexes formed with D-Tyr show more conformational mobility: two out of the three lowest energy conformers observed are tripodal, with three NH…O interactions between the amino acid ammonium and the crown ether oxygens. An additional conformer shows only one NH…O interaction, but has two interactions involving one of the cavity COOH moieties, one with the NH and one with the phenol OH. The increased calculated stability of the complex made from (-) 18C6TA and L-Tyr parallels its higher abundance in the mass spectrum of an isotopically labelled racemic mixture.

Accurate Computation of Aqueous pKas of Biologically Relevant Organic Acids: Overcoming the Challenges Posed by Multiple Conformers, Tautomeric Equilibria, and Disparate Functional Groups with the Fully Black-Box pK-Yay Method

The use of static electronic structure calculations to compute solution-phase pKas offers a great advantage in that a macroscopic bulk property could be computed via microscopic computations involving very few molecules. There are various sources of errors in the quantum chemical calculations though. Overcoming these errors to accurately compute pKas of a plethora of acids is an active area of research in physical chemistry pursued by both computational as well as experimental chemists. We recently developed the pK-Yay method in our attempt to accurately compute aqueous pKas of strong and weak acids. The method is fully black-box, computationally inexpensive, and is very easy for even a nonexpert to use. However, the method was thus far tested on very few molecules (only 16 in all). Herein, in order to assess the future applicability of pK-Yay, we study the effect of multiple conformers, the presence of tautomers under equilibrium, and the impact of a wide variety of functional groups (derivatives of acetic acid with substituents at various positions, dicarboxylic acids, aromatic carboxylic acids, amines and amides, phenols and thiols, and fluorine bearing organic acids). Starting with more than 1000 conformers and tautomers, this study establishes that overall errors of ∼ 1.0 pKa units are routinely obtained for a majority of the molecules. Larger errors are noted in cases where multiple charges, intramolecular hydrogen bonding, and several ionizable functional groups are simultaneously present. An important conclusion to emerge from this work is that, the computed pKas are insensitive (difference <0.5) to whether we consider multiple conformers/tautomers or only choose the most stable conformer/tautomer. Further, pK-Yay captures the stereoelectronic effects arising due to differing axial vs equatorial pattern, and is useful to predict the dominant acid-base equilibrium in a system featuring several equilibria. Overall, pK-Yay may be employed in several chemical applications featuring organic molecules and biomonomers.

How change in chirality prevents β-amyloid type interaction in a protonated cyclic dipeptide dimer

The protonated dimers of the diketopiperazine dipeptide cyclo (LPhe-LHis) and cyclo (LPhe-DHis) are studied by laser spectroscopy combined with mass spectrometry to shed light on the influence of stereochemistry on the clustering propensity of cyclic dipeptides. The marked spectroscopic differences experimentally observed in the hydride stretch region are well accounted for by the results of DFT calculations. Both diastereomeric protonated dimers involve a strong ionic hydrogen bond from the protonated imidazole ring of one monomer to the neutral imidazole nitrogen of the other. While this strong interaction is accompanied by a single NH⋯O hydrogen bond between the amide functions of the two moieties for the protonated dimer of cyclo (LPhe-DHis), that of cyclo (LPhe-LHis) involves two NH⋯O interactions, forming the motif of an antiparallel β sheet. Therefore, a change in chirality of the residue prevents the formation of the β sheet pattern observed in the amyloid type aggregation. These results emphasize the peculiar role of the histidine residue in peptide structure and interaction.

On Aromaticity of the Aromatic α-Amino Acids and Tuning of the NICS Indices to Find the Aromaticity Order

The NICS aromaticity indices of the rings in flexible phenylalanine (Phe), tryptophan (Trp), tyrosine (Tyr), and histidine (His) chiral molecules were analyzed. These molecules have several dozens of conformers, and their rings are slightly non-planar. Therefore, the population-averaged NICS_pav index was defined, and the NICS scans had to be performed with respect to planes found by the least-squares routine. A rule differentiating an obverse and a reverse ring face in aromatic amino acids was formulated. The NICS scan minima corresponding to the obverse and reverse face were unequal, which prompted us to use the term ring face aromaticity/ring face tropicity. It appeared that for Phe, Trp, Tyr, and His, the reverse face has always had higher ring face aromaticity/ring face tropicity than the obverse one. Despite the NICS modifications, uncertainty about the amino acid aromaticity order remained. This motivated us to use the integral INICS index newly proposed by Stanger as well. Then, the following sequence was obtained: Trp(phenyl) > Phe > Trp(pyrrole) > His > Tyr. The juxtaposition of the INICS indices of amino acids with that of some model rings revealed a fair transferability of the values. Finally, analysis of the substituent effect on INICS demonstrated that the aromaticity of Tyr is the lowest due to the strength of the OH group π-electron-donating effect able to perturb enough the ring charge distribution and its magnetic aromaticity. The NICS calculations were executed using the ARONICS program written within the project.

Neutral Peptides in the Gas Phase: Conformation and Aggregation Issues

Combined IR and UV laser spectroscopic techniques in molecular beams merged with theoretical approaches have proven to be an ideal tool to elucidate intrinsic structural properties on a molecular level. It offers the possibility to analyze structural changes, in a controlled molecular environment, when successively adding aggregation partners. By this, it further makes these techniques a valuable starting point for a bottom-up approach in understanding the forces shaping larger molecular systems. This bottom-up approach was successfully applied to neutral amino acids starting around the 1990s. Ever since, experimental and theoretical methods developed further, and investigations could be extended to larger peptide systems. Against this background, the review gives an introduction to secondary structures and experimental methods as well as a summary on theoretical approaches. Vibrational frequencies being characteristic probes of molecular structure and interactions are especially addressed. Archetypal biologically relevant secondary structures investigated by molecular beam spectroscopy are described, and the influences of specific peptide residues on conformational preferences as well as the competition between secondary structures are discussed. Important influences like microsolvation or aggregation behavior are presented. Beyond the linear α-peptides, the main results of structural analysis on cyclic systems as well as on β- and γ-peptides are summarized. Overall, this contribution addresses current aspects of molecular beam spectroscopy on peptides and related species and provides molecular level insights into manifold issues of chemical and biochemical relevance.

Revealing Single-Bond Anomeric Selectivity in Carbohydrate-Protein Interactions

The noncovalent binding of proteins to glycans is amazingly selective to the isoforms of carbohydrates, including α/β anomers that coexist in solution. We isolate in the gas phase and study at the atomic level the simplest model system: noncovalent complexes of monosaccharide α/β-GalNAc and protonated aromatic molecule tyramine. IR/UV cold ion spectroscopy and quantum chemistry calculations jointly solve the structures of the two complexes. Although the onsets of the measured UV absorptions of the complexes differ significantly, the networks of H bonds in both complexes appear identical and do not include the anomeric hydroxyl. The detailed analysis reveals that, through inductive polarization, the α- to β-reorientation of this group nevertheless reduces the length of one remote short intermolecular H-bond by 0.03 Å. Although small, this change substantially strengthens the bond, thus contributing to the anomeric selectivity of the binding.

Conformational Study of the Jet-Cooled Diketopiperazine Peptide Cyclo Tyrosyl-Prolyl

The conformational landscape of the diketopiperazine (DKP) dipeptide built on tyrosine and proline, namely, cyclo Tyr-Pro, is studied by combining resonance-enhanced multiphoton ionization, double resonance infrared ultraviolet (IR-UV) spectroscopy, and quantum chemical calculations. Despite the geometrical constraints due the two aliphatic rings, DKP and proline, cyclo Tyr-Pro is a flexible molecule. For both diastereoisomers, cyclo LTyr-LPro and cyclo LTyr-DPro, two structural families coexist under supersonic jet conditions. In the most stable conformation, the aromatic tyrosine substituent is folded over the DKP ring (g⁺ geometry of the aromatic ring) as it is in the solid state. The other structure is completely extended (g^- geometry of the aromatic ring) and resembles that proposed for the vapor phase. IR-UV results are not sufficient for unambiguous assignment of the observed spectra to either folded or extended conformations and the simulation of the vibronic pattern of the S₀-S₁ transition is necessary. Still, the comparison between IR-UV results and anharmonic calculations allows explanation of the minor structural differences between cyclo LTyr-LPro and cyclo LTyr-DPro in terms of different NH···π and CH···π interactions.

N-trimethyl chitosan nanoparticles and CSKSSDYQC peptide: N-trimethyl chitosan conjugates enhance the oral bioavailability of gemcitabine to treat breast cancer

Gemcitabine is a nucleoside analogue effective against a number of cancers. However, the full potential of this drug has not been realised, in part due to low oral bioavailability and frequent dosing requirements. This study reports the synthesis, in-vitro, ex-vivo and in-vivo evaluation of trimethyl chitosan (TMC) - CSKSSDYQC (CSK) peptide conjugates capable of enhancing the oral bioavailability of gemcitabine due to the ability to target intestinal goblet cells and promote intestinal cellular uptake. TMC was synthesized by a novel two-step methylation method to improve quanternization and yield. The CSK-TMC conjugates were prepared by ionic gelation to achieve particles sized at 173.6 ± 6.8 nm, zeta potential of +18.5 ± 0.2 mV and entrapment efficiency of 66.4 ± 0.1%, capable of sustained drug release. By encapsulating gemcitabine into CSK-TMC conjugates, an increased amount of drug permeated through porcine intestinal epithelial membranes compared with the unconjugated TMC nanoparticles (NPs). The rate of cellular uptake of drug loaded conjugates into HT29-MTX-E12 intestinal goblet cells, was time- and concentration-dependant. The conjugates underwent active transport associated with adsorptive mediated, clathrin and caveolae mediated endocytosis. In cellular transport studies, drug loaded conjugates had greater drug transport capability compared with drug solution and TMC NPs over the co-cultured Caco-2/HT29-MTX-E12 cell monolayer. The drug loaded conjugates exhibited electrostatic interaction with the intestinal epithelial cells. Both P-glycoprotein (P-gp) and multiple resistance protein-2 (MRP2) efflux affected the cellular transport of the conjugates. Importantly, during the pharmacokinetic studies, the orally administrated drug loaded into TMC NPs showed an improved oral bioavailability of 54.0%, compared with gemcitabine solution of 9.9%. Notable, the CSK-TMC conjugates further improved oral bioavailability to 60.1% and reduced the tumour growth rate in a BALB/c nude mouse model, with a 5.1-fold and 3.3-fold reduction compare with the non-treated group and gemcitabine solution group. Furthermore, no major evidence of toxicity was discernible on histologic studies of selected organs. In conclusion, the presented CSK-TMC conjugates and TMC nanoparticles both significantly improve the oral bioavailability of gemcitabine and have the potential for the treatment of breast cancer.

Complexation of halide ions to tyrosine: role of non-covalent interactions evidenced by IRMPD spectroscopy

The binding motifs in the halide adducts with tyrosine ([Tyr + X]^-, X = Cl, Br, I) have been investigated and compared with the analogues with 3-nitrotyrosine (nitroTyr), a biomarker of protein nitration, in a solvent-free environment by mass-selected infrared multiple photon dissociation (IRMPD) spectroscopy over two IR frequency ranges, namely 950-1950 and 2800-3700 cm^-1. Extensive quantum chemical calculations at B3LYP, B3LYP-D3 and MP2 levels of theory have been performed using the 6-311++G(d,p) basis set to determine the geometry, relative energy and vibrational properties of likely isomers and interpret the measured spectra. A diagnostic carbonyl stretching band at ∼1720 cm^-1 from the intact carboxylic group characterizes the IRMPD spectra of both [Tyr + X]^- and [nitroTyr + X]^-, revealing that the canonical isomers (maintaining intact amino and carboxylic functions) are the prevalent structures. The spectroscopic evidence reveals the presence of multiple non-covalent forms. The halide complexes of tyrosine conform to a mixture of plane and phenol isomers. The contribution of phenol-bound isomers is sensitive to anion size, increasing from chloride to iodide, consistent with the decreasing basicity of the halide, with relative amounts depending on the relative energies of the respective structures. The stability of the most favorable phenol isomer with respect to the reference plane geometry is in fact 1.3, -2.1, -6.8 kJ mol^-1, for X = Cl, Br, I, respectively. The change in π-acidity by ring nitration also stabilizes anion-π interactions yielding ring isomers for [nitroTyr + X]^-, where the anion is placed above the face of the aromatic ring.

Stereochemistry-dependent hydrogen bonds stabilise stacked conformations in jet-cooled cyclic dipeptides: (LD) vs. (LL) cyclo tyrosine-tyrosine

Tyrosine-containing cyclic dipeptides based on a diketopiperazine (DKP) ring are studied under jet-cooled conditions using resonance-enhanced multi-photon ionisation (REMPI), conformer-selective IR-UV double resonance vibrational spectroscopy and quantum chemical calculations. The conformational landscape of the dipeptide containing natural L tyrosine (Tyr), namely c-LTyr-LTyr strongly differs from that of its diastereomer c-LTyr-DTyr. A similar family of conformers exists in both systems, with one aromatic ring folded on the dipeptide DKP ring and the other one extended. Weak NHπ and CHπ interactions are observed, which are slightly different in c-LTyr-LTyr and c-LTyr-DTyr. These structures are identical to those of LL and LD cyclo diphenylalanine, which only differ from c-Tyr-Tyr by the absence of hydroxyl on the benzene rings. While this is the only conformation observed for c-LTyr-DTyr, c-LTyr-LTyr exhibits an additional form stabilised by the interaction of the two hydroxyls, in which the two aromatic rings are in a stacked geometry. Stereochemical effects are still visible in the radical cation, for which one structure is observed for c-LTyr-DTyr, while the spectrum of the c-LTyr-LTyr radical cation is explained in terms of two co-existing structures.

Computational study on single molecular spectroscopy of tyrosin-glycine, tryptophane-glycine and glycine-tryptophane

Quantum chemistry calculations play a fundamental role in revealing the molecular structures observed in gas-phase spectroscopic measurements. The supersonic jet cooling widely used in single molecular spectroscopy experiment is a non-equilibrium process and often causes confusion on the theoretical and experimental comparison. A computational approach is proposed here to account for the effect of the non-equilibrium cooling on the experimental spectra and applied to the cases of tyrosin-glycine (YG), tryptophane-glycine (WG) and glycine-tryptophane (GW). The low energy conformers of YG, WG and GW are obtained through thorough conformational searches. The structural features and equilibrium distributions of conformations and the energy barriers for conformer conversions are then determined. Three classes of transition energy barriers, high, medium and low, are found for the conversions among conformers with distinctly different, similar and the same structural types, respectively. The final conformation populations are determined by assuming an initial temperature of about 450 K and allowing for only the conformation conversion with a low energy barrier to occur during the rapid cooling process. The results provide a natural explanation for the numbers of YG, WG and GW conformations observed experimentally. The theoretical conformation assignments are also in good agreement with the experimental IR data.

UV and IR Spectroscopy of Cryogenically Cooled, Lanthanide-Containing Ions in the Gas Phase

We measure UV and IR spectra in the gas phase for EuOH⁺, EuCl⁺, and TbO⁺ ions, which are produced by an electrospray ionization source and cooled to ∼10 K in a cold, 22-pole ion trap. The UV photodissociation (UVPD) spectra of these ions show a number of sharp, well-resolved bands in the 30000-38000 cm^-1 region, although a definite assignment of the spectra is difficult because of a high degree of congestion. We also measure an IR spectrum of the EuOH⁺ ion in the 3500-3800 cm^-1 region by IR-UV double-resonance spectroscopy, which reveals an OH stretching band at 3732 cm^-1. We perform density functional theory (DFT) and time-dependent DFT (TD-DFT) calculations of these ions in order to examine the nature of the transitions. The DFT results indicate that the states of highest-spin multiplicity (octet for EuOH⁺ and EuCl⁺ and septet for TbO⁺) are substantially more stable than other states of lower-spin multiplicity. The TD-DFT calculations suggest that UV absorption of the EuOH⁺ and EuCl⁺ ions arises from Eu(4f) → Eu(5d,6p) transitions, whereas electronic transitions of the TbO⁺ ion are mainly due to the electron promotion of O(2p) → Tb(4f,6s). The UVPD results of the lanthanide-containing ions in this study suggest the possibility of using lanthanide ions as "conformation reporters" for gas-phase spectroscopy for large molecules.

Structural analyses of isolated cyclic tetrapeptides with varying amino acid residues

Structural analyses of isolated cyclic tetrapeptides with varying amino acid residues were performed by applying combined IR/UV spectroscopy in the molecular beam and DFT calculations. The intrinsic structural properties especially with regard to the influence of different amino acid residues are fundamental for optimizing their binding ability.

Structural investigations on a linear isolated depsipeptide: the importance of dispersion interactions

The first molecular beam investigations of an isolated linear depsipeptide are presented. By applying IR/UV spectroscopic methods and DFT calculations three structural arrangements are identified with the most stable structure being only stable by including dispersion interactions.

Isolated neutral peptides

This chapter examines the structural characterisation of isolated neutral amino-acids and peptides. After a presentation of the experimental and theoretical state-of-the-art in the field, a review of the major structures and shaping interactions is presented. Special focus is made on conformationally-resolved studies which enable one to go beyond simple structural characterisation; probing flexibility and excited-state photophysics are given as examples of promising future directions.

A conformation-selective IR-UV study of the dipeptides Ac-Phe-Ser-NH2 and Ac-Phe-Cys-NH2: probing the SH···O and OH···O hydrogen bond interactions

The conformational preferences of peptides are mainly controlled by the stabilizing effect of intramolecular interactions. In peptides with polar side chains, not only the backbone but also the side chain interactions determine the resulting conformations. In this paper, the conformational preferences of the capped dipeptides Ac-Phe-Ser-NH2 (FS) and Ac-Phe-Cys-NH2 (FC) are resolved under laser-desorbed jet cooling conditions using IR-UV ion dip spectroscopy and density functional theory (DFT) quantum chemistry calculations. As serine (Ser) and cysteine (Cys) only differ in an OH (Ser) or SH (Cys) moiety; this subtle alteration allows us to study the effect of the difference in hydrogen bonding for an OH and SH group in detail, and its effect on the secondary structure. IR absorption spectra are recorded in the NH stretching region (3200-3600 cm(-1)). In combination with quantum chemical calculations the spectra provide a direct view of intramolecular interactions. Here, we show that both FS as FC share a singly γ-folded backbone conformation as the most stable conformer. The hydrogen bond strength of OH···O (FS) is stronger than that of SH···O (FC), resulting in a more compact gamma turn structure. A second conformer is found for FC, showing a β turn interaction.

Can far-IR action spectroscopy combined with BOMD simulations be conformation selective?

The combination of conformation selective far-IR/UV double resonance spectroscopy with Born–Oppenheimer molecular dynamics (BOMD) simulations is presented here for the structural characterization of the Ac-Phe-Pro-NH₂ peptide in the far-infrared spectral domain, i.e. for radiation below 800 cm⁻¹.

Isomer-selective detection of hydrogen-bond vibrations in the protonated water hexamer

The properties of hydrogen ions in aqueous solution are governed by the ability of water to incorporate ions in a dynamical hydrogen bond network, characterized by a structural variability that has complicated the development of a consistent molecular level description of H(+)(aq). Isolated protonated water clusters, H(+)(H2O)n, serve as finite model systems for H(+)(aq), which are amenable to highly sensitive and selective gas phase spectroscopic techniques. Here, we isolate and assign the infrared (IR) signatures of the Zundel-type and Eigen-type isomers of H(+)(H2O)6, the smallest protonated water cluster for which both of these characteristic binding motifs coexist, down into the terahertz spectral region. We use isomer-selective double-resonance population labeling spectroscopy on messenger-tagged H(+)(H2O)6·H2 complexes from 260 to 3900 cm(-1). Ab initio molecular dynamics calculations qualitatively recover the IR spectra of the two isomers and allow attributing the increased width of IR bands associated with H-bonded moieties to anharmonicities rather than excited state lifetime broadening. Characteristic hydrogen-bond stretching bands are observed below 400 cm(-1).

Conformationally resolved spectra of acetaminophen by UV-UV hole burning and IR dip spectroscopy in the gas phase

Electronic and vibrational spectra of acetaminophen were measured by using UV-UV hole burning (HB) and IR dip spectroscopy. HB spectra show the coexistence of 4 different species, which include two new ones. Low-frequency transitions in the spectra are reproduced by a one-dimensional periodic potential with a free-rotor basis set for the methyl group. From the analysis, we concluded that acetaminophen has two conformers and each conformer gives two independent transitions starting from the most stable 0a(1) and the hot 1e internal rotational levels. It is also found that the HB spectrum of the trans-conformer in the previous report is that from the 1e excited level, while the HB spectrum of the cis-conformer is contaminated by the transitions of the trans-conformer. Potential curves of the methyl rotational motion are determined both in S(0) and S(1).

Isomer-Specific IR-IR Double Resonance Spectroscopy of D2-Tagged Protonated Dipeptides Prepared in a Cryogenic Ion Trap

Isomer-specific vibrational predissociation spectra are reported for the gas-phase GlySarH(+) and SarSarH(+) [Gly = glycine; Sar = sarcosine] ions prepared by electrospray ionization and tagged with weakly bound D2 adducts using a cryogenic ion trap. The contributions of individual isomers to the overlapping vibrational band patterns are completely isolated using a pump-probe photochemical hole-burning scheme involving two tunable infrared lasers and two stages of mass selection (hence IR(2)MS(2)). These patterns are then assigned by comparison with harmonic (MP2/6-311+G(d,p)) spectra for various possible conformers. Both systems occur in two conformations based on cis and trans configurations with respect to the amide bond. In addition to the usual single intramolecular hydrogen bond motif between the protonated amine and the nearby amide oxygen atom, cis-SarSarH(+) adopts a previous unreported conformation in which both amino NH's act as H-bond donors. The correlated red shifts in the NH donor and C═O acceptor components of the NH···O═C linkage to the acid group are unambiguously assigned in the double H-bonded conformer.

On the perturbation of the H-bonding interaction in ethylene glycol clusters upon hydration

Ab initio and density functional methods have been employed to study the structure, stability, and spectral properties of various ethylene glycol (EG(m)) and ethylene glycol-water (EG(m)W(n)) (m = 1-3, n = 1-4) clusters. The effective fragment potential (EFP) approach was used to explore various possible EG(m)W(n) clusters. Calculated interaction energies of EG(m)W(n) clusters confirm that the hydrogen-bonding interaction between EG molecules is perturbed by the presence of water molecules and vice versa. Further, energy decomposition analysis shows that both electrostatic and polarization interactions predominantly contribute to the stability of these clusters. It was found from the same analysis that ethylene glycol-water interaction is predominant over the ethylene glycol-ethylene glycol and water-water interactions. Overall, the results clearly illustrate that the presence of water disrupts the ethylene glycol-ethylene glycol hydrogen bonds.

Characterizing the intramolecular H-bond and secondary structure in methylated GlyGlyH+ with H2 predissociation spectroscopy

We report vibrational predissociation spectra of the four protonated dipeptides derived from glycine and sarcosine, GlyGlyH(+)•(H(2))(1,2), GlySarH(+)•(D(2))(2), SarGlyH(+)•(H(2))(2), and SarSarH(+)•(D(2))(2), generated in a cryogenic ion trap. Sharp bands were recovered by monitoring photoevaporation of the weakly bound H(2) (D(2)) molecules in a linear action regime throughout the 700-4200 cm(-1) range using a table-top laser system. The spectral patterns were analyzed in the context of the low energy structures obtained from electronic structure calculations. These results indicate that all four species are protonated on the N-terminus, and feature an intramolecular H-bond involving the amino group. The large blue-shift in the H-bonded N-H fundamental upon incorporation of a methyl group at the N-terminus indicates that this modification significantly lowers the strength of the intramolecular H-bond. Methylation at the amide nitrogen, on the other hand, induces a significant rotation (~110°) about the peptide backbone.