Gas-phase folding of a two-residue model peptide chain: on the importance of an interplay between experiment and theory

Length-Dependent Transition from Extended to Folded Shapes in Short Oligomers of an Azetidine-Based α-Amino Acid: The Critical Role of NH···N H-Bonds

Hydrogen bonds (H-bonds) are ubiquitous in peptides and proteins and are central to the stabilization of their structures. Inter-residue H-bonds between non-adjacent backbone amide NH and C=O motifs lead to the well-known secondary structures of helices, turns and sheets, but it is recognized that other H-bonding modes may be significant, including the weak intra-residue H-bond (called a C5 H-bond) that implicates the NH and C=O motifs of the same amino acid residue. Peptide model compounds that adopt stable C5 H-bonds are not readily available and the so-called 2.0₅-helix, formed by successive C5 H-bonds, is an elusive secondary structure. Using a combination of theoretical chemistry and spectroscopic studies in both the gas phase and solution phase, we have demonstrated that derivatives of 3-amino-1-methylazetidine-3-carboxylic acid, Aatc(Me) can form sidechain-backbone N-H···N C6γ H-bonds that accompany-and thereby stabilize-C5 H-bonds. In the capped trimer of Aatc(Me), extended C5/C6γ motifs are sufficiently robust to challenge classical 3₁₀-helix formation in solution and the fully-extended 2.0₅-helix conformer has been characterized in the gas phase. Concurrent H-bonding support for successive C5 motifs is a new axiom for stabilizing the extended backbone secondary structure in short peptides.

Sequence dependent folding motifs of the secondary structures of Gly-Pro and Pro-Gly containing oligopeptides

Folding motifs of the secondary structures of peptides and proteins are primarily based on the hydrogen bonding interactions in the backbone as well as the sequence of the amino acid residues present. For instance, the β-turn structure directed by the Pro-Gly sequence is the key to the β-hairpin structure of peptides/proteins as well as a selective site for the enzymatic hydroxylation of pro-collagen. Herein, we have investigated the sequence dependent folding motifs of end-protected Gly-Pro and Pro-Gly dipeptides using a combination of gas phase laser spectroscopy, quantum chemistry calculations, solution phase IR and NMR spectroscopy and single crystal X-Ray diffraction (XRD). All three observed conformers of the Gly-Pro peptide in the gas phase have been found to have extended β-strand or polyproline-II (PP-II) structures with C5-C7 hydrogen bonding interactions, which correlates well with the structure obtained from solution phase spectroscopy and XRD. On the other hand, we have found that the Pro-Gly peptide has a C10/β-turn structure in the solution phase in contrast to the C7-C7 (i.e. 2₇-ribbon) structure observed in the gas phase. Although the lowest energy structure in the gas phase is not C10, we find that C7-C7 is an abundantly found structural motif of Pro-Gly containing peptides in the Cambridge Structural Database, indicating that the gas phase conformers are not sampling any unusual forms. We surmise that the role of the solvent could be crucial in dictating the preferential stabilization of the C10 structure in the solution phase. The present investigation provides a comprehensive picture of the folding motifs of the Gly-Pro and Pro-Gly peptides observed in the gas phase and condensed phase weaving a fine interplay of the intrinsic conformational properties, solvation, and crystal packing of the peptides.

Selenium in Proteins: Conformational Changes Induced by Se Substitution on Methionine, as Studied in Isolated Model Peptides by Optical Spectroscopy and Quantum Chemistry

The side-chain of methionine residues is long enough to establish NH⋯S H-bonds with neighboring carbonyl groups of the backbone, giving rise to so-called intra-residue 6^δ and inter-residue 7^δ H-bonds. The aim of the present article is to document how the substitution of sulfur with a selenium atom affects the H-bonding of the Met system. This was investigated both experimentally and theoretically by conformation-resolved optical spectroscopy, following an isolated molecule approach. The present work emphasizes the similarities of the Met and Sem residues in terms of conformational structures, energetics, NH⋯Se/S H-bond strength and NH stretch spectral shifts, but also reveals subtle behavior differences between them. It provides evidence for the sensitivity of the H-bonding network with the folding type of the Sem/Met side-chains, where a simple flip of the terminal part of the side-chain can induce an extra 50 cm⁻¹ spectral shift of the NH stretch engaged in a 7^δ NH⋯S/Se bond.

Structural Properties of Phenylalanine-Based Dimers Revealed Using IR Action Spectroscopy

