Computational study of conformational preferences of thioamide-containing azaglycine peptides

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.

Induction of chirality: experimental evidence of atropisomerism in azapeptides

Methylation of the peptide bond in model azadipeptides leads to the E configuration and hence to atropisomerism due to a restricted rotation around the N-N axis.

The β-turn scaffold of tripeptide containing an azaphenylalanine residue

The conformational preferences of azaphenylalanine-containing peptide were investigated using a model compound, Ac-azaPhe-NHMe with ab initio method at the HF/3-21G and HF/6-31G(*) levels, and the seven minimum energy conformations with trans orientation of acetyl group and the 4 minimum energy conformations with cis orientation of acetyl group were found at the HF/6-31G(*) level if their mirror images were not considered. An average backbone dihedral angle of the 11 minimum energy conformations is phi=+/-91 degrees +/-24 degrees , psi =+/-18 degrees +/-10 degrees (or +/-169 degrees +/-8 degrees ), corresponding to the i+2 position of beta-turn (delta(R)) or polyproline II (beta(P)) structure, respectively. The chi(1) angle in the aromatic side chain of azaPhe residue adopts preferentially between +/-60 degrees and +/-130 degrees, which reflect a steric hindrance between the N-terminal carbonyl group or the C-terminal amide group and the aromatic side chain with respect to the configuration of the acetyl group. These conformational preferences of Ac-azaPhe-NHMe predicted theoretically were compared with those of For-Phe-NHMe to characterize the structural role of azaPhe residue. Four tripeptides containing azaPhe residue, Boc-Xaa-azaPhe-Ala-OMe [Xaa=Gly(1), Ala(2), Phe(3), Asn(4)] were designed and synthesized to verify whether the backbone torsion angles of azaPhe reside are still the same as compared with theoretical conformations and how the preceding amino acids of azaPhe residue perturb the beta-turn skeleton in solution. The solution conformations of these tripeptide models containing azaPhe residue were determined in CDCl(3) and DMSO solvents using NMR and molecular modeling techniques. The characteristic NOE patterns, the temperature coefficients of amide protons and small solvent accessibility for the azapeptides 1-4 reveal to adopt the beta-turn structure. The structures of azapeptides containing azaPhe residue from a restrained molecular dynamics simulation indicated that average dihedral angles [(phi(1), psi(1)), (phi(2), psi(2))] of Xaa-azaPhe fragment in azapeptide, Boc-Xaa-azaPhe-Ala-OMe were [(-68 degrees, 135 degrees ), (116 degrees, -1 degrees )], and this implies that the intercalation of an azaPhe residue in tripeptide induces the betaII-turn conformation, and the volume change of a preceding amino acid of azaPhe residue in tripeptides would not perturb seriously the backbone dihedral angle of beta-turn conformation. We believe such information could be critical in designing useful molecules containing azaPhe residue for drug discovery and peptide engineering.

Origin of Rotational Barriers of the N−N Bond in Hydrazine: NBO Analysis

Hydrazine passes through two transition states, TS1 (phi = 0 degrees ) and TS2 (phi = 180 degrees ), in the course of internal rotation around its N-N bond. The origin of the corresponding rotational barriers in hydrazine has been extensively studied by experimental and theoretical methods. Here, we used natural bond orbital (NBO) analysis and energy decomposition of rotational barrier energy (DeltaE(barrier)) to understand the origin of the torsional potential energy profile of this molecule. DeltaE(barrier) was dissected into structural (DeltaE(struc)), steric exchange (DeltaE(steric)), and hyperconjugative (DeltaE(deloc)) energy contributions. In both transition states, the major barrier-forming contribution is DeltaE(deloc). The TS2 barrier is lowered by pyramidalization of nitrogen atoms through lowering DeltaE(struc), not by N-N bond lengthening through lowering DeltaE(steric). Higher pyramidality of nitrogen atoms of TS2 than that of TS1 explains well why the N-N bond of TS2 is longer than that of TS1. Finally, the steric repulsion between nitrogen lone pairs does not determine the rotational barrier; nuclear-nuclear Coulombic repulsion between outer H/H atoms in TS1 plays an important role in increasing DeltaE(struc). Taken together, we explain the reason for the different TS1 and TS2 barriers. We show that NBO analysis is a useful tool for understanding structures and potential energy surfaces of compounds containing the N-N bond.

A new scheme for determining the intramolecular seven-membered ring N–H⋯OC hydrogen-bonding energies of glycine and alanine peptides

In this paper a new scheme was proposed to calculate the intramolecular hydrogen-bonding energies in peptides and was applied to calculate the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies of the glycine and alanine peptides. The density-functional theory B3LYP6-31G(d) and B3LYP6-311G(d,p) methods and the second-order Moller-Plesset perturbation theory MP26-31G(d) method were used to calculate the optimal geometries and frequencies of glycine and alanine peptides and related structures. MP26-311++G(d,p), MP26-311++G(3df,2p), and MP2/aug-cc-pVTZ methods were then used to evaluate the single-point energies. It was found that the B3LYP6-31G(d), MP26-31G(d), and B3LYP6-311G(d,p) methods yield almost similar structural parameters for the conformers of the glycine and alanine dipeptides. MP2/aug-cc-pVTZ predicts that the intramolecular seven-membered ring N-H...O=C hydrogen-bonding strength has a value of 5.54 kcal/mol in glycine dipeptide and 5.73 and 5.19 kcal/mol in alanine dipeptides, while the steric repulsive interactions of the seven-membered ring conformers are 4.13 kcal/mol in glycine dipeptide and 6.62 and 3.71 kcal/mol in alanine dipeptides. It was also found that MP26-311++G(3df,2p) gives as accurate intramolecular N-H...O=C hydrogen-bonding energies and steric repulsive interactions as the much more costly MP2/aug-cc-pVTZ does.