Isolated neutral peptides

Experimental/Computational Study on the Impact of Fluorine on the Structure and Noncovalent Interactions in the Monohydrated Cluster of ortho-Fluorinated 2-Phenylethylamine

Fluorine-containing medicinal compounds frequently allow modulation of physical-chemical properties. Here, we address the effect of fluorine, near the ethylamino side chain, on conformational flexibility and noncovalent interactions (NCIs) of the selected jet-cooled monohydrated cluster of 2-(2-fluoro-phenyl)-ethylamine (2-FPEA) by mass-selected resonance-enhanced two-photon ionization and ionization-loss stimulated Raman spectroscopies. Our results show that Raman spectral signatures of the 2-FPEA-H₂O cluster match the scaled harmonic vibrational Raman frequencies, resulting from density functional theory calculations of the most stable 2-FPEA gauche conformer hydrogen-bonded (HB) to water, confirming the three-dimensional cluster structure. This predicted electronic structure, together with NCI analysis, allows visualization and assessment of the attractive and repulsive interactions. The comparison of the NCIs and revealed red (O-H and N-H stretches) and blue shifts (C-H stretches and CH₂ out-of-plane bends) of the cluster to other class members confirm O-H···N, N-H···π, C-H···O, and C-H···F HB formation and their contribution to structure stabilization, uncovering the potential of the approach.

Revealing the Structure and Noncovalent Interactions of Isolated Molecules by Laser-Desorption/Ionization-Loss Stimulated Raman Spectroscopy and Quantum Calculations

The structural and dynamical characteristics of isolated molecules are essential, yet obtaining this information is difficult. We demonstrate laser-desorption jet-cooling/ionization-loss stimulated Raman spectroscopy to obtain Raman spectral signatures of nonvolatile molecules in the gas phase. The vibrational features of a test substance, the most abundant conformer of tryptamine, are compared and found to match those resulting from the scaled harmonic Raman spectrum obtained by density functional theory calculations. The vibrational signatures serve to identify the most prominent gauche conformer and evaluate its predicted electronic structure. These findings, together with noncovalent interaction (NCI) analysis, provide new insights into electron densities and reduced density gradients, assessing the hydrogen bonds (N-H···π and C-H···H-C) and interplay between attractive and repulsive NCIs affecting the structure. This approach accesses vibrational signatures of isolated nonvolatile molecules by tabletop lasers at uniform resolution and in a broad frequency range, promising great benefit to future studies.

Some opinions on MD-based vibrational spectroscopy of gas phase molecules and their assembly: An overview of what has been achieved and where to go

We hereby review molecular dynamics simulations for anharmonic gas phase spectroscopy and provide some of our opinions of where the field is heading. With these new directions, the theoretical IR/Raman spectroscopy of large (bio)-molecular systems will be more easily achievable over longer time-scale MD trajectories for an increase in accuracy of the MD-IR and MD-Raman calculated spectra. With the new directions presented here, the high throughput 'decoding' of experimental IR/Raman spectra into 3D-structures should thus be possible, hence advancing e.g. the field of MS-IR for structural characterization by spectroscopy. We also review the assignment of vibrational spectra in terms of anharmonic molecular modes from the MD trajectories, and especially introduce our recent developments based on Graph Theory algorithms. Graph Theory algorithmic is also introduced in this review for the identification of the molecular 3D-structures sampled over MD trajectories.

A theoretical and experimental case study of the hydrogen bonding predilection of S-methylcysteine

S-containing amino acids can lead to two types of local NH···S interactions which bridge backbone NH sites to the side chain to form either intra- or inter-residue H-bonds. The present work reports on the conformational preferences of S-methyl-L-cysteine, Cys(Me), using a variety of investigating tools, ranging from quantum chemistry simulations, gas-phase UV and IR laser spectroscopy, and solution state IR and NMR spectroscopies, on model compounds comprising one or two Cys(Me) residues. We demonstrate that in gas phase and in low polarity solution, the C- and N-capped model compound for one Cys(Me) residue adopts a preferred C5-C6γ conformation which combines an intra-residue N-H···O=C backbone interaction (C5) and an inter-residue N-H···S interaction implicating the side-chain sulfur atom (C6γ). In contrast, the dominant conformation of the C- and N-capped model compound featuring two consecutive Cys(Me) residues is a regular type I β-turn. This structure is incompatible with concomitant C6γ interactions, which are no longer in evidence. Instead, C5γ interactions occur, that are fully consistent with the turn geometry and additionally stabilize the structure. Comparison with the thietane amino acid Attc, which exhibits a rigid cyclic side chain, pinpoints the significance of side chain flexibility for the specific conformational behavior of Cys(Me).