Peptide segments with phenylalanine residues are commonly found in proteins that are related to neurodegenerative diseases. However, the self-assembly of phenylalanine-based peptides can be also functional. Peptides containing phenylalanine residues with different side caps, composition, and chemical alteration can form different types of nanostructures that find many applications in technology and medicine. Various studies have been performed in order to explain the remarkable stability of the resulting nanostructures. Here, we study the early stages of self-assembly of two phenylalanine derived peptides in the gas phase using IR action spectroscopy. Our focus lies on the identification of the key intra- and intermolecular interactions that govern the formation of the dimers. The far-IR region allowed us to distinguish between structural families and to assign the 2-(2-amino-2-phenylacetamido)-2-phenylacetic acid (PhgPhg) dimer to a very symmetric structure with two intermolecular hydrogen bonds and its aromatic rings folded away from the backbone. By comparison with the phenylalanine-based peptide cyclic L-phenylalanyl-L-phenylalanine (cyclo-FF), we found that the linear FF dimer likely adopts a less ordered structure. However, when one more phenylalanine residue is added (FFF), a more structurally organized dimer is formed with several intermolecular hydrogen bonds.

Double Ion Trap Laser Spectroscopy of Alkali Metal Ion Complexes with a Partial Peptide of the Selectivity Filter in K+ Channels─Temperature Effect and Barrier for Conformational Conversions

Potassium ion channels selectively permeate K⁺, as well as Rb⁺ and Cs⁺ to some degree, while excluding Na⁺ and Li⁺. Conformations of alkali metal complexes of Ac-Tyr-NHMe, a model peptide of the selectivity filter in a K⁺ channel, were previously found to correlate with the permeability of alkali metal ions to a K⁺ channel by cold ion trap infrared spectroscopy. With an additional temperature-controlled ion trap, we examined the conformations of the alkali metal complexes, allowing the ions to collide with a He buffer gas at different temperatures, prior to spectroscopic investigation. The conformational distribution of the K⁺-peptide complex shows the most significant variation with temperature, which suggests that this complex has more flexibility when complexed with K⁺ and suggests lower barrier heights than other metal-peptide complexes. The variability of the conformational distribution with temperature for the ions follows the same order of ion permeability of a K⁺ channel. This work demonstrates that the additional temperature-controlled ion trap is a powerful tool to explore the conformational landscape of flexible molecular systems.

Competition between inter and intramolecular hydrogen bond evidenced by vibrational circular dichroism spectroscopy: The case of (1S,2R)-(-)-cis-1-amino-2-indanol

The infrared (IR) absorption and vibrational circular dichroism (VCD) spectra of an intramolecularly hydrogen-bonded chiral amino-alcohol, (1S,2R)-(-)-cis-1-amino-2-indanol, are studied in DMSO-d₆ . The spectra are simulated at the density functional theory (DFT) level within the frame of the cluster-in-the-liquid model. Both IR and VCD spectra show a clear signature of the formation of intermolecular hydrogen bonds at the detriment of the intramolecular OH … N interaction present in the isolated molecule. Two solvent molecules are necessary to reproduce the experimental spectra. Whereas the first DMSO molecule captures the main spectral modifications due to hydrogen bond formation between the solute and the solvent, the second DMSO molecule is necessary for a good description of the Boltzmann contribution of the different complexes, based on their Gibbs free energy.

Side-Chain Polarity Modulates the Intrinsic Conformational Landscape of Model Dipeptides

The intrinsic conformational preferences of small peptides may provide additional insight into the thermodynamics and kinetics of protein folding. In this study, we explore the underlying energy landscapes of two model peptides, namely, Ac-Ala-NH₂ and Ac-Ser-NH₂, using geometry-optimization-based tools developed within the context of energy landscape theory. We analyze not only how side-chain polarity influences the structural preferences of the dipeptides, but also other emergent properties of the landscape, including heat capacity profiles, and kinetics of conformational rearrangements. The contrasting topographies of the free energy landscape agree with recent results from Fourier transform microwave spectroscopy experiments, where Ac-Ala-NH₂ was found to exist as a mixture of two conformers, while Ac-Ser-NH₂ remained structurally locked, despite exhibiting an apparently rich conformational landscape.

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.

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.

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.

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.

Molecular Vibrations of an Oxygen-Evolving Complex and Its Synthetic Mimic

