Ramachandran Plot
From Proteopedia
This page, as it appeared on November 30, 2010, was featured in this article in the journal Biochemistry and Molecular Biology Education.
The Ramachandran plot is a plot of the torsional angles - phi (φ)and psi (ψ) - of the residues (amino acids) contained in a peptide. In sequence order, φ is the C(i-1),N(i),Ca(i),C(i) torsion angle and ψ is the N(i),Ca(i),C(i),N(i+1) torsion angle. The plot was developed in 1963 by G. N. Ramachandran, et. al.[1] by plotting the φ values on the x-axis and the ψ values on the y-axis, as for the image at left[2]. 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 regionsSecondary 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. α-helixThe 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. Planes 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 Ramachandran plot of this peptide has points clustered about the values of φ= -57o and ψ= -47o which are the average values for α-helices. Adding the values 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. β-sheetsDisplay a two segment twisted β-sheet. Draw planes of the peptide bonds. Most β-sheets in globular proteins are twisted sheets which do not have flat parallel pleats. Closer view of β-sheets. The Ramachandran plot of this twisted sheet has points clustered about the values of φ= -130o and ψ= +140o which are the average values for twisted sheets. Adding the values of three other sheet segments more clearly defines the area in which values for twisted sheets are located. Plot regions limited by steric hindranceMost combinations of φ and ψ are sterically forbidden, as illustrated in the tripeptide, Glu-Ser-Ala. With Ser having values of φ = 55o and ψ = -116o 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[3] were developed 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 55o, -116o 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 -57o and -47o which are in the core region in the lower left, the Ser side chain 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.
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. 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 tripeptide Glu-Gly-Ala with Gly having φ and ψ values of +55o and -116o 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[4][5] use separate functions for subsets of the amino acids that have different local steric-hindrance properties, and may even consider the effects of sequence neighbors[6]. 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).Functionally relevant residuesFunctionally relevant residues are more likely than others to have torsion angles that plot to the allowed but disfavored regions of a Ramachandran plot. The specific geometry of these functionally relevant residues, while somewhat energetically unfavorable, may be important for the protein's function, catalytic or otherwise. Such conformations need to be stabilized by the protein using H-bonds, steric packing, or other means, and should very seldom occur for highly solvent-exposed residues. Plots of proteinsView of structure : Secondary structure consists of α-helix, loops and ordered, nonrepetitive structures. Ramachandran plot: Red data points outside of the area expected for α-helix most likely involve residues at the end of the α-helix because often these have angle values that are not typical for α-helix. White points are those for loops and ordered, nonrepetitive structures. The few residues that map to the disallowed region are Gly. Return to model after viewing plot
View of structure: Twisted β-sheet with small segments of α-helix. Ramachandran plot: Most of the yellow points are located in the area for twisted β-sheets where one would expect them, and again the points in the disallowed region are Gly. Return to model after viewing plot.
Ramachandran Plot: Close to equal amounts of α-helix, β-sheet, and ordered, nonrepetitive structures. One important exception to Gly in the disallowed region is Ser:200. Locate this residue that is located in a disallowed region (lower right quadrant). An interesting aspect concerning Ser:200 is that it is one of a triad of residues that are part of the catalytic site and are involved in the catalytic action of this enzyme. The unique φ and ψ values for Ser:200 is the major factor in positioning the side chain so that it can participate in the catalysis. Return to model after viewing plot.
|
|
Other entries in Proteopedia
Interactive Ramachandran plots can be generated for any entry in Proteopedia with the use of a typed Jmol command[7]:
- For example, in a new browser window open the entry in Proteopedia for 1bhl
- If the JSmol panel shows a "Displaying simplified model" message, click on the "load full" orange button below it. Once the model is reloaded,
- Right-click on an empty space of the JSmol panel showing the 3D structure on the page, or click on the JSmol logo (or frank) in the bottom right corner.
- When the menu comes up, select
Console
- Click in the lower text panel of the console that comes up and type the command
Ramachandran
, followed by the return key. - After some processing the Ramachandran plot will be visible and you can hover over and click on the points in the plot just as you can with atoms in a Jmol scene window. (To return to the model, an easy solution is to reload the page or open a new browser instance of that page, or enter into the console
model 1.1
.) With the console window open, the values will be listed as you click on the spheres. - To limit the plot to displaying certain residues or portions of the structure, you can issue commands in the console, such as
display helix
ordisplay gly
. The latter command will limit the plotted display to just glycine residues. In order to return to showing all values on the plot, issue the commanddisplay all
in the Jmol console.
This method to generate interactive Ramachandran plots will also work for other instances of the Jmol applet / JSmol object elsewhere on the web as long as the version of the J(S)mol is 11.4 or greater.
If you just need to report φ and ψ values for a few residues, use the Scene Authoring Tools to select the residues of interest and enter the command draw RAMACHANDRAN
in the console.
Notes
- ↑ Ramachandran GN, Ramakrishnan C, Sasisekharan V (July 1963). "Stereochemistry of polypeptide chain configurations". J. Mol. Biol. 7: 95–9. PMID 13990617
- ↑ Lovell SC, Davis IW, Arendall WB 3rd, de Bakker PI, Word JM, Prisant MG, Richardson JS, Richardson DC. Structure validation by Calpha geometry: phi,psi and Cbeta deviation. Proteins. 2003 Feb 15;50(3):437-50. PMID:12557186 doi:10.1002/prot.10286
- ↑ Laskowski,RA, MacArthur,MW, Moss,DS and Thornton,JM (1993) PROCHECK - a program to check the stereochemical quality of protein structures. J. Appl. Cryst., 26, 283–291
- ↑ Lovell SC, Davis IW, Arendall WB 3rd, de Bakker PI, Word JM, Prisant MG, Richardson JS, Richardson DC. Structure validation by Calpha geometry: phi,psi and Cbeta deviation. Proteins. 2003 Feb 15;50(3):437-50. PMID:12557186 doi:10.1002/prot.10286
- ↑ Read RJ, Adams PD, Arendall WB 3rd, Brunger AT, Emsley P, Joosten RP, Kleywegt GJ, Krissinel EB, Lutteke T, Otwinowski Z, Perrakis A, Richardson JS, Sheffler WH, Smith JL, Tickle IJ, Vriend G, Zwart PH. A new generation of crystallographic validation tools for the protein data bank. Structure. 2011 Oct 12;19(10):1395-412. PMID:22000512 doi:10.1016/j.str.2011.08.006
- ↑ Ting D, Wang G, Shapovalov M, Mitra R, Jordan MI, Dunbrack RL Jr. Neighbor-dependent Ramachandran probability distributions of amino acids developed from a hierarchical Dirichlet process model. PLoS Comput Biol. 2010 Apr 29;6(4):e1000763. PMID:20442867 doi:10.1371/journal.pcbi.1000763
- ↑ Command defined at site for official Jmol documentation
See Also
- Dihedral/Index which has links to many tutorials and resources about the Ramachandran principle, and phi and psi angles in proteins.
- Ramachandran Plot in Wikipedia
- Another example of a Ramachandran Plot showing the different regions. at the European Bioinformatics Institute (EBI)
- Basics of Protein Structure
Proteopedia Page Contributors and Editors (what is this?)
Karl Oberholser, Wayne Decatur, Eran Hodis, Jane S. Richardson, Jaime Prilusky, Alexander Berchansky, Angel Herraez, Norbert Sträter, Joel L. Sussman, Shelly Livne, Eric Martz
DOI: https://dx.doi.org/10.14576/381225.1772128 (?)Citation: Oberholser K, Sussman J L, Hodis E, Decatur W, Livne S, Prilusky J, Richardson J S, Berchansky A, 2013, "Ramachandran Plot",