{"id":3802,"date":"2011-04-09T16:26:55","date_gmt":"2011-04-09T16:26:55","guid":{"rendered":"http:\/\/www.ch.imperial.ac.uk\/rzepa\/blog\/?p=3802"},"modified":"2011-04-09T16:26:55","modified_gmt":"2011-04-09T16:26:55","slug":"why-are-%ce%b1-helices-in-proteins-mostly-right-handed","status":"publish","type":"post","link":"https:\/\/rzepa.net\/blog\/2011\/04\/09\/why-are-%ce%b1-helices-in-proteins-mostly-right-handed\/","title":{"rendered":"Why are \u03b1-helices in proteins mostly right handed?"},"content":{"rendered":"<p>Understanding why and how proteins fold continues to be a <em>grand challenge<\/em> in science. I have <a href=\"http:\/\/www.ch.imperial.ac.uk\/rzepa\/blog\/?p=3746\" target=\"_blank\" rel=\"noopener noreferrer\">described<\/a> how Wrinch in 1936 made a bold proposal for the mechanism, which however flew in the face of much of then known chemistry. Linus Pauling took most of the credit (and a Nobel prize) when in a famous paper[cite]10.1073\/pnas.37.4.205[\/cite] in 1951 he suggested a mechanism that involved (<em>inter alia<\/em>) the formation of what he termed \u03b1-helices. Jack Dunitz in 2001[cite]10.1002\/1521-3773(20011119)40:22%3C4167::AID-ANIE4167%3E3.0.CO;2-Q[\/cite] wrote a must-read article[cite]10.fgkwqb[\/cite] on the topic of &#8220;Pauling&#8217;s Left-handed \u03b1-helix&#8221; (it is now known to be right handed).\u00a0I thought I would revisit this famous example with a calculation of my own and here I have used the\u00a0\u03c9B97XD\/6-311G(d,p) DFT procedure[cite]10.1021\/ct100469b[\/cite] to calculate some of the energy components of a small helix comprising (ala)<sub>6 <\/sub>in both left and right handed form.<\/p>\n<p>Firstly, it is important to note that Pauling was apparently not aware of the absolute handedness of amino acids (which are (S) in CIP terminology). This had in fact only been <a href=\"http:\/\/www.ch.imperial.ac.uk\/rzepa\/blog\/?p=3235\" target=\"_blank\" rel=\"noopener noreferrer\">established<\/a> a few months before Pauling&#8217;s publication by Bijvoet[cite]10.1038\/168271a0[\/cite],\u00a0and news of this might not have reached Pauling. So Pauling guessed (or perhaps, he had already built his models, and did not have time to reconstruct them) and his famous\u00a0\u03b1-helix diagram[cite]10.1073\/pnas.37.4.205[\/cite] turned out to be the enantiomer of the real McCoy. As with <a href=\"http:\/\/www.ch.imperial.ac.uk\/rzepa\/blog\/?p=3235\" target=\"_blank\" rel=\"noopener noreferrer\">DNA itself<\/a>, the helix bears a diastereomeric relationship to the chirality of the amino acids; both have to be inverted to get the proper enantiomer (which is what Pauling did). The secret that Pauling had discovered was hydrogen bonding, and particular, weak N-H&#8230;O=C interactions (Wrinch had thought it was strong covalent N-C-OH bonding instead). Of course, there are other effects at work, which include van der Waals or dispersion interactions, electrostatic effects resulting from the large dipoles in peptides (not least due to the zwitterionic character), the <a href=\"http:\/\/www.ch.imperial.ac.uk\/rzepa\/blog\/?p=1174\" target=\"_blank\" rel=\"noopener noreferrer\">planarity of the peptide bond<\/a> itself, the potential for other types of hydrogen bond (<em>e.g.<\/em> C-H&#8230;O) and entropic effects. I have split some of these down for left and right handed forms of DNA in <a href=\"http:\/\/www.ch.imperial.ac.uk\/rzepa\/blog\/?p=3621\" target=\"_blank\" rel=\"noopener noreferrer\">another post<\/a>.<\/p>\n<p>It turns out calculating most of these effects on an even-handed basis is not that easy. Only the recent advent of dispersion-corrected DFT procedures, together with solvation algorithms that allow for accurate geometry optimisation and subsequent evaluation of free energies allows such a calculation to be performed. Hitherto, it has been mostly molecular mechanics that has been used (which itself relies on many parameters from quantum mechanics, such as atom charges, and explicitly\u00a0identifying interactions for hydrogen bonding). By returning to a quantum-mechanical model, some of these assumptions inherent in the mechanics method need not be made.<\/p>\n<p>We\u00a0showed in 1991[cite]10.1039\/P29910000531[\/cite] that an effective solvation treatment required for the zwitterionic form of amino acids in aqueous solutions would ideally comprise not only a self-consistent-reaction-field, but also explicit water molecules as solvent. Here only the former solvation term is included, but expanding the model to include water is certainly possible. Both the zwitterionic and the neutral forms of (ala)<sub>6<\/sub> are included below, so that the effect of a large dipole on the structure and relative helical stability can be estimated. One notes that (even in a dielectric cavity corresponding to water), the extended zwitterions are high energy species. \u00a0In a protein, they of course would be stabilized by the immediate environment of the ions. The right-handed helix clearly comes out as more stable (by about 1 kcal\/mol per residue, see also[cite]10.1021\/ja960665u[\/cite]\u00a0but this is not really due to either dispersion effects or entropy and must therefore arise largely from the hydrogen-bond like interactions. Ionizing the termini to form a zwitterion increases the propensity for a right handed helix slightly.<\/p>\n<table summary=\"Thermodynamic relative energies\" border=\"1\">\n<thead>\n<tr>\n<th colspan=\"6\">Relative thermodynamic energies (kcal mol<sup>-1<\/sup>) of (ala)<sub>6<\/sub> \u03b1-helices<\/th>\n<\/tr>\n<\/thead>\n<tbody>\n<tr>\n<th>System<\/th>\n<th>Total energy<\/th>\n<th>Dispersion<\/th>\n<th>\u0394\u0394H<sub>298<\/sub><\/th>\n<th>\u0394(T.\u0394S<sub>298<\/sub>)<\/th>\n<th>\u0394\u0394G<sub>298<\/sub><\/th>\n<\/tr>\n<tr>\n<td>Left, neutral<\/td>\n<td>0.0<br \/> <!-- -10061.0273200 --><\/td>\n<td>0.0<\/td>\n<td>0.0<\/td>\n<td>0.0<\/td>\n<td><a href=\"http:\/\/hdl.handle.net\/10042\/to-8249\" target=\"_blank\" rel=\"noopener noreferrer\">0.0<\/a><\/td>\n<\/tr>\n<tr>\n<td>Right, neutral<\/td>\n<td>-4.0<br \/> <!-- -10061.0174426 --><\/td>\n<td>+0.2<\/td>\n<td>-4.0<\/td>\n<td>0.9<\/td>\n<td><a href=\"http:\/\/hdl.handle.net\/10042\/to-8250\" target=\"_blank\" rel=\"noopener noreferrer\">-4.9<\/a><\/td>\n<\/tr>\n<tr>\n<td>Left, zwitterion<\/td>\n<td>0.0<\/td>\n<td>0.0<\/td>\n<td>0.0<\/td>\n<td>0.0<\/td>\n<td><a href=\"http:\/\/hdl.handle.net\/10042\/to-8252\" target=\"_blank\" rel=\"noopener noreferrer\">0.0<\/a><\/td>\n<\/tr>\n<tr>\n<td>Right, zwitterion<\/td>\n<td>-7.1<\/td>\n<td>0.1<\/td>\n<td>-6.3<\/td>\n<td>1.7<\/td>\n<td><a href=\"http:\/\/hdl.handle.net\/10042\/to-8251\" target=\"_blank\" rel=\"noopener noreferrer\">-8.0<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>Shown below are the calculated structures. The chains have (inter alia) unusual bifurcated hydrogen-bonding interactions, between one carbonyl group and two N-H groups (show as atom with halo). These are not quite the models that Linus Pauling built!<\/p>\n<table border=\"0\">\n<tbody>\n<tr>\n<td>\n<figure id=\"attachment_3821\" aria-describedby=\"caption-attachment-3821\" style=\"width: 200px\" class=\"wp-caption aligncenter\"><span style=\"line-height: 17px;\"><img loading=\"lazy\" decoding=\"async\" class=\"size-full wp-image-3872\" title=\"left-n\" onclick=\"jmolInitialize('..\/Jmol\/',true);jmolSetAppletColor('yellow');jmolApplet([600,600],'load wp-content\/uploads\/2011\/04\/left-n.xml;zoom 90;select atomno=30;halo on;measure 30 9;measure 25 5;measure 30 14;measure 16 36;');\" src=\"http:\/\/www.ch.imperial.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2011\/04\/left-n.jpg\" alt=\"\" width=\"200\" height=\"145\" \/><\/span><figcaption id=\"caption-attachment-3821\" class=\"wp-caption-text\">Left handed. Click for 3D<\/figcaption><\/figure><\/td>\n<td>\n<figure id=\"attachment_3822\" aria-describedby=\"caption-attachment-3822\" style=\"width: 200px\" class=\"wp-caption aligncenter\"><img loading=\"lazy\" decoding=\"async\" class=\"size-full wp-image-3822\" title=\"right-zw\" onclick=\"jmolInitialize('..\/Jmol\/',true);jmolSetAppletColor('yellow');jmolApplet([600,600],'load wp-content\/uploads\/2011\/04\/right-n.xml;zoom 90;select atomno=50;halo on;measure 61 21;measure 7 31;measure 19 50;measure 29 50;measure 57 48;measure 10 40;');\" src=\"http:\/\/www.ch.imperial.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2011\/04\/right-n.jpg\" alt=\"\" width=\"200\" height=\"137\" \/><figcaption id=\"caption-attachment-3822\" class=\"wp-caption-text\">Right handed. Click for 3D<\/figcaption><\/figure><\/td>\n<\/tr>\n<tr>\n<td>\n<figure id=\"attachment_3868\" aria-describedby=\"caption-attachment-3868\" style=\"width: 200px\" class=\"wp-caption aligncenter\"><img loading=\"lazy\" decoding=\"async\" class=\"size-full wp-image-3868\" title=\"left-zw\" onclick=\"jmolInitialize('..\/Jmol\/',true);jmolSetAppletColor('yellow');jmolApplet([600,600],'load wp-content\/uploads\/2011\/04\/left-zw.xml;zoom 90;select atomno=25;halo on;measure 5 25;measure 9 30;measure 9 25;');\" src=\"http:\/\/www.ch.imperial.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2011\/04\/left-zw.jpg\" alt=\"\" width=\"200\" height=\"122\" \/><figcaption id=\"caption-attachment-3868\" class=\"wp-caption-text\">Left  handed zwitterion. Click for 3D<\/figcaption><\/figure><\/td>\n<td>\n<figure id=\"attachment_3869\" aria-describedby=\"caption-attachment-3869\" style=\"width: 200px\" class=\"wp-caption aligncenter\"><img loading=\"lazy\" decoding=\"async\" class=\"size-full wp-image-3869\" title=\"right-zw\" onclick=\"jmolInitialize('..\/Jmol\/',true);jmolSetAppletColor('yellow');jmolApplet([600,600],'load wp-content\/uploads\/2011\/04\/right-zw.xml;zoom 90;select atomno=50;halo on;measure 40 10;measure 50 29;measure 50 19;');\" src=\"http:\/\/www.ch.imperial.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2011\/04\/right-zw1.jpg\" alt=\"\" width=\"200\" height=\"145\" \/><figcaption id=\"caption-attachment-3869\" class=\"wp-caption-text\">Right  handed zwitterion. Click for 3D<\/figcaption><\/figure><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>For a more objective analysis of the interactions within the system, a QTAIM analysis is shown below.<\/p>\n<table border=\"0\">\n<tbody>\n<tr>\n<td>\n<figure id=\"attachment_3890\" aria-describedby=\"caption-attachment-3890\" style=\"width: 200px\" class=\"wp-caption aligncenter\"><img loading=\"lazy\" decoding=\"async\" class=\"size-full wp-image-3890\" title=\"left-n-aim\" onclick=\"jmolInitialize('..\/Jmol\/',true);jmolSetAppletColor('yellow');jmolApplet([600,600],'load wp-content\/uploads\/2011\/04\/left-n-aim.mol2;zoom 90;');\" src=\"http:\/\/www.ch.imperial.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2011\/04\/left-n-aim.jpg\" alt=\"\" width=\"200\" height=\"138\" \/><figcaption id=\"caption-attachment-3890\" class=\"wp-caption-text\">Left helix. Bond critical points in green. Click for 3D.<\/figcaption><\/figure><\/td>\n<td>\n<figure id=\"attachment_3891\" aria-describedby=\"caption-attachment-3891\" style=\"width: 200px\" class=\"wp-caption aligncenter\"><img loading=\"lazy\" decoding=\"async\" class=\"size-full wp-image-3891\" title=\"right-n-aim\" onclick=\"jmolInitialize('..\/Jmol\/',true);jmolSetAppletColor('yellow');jmolApplet([600,600],'load wp-content\/uploads\/2011\/04\/right-n-aim.mol2;zoom 90;');\" src=\"http:\/\/www.ch.imperial.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2011\/04\/right-n-aim.jpg\" alt=\"\" width=\"200\" height=\"135\" \/><figcaption id=\"caption-attachment-3891\" class=\"wp-caption-text\">Right helix. Click for 3D<\/figcaption><\/figure><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>Whilst the overall conclusion is that theory agrees well with the experimental observation that peptide sequences tend to coil into right rather than left handed helices,\u00a0the reasons they do so is a little more subtle than simple model building alone can reveal.\u00a0As the AIM shows, a plethora of unusual and weaker interactions occur within these helices, a full analysis of which must await presentation elsewhere.<\/p>\n<p>An NCI analysis reveals strong hydrogen bonds as blue-shaded surfaces.<\/p>\n<figure id=\"attachment_4214\" aria-describedby=\"caption-attachment-4214\" style=\"width: 252px\" class=\"wp-caption aligncenter\"><img loading=\"lazy\" decoding=\"async\" class=\"size-full wp-image-4214\" title=\"a-helix\" onclick=\"jmolInitialize('..\/Jmol\/',true);jmolSetAppletColor('white');jmolApplet([600,600],'load wp-content\/uploads\/2011\/04\/right-zw.xyz;isosurface wp-content\/uploads\/2011\/04\/right-zw.jvxl colorscheme translucent bgyor;');\" src=\"http:\/\/www.ch.imperial.ac.uk\/rzepa\/blog\/wp-content\/uploads\/2011\/04\/a-helix.jpg\" alt=\"\" width=\"252\" height=\"172\" \/><figcaption id=\"caption-attachment-4214\" class=\"wp-caption-text\">NCI surface. Click for 3D.<\/figcaption><\/figure>\n","protected":false},"excerpt":{"rendered":"<p>Understanding why and how proteins fold continues to be a grand challenge in science. I have described how Wrinch in 1936 made a bold proposal for the mechanism, which however flew in the face of much of then known chemistry. Linus Pauling took most of the credit (and a Nobel prize) when in a famous [&hellip;]<\/p>\n","protected":false},"author":2,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[11],"tags":[132,198,552,647,788,902,906,1179,1188,1346,1381,1995,2280,2489,2562],"class_list":["post-3802","post","type-post","status-publish","format-standard","hentry","category-interesting-chemistry","tag-alpha-helix","tag-aqueous-solutions","tag-chiroptical","tag-conformational-analysis","tag-dielectric","tag-energy","tag-energy-components","tag-high-energy-species","tag-historical","tag-jack-dunitz","tag-julia-contreras-garcia","tag-protein","tag-solvation-algorithms","tag-tutorial-material","tag-watoc11"],"_links":{"self":[{"href":"https:\/\/rzepa.net\/blog\/wp-json\/wp\/v2\/posts\/3802","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/rzepa.net\/blog\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/rzepa.net\/blog\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/rzepa.net\/blog\/wp-json\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/rzepa.net\/blog\/wp-json\/wp\/v2\/comments?post=3802"}],"version-history":[{"count":0,"href":"https:\/\/rzepa.net\/blog\/wp-json\/wp\/v2\/posts\/3802\/revisions"}],"wp:attachment":[{"href":"https:\/\/rzepa.net\/blog\/wp-json\/wp\/v2\/media?parent=3802"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/rzepa.net\/blog\/wp-json\/wp\/v2\/categories?post=3802"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/rzepa.net\/blog\/wp-json\/wp\/v2\/tags?post=3802"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}