N-H⋯X interactions stabilize intra-residue C5 hydrogen bonded conformations in heterocyclic α-amino acid derivatives

Kinetic Isotope Effect in the Unfolding of a Protein Secondary Structure: Calculations for Beta-Sheet Polyglycine Dimers as a Model

In the present work, we performed calculations of the kinetic isotope effect (KIE) on H/D, 14N/15N, 16O/18O, and 12C/13C isotopic substitution in the dissociation of beta-sheet polyglycine dimers of different lengths into two monomer chains. This dissociation reaction, proceeding via breaking of the interchain hydrogen bonds (H-bonds), is considered to be a model of unfolding of the secondary structure of proteins. The calculated strengthening of the interchain hydrogen bonds N-H⋯O=C due to heavy isotope substitution decreases in the row H/D >> 14N/15N > 16O/18O > 12C/13C. The KIE for H/D substitution, defined as the ratio of the rate constants k(H)k(D), was calculated with the use of a "completely loose" transition state model. The results of the calculations show that a very high H/D isotope effect can be achieved for proteins even with moderately long chains connected by dozens of interchain H-bonds. The results obtained also indicate that the heavy isotope substitution in the internal (interchain) and external H-bonds, located on the periphery of a dimer, can provide comparable effects on secondary structure stabilization.

Effects of sulfoxide and sulfone sidechain-backbone hydrogen bonding on local conformations in peptide models

We examine peptide model systems designed to probe short-range N-H⋯OS sidechain-backbone hydrogen bonding involving amino acid residues with sidechain sulfoxide or sulfone functional groups and its effects on local conformations. A strong 7-membered ring hydrogen bond of this type accompanies an intra-residue N-H⋯OC interaction and stabilizes an extended backbone conformation in preference to classical folded structures.

Unveiling the Configurational Landscape of Carbamate: Paving the Way for Designing Functional Sequence-Defined Polymers

Carbamate is an emerging class of a polymer backbone for constructing sequence-defined, abiotic polymers. It is expected that new functional materials can be de novo designed by controlling the primary polycarbamate sequence. While amino acids have been actively studied as building blocks for protein folding and peptide self-assembly, carbamates have not been widely investigated from this perspective. Here, we combined infrared (IR), vibrational circular dichroism (VCD), and nuclear magnetic resonance (NMR) spectroscopy with density functional theory (DFT) calculations to understand the conformation of carbamate monomer units in a nonpolar, aprotic environment (chloroform). Compared with amino acid building blocks, carbamates are more rigid, presumably due to the extended delocalization of π-electrons on the backbones. Cis configurations of the amide bond can be energetically stable in carbamates, whereas peptides often assume trans configurations at low energies. This study lays an essential foundation for future developments of carbamate-based sequence-defined polymer material design.

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.

The Realm of Unconventional Noncovalent Interactions in Proteins: Their Significance in Structure and Function

Proteins and their assemblies are fundamental for living cells to function. Their complex three-dimensional architecture and its stability are attributed to the combined effect of various noncovalent interactions. It is critical to scrutinize these noncovalent interactions to understand their role in the energy landscape in folding, catalysis, and molecular recognition. This Review presents a comprehensive summary of unconventional noncovalent interactions, beyond conventional hydrogen bonds and hydrophobic interactions, which have gained prominence over the past decade. The noncovalent interactions discussed include low-barrier hydrogen bonds, C5 hydrogen bonds, C-H···π interactions, sulfur-mediated hydrogen bonds, n → π* interactions, London dispersion interactions, halogen bonds, chalcogen bonds, and tetrel bonds. This Review focuses on their chemical nature, interaction strength, and geometrical parameters obtained from X-ray crystallography, spectroscopy, bioinformatics, and computational chemistry. Also highlighted are their occurrence in proteins or their complexes and recent advances made toward understanding their role in biomolecular structure and function. Probing the chemical diversity of these interactions, we determined that the variable frequency of occurrence in proteins and the ability to synergize with one another are important not only for ab initio structure prediction but also to design proteins with new functionalities. A better understanding of these interactions will promote their utilization in designing and engineering ligands with potential therapeutic value.

Non-covalent interactions reveal the protein chain δ conformation in a flexible single-residue model