An intraresidue H-bonding motif in selenocysteine and cysteine, revealed by gas phase laser spectroscopy and quantum chemistry calculations

Models of protein chains containing a seleno-cysteine (Sec) residue have been investigated by gas phase laser spectroscopy in order to document the effect of the H-bonding properties of the SeH group in the folding of the Sec side chain, by comparison with recent data on Ser- and Cys-containing sequences. Experimental data, complemented by quantum chemistry calculations and natural bonding orbital (NBO) analyses, are interpreted in terms of the formation of a so-called 5γ intra-residue motif, which bridges the acceptor chalcogen atom of the side chain to the NH bond of the same residue. This local structure, in which the O/S/Se atom is close to the plane of the N-terminal side amide, is constrained by local backbone-side chain hyperconjugation effects involving the S and Se atoms. Theoretical investigations of the Cys/Sec side chain show that (i) this 5γ motif is an intrinsic feature of these residues, (ii) the corresponding H-bond is strongly non-linear and intrinsically weak, (iii) but enhanced by γ- and β-turn secondary structures, which promote a more favorable 5γ H-bonding approach and distance. The resulting H-bonds are slightly stronger in selenocysteine than in cysteine, but nearly inexistent in serine, whose side chain in contrast behaves as a H-bonding donor. The modest spectral shifts of the Cys/Sec NH stretches measured experimentally reflect the moderate strength of the 5γ H-bonding, in agreement with the correlation obtained with a NBO-based H-bond strength indicator. The evolution along the Ser, Cys and Sec series emphasizes the compromise between the several factors that control the H-bonding in a hyperconjugation-constrained geometry, among them the chalcogen van der Waals and covalent radii. It also illustrates the 5γ H-bond enhancements with the Sec and Cys residues favoured by the constraints imposed by the γ- and β-turn structures of the peptide chain.

Conformational analysis by UV spectroscopy: the decisive contribution of environment-induced electronic Stark effects

UV chromophores are frequently used as probes of the molecular structure. In particular, they are sensitive to the electric field generated by the molecular environment, resulting in the observation of Stark effects on UV spectra. While these environment-induced electronic Stark effects (EI-ESE) are already used for conformational analysis in the condensed phase, this work explores the potential of such an approach when performed at much higher conformational resolution in the gas phase. By investigating model alkali benzylacetate and 4-phenylbutyrate ion pairs, where the electric field applied to the phenyl ring is chemically tuned by changing the nature of the alkali cation, this work demonstrates that precise conformational assignments can be proposed based on the correlation between the conformation-dependent calculated electric fields and the frequency of the electronic transitions observed in the experimental UV spectra. Remarkably, the sole analysis of Stark effects and fragmentation patterns in mass-selected UV spectra provided an accurate and complete conformational analysis, where spectral differences as small as a few cm-1 between electronic transitions were rationalized. This case study illustrates that the identification of EI-ESE together with their interpretation at the modest cost of a ground state electric field calculation qualify UV spectroscopy as a powerful tool for conformational analysis.

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.

Far-IR Absorption of Neutral Polycyclic Aromatic Hydrocarbons (PAHs): Light on the Mechanism of IR-UV Ion Dip Spectroscopy

Gas-phase IR-UV double-resonance laser spectroscopy is an IR absorption technique that bridges the gap between experimental IR spectroscopy and theory. The IR experiments are used to directly evaluate predicted frequencies and potential energy surfaces as well as to probe the structure of isolated molecules. However, a detailed understanding of the underlying mechanisms is, especially in the far-IR regime, still far from complete, even though this is crucial for properly interpreting the recorded IR absorption spectra. Here, events occurring upon excitation to vibrational levels of polycyclic aromatic hydrocarbons by far-IR radiation from the FELIX free electron laser are followed using resonance-enhanced multiphoton ionization spectroscopy. These studies provide detailed insight into how ladder climbing and anharmonicity influence IR-UV spectroscopy and therefore the resulting IR signatures in the far-IR region. Moreover, the potential energy surfaces of these low-frequency delocalized modes are investigated and shown to have a strong harmonic character.

