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.

Conformational preferences of Ac-Pro-azaXaa-NHMe (Xaa = Asn, Asp, Ala) and the effect of intramolecular hydrogen bonds on their stability in gas phase and solution.

The conformational preferences of three azadipeptides Ac-Pro-azaXaa-NHMe [Xaa = Asn (1), Asp (2), Ala (3)] have been carried out in gas phase and solution (water) using the density functional method B3LYP/6-311 + + G(d,p) to explore the effect of the change of side chain of azaamino acids at the i + 2 position on the stability of these components. The most stable conformations of compounds (1), (2), and (3) have an amid bond oriented trans, trans, and cis, respectively, in gas phase, whereas the orientation of amid bond in water solvent of compounds (2) and (3) has changed to cis and trans, respectively. We have also noticed the importance of backbone-side chain hydrogen bonds in the stabilization of the ß turn motif in gas phase since this motif is more stable in the case of compounds (1) and (2) and less stable in the case of compound (3) in which these hydrogen bonds are absent. Furthermore, the ßII(ßII') turn structure is more stable than ßI turn for all conformations of the three compounds in gas phase, while it is not true in the case of some conformations in solution. Moreover, the stability of ß turn increases from azaAsn to azaAsp which could be due to the side chain's basic nature of azaAsn. In general, hydrogen bonds were found to play a key role in the stabilization of these compounds since most of conformers are lower in energy when they have more than two hydrogen bond interactions while conformations with one or no hydrogen bonds are higher in energy and thus less stable.