The δ conformation is a local secondary structure in proteins that implicates a π_amide N-H⋯N interaction between a backbone N atom and the NH of the following residue. Small-molecule models thereof have been limited so far to rigid proline-type compounds. We show here that in derivatives of a cyclic amino acid with a sulphur atom in the γ-position, specific side-chain/backbone N-H⋯S interactions stabilize the δ conformation sufficiently to allow it to compete with classical C5 and C7 H-bonded conformers.

Upon DFT-D3 dispersion correction and ECD spectral confirmation, only several conformers can stably coexist for three fungal cycloaspeptides (A, D, G)

Dispersion correction in theoretical determination of cyclopeptide conformations is emphasized. Whether in gas approximation or in solvation simulation, the density functional theory with London dispersion correction (DFT-D3) demonstrates that only 2-3 conformers can stably coexist for cycloaspeptides (A, D, G) at B3LYP-D3 and CAM-B3LYP-D3. Conformational rationality is confirmed by electronic circular dichroism (ECD). Whether for Cotton effect or for excitation energy, TD-B3LYP-D3 has better performances than TD-CAM-B3LYP-D3 because the former can better reproduce the experiment. A molecular orbital analysis is used to interpret ECD, where two energy bands observed in experiment originates from the ππ* transitions other than the σπ* transitions. Long-range correction and solvent effect make H-bonds shorten, and dispersion correction makes them further shorten.

Intermolecular noncovalent interactions with carbon in solution

One of the most familiar carbon-centered noncovalent interactions (NCIs) involving an antibonding π*-orbital situated at the Bürgi-Dunitz angle from the electron donor, mostly lone pairs of electrons, is known as n → π* interactions, and if it involves a σ* orbital in a linear fashion, then it is known as the carbon bond. These NCIs can be intra- or inter-molecular and are usually weak in strength but have a paramount effect on the structure and function of small-molecular crystals and proteins. Surprisingly, the experimental evidence of such interactions in the solution phase is scarce. It is even difficult to determine the interaction energy in the solution. Using NMR spectroscopy aided with molecular dynamics (MD) simulation and high-level quantum mechanical calculations, herein we provide the experimental evidence of intermolecular carbon-centered NCIs in solution. The challenge was to find appropriate heterodimers that could sustain room temperature thermal energy and collisions from the solvent molecules. However, after several trial model compounds, the pyridine-N-oxide:dimethyltetracyanocyclopropane (PNO-DMTCCP) complex was found to be a good candidate for the investigation. NBO analyses show that the PNO:DMTCCP complex is stabilized mainly by intermolecular n → π* interaction when a weaker carbon bond gives extra stability to the complex. From the NMR study, it is observed that the NCIs between DMTCCP and PNO are enthalpy driven with an enthalpy change of -28.12 kJ mol-1 and dimerization energy of ∼-38 kJ mol-1 is comparable to the binding energies of a conventional hydrogen-bonded dimer. This study opens up a new strategy to investigate weak intermolecular interactions such as n → π* interaction and carbon bonds in the solution phase.

Detecting the Hydrogen Bond Cooperativity in a Protein β-Sheet by H/D Exchange

The hydrogen bond (H-bond) cooperativity in the β-sheet of GB3 is investigated by a NMR hydrogen/deuterium (H/D) exchange method. It is shown that the weakening of one backbone N-H…O=C H-bond between two β-strands, β1 and β2, due to the exchange of NH to ND of the H-bond donor in β1, perturbs the chemical shift of ¹³C_α, ¹³C_β, ¹H_α, ¹H_N, and ¹⁵N of the H-bond acceptor and its following residue in β2. Quantum mechanical calculations suggest that the -H-bond chemical shift isotope effect is caused by the structural reorganization in response to the H-bond weakening. This structural reorganization perturbs four neighboring H-bonds, with three being weaker and one being stronger, indicating that three H-bonds are cooperative and one is anticooperative with the perturbed H-bond. The sign of the cooperativity depends on the relative position of the H-bonds. This H-bond cooperativity, which contributes to β-sheet stability overall, can be important for conformational coupling across the β-sheet.

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.

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.

Newly identified C–H⋯O hydrogen bond in histidine

Histidine C–H bonds observed in protein structures include (clockwise from top left): myoglobin, β-lactamase, and photoactive yellow protein; calculations indicate that tautomeric/protonation state influences H-bonding ability (bottom left).