Probing the conformational landscape and thermochemistry of DNA dinucleotide anions via helium nanodroplet infrared action spectroscopy

Isolation of biomolecules in vacuum facilitates characterization of the intramolecular interactions that determine three-dimensional structure, but experimental quantification of conformer thermochemistry remains challenging. Infrared spectroscopy of molecules trapped in helium nanodroplets is a promising methodology for the measurement of thermochemical parameters. When molecules are captured in a helium nanodroplet, the rate of cooling to an equilibrium temperature of ca. 0.4 K is generally faster than the rate of isomerization, resulting in "shock-freezing" that kinetically traps molecules in local conformational minima. This unique property enables the study of temperature-dependent conformational equilibria via infrared spectroscopy at 0.4 K, thereby avoiding the deleterious effects of spectral broadening at higher temperatures. Herein, we demonstrate the first application of this approach to ionic species by coupling electrospray ionization mass spectrometry (ESI-MS) with helium nanodroplet infrared action spectroscopy to probe the structure and thermochemistry of deprotonated DNA dinucleotides. Dinucleotide anions were generated by ESI, confined in an ion trap at temperatures between 90 and 350 K, and entrained in traversing helium nanodroplets. The infrared action spectra of the entrained ions show a strong dependence on pre-pickup ion temperature, consistent with the preservation of conformer population upon cooling to 0.4 K. Non-negative matrix factorization was utilized to identify component conformer infrared spectra and determine temperature-dependent conformer populations. Relative enthalpies and entropies of conformers were subsequently obtained from a van't Hoff analysis. IR spectra and conformer thermochemistry are compared to results from ion mobility spectrometry (IMS) and electronic structure methods. The implementation of ESI-MS as a source of dopant molecules expands the diversity of molecules accessible for thermochemical measurements, enabling the study of larger, non-volatile species.

Conformation control through concurrent N-H⋯S and N-H⋯O[double bond, length as m-dash]C hydrogen bonding and hyperconjugation effects

In addition to the classical N-H⋯O[double bond, length as m-dash]C non-covalent interaction, less conventional types of hydrogen bonding, such as N-H⋯S, may play a key role in determining the molecular structure. In this work, using theoretical calculations in combination with spectroscopic analysis in both gas phase and solution phase, we demonstrate that both these H-bonding modes exist simultaneously in low-energy conformers of capped derivatives of Attc, a thietane α-amino acid. 6-Membered ring inter-residue N-H⋯S interactions (C6^γ), assisted by hyperconjugation between the thietane ring and the backbone, combine with 5-membered ring intra-residue backbone N-H⋯O[double bond, length as m-dash]C interactions (C5) to provide a C5-C6^γ feature that stabilizes a planar geometry in the monomer unit. Two contiguous C5-C6^γ features in the planar dimer implicate an unprecedented three-centre H-bond of the type C[double bond, length as m-dash]O⋯H(N)⋯SR₂, while the trimer adopts two C5-C6^γ features separated by a Ramachandran α-type backbone configuration. These low-energy conformers are fully characterized in the gas phase and support is presented for their existence in solution state.

Gas-Phase Infrared Spectroscopy of Neutral Peptides: Insights from the Far-IR and THz Domain

Gas-phase, double resonance IR spectroscopy has proven to be an excellent approach to obtain structural information on peptides ranging from single amino acids to large peptides and peptide clusters. In this review, we discuss the state-of-the-art of infrared action spectroscopy of peptides in the far-IR and THz regime. An introduction to the field of far-IR spectroscopy is given, thereby highlighting the opportunities that are provided for gas-phase research on neutral peptides. Current experimental methods, including spectroscopic schemes, have been reviewed. Structural information from the experimental far-IR spectra can be obtained with the help of suitable theoretical approaches such as dynamical DFT techniques and the recently developed Graph Theory. The aim of this review is to underline how the synergy between far-IR spectroscopy and theory can provide an unprecedented picture of the structure of neutral biomolecules in the gas phase. The far-IR signatures of the discussed studies are summarized in a far-IR map, in order to gain insight into the origin of the far-IR localized and delocalized motions present in peptides and where they can be found in the electromagnetic spectrum.

Rationalizing the diversity of amide-amide H-bonding in peptides using the natural bond orbital method

