In this theoretical study, the role of the side chain moiety of C-terminal residue in influencing the structural and molecular properties of dipeptides is analyzed by considering a series of seven dipeptides. The C-terminal positions of the dipeptides are varied with seven different amino acid residues, namely. Val, Leu, Asp, Ser, Gln, His, and Pyl while their N-terminal positions are kept constant with Sec residues. Full geometry optimization and vibrational frequency calculations are carried out at B3LYP/6-311++G(d,p) level in gas and aqueous phase. The stereo-electronic effects of the side chain moieties of C-terminal residues are found to influence the values of and dihedrals, planarity of the peptide planes, and geometry around the C7？？ -carbon atoms of the dipeptides. The gas phase intramolecular H-bond combinations of the dipeptides are similar to those in aqueous phase. The theoretical vibrational spectra of the dipeptides reflect the nature of intramolecular H-bonds existing in the dipeptide structures. Solvation effects of aqueous environment are evident on the geometrical parameters related to the amide planes, dipole moments, HOMOLUMO energy gaps as well as thermodynamic stability of the dipeptides. 1. Introduction Generally, twenty canonical amino acid residues adequately build up the proteins and enzymes necessary to support most of the cellular functions in all the three domains of life on earth. Selenocysteine (Sec) and pyrrolysine (Pyl) are the two rarely occurring genetically encoded amino acids whose presence in the active sites of some enzymes enables them to sustain life in some extraordinarily unique ways [1–7]. The chemical structures of Sec and Pyl are portrayed in Figure 1. Although the distribution of Pyl is limited to methanogenic archea and certain bacteria, Sec is commonly found in eubacteria, archaea, and eukarya . Sec and Pyl are cotranslationally inserted into proteins corresponding to UGA (opal codon) and UAG codons (canonical stop codon) respectively, which are generally responsible for terminating the process of protein biosynthesis. Figure 1: Chemical structures of selenocysteine (Sec) and pyrrolysine (Pyl). The dynamic properties and functional specificity of the proteins and polypeptides are known to depend primarily on the linear sequence of amino acid residues . Therefore, over the last few decades, small amino acid sequences like di- or tripeptides have been used extensively as model systems in the experimental and theoretical studies concerning the structure of proteins and energetics of protein folding.
I. Chambers, J. Frampton, P. Goldfarb, N. Affara, W. McBain, and P. R. Harrison, “The structure of the mouse glutathione peroxidase gene: the selenocysteine in the active site is encoded by the ‘termination’ codon, TGA,” The EMBO Journal, vol. 5, no. 6, pp. 1221–1227, 1986.
F. Zinoni, A. Birkmann, T. C. Stadtman, and A. Bock, “Nucleotide sequence and expression of the selenocysteine-containing polypeptide of formate dehydrogenase (formate-hydrogen-lyase-linked) from Escherichia coli,” Proceedings of the National Academy of Sciences of the United States of America, vol. 83, no. 13, pp. 4650–4654, 1986.
M. R. Wormald, A. J. Petrescu, Y.-L. Pao, A. Glithero, T. Elliott, and R. A. Dwek, “Conformational studies of oligosaccharides and glycopeptides: complementarity of NMR, X-ray crystallography, and molecular modelling,” Chemical Reviews, vol. 102, no. 2, pp. 371–386, 2002.
N. Foloppe, B. Hartmann, L. Nilsson, and A. D. MacKerell Jr., “Intrinsic conformational energetics associated with the glycosyl torsion in DNA: a quantum mechanical study,” Biophysical Journal, vol. 82, no. 3, pp. 1554–1569, 2002.
J. Sponer, M. Zgarbova, P. Jurecka, K. E. Riley, J. E. Sponer, and P. Hobza, “Reference quantum chemical calculations on RNA base pairs directly involving the 2′-OH group of ribose,” Journal of Chemical Theory and Computation, vol. 5, pp. 1166–1179, 1999.
D. J. Tobiast and C. L. Brooks III, “Conformational equilibrium in the alanine dipeptide in the gas phase and aqueous solution: a comparison of theoretical results,” Journal of Physical Chemistry A, vol. 96, pp. 3864–3874, 1992.
Z.-X. Wang and Y. Duan, “Solvation effects on alanine dipeptide: a MP2/cc-pVTZ//MP2/6-31G** study of (Φ, Ψ) energy maps and conformers in the gas phase, ether, and water,” Journal of Computational Chemistry, vol. 25, no. 14, pp. 1699–1716, 2004.
F. F. García-Prieto, I. F. Galvan, M. A. Aguilar, and M. E. Martín, “Study on the conformational equilibrium of the alanine dipeptide in water solution by using the averaged solvent electrostatic potential from molecular dynamics methodology,” Journal of Chemical Physics, vol. 135, pp. 194502–194509, 2011.
I. R. Gould, W. D. Cornell, and I. H. Hillier, “A quantum mechanical investigation of the conformational energetics of the alanine and glycine dipeptides in the gas phase and in aqueous solution,” Journal of the American Chemical Society, vol. 116, no. 20, pp. 9250–9256, 1994.
T. Head-Gordon, M. Head-Gordon, M. J. Frisch, C. L. Brooks III, and J. A. Pople, “Theoretical study of blocked glycine and alanine peptide analogues,” Journal of the American Chemical Society, vol. 113, no. 16, pp. 5989–5997, 1991.
C. Adamo, V. Dillet, and V. Barone, “Solvent effects on the conformational behavior of model peptides. A comparison between different continuum models,” Chemical Physics Letters, vol. 263, no. 1-2, pp. 113–118, 1996.
G. Das, “Investigations of dipeptide structures containing pyrrolysine as N-terminal residues: a DFT study in gas and aqueous phase,” Journal of Molecular Modeling, vol. 19, no. 4, pp. 1901–1911, 2013.
S. Miertus, E. Scrocco, and J. Tomasi, “Electrostatic interaction of a solute with a continuum. A direct utilizaion of AB initio molecular potentials for the prevision of solvent effects,” Chemical Physics, vol. 55, no. 1, pp. 117–129, 1981.
M. P. Andersson and P. Uvdal, “New scale factors for harmonic vibrational frequencies using the B3LYP density functional method with the triple-ζ basis Set 6-311+G(d,p),” Journal of Physical Chemistry A, vol. 109, no. 12, pp. 2937–2941, 2005.
S. G. Stepanian, I. D. Reva, E. D. Radchenko et al., “Matrix-isolation infrared and theoretical studies of the glycine conformers,” Journal of Physical Chemistry A, vol. 102, no. 6, pp. 1041–1054, 1998.
S. G. Stepanian, I. D. Reva, E. D. Radchenko, and L. Adamowicz, “Conformational behavior of α-alanine. Matrix-isolation infrared and theoretical DFT and ab initio study,” Journal of Physical Chemistry A, vol. 102, no. 24, pp. 4623–4629, 1998.
M. T. Baei and S. Zahra Sayyed-Alangi, “Effects of zinc binding on the structure and stability of glycylglycine dipeptide: a computational study,” E-Journal of Chemistry, vol. 9, no. 3, pp. 1244–1250, 2012.
L. Padmaja, C. Ravikumar, D. Sajan et al., “Density functional study on the structural conformations and intramolecular charge transfer from the vibrational spectra of the anticancer drug combretastatin-A2,” Journal of Raman Spectroscopy, vol. 40, no. 4, pp. 419–428, 2009.
C. Ravikumar, J. I. Hubert, and V. S. Jayakumar, “Charge transfer interactions and nonlinear optical properties of push-pull chromophore benzaldehyde phenylhydrazone: a vibrational approach,” Chemical Physics Letters, vol. 460, pp. 552–558, 2008.