Conformational preference of azaglycine-containing dipeptides studied by PCM and IPCM methods

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.

Peptide-based capsules with chirality-controlled functionalized interiors - rational design and amplification from dynamic combinatorial libraries

Self-assembled capsules were designed by two approaches and synthesized using azapeptides of various chiralities attached to resorcin[4]arene by semicarbazone linkers.

Peptides are commonly perceived as inapplicable components for construction of porous structures. Due to their flexibility the design is difficult and shape persistence of such putative structures is diminished. Notwithstanding these limitations, the advantages of peptides as building blocks are numerous: they are functional and functionalizable, widely available, diverse and biocompatible. We aimed at the construction of discrete porous structures that exploit the inherent functionality of peptides by an approach that is inspired by nature: structural pockets are defined by the backbones of peptides while functionality is introduced by their side chains. In this work peptide ribbons were preorganized on a macrocyclic scaffold using azapeptide–aldehyde reactions. The resulting cavitands with semicarbazone linkers arrange the peptide backbones at positions that are suitable for self-assembly of dimeric capsules by formation of binding motifs that resemble eight-stranded β-barrels. Self-assembly properties and inside/outside positions of the side chains depend crucially on the chirality of peptides. By rational optimization of successive generations of capsules we have found that azapeptides containing three amino acids in a (l, d, d) sequence give well-defined dimeric capsules with side chains inside their cavities. Taking advantage of the reversibility of the reaction of semicarbazone formation we have also employed the dynamic covalent chemistry (DCC) for a combinatorial discovery of capsules that could not be rationally designed. Indeed, the results show that stable capsules with side chains positioned internally can be obtained even for shorter sequences but only for combination peptides of (l, l) and (d, l) chirality. The hybrid (l, l)(d, l) capsule is amplified directly from a reaction mixture containing two different peptides. All capsules gain substantial ordering upon self-assembly, which is manifested by a two orders of magnitude increase of the intensity of CD spectra of capsules compared with non-assembled analogs. Temperature-dependent CD measurements indicate that the capsules remain stable over the entire temperature range tested (20–100 °C). Circular dichroism coupled with TD DFT calculations, DOSY measurements and X-ray crystallography allow for elucidation of the structures in the solid state and in solution and guide their iterative evolution for the current goals.

β-Turn structure in glycinylphenylalanine dipeptide based N-amidothioureas

Transforming the C-terminal amide of a glycinylphenylalanine dipeptide into N-amidothiourea affords a β-turn structure in the formed dipeptide based N-amidothioureas, which can be readily identified by an induced CD signal from the achiral phenylthiourea chromophore.

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.

Azapeptides and their therapeutic potential

Azapeptides are peptide analogs in which one or more of the amino residues is replaced by a semicarbazide. This substitution of a nitrogen for the α-carbon center results in conformational restrictions, which bend the peptide about the aza-amino acid residue away from a linear geometry. The resulting azapeptide turn conformations have been observed by x-ray crystallography and spectroscopy, as well as predicted based on computational models. In biologically active peptide analogs, the aza-substitution has led to enhanced activity and selectivity as well as improved properties, such as prolonged duration of action and metabolic stability. In light of these characteristics, azapeptides have found important uses as receptor ligands, enzyme inhibitors, drugs, pro-drugs, probes and imaging agents. Recent improvements in synthetic methods for their procurement have ushered in a new era of azapeptide chemistry. This review aims to provide a historical look at the development of azapeptide science along with a focus on recent developments and perspectives on the future of this useful tool for medicinal chemistry.

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.