Natural bond orbital (NBO) analysis of electron delocalization in a series of capped isolated peptides is used to diagnose amide-amide H-bonding and backbone-induced hyperconjugative interactions, and to rationalize their spectral effects. The sum of the stabilization energies corresponding to the interactions between NBOs that are involved in the H-bonding is demonstrated as an insightful indicator for the H-bond strength. It is then used to decouple the effect of the H-bond distance from that, intrinsic, of the donor/acceptor relative orientation, i.e., the geometrical approach. The diversity of the approaches given by the series of peptides studied enables us to illustrate the crucial importance of the approach when the acceptor is a carbonyl group, and emphasizes that efficient approaches can be achieved despite not matching the usual picture of a proton donor directly facing a lone pair of the proton acceptor, i.e., that encountered in intermolecular H-bonds. The study also illustrates the role of backbone flexibility, partly controlled by backbone-amide hyperconjugative interactions, in influencing the equilibrium structures, in particular by frustrating or enhancing the HB for a given geometrical approach. Finally, the presently used NBO-based HB strength indicator enables a fair prediction of the frequency of the proton donor amide NH stretching mode, but this simple picture is blurred by ubiquitous hyperconjugative effects between the backbone and amide groups, whose magnitude can be comparable to that of the weakest H-bonds.

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.

Can the Partial Peptide SIVSF of the β2-Adrenergic Receptor Recognize Chirality of the Epinephrine Neurotransmitter?

Chirality plays an essential role in biological molecular recognition, such as neurotransmission. Here, we applied electrospray-cold ion trap spectroscopy to complexes of a partial binding motif SIVSF of a β₂-adrenergic receptor pocket with L- and D-epinephrine AdH⁺. The ultraviolet spectrum of the SIVSF-AdH⁺ complex changed drastically when L-AdH⁺ was replaced by its enantiomer. The isomer-selected infrared spectra revealed that D-AdH⁺ was bound to SIVSF by its protonated amino-group or a single catechol OH and induced nonhelical secondary structures of SIVSF. This is in sharp contrast to the helical SIVSF complex with L-AdH⁺, which is close to the natural binding structure with two catechol OHs binding in the receptor. This shows that a short pentapeptide SIVSF can distinguish the chirality of the ligand AdH⁺ as well as the receptor. This stereoselectivity is suggested to arise from an additional interaction involving the hydroxyl group on the chiral carbon.

Interactions of aggregating peptides probed by IR-UV action spectroscopy

The interplay between intramolecular and formed inter-sheet hydrogen bonds and the effect of dispersion interactions on the formation of peptide dimers is studied using IR-UV action spectroscopy.

Conformational assignment of gas phase peptides and their H-bonded complexes using far-IR/THz: IR-UV ion dip experiment, DFT-MD spectroscopy, and graph theory for mode assignment

Graph theory based vibrational modes as new entities for vibrational THz spectroscopy.

Editorial of the PCCP themed issue on "Physical Chemistry for Life Sciences"

This themed issue features high quality and original research on the theme of Physical Chemistry for Life Sciences.

On the turn-inducing properties of asparagine: the structuring role of the amide side chain, from isolated model peptides to crystallized proteins

The anchoring properties of an asparagine (Asn) residue to its local backbone environment in turn model peptides is characterized using gas phase laser spectroscopy and compared to crystallized protein 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.

Unifying the microscopic picture of His-containing turns: from gas phase model peptides to crystallized proteins

The presence in crystallized proteins of a local anchoring between the side chain of a His residue, located in the central position of a γ- or β-turn, and its local main chain environment, is assessed by the comparison of protein structures with relevant isolated model peptides.

Mapping gas phase dipeptide motions in the far-infrared and terahertz domain

Vibrational signatures of Ac-Phe-AA-NH₂ dipeptides are recorded and analysed in the far IR/THz spectral domain (100–800 cm⁻¹, 3–24 THz), with the ‘AA’ amino acid chosen within the series ‘AA’ = Gly, Ala, Pro, Cys, Ser, Val. Phe stands for phenylalanine.

Molecular gels in the gas phase? Gelator-gelator and gelator-solvent interactions probed by vibrational spectroscopy

