Conformational Analysis of n→π* Interactions in Collagen Triple Helix Models

Neighbor effect on conformational spaces of alanine residue in azapeptides

The conformational properties of Alanine (Ala) residue have been investigated to understand protein folding and develop force fields. In this work, we examined the neighbor effect on the conformational spaces of Ala residue using model azapeptides, Ac-Ala-azaGly-NHMe (3, AaG), and Ac-azaGly-Ala-NHMe (4, aGA1). Ramachandran energy maps were generated by scanning (φ, ψ) dihedral angles of the Ala residues in models with the fixed dihedral angles (φ = ±90°, ψ = ±0° or ±180°) of azaGly residue using LCgau-BOP and LCgau-BOP + LRD functionals in the gas and water phases. The integral-equation-formalism polarizable continuum model (IEF-PCM) and a solvation model density (SMD) were employed to mimic the solvation effect. The most favorable conformation of Ala residue in azapeptide models is found as the polyproline II (βP), inverse γ-turn (γ'), β-sheet (βS), right-handed helix (αR), or left-handed helix (αL) depending on the conformation of neighbor azaGly residue in isolated form. Solvation methods exhibit that the Ala residue favors the βP, δR, and αR conformations regardless of its position in azapeptides 3 and 4 in water. Azapeptide 5, Ac-azaGly-Ala-NH2 (aGA2), was synthesized to evaluate the theoretical results. The X-ray structure showed that azaGly residue adopts the polyproline II (βP) and Ala residue adopts the right-handed helical (αR) structure in aGA2. The conformational preferences of aGA2 and the dimer structure of aGA2 based on the X-ray structure were examined to assess the performance of DFT functionals. In addition, the local minima of azapeptide 6, Ac-Phe-azaGly-NH2 (FaG), were compared with the previous experimental results. SMD/LCgau-BOP + LRD methods agreed well with the reported experimental results. The results suggest the importance of weak dispersion interactions, neighbor effect, and solvent influence in the conformational preferences of Ala residue in model azapeptides.

Computational Investigation of Conformational Properties of Short Azapeptides: Insights from DFT Study and NBO Analysis

Azapeptides have gained much attention due to their ability to enhance the stability and bioavailability of peptide drugs. Their structural preferences, essential to understanding their function and potential application in the peptide drug design, remain largely unknown. In this work, we systematically investigated the conformational preferences of three azaamino acid residues in tripeptide models, Ac-azaXaa-Pro-NHMe [Xaa = Asn (4), Asp (5), Ala (6)], using the popular DFT functionals, B3LYP and B3LYP-D3. A solvation model density (SMD) was used to mimic the solvation effect on the conformational behaviors of azapeptides in water. During the calculation, we considered the impact of the amide bond in the azapeptide models on the conformational preferences of models 4-6. We analyzed the effect of the HB between the side-chain main chain and main-chain main-chain on the conformational behaviors of azapeptides 4-6. We found that the predicted lowest energy conformation for the three models differs depending on the calculation methods. In the gas phase, B3LYP functional indicates that the conformers tttANP-1 and tttADP-1 of azapeptides 4 and 5 correspond to the type I of β-turn, the lowest energy conformation with all-trans amide bonds. Considering the dispersion correction, B3LYP-D3 functional predicts the conformers tctANP-2 and tctADP-3 of azapeptide 4 and 5, which contain the cis amide bond preceding the Pro residue, as the lowest energy conformation in the gas phase. The results imply that azaAsx and Pro residues may involve cis-trans isomerization in the gas phase. In water, the predicted lowest energy conformer of azapeptides 4 and 5 differs from the gas phase results and depends on the calculational method. For azapeptide 6, regardless of calculation methods and phases, tttAAP-1 (β-I turn) is predicted as the lowest energy conformer. The results imply that the effect of the side chain that can form HBs on the conformational preferences of azapeptides 4 and 5 may not be negligible. We compared the theoretical results of azaXaa-Pro models with those of Pro-azaXaa models, showing that incorporating azaamino acid residue in peptides at different positions can significantly impact the folding patterns and stability of azapeptides.

A fluoro-alkene mimic of Gly-trans-Pro produces a stable collagen triple helix

We report the first experimental evidence for a fluoro-alkene amide isostere participating in n→π* donation, which stabilizes the collagen triple helix. Of the three amide positions in canonical collagen-like peptides, Gly-Pro, Pro-Hyp, and Hyp-Gly, triple helix stability stands to benefit from substitution of only the isomerizable 3° Gly-Pro amide bond with a trans-locked fluoro-alkene. A (Z)-fluoro-alkene isostere of Gly-trans-Pro was synthesized, and its effect on the thermostability of a collagen-like peptide triple helix was measured. The mixture of enantiomers, Boc-Gly-Ψ[(Z)CFC]-L/D-Pro-OH, was synthesized in 8 steps with 27% overall yield, and the Fmoc-Gly-Ψ[(Z)CFC]-L/D-Pro-Hyp-OBn diastereomers were separated. The Gly-Ψ[(Z)CFC]-Pro isostere installed in a collagen-like peptide forms a stable triple helix. By CD, the thermal melting (T_m) value of the fluoro-alkene peptide was +42.2 ± 0.4 °C, and the T_m value of the control peptide was +48.4 ± 0.5 °C, a difference in stability of ΔT_m -6.2 °C. Deshielding of the fluorine nucleus in the ¹⁹F NMR spectra is evidence of a stabilizing n→π* electronic interaction.

Conformational Analysis of Fluoro-, Chloro-, and Proteo-Alkene Gly-Pro and Pro-Pro Isosteres to Mimic Collagen

Collagen is the most abundant human protein, with the canonical sequence (Gly-Pro-Hyp)_n in its triple helix region. Cis-trans isomerization of the Xaa-Pro amide has made two of these amide bonds the target of alkene replacement: the Gly-Pro and the Pro-Hyp positions. The conformations of Gly-Pro and Pro-Pro (as a Pro-Hyp model) fluoro-, chloro-, and proteo-alkene mimic models were investigated computationally to determine whether these alkenes can stabilize the polyproline type II (PPII) conformation of collagen. Second-order Møller-Plesset (MP2) calculations with various basis sets were used to perform the conformational analyses and locate stationary points. The calculation results predict that fluoro- and chloro-alkene mimics of Gly-Pro and Pro-Pro can participate in n→π* donation to stabilize PPII conformations, yet they are poor n→π* acceptors, shifting the global minima away from PPII conformations. For the proteo-alkene mimics, the lack of significant n→π* interactions and unstable PPII-like geometries explains their known destabilization of the triple helix in collagen-like peptides.

Deciphering the Backbone Noncovalent Interactions that Stabilize Polyproline II Conformation and Reduce cis Proline Abundance in Polyproline Tracts

Proline (Pro) has a higher propensity to adopt cis amide geometry than the other natural amino acids, and a poly-Pro (poly-P) tract can adopt either a polyproline I (PPI, all cis amide) or a polyproline II (PPII, all trans amide) helical conformation. Recent studies have revealed a reduced abundance of cis amide geometry among the inner Pro residues of a poly-P tract. However, the forces that stabilize the polyproline helices and the reason for the higher trans amide propensity of the inner Pro residues of a poly-P tract are poorly understood. Herein, we have studied both Pro and non-Pro PPII helical sequences and identified the backbone noncovalent interactions that are crucial to the higher stability of the trans Pro-amide geometry and the preference for a PPII helical conformation. We show the presence of reciprocal CO···CO interactions that extend over the whole PPII helical region. Interestingly, the CO···CO interactions strengthen with the increase in the PPII helical chain length and the inner CO groups possess stronger CO···CO interactions, which could explain the reduced cis abundance of the inner Pro residues of a poly-P tract. We also identified a much stronger (∼0.9 kcal·mol^-1) n_O → σ*_Cα-Cβ interaction between the N-terminal CO oxygen lone pair and the antibonding orbital (σ*) of their C_α-C_β bonds. As the n_O → σ*_Cα-Cβ interaction is possible only in the trans isomers of Pro, this interaction should be crucial for the stabilization of a PPII helix. Finally, an unusual n_N(amide) → σ*_C-N interaction (∼0.3 kcal·mol^-1) was observed between the peptidic nitrogen lone pair (n_N) and the antibonding orbital (σ*_C-N) of the subsequent C-terminal peptide C-N bond. We propose a cumulative effect of these interactions in the stabilization of a PPII helix.

Design of Pseudodiproline Dimers as Mimetics of Pro-Pro Units: Stereocontrolled Synthesis, Configurational Relevance, and Structural Properties

Stereocontrolled methods are described for the synthesis of hitherto unreported pseudodiproline dimers in which a cyclopentane carboxylic acid is linked to a pyrrolidine residue by a stereochemically defined hydroxymethylene tether. These proline-cyclopentane (Pro-Cyp) dimers have interesting structural characteristics as seen in their X-ray crystal structures as well as their nuclear magnetic resonance (NMR) spectra in CDCl₃. They can be considered to be novel Pro-Pro mimetics, which can be used to replace natural diproline sequences with potential applications in medicinal chemistry. They also represent a new concept in the peptidomimetic design of chimeric proline-based amino acids as carbocyclic hydroxyethylene isosteres of inhibitor molecules, in which the stereodefined bridging hydroxyl group can simulate a tetrahedral intermediate in an enzyme complex.

Conformational control of N-methyl-N,N′-diacylhydrazines by noncovalent carbon bonding in solution

Unusual N(amide)⋯C–X noncovalent carbon bonding interactions stabilize the trans–cis (t–c) amide bond rotamers of N-methyl-N,N′-diacylhydrazines over the expected trans–trans (t–t) rotamers in solution.