Competition between Amide Stacking and Intramolecular H Bonds in γ-Peptide Derivatives: Controlling Nearest-Neighbor Preferences

Conformation Controlled Hydrogelation of Minimalistic α, γ Hybrid Peptide

A majority of short peptide (≤7 amino acids) hydrogels are primarily assembled via cross β-structure formation. In contrast to the natural trend, herein, we report the formation of supramolecular hydrogel from the ultrashort hybrid folded peptide composed of canonical α-amino acid and noncanonical γ-amino acid, Fmoc-γPhe-Phe-OH. The designed hybrid peptide hydrogel is composed of entangled fibers, has viscoelastic properties, exhibits proteolytic stability, and exhibits cytocompatibility with L929 fibroblast cells. Mutating the peptide sequence by altering the position of γPhe from the N-termini to C-termini transforms the self-assembly into crystalline aggregates. Combining FTIR, 2D NMR, and DFT calculations revealed that the hydrogel-forming peptide adopts a C9 H-bonded conformation, resembling the well-known γ-turn. However, the isomeric hybrid peptide adopts an extended structure. The present study highlights the importance of secondary structure in the higher order assembly of minimalist hybrid peptides and broadens the range of secondary structures to design short peptide-based hydrogels.

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.

Conformer-Specific Spectroscopy and IR-Induced Isomerization of a Model γ-Peptide: Ac-γ4-Phe-NHMe

Single-conformation IR and UV spectroscopy of the prototypical capped γ-peptide Ac-γ⁴-Phe-NHMe (γ⁴F) was carried out under jet-cooled conditions in the gas phase in order to understand its innate conformational preferences in the absence of a solvent. We obtained conformer-specific IR and UV spectra and compared the results with calculations to make assignments and explore the differences between the γ²- and γ⁴-substituted molecules. We found four conformers of γ⁴F in our experiment. Three conformers form nine-membered hydrogen-bonded rings (C9) enclosed by an NH···O═C H-bond but differing in their phenyl ring positions (a, g+, and g-). The fourth conformer forms a strained seven-membered hydrogen-bonded ring in which the amide groups lie in a nominally anti-parallel arrangement stacked on top of one another (labeled S7). This conformer is a close analogue of the amide-stacked conformer (S) found previously in γ²F, in which the Phe side chain is substituted at the γ² position, Ac-γ2-Phe-NHMe (J. Am. Chem. Soc. 2009, 131, 14243-14245). IR population transfer spectroscopy was used to determine the fractional abundances of the γ⁴F conformers in the expansion. A combination of force field and density functional theory calculations is used to map out the conformational potential energy surfaces for γ⁴F and compare it with its γ²F counterpart. Based on this analysis, the phenyl ring prefers to take up structures that facilitate NH···π interactions in γ⁴F or avoid phenyl interactions with the C═O group in γ²F. The disconnectivity graph for γ⁴F reveals separate basins associated with the C9 and amide-stacked conformational families, which are separated by a barrier of about 42 kJ/mol. The overall shape of the potential energy surface bears a resemblance to peptides and proteins that have a misfolding pathway that competes with the formation of the native structure.

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.

Urea-aromatic interactions in biology

Noncovalent interactions are key determinants in both chemical and biological processes. Among such processes, the hydrophobic interactions play an eminent role in folding of proteins, nucleic acids, formation of membranes, protein-ligand recognition, etc.. Though this interaction is mediated through the aqueous solvent, the stability of the above biomolecules can be highly sensitive to any small external perturbations, such as temperature, pressure, pH, or even cosolvent additives, like, urea-a highly soluble small organic molecule utilized by various living organisms to regulate osmotic pressure. A plethora of detailed studies exist covering both experimental and theoretical regimes, to understand how urea modulates the stability of biological macromolecules. While experimentalists have been primarily focusing on the thermodynamic and kinetic aspects, theoretical modeling predominantly involves mechanistic information at the molecular level, calculating atomistic details applying the force field approach to the high level electronic details using the quantum mechanical methods. The review focuses mainly on examples with biological relevance, such as (1) urea-assisted protein unfolding, (2) urea-assisted RNA unfolding, (3) urea lesion interaction within damaged DNA, (4) urea conduction through membrane proteins, and (5) protein-ligand interactions those explicitly address the vitality of hydrophobic interactions involving exclusively the urea-aromatic moiety.

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.

Single Conformation Spectroscopy of Suberoylanilide Hydroxamic Acid: A Molecule Bites Its Tail

