Calculate structure
From Proteopedia
An important part of protein structure is the secondary structure which is made up of helices, sheets and turns, and Jmol is capable of determining and displaying these three types of structures with limitations as described in How Jmol Determines Secondary Structure . The calculate structure[1] is a command which does a more fundamental identification of these secondary structures by re-calculating the secondary structure, but it is not available in Jmol 11.8 which is used in Proteopedia as of June 2011. It is available in Jmol ver. 12. Calculate hbonds structure is also available in ver. 12, and it identifies and displays the hbonds involved in these three types of secondary structures[1].
Any page of Proteopedia can be run in the signed ver. 12 by appending "?JMOLJAR=http://chemapps.stolaf.edu/jmol/docs/examples-12/JmolAppletSigned0.jar" to the url of the page and reloading the page. The user must give permission for the signed version of Jmol to open, and when it does it has a red frank, whereas in the unsigned version it is grey. Click on the Jmol frank, in the main menu click on Console, in the bottom box of the console enter the commands:select protein; calculate structure; cartoon; color structure; calculate hbonds structure and then click Run.
The objectives of this article are:
- Describe briefly how calculate structure identifies secondary structures, with a focus on turns, and relate its turn identification to β and γ-turns.
- Summarize the results of using calculate structure to identify turns in two proteins.
- Show details of the above identifications.
Basis of Secondary Structure DeterminationCalculate structure is based on Defined Secondary Structure of Protein (DSSP), a program written in Pascal.[2] The secondary structure recognition algorithms used in DSSP are based mainly on hydrogen-bonding patterns along with some geometric structures , such as bends. There are two different hydrogen-bonding patterns which are recognized. The one determines the value of n in the expression i + n (i is a residue that forms a hydrogen bond with a residue n residues removed from residue i.) where n = 3, 4 or 5. These values define three types of turns. A peptide segment that has repeating turns of the same type are called 310-helix, α-helix, or Π-helix, respectively. If the turn is isolate, it is simply called an n-turn. The other recognized pattern is a hydrogen bond which is between residues which are not close together in sequence. This type of hydrogen bond is called a bridge. Kabsch & Sanders define a ladder as a "set of one or more consecutive bridges of identical type" and a sheet as a "set of one or more ladders connected by shared residues"[2]. Bends are peptide segments with high curvature, and the determination of curvature involves angles of the Cα. Bends can overlap with helices and turns. After its completion the results of calculate structure computation are printed in the upper box of the console. One part of that output is a summary of the different types of secondary structure with each type having a one letter identifier. (The summary for myohemerytherin (2mhr) is given below with a key for the one letter identifier.) It is possible for a residue or a segment of residues to be assigned more than one structural type, and for this reason the key for the identifiers is rank ordered in decreasing priority of assignment. With bend having the lowest priority in assignment a structure is identified as a bend only if it is not assigned any other structural type. The helices and sheets which are identified on the summary can easily be associated with the corresponding structures in the applet, but the turns need some additional explanation. Relationship to β and γ TurnsThe DSSP determination of helices and β-sheets is in agreement with the generally accepted view of these two structures, but the DSSP determination of turns is not as specific as the generally accepted definition of turns. As described above DSSP identifies turns that have 4, 5, or 6 residues with a backbone hbond being present between the first and the last residues. The presence of the hbond is a requirement to be classified as a turn. Phi and psi torsional angles of the Cα are not used by the DSSP procedure, but the generally accepted definitions of β and γ turns involve these angles. β-turns contain four residues and therefore are 3-turns found by DSSP. The classes of β-turns are defined by the range of psi and phi values for the second and third residues.[3] β-turns often have a hbond between residues one and four (i + 3) of β-turns, but there is not an absolute requirement for one. In three classes (VIa1, VIa2, VIb) a Pro in the third position has the cis configuration which does not permit the formation of a hbond (View display of structure.). The turns in these three classes are not detected by DSSP since they do not contain a hbond. γ-turns contain three residues having a hbond between residues i and i + 2 and therefore are not included among the turns found by DSSP. The classic γ-turns have phi and psi values at residue i + 1 of +75.0 and -64, respectively, and the inverse γ-turns have phi and psi values at residue i + 1 of -79 and +69, respectively.[4] Summary of the use of Calculate structure to identify turnsTwo proteins is a small sample, but it does give some indication of the nature of the T: segments (turns) reported in the summary and of the trace segments colored blue.
IllustrationsThe user is urged to use the above directions to open Jmol version 12 and to run the calculate structure and the accompanying commands so that the resulting display can be compared with the summary below. Without displaying the images generated by calculate structure and calculate hbonds structure the activities and comparisons described below can not be performed. Unless a green link is designed to change the color and structural representation these two display characteristics will not change after they have been set by calculate structure, but all hbonds are deleted by clicking a green link so calculate hbonds structure has to be run from the console after every green link click in order to display hbonds. Myohemerytherin ()
SUMMARY for Myohemerytherin: Domain 2 of chain A Glycogen Phosphorylase () - If the applet is not running the signed ver. 12 of Jmol, connect with it as you did above, and then click on the above green link.
SUMMARY of T's for Domain 2 of Chain A Glycogen Phosphorylase:(All other segments deleted.) 486-490,493-497,726-730,771-775,776-779,804-809; 524-527,593-596,610-613,633-636,636-639,668-671,675-678,682-685,693-696,748-751,751-754,820-826 Detection of Known γ-turns Out of the eleven classic γ-turns identified by Miner-White et. al.[4] Calculate structure only identifies and displays the hbond of thermolysin. Miner-White et. al. indicate that two hbonds are present in the classic γ-turn in thermolysin (They also state that the turn in thermolysin is the only one of the eleven that has two hbonds.), but Calculate structure only displays one and the T segment of the turn is not listed in the summary displayed in the console.
The hbonds were not identified by Calculate structure in the other ten classic turns and the T segments were not included in the summary displayed in the console. Several illustrations are given.
|
References
- ↑ 1.0 1.1 A detailed description is at [1].
- ↑ 2.0 2.1 W. Kabsch & C. Sanders, Biopolymers, 22, 2577-2636, 1983.
- ↑ Characteristics of β-turn classes
- ↑ 4.0 4.1 Miner-White, EJ, et. al. One type of gamma turn, rather than the other, gives rise to chain reversal in proteins. J. Mol. Bio. 204, 1983, pp. 777-782.