Benzylidene glucose (BzGlc) is a member of the benzylidene glycoside family. These molecules have the ability to form molecular physical gels. These materials are formed when gelator molecules create a non-covalently bound frame where solvent molecules are trapped. Since the gel formation process and its properties are determined by the subtle balance between non-covalent forces, it is difficult to anticipate them. Quantitative and qualitative understanding of the gelator-gelator and gelator-solvent interactions is needed to better control these materials for important potential applications. We have used gas phase vibrational spectroscopy and theoretical chemistry to study the conformational choices of BzGlc, its dimer and the complexes it forms with water or toluene. To interpret the vibrational spectra we have used the dispersion corrected functional B97D which we have calibrated for the calculation of OH stretching frequencies. Even at the most basic molecular level, it is possible to interrogate a large range of non-covalent interactions ranging from OH → OH hydrogen bonding, to OH → π, and CH → π, all being at the center of gel properties at the macroscopic level.

Far-Infrared Signatures of Hydrogen Bonding in Phenol Derivatives

One of the most direct ways to study the intrinsic properties of the hydrogen-bond interaction is by gas-phase far-infrared (far-IR) spectroscopy because the modes involving hydrogen-bond deformation are excited in this spectral region; however, the far-IR regime is often ignored in molecular structure identification due to the absence of strong far-IR light sources and difficulty in assigning the observed modes by quantum chemical calculations. Far-IR/UV ion-dip spectroscopy using the free electron laser FELIX was applied to directly probe the intramolecular hydrogen-bond interaction in a family of phenol derivatives. Three vibrational modes have been identified, which are expected to be diagnostic for the hydrogen-bond strength: hydrogen-bond stretching and hydrogen-bond-donating and -accepting OH torsion vibrations. Their position is evaluated with respect to the hydrogen bond strength, that is, the length of the hydrogen-bonded OH length. This shows that the hydrogen bond stretching frequency is diagnostic for the size of the ring that is closed by the hydrogen bond, while the strength of the hydrogen bond can be determined from the hydrogen-bond-donating OH torsion frequency. The combination of these two normal modes allows the direct probing of intramolecular hydrogen-bond characteristics using conformation-selective far-IR vibrational spectroscopy.

Gas-Phase Spectroscopic Signatures of Carboxylate-Li(+) Contact Ion Pairs: New Benchmarks For Characterizing Ion Pairing in Solution

The coexistence of several types of ion pairs in solution together with their elusive nature hampers their experimental characterization, which relies in practice on theoretical models resorting to numerous approximations. In this context, a series of isolated contact ion pairs between a lithium cation and phenyl-tagged carboxylate anions of various lengths (Ph-(CH2)n-COO(-), n = 1-3) has been investigated in a conformer-selective manner by IR and UV laser spectroscopy, in conjunction with quantum chemistry calculations. The typical gas-phase IR signature of the bidentate structure formed between the carboxylate moiety and Li(+) has thus been obtained in the CO2(-) stretch region. In addition to the cation-anion interaction, a cation-π interaction occurs simultaneously in the largest system investigated (n = 3). The resulting distorted ion pair structure has been evidenced from both the IR signature of the CO2(-) stretches and the unique vibrationally resolved UV spectroscopy of a phenyl ring interacting with a cation. Such specific spectroscopic signatures of contact ion pairs provide experimental benchmarks, alternative to theoretical predictions, that can assist the assignment of vibrational spectra in solution.

Local NH–π interactions involving aromatic residues of proteins: influence of backbone conformation and ππ* excitation on the π H-bond strength, as revealed from studies of isolated model peptides

Gas phase conformer-selective IR spectroscopy combined and relevant quantum chemistry methods document the NH–π interactions in Phe residues.

Intra-residue interactions in proteins: interplay between serine or cysteine side chains and backbone conformations, revealed by laser spectroscopy of isolated model peptides

Intra-residue interactions play an important role in proteins by influencing local folding of the backbone. Taking advantage of the capability of gas phase experiments to provide relevant information on the intrinsic H-bonding pattern of isolated peptide chains, the intra-residue interactions of serine and cysteine residues, i.e., OH/SH···OC(i) C6 and NH(i···)O/S C5 interactions in Ser/Cys residues, are probed by laser spectroscopy of isolated peptides. The strength of these local side chain-main chain interactions, elegantly documented from their IR spectral features for well-defined conformations of the main chain, demonstrates that a subtle competition exists between the two types of intra-residue bond: the C6 H-bond is the major interaction with Ser, in contrast to Cys where C5 interaction takes over. The restricted number of conformers observed in the gas phase experiment with Ser compared to Cys (where both extended and folded forms are observed) also suggests a significant mediation role of these intra-residue interactions on the competition between the several main chain folding patterns.