Suberoylanilide hydroxamic acid (SAHA) is a histone deacetylase inhibitor that causes growth arrest and differentiation of many tumor types and is an approved drug for the treatment of cancer. The chemical structure of SAHA consists of formanilide "head" and a hydroxamic acid "tail" separated by an n-hexyl chain, C₆H₅NH(C═O)-(CH₂)₆-(C═O)NHOH. The alkyl chain's preference for extended structures is in competition with tail-to-head (T-H) or head-to-tail (H-T) hydrogen bonds between the amide and hydroxamic acid groups. Laser desorption was used to bring SAHA into the gas phase and cool it in a supersonic expansion before interrogation with mass-resolved resonant two-photon ionization spectroscopy. Single conformation UV spectra in the S₀-S₁ region and infrared spectra in the hydride stretch and mid-IR regions were recorded using IR-UV hole-burning and resonant ion-dip infrared spectroscopy, respectively. Three conformers of SAHA were distinguished and spectroscopically characterized. Comparison of the experimental IR spectra with the predictions of density functional theory calculations (DFT, B3LYP D3BJ/6-31+G(d)) leads to assignments for the three conformers, all of which possess tightly folded alkyl chains that enable formation of a T-H (conformer A) or H-T (conformers B and C) hydrogen bonds. A modified version of the generalized Amber force field was developed to more accurately describe the hydroxamic acid OH internal rotor potential, leading to predictions for the relative energies in reasonable agreement with experiment. This force field was used to generate a disconnectivity graph for the low-energy portion of the potential energy landscape of SAHA. This disconnectivity graph contains more than one hundred minima and maps out the lowest-energy pathways between them, which could then be characterized via DFT calculations. This combination of force field and DFT calculations provides insight into the potential energy landscape and how population was funneled into the three observed conformers.

Natural and non-natural amino-acid side-chain substitutions: affinity and diffraction studies of meditope-Fab complexes

Herein, multiple crystal structures of meditope peptide derivatives incorporating natural and unnatural amino acids bound to the cetuximab Fab domain are presented. The affinity of each derivative was determined by surface plasmon resonance and correlated to the atomic structure. Overall, it was observed that the hydrophobic residues in the meditope peptide, Phe3, Leu5 and Leu10, could accommodate a number of moderate substitutions, but these invariably reduced the overall affinity and half-life of the interaction. In one case, the substitution of Phe3 by histidine led to a change in the rotamer conformation, in which the imidazole ring flipped to a solvent-exposed position. Based on this observation, Phe3 was substituted by diphenylalanine and it was found that the phenyl rings in this variant mimic the superposition of the Phe3 and His3 structures, producing a moderate increase, of 1.4-fold, in the half-life of the complex. In addition, it was observed that substitution of Leu5 by tyrosine and glutamate strongly reduced the affinity, whereas the substitution of Leu5 by diphenylalanine moderately reduced the half-life (by approximately fivefold). Finally, it was observed that substitution of Arg8 and Arg9 by citrulline dramatically reduced the overall affinity, presumably owing to lost electrostatic interactions. Taken together, these studies provide insight into the meditope-cetuximab interaction at the atomic level.

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.

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.

Synthesis, structural and spectroscopic evaluations and nonlinear optical properties of 3,5-bis(4-methoxyphenyl)-4,5-dihydro-1H-pyrazole-1-carbothioic O-acid

In this study, we report a combined experimental and theoretical study on 3,5-bis(4-methoxyphenyl)-4,5-dihydro-1H-pyrazole-1-carbothioic O-acid (C18H18N2O3S) molecule. The compound crystallizes in the trigonal space group R-3 with a=b=27.7151(12) Å, c=12.4866(6) Å, α=β=90.0°, γ=120.0° and Z=18. The crystal packing is stabilized by O-H⋯O and O-H⋯S intermolecular hydrogen bonds. These hydrogen bond interactions are also proved by NBO analysis. A detailed spectroscopic investigation is performed by the application of FT-IR and FT-NMR in addition to the theoretical approaches. Small energy gap between the frontier molecular orbitals is responsible for the nonlinear optical activity of the title molecule.

Robustness of thioamide dimer synthon, carbon bonding and thioamide–thioamide stacking in ferrocene-based thiosemicarbazones

The role of thioureas in crystal engineering as robust supramolecular synthons is now recognized, but their analogs, namely thiosemicarbazones/N-iminothioureas, have not received the attention they deserve.

Effect of intermolecular hydrogen bonding, vibrational analysis and molecular structure of a biomolecule: 5-Hydroxymethyluracil

