Conformational Preferences of the Asparagine Residue. Gas-Phase, Aqueous Solution, and Chloroform Solution Calculations on the Model Dipeptide

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.

Intrinsic conformational characteristics of alpha,alpha-diphenylglycine

Quantum mechanical calculations at the B3LYP/6-31+G(d,p) level have been used to investigate the intrinsic conformational preferences of alpha,alpha-diphenylglycine, a simple alpha,alpha-dialkylated amino acid bearing two phenyl substituents on the alpha-carbon, in both the gas phase and aqueous solution. Nine minimum energy conformations have been characterized for the N-acetyl-N'-methylamide derivative within a relative energy range of about 9 kcal/mol. The relative stability of these structures is largely influenced by specific backbone...side chain and side chain...side chain interactions that can be attractive (N-H...pi and C-H...pi) or repulsive (C=O...pi). On the other hand, comparison with the minimum energy conformations calculated for alpha-aminoisobutyric acid, in which the two phenyl substituents are replaced by methyl groups, revealed that the bulky aromatic rings of alpha,alpha-diphenylglycine induce strain in the internal geometry of the peptide. Finally, a set of force-field parameters for classical Molecular Mechanics calculations was developed for the investigated amino acid. Molecular Dynamics simulations in aqueous solutions have been carried out to validate the parameters obtained.

Toward direct determination of conformations of protein building units from multidimensional NMR experiments VI. Chemical shift analysis of his to gain 3D structure and protonation state information

NMR--chemical shift structure correlations were investigated by using GIAO-RB3LYP/6-311++G(2d,2p) formalism. Geometries and chemical shifts (CSI values) of 103 different conformers of N'-formyl-L-histidinamide were determined including both neutral and charged protonation forms. Correlations between amino acid torsional angle values and chemical shifts were investigated for the first time for an aromatic and polar amino acid residue whose side chain may carry different charges. Linear correlation coefficients of a significant level were determined between chemical shifts and dihedral angles for CSI[(1)H(alpha)]/phi, CSI[(13)C(alpha)]/phi, and CSI[(13)C(alpha)]/psi. Protonation of the imidazole ring induces the upfield shift of CSI[(13)C(alpha)] for positively charged histidines and an opposite effect for the negative residue. We investigated the correspondence of theoretical and experimental (13)C(alpha), (13)C(beta), and (1)H(alpha) chemical shifts and the nine basic conformational building units characteristic for proteins. These three chemical shift values allow the identification of conformational building units at 80% accuracy. These results enable the prediction of additional regular secondary structural elements (e.g., polyProlineII, inverse gamma-turns) and loops beyond the assignment of chemical shifts to alpha-helices and beta-pleated sheets. Moreover, the location of the His residue can be further specified in a beta-sheet. It is possible to determine whether the appropriate residue is located at the middle or in a first/last beta-strand within a beta-sheet based on calculated CSI values. Thus, the attractive idea of establishing local residue specific backbone folding parameters in peptides and proteins by employing chemical shift information (e.g., (1)H(alpha) and (13)C(alpha)) obtained from selected heteronuclear correlation NMR experiments (e.g., 2D-HSQC) is reinforced.

A Dynamic Coupling Model for Sum Frequency Chiral Response from Liquids Composed of Molecules with a Chiral Side Chain and an Achiral Chromophore

A theoretical formulation for optically active sum frequency generation (OA-SFG) from isotropic chiral solutions was proposed for molecules with a chiral side chain and an intrinsically achiral chromophore. Adapting an electron correlation model first proposed by Höhn and Weigang for linear optical activity, we presented a dynamic coupling model for OA-SFG near the electronic resonance of the achiral chromophore. As a demonstration, we used this model to explain the observed OA-SFG spectra of a series of amino acids near the electronic resonance of the intrinsically achiral carboxyl group. Our model shows that the nonlinear chiroptical response comes about by the through-space correlative electronic interactions between the chiral side chain and the achiral chromophore, and its magnitude is determined by the position and orientation of the bonds that make up the chiral side chain. Using the bond polarizability values in the literature and the conformations of amino acids obtained from calculation, we were able to reproduce the relative OA-SFG strength from a series of amino acids.

N-acetyl-N'-methylamide derivative of (2S,3S)-1-amino-2,3-diphenylcyclopropanecarboxylic acid: theoretical analysis of the conformational impact produced by the incorporation of the second phenyl group to the cyclopropane analogue of phenylalanine

The intrinsic conformational preferences of (2S,3S)-1-amino-2,3-diphenylcyclopropanecarboxylic acid, a phenylalanine cyclopropane analogue bearing two phenyl substituents, have been examined theoretically. For this purpose, its N-acetyl-N'-methylamide derivative, Ac-(2S,3S)c(3)diPhe-NHMe, has been investigated by using ab initio HF and DFT methods. Results have been compared with those previously reported for other cyclopropane analogues of phenylalanine, and with experimental data available for c(3)diPhe-containing peptides.

Toward direct determination of conformations of protein building units from multidimensional NMR experiments. V. NMR chemical shielding analysis of N-formyl-serinamide, a model for polar side-chain containing peptides

Knowledge of chemical shift-structure relationships could greatly facilitate the NMR chemical shift assignment and structure refinement processes that occur during peptide/protein structure determination via NMR spectroscopy. To determine whether such correlations exist for polar side chain containing amino acid residues the serine dipeptide model, For-L-Ser-NH(2), was studied. Using the GIAO-RHF/6-31+G(d) and GIAO-RHF/TZ2P levels of theory the NMR chemical shifts of all hydrogen ((1)H(N), (1)H(alpha), (1)H(beta1), (1)H(beta2)), carbon ((13)C(alpha), (13)C(beta), (13)C') and nitrogen ((15)N) atoms have been computed for all 44 stable conformers of For-L-Ser-NH(2). An attempt was made to establish correlation between chemical shift of each nucleus and the major conformational variables (omega(0), phi, psi, omega(1), chi,(1) and chi(2)). At both levels of theory a linear correlation can be observed between (1)H(alpha)/phi, (13)C(alpha)/phi, and (13)C(alpha)/psi. These results indicate that the backbone and side-chain structures of For-L-Ser-NH(2) have a strong influence on its chemical shifts.