Ramachandran Plot: Difference between revisions
try change |
phi,psi order; 2 factual corrections; refs |
||
| Line 4: | Line 4: | ||
}} | }} | ||
<applet load='1rnh.pdb' size=' | <applet load='1rnh.pdb' size='500' frame='true' align='right' scene ='Ramachandran_Plots/Plot_1rnh/2' caption='Ramachandran plot of [[1rnh|Ribonuclease H]] ' /> | ||
[[Image:Ramachandran plot general 100K.jpg|thumb|left|240px|Ramachandran plot and contours from 100,000 high-quality general-case datapoints]] | [[Image:Ramachandran plot general 100K.jpg|thumb|left|240px|Ramachandran plot and contours from 100,000 high-quality general-case datapoints]] | ||
The Ramachandran plot is a plot of the torsional angles - [[Psi_and_Phi_Angles|phi (φ)and psi (ψ)]] - of the residues (amino acids) contained in a peptide. In sequence order, φ is the N(i-1),C(i),Ca(i),N(i) torsion angle and ψ is the C(i),Ca(i),N(i),C(i+1) torsion angle. The plot was developed in 1963 by G. N. Ramachandran, et. al.<ref>RAMACHANDRAN GN, RAMAKRISHNAN C, SASISEKHARAN V (July 1963). "Stereochemistry of polypeptide chain configurations". J. Mol. Biol. 7: 95–9. PMID 13990617</ref> by plotting the φ values on the x-axis and the ψ values on the y-axis, as for the image at left<ref>doi:10.1002/prot.10286</ref>. Plotting the torsional angles in this way graphically shows which combination of angles are possible. The torsional angles of each residue in a peptide define the geometry of its attachment to its two adjacent residues by positioning its planar peptide bond relative to the two adjacent planar peptide bonds, thereby the torsional angles determine the conformation of the residues and the peptide. Many of the angle combinations, and therefore the conformations of residues, are not possible because of steric hindrance. By making a Ramachandran plot, protein structural scientists can determine which torsional angles are permitted and can obtain insight into the structure of peptides. The scene on the right is the Ramachandran plot of ribonuclease H. | |||
==Secondary structure plot regions== | ==Secondary structure plot regions== | ||
Secondary structures of a peptide are segments of the peptide that have ordered and repetitive structure, and the repetitive structure is due to a repetitive conformation of the residues and, ultimately, repetitive values of & | Secondary structures of a peptide are segments of the peptide that have ordered and repetitive structure, and the repetitive structure is due to a repetitive conformation of the residues and, ultimately, repetitive values of φ and ψ. The different secondary structures can be distinguished by their range of φ and ψ values with the values of different secondary structures mapping to different regions of the Ramachandran plot. Two common examples of secondary structure are illustrated below. | ||
=== α-helix === | === α-helix === | ||
<table width='500' align='right' cellpadding='5'><tr><td rowspan='2'> </td><td bgcolor='#eeeeee'><applet load='Ahelix_2.pdb' size='490' frame='true' align='right' scene =' Ramachandran_Plots/Helix_first/1' caption='' /></td></tr><tr><td bgcolor='#eeeeee'><center>'''Peptide Segments Illustrating Secondary Structures'''<scene name='Ramachandran_Plots/Helix_first/1'> (Initial scene)</scene></center></td></tr></table> | <table width='500' align='right' cellpadding='5'><tr><td rowspan='2'> </td><td bgcolor='#eeeeee'><applet load='Ahelix_2.pdb' size='490' frame='true' align='right' scene =' Ramachandran_Plots/Helix_first/1' caption='' /></td></tr><tr><td bgcolor='#eeeeee'><center>'''Peptide Segments Illustrating Secondary Structures'''<scene name='Ramachandran_Plots/Helix_first/1'> (Initial scene)</scene></center></td></tr></table> | ||
The scene on the right shows the axis of the α-helix rotating in the y-plane. When viewing the helix on end, observe the open center of the helix. <scene name='Ramachandran_Plots/Helix_planes/2'>Planes</scene> are drawn on some of the peptide bonds to emphasize that in an α-helix the planar peptide bonds rotate about the axis of the helix. The <scene name='User:Karl_Oberholser/Ramachandran_Plots/Plot1/3'>Ramachandran plot</scene> of this peptide has points clustered about the values of & | The scene on the right shows the axis of the α-helix rotating in the y-plane. When viewing the helix on end, observe the open center of the helix. <scene name='Ramachandran_Plots/Helix_planes/2'>Planes</scene> are drawn on some of the peptide bonds to emphasize that in an α-helix the planar peptide bonds rotate about the axis of the helix. The <scene name='User:Karl_Oberholser/Ramachandran_Plots/Plot1/3'>Ramachandran plot</scene> of this peptide has points clustered about the values of φ= -57<sup>o</sup> and ψ= -47<sup>o</sup> which are the average values for α-helices. <scene name='User:Karl_Oberholser/Ramachandran_Plots/Plot2/3'>Adding the values</scene> of two other helical segments demonstrates that data from all three appear in one large cluster and that the helical segments can not be distinguished by the differences in their φ and ψ values. | ||
=== β-sheets === | === β-sheets === | ||
<scene name='Ramachandran_Plots/Sheet_first/1'>Display</scene> a two segment twisted β- | <scene name='Ramachandran_Plots/Sheet_first/1'>Display</scene> a two segment twisted β-sheet. <scene name='User:Karl_Oberholser/Ramachandran_Plots/Opening_sheet/7'>Draw planes</scene> of the peptide bonds. Most β-sheets in globular proteins are twisted sheets which do not have flat parallel pleats. <scene name='User:Karl_Oberholser/Ramachandran_Plots/Opening_sheet2/1'>Closer view</scene> of β-sheets. The <scene name='User:Karl_Oberholser/Ramachandran_Plots/Plot3/1'>Ramachandran plot</scene> of this twisted sheet has points clustered about the values of φ= -130<sup>o</sup> and ψ= +140<sup>o</sup> which are the average values for twisted sheets. <scene name='User:Karl_Oberholser/Ramachandran_Plots/Plot4/1'>Adding the values</scene> of three other sheet segments more clearly defines the area in which values for twisted sheets are located.<br> | ||
{{clear}} | {{clear}} | ||
== Plot regions limited by steric hindrance == | == Plot regions limited by steric hindrance == | ||
<table width='400' align='right' cellpadding='5'><tr><td rowspan='2'> </td><td bgcolor='#eeeeee'><applet load='1MBO' size='390' frame='true' align='right' scene ='Ramachandran_Plots/Tripeptide_disallowed/1' caption='' /></td></tr><tr><td bgcolor='#eeeeee'><center>'''Tripeptides Illustrating the Presence or Absence of Steric Hindrance'''<scene name='Ramachandran_Plots/Tripeptide_disallowed/1'> (Initial scene)</scene></center></td></tr></table> | <table width='400' align='right' cellpadding='5'><tr><td rowspan='2'> </td><td bgcolor='#eeeeee'><applet load='1MBO' size='390' frame='true' align='right' scene ='Ramachandran_Plots/Tripeptide_disallowed/1' caption='' /></td></tr><tr><td bgcolor='#eeeeee'><center>'''Tripeptides Illustrating the Presence or Absence of Steric Hindrance'''<scene name='Ramachandran_Plots/Tripeptide_disallowed/1'> (Initial scene)</scene></center></td></tr></table> | ||
Most combinations of & | Most combinations of φ and ψ are sterically forbidden, as illustrated in the tripeptide, '''<font color=limegreen>Glu</font>-Ser-<font color=blue>Ala</font>'''. With Ser having values of φ = 55<sup>o</sup> and ψ = -116<sup>o</sup> the Ser side chain is in contact with the Ala, colored blue. In plots of native peptides the data points will form clusters in the several areas in which steric hindrance does not occur. These regions are illustrated in the Figure on the left, according to our understanding in the early 1990's when the first validation programs such as ProCheck were developed [[Image:Ramaplot.png|thumb|left|Figure: Ramachandran Plot showing two core regions (blue) and three allowed regions (green).]] [[Image:Ramachandran_general-case_data_and_contours_T8000_small.jpg|thumb|left|240px| General-case Ramachandran plot for >1,000,000 high-quality datapoints]] The core regions (blue in the Figure) contain the most favorable combinations of φ and ψ and contain the greatest number of points. Plots of some proteins contain a small third core region in the upper right quadrant. The allowed regions (green in the Figure) can be located around the core regions or can be unassociated with a core region, but they contain fewer data points than the core regions. The generous regions (not shown in the Figure) extend beyond the allowed regions. The remaining areas are considered disallowed. Observe that the data point 55<sup>o</sup>, -116<sup>o</sup> for Ser of the above tripeptide would fall in the lower right quadrant which contains only a disallowed region in ProCheck. If the Ser in the tripeptide has φ and ψ values of -57<sup>o</sup> and -47<sup>o</sup> which are in the core region in the lower left, the <scene name='Ramachandran_Plots/Tripeptide_allowed/1'>Ser side chain</scene> is rotated away from the Ala and is not in contact with the Ala. Most, if not all, of the points in the above plots for α-helix and β-sheet are in one of the core areas.<p> | ||
Since the 1990's, great expansion in the number and quality of macromolecular crystal structures and advances in methodology have greatly improved understanding of the energetically favored, allowed, and truly disallowed conformations of proteins and nucleic acids. The lower figure plots Ramachandran values for over a million general-case residues with resolution <2.0Å and backbone B-factor <30. Over half the plot is entirely empty, but there are further clusters evident that have moderate population and presumably have somewhat unfavorable energy but are possible; this actually includes the region near +55, -116, which is found in one type of beta turn. | Since the 1990's, great expansion in the number and quality of macromolecular crystal structures and advances in methodology have greatly improved understanding of the energetically favored, allowed, and truly disallowed conformations of proteins and nucleic acids. The lower figure plots Ramachandran values for over a million general-case residues with resolution <2.0Å and backbone B-factor <30. Over half the plot is entirely empty, but there are further clusters evident that have moderate population and presumably have somewhat unfavorable energy but are possible; this actually includes the region near +55, -116, which is found in one type of beta turn. | ||
| Line 25: | Line 25: | ||
=== Glycine, Proline, etc.=== | === Glycine, Proline, etc.=== | ||
Since Gly has only a hydrogen as a side chain, steric hindrance is not as likely to occur as φ and ψ are rotated through a series of values. The <scene name='Ramachandran_Plots/Tripeptide_gly/1'>tripeptide Glu-Gly-Ala</scene> with Gly having φ and ψ values of +55<sup>o</sup> and -116<sup>o</sup> respectively, does not show the steric hindrance that the Glu-Ser-Ala had. For that reason Gly will frequently plot in the disallowed region of a general-case Ramachandran plot. Nearly all of the data points in the disallowed region in the above Figure are Gly points. Therefore modern Ramachandran criteria use separate functions for subsets of the amino acids that have different local steric-hindrance properties. Proline, with the sidechain covalently linked to the preceding backbone N, is more tightly constrained than general-case residues. Residues just before a Pro (called "prePro") have some limits from steric interaction with the proline ring. The branched beta-carbons of Ile and Val also give them a distinct shape of disallowed Ramachandran regions. The other 16 amino-acid types vary in their preference for the very favorable regions, but their outer contours that separate allowed from outlier regions all agree very closely, so they are all grouped together into the "general-case" distribution (as shown above) for structure validation purposes. The figure below shows the plots for Gly and for trans Pro (cis Pro is slightly shifted up and left, and lacks the small central cluster).</p> | Since Gly has only a hydrogen as a side chain, steric hindrance is not as likely to occur as φ and ψ are rotated through a series of values. The <scene name='Ramachandran_Plots/Tripeptide_gly/1'>tripeptide Glu-Gly-Ala</scene> with Gly having φ and ψ values of +55<sup>o</sup> and -116<sup>o</sup> respectively, does not show the steric hindrance that the Glu-Ser-Ala had. For that reason Gly will frequently plot in the disallowed region of a general-case Ramachandran plot. Nearly all of the data points in the disallowed region in the above Figure are Gly points. Therefore modern Ramachandran criteria<ref>PMID:12557186</ref><ref>doi: 10.1016/j.str.2011.08.006</ref> use separate functions for subsets of the amino acids that have different local steric-hindrance properties. Proline, with the sidechain covalently linked to the preceding backbone N, is more tightly constrained than general-case residues. Residues just before a Pro (called "prePro") have some limits from steric interaction with the proline ring. The branched beta-carbons of Ile and Val also give them a distinct shape of disallowed Ramachandran regions. The other 16 amino-acid types vary in their preference for the very favorable regions, but their outer contours that separate allowed from outlier regions all agree very closely, so they are all grouped together into the "general-case" distribution (as shown above) for structure validation purposes. The figure below shows the plots for Gly and for trans Pro (cis Pro is slightly shifted up and left, and lacks the small central cluster).</p> | ||
[[Image:Ramachandran_Gly_Pro_data_and_contours_T8000_small.jpg|thumb|left|480px|Ramachandran plots for Gly and for trans Pro]] | [[Image:Ramachandran_Gly_Pro_data_and_contours_T8000_small.jpg|thumb|left|480px|Ramachandran plots for Gly and for trans Pro]] | ||