In the present work, the experimental and theoretical vibrational spectra of 5-hydroxymethyluracil were investigated. The FT-IR (4000-400cm(-1)) spectrum of the molecule in the solid phase was recorded. The geometric parameters (bond lengths and bond angles), vibrational frequencies, Infrared intensities of the title molecule in the ground state were calculated using density functional B3LYP and M06-2X methods with the 6-311++G(d,p) basis set for the first time. The optimized geometric parameters and theoretical vibrational frequencies were found to be in good agreement with the corresponding experimental data, and with the results found in the literature. The vibrational frequencies were assigned based on the potential energy distribution using the VEDA 4 program. The dimeric form of 5-hydroxymethyluracil molecule was also simulated to evaluate the effect of intermolecular hydrogen bonding on its vibrational frequencies. It was observed that the NH stretching modes shifted to lower frequencies, while its in-plane and out-of-plane bending modes shifted to higher frequencies due to the intermolecular NH⋯O hydrogen bond. Also, the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energies and diagrams were presented.

Cyclic constraints on conformational flexibility in γ-peptides: conformation specific IR and UV spectroscopy

Single-conformation spectroscopy has been used to study two cyclically constrained and capped γ-peptides: Ac-γACHC-NHBn (hereafter γACHC, Figure 1a), and Ac-γACHC-γACHC-NHBn (γγACHC, Figure 1b), under jet-cooled conditions in the gas phase. The γ-peptide backbone in both molecules contains a cyclohexane ring incorporated across each Cβ-Cγ bond and an ethyl group at each Cα. This substitution pattern was designed to stabilize a (g+, g+) torsion angle sequence across the Cα-Cβ-Cγ segment of each γ-amino acid residue. Resonant two-photon ionization (R2PI), infrared-ultraviolet hole-burning (IR-UV HB), and resonant ion-dip infrared (RIDIR) spectroscopy have been used to probe the single-conformation spectroscopy of these molecules. In both γACHC and γγACHC, all population is funneled into a single conformation. With RIDIR spectra in the NH stretch (3200-3500 cm(-1)) and amide I/II regions (1400-1800 cm(-1)), in conjunction with theoretical predictions, assignments have been made for the conformations observed in the molecular beam. γACHC forms a single nearest-neighbor C9 hydrogen-bonded ring whereas γγACHC takes up a next-nearest-neighbor C14 hydrogen-bonded structure. The gas-phase C14 conformation represents the beginning of a 2.614-helix, suggesting that the constraints imposed on the γ-peptide backbone by the ACHC and ethyl groups already impose this preference in the gas-phase di-γ-peptide, in which only a single C14 H-bond is possible, constituting one full turn of the helix. A similar conformational preference was previously documented in crystal structures and NMR analysis of longer γ-peptide oligomers containing the γACHC subunit [Guo, L., et al. Angew. Chem. Int. Ed. 2011, 50, 5843-5846]. In the gas phase, the γACHC-H2O complex was also observed and spectroscopically interrogated in the molecular beam. Here, the monosolvated γACHC retains the C9 hydrogen bond observed in the bare molecule, with the water acting as a bridge between the C-terminal carbonyl and the π-cloud of the UV chromophore. This is in contrast to the unconstrained γ-peptide-H2O complex, which incorporates H2O into both C9 and amide-stacked conformations.

Vibrational analysis of 4-chloro-3-nitrobenzonitrile by quantum chemical calculations

In the present study, the experimental and theoretical harmonic and anharmonic vibrational frequencies of 4-chloro-3-nitrobenzonitrile were investigated. The experimental FT-IR (400-4000 cm(-1)) and μ-Raman spectra (100-4000 cm(-1)) of the molecule in the solid phase were recorded. Theoretical vibrational frequencies and geometric parameters (bond lengths and bond angles) were calculated using ab initio Hartree Fock (HF), density functional B3LYP and M06-2X methods with 6-311++G(d,p) basis set by Gaussian 09 W program, for the first time. The assignments of the vibrational frequencies were performed by potential energy distribution (PED) analysis by using VEDA 4 program. The theoretical optimized geometric parameters and vibrational frequencies were compared with the corresponding experimental data, and they were seen to be in a good agreement with each other. Also, the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energies were found.

Efficient stacking on protein amide fragments

The less polar π-surface of protein amide groups is exposed in many receptor binding sites, either as part of the backbone or in Gln/Asn side chains. Using quantum chemical calculations and Protein Data Bank (PDB) searches on model systems, we investigate the energetics and geometric preferences for the stacking on amide groups of a large number of heteroarenes that are relevant to medicinal chemistry. From this study, we discern that the stacking energy of an aromatic ligand substituent can be improved by: 1) orienting the fragment dipole vector such that it is aligned in an antiparallel fashion with the dipole of the interacting protein amide group, 2) increasing its dipole moment, and 3) decreasing its π-electron density. These guidelines should be helpful to more rationally exploit this interaction type in future structure-based drug design.