Anales de la Asociación Química Argentina
Print version ISSN 0365-0375
An. Asoc. Quím. Argent. vol.94 no.1-3 Buenos Aires Jan./July 2006
Electronic factors favouring the cis conformation in proline peptidic bonds
Rivail,* J. L. ; Bouchy, A.; Loos, P. F.
Equipe de Chimie et Biochimie théoriques, UMR 7565 CNRS-Université Henri Poincaré B.P.239, 54506 Nancy-Vandoeuvre, France.
Fax : +33 383 68 43 71, E-Mail: Jean-Louis.Rivail@cbt.uhp-nancy.fr
Received March 22nd, 2006. In final form June 4th, 2006
Dedicated to Prof. Imre G. Csizmadia on the occasion of his 75th birthday
A quantum chemical study based on B3LYP, 6-31+G* computations offered the possibility of comparing three dipeptides containing as a central residue proline, N ethyl alanine and an alanine respectively in their cis and trans conformations. The presence of the electron donating alkyl substituent on the nitrogen atom is an important factor favouring the cis conformation. The stabilisation of the cis conformer is interpreted as being due to electrostatic and induction interactions between the carbonyl group and the substituent on the nitrogen atom.
Un estudio químico cuántico basado en cálculos computacionales B3LYP, 6.31+G* permitió la comparación de tres dipéptidos que contienen como residuo central prolina, N etil alanina y alanina, respectivamente, en sus conformaciones cis y trans. La presencia de un sustituyente alquilo donador de electrones sobre el átomo de nitrógeno es un factor importante que favorece la conformación cis. La estabilización del confórmero cis es debida a interacciones electrostáticas y de inducción entre el grupo carbonilo y el sustituyente sobre el átomo de nitrógeno.
Protein folding is one of the major problems that chemistry and biochemistry have to face at the beginning of the 21st century. Among the information to be gained to rationalize the phenomenon, peptidic conformation plays a central role, and I.G. Csizmadia and coworkers made an important contribution to the understanding of what may influence the three dimensional structure of a protein at the level of the peptidic bonds . Although the great majority of the peptidic bonds in proteins adopt the trans conformation, it happens that a few of them take the other possible structure: the cis conformation . This occurrence is particularly frequent when the bond involves a proline residue bonded at its amino moiety . This remark suggests that, in addition to the many factors which may be responsible for these unusual structures, electrostatic or steric interactions with the surrounding parts of the protein, proline may have some intrinsic properties which favour such conformations. Some theoretical investigations have already been done on the conformations of the proline peptidic bond [4-6], which confirm the observed tendency, but to our knowledge, no convincing explanation has been proposed to rationalise this result. The aim of this study is essentially a qualitative investigation of the electronic effects which may help understanding the special role of proline in proteins. Several effects may be invoked to explain this particular behaviour compared to the other aminoacids. One immediately thinks of a possible modification of the electronic properties of the nitrogen atom induced by the five membered ring which constrains the bond angles and also of a possible modification of the local electronic distribution due to the electron donating properties of the (cyclo) alkyl group bonded to the nitrogen atom. Interactions between non bonded atoms may also occur. The flexibility of the five membered ring does not seem to be a major factor since there are several examples of cis-proline bonds in which the conformation of the ring is changed. This is the case, for instance, of the cis peptidic bonds between Tyr 92 and Pro 93 or Asn 113 and Pro 114 in ribonuclease A .
In order to gain some insight in the cis-proline conformation, we performed a series of computations on a model molecule with a central glycine-proline peptidic bond and compared the results with those obtained when proline is replaced by alanine and by the hypothetical N ethyl alanine which mimics the influence of the alkyl chain free of ring constrains. The reference model system is the N methyl glycine–proline methylamide (scheme 1).
Scheme 1: Structure of the dipeptides (trans conformers)
The computations have been performed with the Gaussian 03 set of programs . Owing to the aim of this study, an intermediate computational level has been adopted: DFT formalism with the B3LYP exchange-correlation functional  and the 6-31+G* standard basis set. The geometries are fully optimized and the nature of the stationary points assessed by computing the vibrational frequencies. The electronic population is analysed by means of Natural Bond Orbitals (NBO)  and Mayer bond orders .
Conformations and energies
The global minimum of the proline dipeptide is the trans conformation of the peptidic bond, in agreement with a previous result on a dipeptide having a slightly different composition . But the cis conformer, represented on figure 1, is less than 1 kcal/mol above the global minimum (figure 2). This is also true when proline is replaced by N ethyl alanine. On the contrary, the difference in favour of the trans conformer is noticeably greater in the case of alanine (more than 5 kcal/mol). The conformation of the chain is characterized by the angles w, j, y given on table 1. For both conformations, the backbones of the three peptides are rather similar. The only noticeable difference in the case of the cis proline is the orientation of the CO bond of the second carbonyl group, not belonging to the peptidic bond of interest, which is characterized by the y angle, negative in the case of cis proline and positive in both other cis conformers.
Figure 1: Three dimensional structure of the cis-proline dipeptide
Figure 2: Three dimensional structure of the trans-proline dipeptide
Table 1: Main geometrical features of the dipeptides in cis and trans conformations
This twist of ca. 60° may be due to constrains introduced by the five membered cycle. All the attempts to localize another minimum with a negative y angle failed, and similarly, the attempts to find another conformation of the other cis conformers with a positive y angle were unsuccessful. The j angle indicates a flattening of the cycle in both conformations of proline: 21.52° and 8.24° for the cis and trans conformations respectively. For comparison, in tetrahydropyrole (the cycle alone), the corresponding value is 44.23°. Finally, intramolecular hydrogen bonds occur in both conformations: between the hydrogen atom bonded to nitrogen 10 and N1 in the cis conformations and between this hydrogen atom and O4 in the trans conformation. It is worth noticing that the latter hydrogen bond is shorter in the trans proline dipeptide than in both other trans molecules. This hydrogen bond, obviously stronger in the trans than in the cis conformation may contribute to the relative stability of the trans proline conformer. In proteins, these hydrogen bonds sometimes cannot be formed and this is probably an extra factor favouring the cis conformation.
The most significant data regarding the Mayer bond indices and the charges are collected in tables 2 to 4, for both conformations of the three dipeptides. In each case, the changes are not very important from one conformation to the other. The most characteristic feature is the charge on the methylene group bonded to the nitrogen atom in the case of proline and N ethyl alanine (position 11) which is strongly positive, especially in the case of cis proline (0.245), and greater than that for tetrahydropyrole in which it is only 0.163. A consequence of this intramolecular electron transfer is an increase of the negative charge on the oxygen atom.
Tabla 2: Main characters of the electronic structure of the proline dipeptide conformers
Table 3: Main characteristics of the electronic structure of the N ethyl alanine dipeptide conformers
Table 4: Main characteristics of the electronic structure of the alanine dipeptide conformers
The above results clearly show that the properties of the backbone do not vary substantially neither from one dipeptide to another nor from the cis to the trans conformations. The presence of an electron donor substituent on the nitrogen atom favours the cis conformer, but the effect of this substituent is not sufficient to make the cis conformer more stable than the trans form. The comparison of the proline dipeptide and the N ethyl alanine one shows that the positive charge on the methylene group at position 11 is greater for proline than for N ethyl alanine in both conformations. This is probably linked to the fact that the cycle is a better electron donor than the ethyl group. One also notices that, in the case of proline, this positive charge is greater in the cis conformation than in the trans form. If now one consider the sum of the charges of the cycle: N plus (CH)6 plus the three methylene groups of the cycle, one finds -0.043 for the cis conformation and -0.042 for the trans form, so that, although the charge of the methylene group at position 11 is greater in the cis conformation than in the trans position, the electron transfer does not vary noticeably from one conformation to the other. One is led to explain the main difference between these charge distributions by and induced polarization which is greater for the cis conformation due to the neighbourhood of the negative oxygen atom. This phenomenon also occurs in the conformers of the N ethyl alanine dipeptide, but there may be some additional effects due to the flexibility of the ethyl group. For instance, the distance O4-C11 is 2.729 Å in the cis proline dipeptide and 2.685 Å in the parent N ethyl alanine compound. Nevertheless the polarization is larger in the first case and obviously contributes more efficiently to stabilize the cis conformer. Finally, the main factor which favours the cis conformer in the proline derivative appears to be the sum of the electrostatic and induction interactions.
An analysis of the structural features of the cis and trans conformations of the three dipeptides of interest clearly shows the major role of the alkyl substituent of the nitrogen atom which reduces the energy difference between both conformer to such an extent that the cis conformer becomes more stable than the trans position. This effect is obviously linked to the electron donating properties of this substituent. The consequence is that the methylene group linked to the nitrogen atom is positively charged. In the cis conformation, this group is close to the carbonyl oxygen atom so that an electrostatic interaction is a factor which stabilizes this conformation. This effect is enhanced by an additional induction interaction. Altogether these effects stabilize the cis conformation by ca. 5 kcal/mol with respect to the trans form, if one refers to a similar compound built with an alanine residue. This is not negligible and is obviously one factor explaining the rather frequent occurrence of cis proline conformations in proteins. But it may not be the only explanation and probably, during protein synthesis, some additional effects, such as steric constrains, may occur. The importance of a substitution on the nitrogen atom may help to design artificial amino acids for synthesising peptides with specific properties.
 Lang, A.; Csizmadia I.G. and Perczel, A., Proteins: Structure function and genetics, 2005, 58, 571. [ Links ]
 See for instance: Sawaya M.R.; Loop J., Biochemistry 1997, 36, 586. [ Links ]
 Wedemeyer, W.J.; Welker E.; Scheraga H.A., Biochemistry, 2002, 41, 14637. [ Links ]
 Sahai, M.A.; Kehoe, T.A.K.; Koo, J.C.P.; Setiadi, D.H.; Chass, G.A.; Viskolcz, B.; Penke, B.; Pai E.F.; Csizmadia, I.G., J. Phys. Chem. A 2005, 109, 2660. [ Links ]
 Kang Y.K., J. Phys. Chem. B 2002, 106, 2074. [ Links ]
 Kang Y.K., J. Mol. Struct. (Theochem) 2004, 675, 37. [ Links ]
 Berisio, R.; Sica, F.; Lamzin, V.S.; Wilson, K.S.; Zagari A.; Mazzarella L., Acta Crystallogr., Sect.D 2002, 58, 441. [ Links ]
 Gaussian 03, Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Montgomery, J. A., Jr.; Vreven, T.; Kudin, K. N.; Burant, J. C.; Millam, J. M.; Iyengar, S. S.; Tomasi, J.; Barone, V.; Mennucci, B.; Cossi, M.; Scalmani, G.; Rega, N.; Petersson, G. A.; Nakatsuji, H.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Klene, M.; Li, X.; Knox, J. E.; Hratchian, H. P.; Cross, J. B.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Ayala, P. Y.; Morokuma, K.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Zakrzewski, V. G.; Dapprich, S.; Daniels, A. D.; Strain, M. C.; Farkas, O.; Malick, D. K.; Rabuck, A. D.; Raghavachari, K.; Foresman, J. B.; Ortiz, J. V.; Cui, Q.; Baboul, A. G.; Clifford, S.; Cioslowski, J.; Stefanov, B. B.; Liu, G.; Liashenko, A.; Piskorz, P.; Komaromi, I.; Martin, R. L.; Fox, D. J.; Keith, T.; Al-Laham, M. A.; Peng, C. Y.; Nanayakkara, A.; Challacombe, M.; Gill, P. M. W.; Johnson, B.; Chen, W.; Wong, M. W.; Gonzalez, C.; Pople, J. A. Gaussian 03, Revision B.05 2003, Gaussian, Inc., Pittsburgh PA. [ Links ]
 Carpenter J.E.; Weinhold, F., J. Mol. Struct. (Theochem) 1988, 46, 41. [ Links ]; Foster J. P.; Weinhold, F., J. Am. Chem. Soc. 1980, 102, 7211-8; [ Links ] Reed A. E.; Weinhold, F., J. Chem. Phys. 1983, 78, 4066-73; [ Links ] Reed A.E.; Weinhold, F., J. Chem. Phys. 1985, 83,1736-40; [ Links ] Reed, A.E.; Weinhold R.B.; Weinhold, F., J. Chem. Phys 1985, 83,735-46; [ Links ] Reed, A. E.; Curtiss, L. A.; Weinhold, F. Chem. Rev. 1988, 88, 899. [ Links ]
 Mayer, I., Chem. Phys. Lett. 1983, 97, 270. ; Int. J. Quantum Chem. 1986, 29, 73, 477. [ Links ]