Bio-inspired catalysis for artificial photosynthesis has been widely studied for decades, in particular, with the purpose of using bio-disposable and non-toxic metals as building blocks. The characterisation of such catalysts has been achieved by using different kinds of spectroscopic methods, from X-ray crystallography to NMR spectroscopy. An artificial Mn₄ CaO₄ cubane cluster with dangling Mn₄ was synthesised in 2015 [Zhang et al. Science 2015, 348, 690-693]; this cluster showed many structural similarities to that of the natural oxygen-evolving complex. An accurate structural and spectroscopic comparison between the natural and artificial systems is highly relevant to understand the catalytic mechanism. Among data from different techniques, the differential FTIR spectra (S_n+1 -S_n ) of photosystem II are still lacking a complete interpretation. The availability of IR data of the artificial cluster offers a unique opportunity to assign absolute absorption spectra on a well-defined and easier to interpret analogous moiety. The present work aims to investigate the novel inorganic compound as a model system for an oxygen-evolving complex through measurement of its spectroscopic properties. The experimental results are compared with calculations by using a variety of theoretical methods (normal mode analysis, effective normal mode analysis) in the S₁ state. We underline the similarities and the differences in the computational spectra based on atomistic models of Mn₄ CaO₅ and Mn₄ CaO₄ complexes.

Formation of Neutral Peptide Aggregates as Studied by Mass-Selective IR Action Spectroscopy

The spontaneous aggregation of proteins and peptides is widely studied owing to its relation to neurodegenerative diseases. To understand the underlying principles of peptide aggregation, elucidation of structure and structural changes upon their formation is key. This level of detail can be obtained by studying the peptide self‐assembly in the gas phase. Structural characterization of aggregates is mainly done on charged species, as adding charges is an intrinsic part of the technique to bring molecules into the gas phase. Studying neutral peptide aggregates will complement the existing picture. These studies are restricted to dimers due to experimental limitations. Herein, we present advances in laser desorption molecular beam spectroscopy to form neutral peptide aggregates consisting of up to 14 monomeric peptides in the gas phase. The combination of this technique with IR–UV spectroscopy allowed us to select each aggregate by size and subsequently characterize its structure.

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.

Competition between folded and extended structures of alanylalanine (Ala-Ala) in a molecular beam

Structural analysis of the neutral dipeptide Ala-Ala by action spectroscopy using IRMPD-VUV spectroscopy reveals predominance of extended structure.

The axial/equatorial conformational landscape and intramolecular dispersion: new insights from the rotational spectra of monoterpenoids

Using rotational spectroscopy and quantum chemistry calculations, we show that intramolecular dispersion stabilises the axial conformers of monoterpenoids, and that an accurate account of these interactions is challenging for theoretical methods.

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.

N-Protonated Isomers and Coulombic Barriers to Dissociation of Doubly Protonated Ala8Arg

Collision-induced dissociation (or tandem mass spectrometry, MS/MS) of a protonated peptide results in a spectrum of fragment ions that is useful for inferring amino acid sequence. This is now commonplace and a foundation of proteomics. The underlying chemical and physical processes are believed to be those familiar from physical organic chemistry and chemical kinetics. However, first-principles predictions remain intractable because of the conflicting necessities for high accuracy (to achieve qualitatively correct kinetics) and computational speed (to compensate for the high cost of reliable calculations on such large molecules). To make progress, shortcuts are needed. Inspired by the popular mobile proton model, we have previously proposed a simplified theoretical model in which the gas-phase fragmentation pattern of protonated peptides reflects the relative stabilities of N-protonated isomers, thus avoiding the need for transition-state information. For singly protonated Ala _n (n = 3-11), the resulting predictions were in qualitative agreement with the results from low-energy MS/MS experiments. Here, the comparison is extended to a model tryptic peptide, doubly protonated Ala₈Arg. This is of interest because doubly protonated tryptic peptides are the most important in proteomics. In comparison with experimental results, our model seriously overpredicts the degree of backbone fragmentation at N₉. We offer an improved model that corrects this deficiency. The principal change is to include Coulombic barriers, which hinder the separation of the product cations from each other. Coulombic barriers may be equally important in MS/MS of all multiply charged peptide ions. Graphical Abstract ᅟ.

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.

Structure and stability of clusters of β-alanine in the gas phase: importance of the nature of intermolecular interactions

We present a theoretical study of neutral clusters of β-alanine molecules in the gas phase, (β-ala)_nn ≤ 5.

The role of amino acid side chains in stabilizing dipeptides: the laser ablation Fourier transform microwave spectrum of Ac-Val-NH2

The steric effects imposed by the isopropyl group of valine in the conformational stabilization of the capped dipeptide N-acetyl-l-valinamide (Ac-Val-NH₂) have been studied by laser ablation molecular beam Fourier transform microwave (LA-MB-FTMW) spectroscopy.

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.

Potential energy surface of fluoroxene: experiment and theory

The potential energy surface (PES) of the general anesthetic fluoroxene was probed in a supersonic jet expansion using broadband CP-FTMW spectroscopy and theoretical calculations.

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.

Conformation-specific spectroscopy of capped glutamine-containing peptides: role of a single glutamine residue on peptide backbone preferences

The conformational preferences of a series of short, aromatic-capped, glutamine-containing peptides have been studied under jet-cooled conditions in the gas phase.

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.

Unraveling non-covalent interactions within flexible biomolecules: from electron density topology to gas phase spectroscopy

The NCI (Non-Covalent Interactions) method, a recently-developed theoretical strategy to visualize weak non-covalent interactions from the topological analysis of the electron density and of its reduced gradient, is applied in the present paper to document intra- and inter-molecular interactions in flexible molecules and systems of biological interest in combination with IR spectroscopy. We first describe the conditions of application of the NCI method to the specific case of intramolecular interactions. Then we apply it to a series of stable conformations of isolated molecules as an interpretative technique to decipher the different physical interactions at play in these systems. Examples are chosen among neutral molecular systems exhibiting a large diversity of interactions, for which an extensive spectroscopic characterization under gas-phase isolation conditions has been obtained using state-of-the-art conformer-specific experimental techniques. The interactions presently documented range from weak intra-molecular H-bonds in simple amino-alcohols, to more complex patterns, with interactions of various strengths in model peptides, as well as in chiral bimolecular systems, where invaluable hints for the understanding of chiral recognition are revealed. We also provide a detailed technical appendix, which discusses the choices of cut-offs as well as the applicability of the NCI analysis to specific constrained systems, where local effects require attention. Finally, the NCI technique provides IR spectroscopists with an elegant visualization of the interactions that potentially impact their vibrational probes, namely the OH and NH stretching motions. This contribution illustrates the power and the conditions of use of the NCI technique, with the aim of providing an easy tool for all chemists, experimentalists and theoreticians, for the visualization and characterization of the interactions shaping complex molecular systems.

Conformational analysis of a polyconjugated protein-binding ligand by joint quantum chemistry and polarizable molecular mechanics. Addressing the issues of anisotropy, conjugation, polarization, and multipole transferability

We investigate the conformational properties of a potent inhibitor of neuropilin-1, a protein involved in cancer processes and macular degeneration. This inhibitor consists of four aromatic/conjugated fragments: a benzimidazole, a methylbenzene, a carboxythiourea, and a benzene-linker dioxane, and these fragments are all linked together by conjugated bonds. The calculations use the SIBFA polarizable molecular mechanics procedure. Prior to docking simulations, it is essential to ensure that variations in the ligand conformational energy upon rotations around its six main-chain torsional bonds are correctly represented (as compared to high-level ab initio quantum chemistry, QC). This is done in two successive calibration stages and one validation stage. In the latter, the minima identified following independent stepwise variations of each of the six main-chain torsion angles are used as starting points for energy minimization of all the torsion angles simultaneously. Single-point QC calculations of the minimized structures are then done to compare their relative energies ΔE conf to the SIBFA ones. We compare three different methods of deriving the multipoles and polarizabilities of the central, most critical moiety of the inhibitor: carboxythiourea (CTU). The representation that gives the best agreement with QC is the one that includes the effects of the mutual polarization energy E pol between the amide and thioamide moieties. This again highlights the critical role of this contribution. The implications and perspectives of these findings are discussed.

Accurate Characterization of the Peptide Linkage in the Gas Phase: A Joint Quantum-Chemical and Rotational Spectroscopy Study of the Glycine Dipeptide Analogue

Accurate structures of aminoacids in the gas phase have been obtained by joint microwave and quantum-chemical investigations. However, the structure and conformational behavior of α-aminoacids once incorporated into peptide chains are completely different and have not yet been characterized with the same accuracy. To fill this gap, we present here an accurate characterization of the simplest dipeptide analogue (N-acetyl-glycinamide) involving peptidic bonds. State-of-the-art quantum-chemical computations are complemented by a comprehensive study of the rotational spectrum using a combination of Fourier transform microwave spectroscopy with laser ablation. The coexistence of the C7 and C5 conformers has been proved and energetically as well as spectroscopically characterized. This joint theoretical-experimental investigation demonstrated the feasibility of obtaining accurate structures for flexible small biomolecules, thus paving the route to the elucidation of the inherent behavior of peptides.