View source for Protein Community

{{Jmol Community Sections}}

This subcommunity is for Jmol users that use Jmol for display and analysis of protein structures. Please add on this page your favorite scripts etc. for displaying proteins in Jmol.
* [[User:AngelHerraez|Angel Herraez]] page
* [[User:EricMartz|Eric Martz]]: [http://bioinformatics.org/firstglance FirstGlance in Jmol] is a simple tool for exploring the structure of any on-line macromolecular model. [http://proteopedia.org Proteopedia.Org] is a Jmol-based wiki with easy Scene-Authoring Tools for creating and sharing customized molecular scenes without learning Jmol language.

== Molecular surfaces ==
[http://biomodel.uah.es/Jmol/surfaces/ documentation and testing page]

== Secondary structure ==
Starting with versions '''12.0.18 and 12.1.15''', Jmol implements the DSSP algorithm for determination of secondary structure in proteins. This is accompanied by some changes in the defaults and some new commands.

=== When reading files ===
# '''If the {{file|pdb}} or {{file|mmcif}} file contains HELIX/SHEET/TURN information''': that information is respected in '''all''' Jmol versions. In '''new''' Jmol versions, {{file|pdb}} and {{file|mmcif}} readers also read helix types alpha, pi, and 3<sub>10</sub>, and color them slightly differently (see colors below).
# '''If the {{file|pdb}} or {{file|mmcif}} file des not contain HELIX/SHEET/TURN information''': 
#* '''Old''' versions of Jmol used Ramachandran-angle based calculation. 
#* '''New''' versions of Jmol use DSSP calculation. 
# In the case of '''alpha-carbon-only chains''', '''all''' versions of Jmol use the method of Levitt and Greer ([http://dx.doi.org/10.1016/0022-2836(77)90207-8 doi]).
# In '''new''' versions of Jmol the calculation based on Ramachandran angles is still available (see below).

=== Forced recalculation of structure ===
==== <code>calculate structure</code>  ====
This command overwrites any secondary structure assignment with a new one:
* In '''old''' Jmol versions, calculated using the Ramachandran method.
* In '''new''' Jmol versions, calculated using the DSSP method.
Note that, intentionally, Jmol turns off all bioshapes when the <code>calculate structure</code> command is given. So, you will usually want to follow it with <code>cartoons on</code> or something of the sort.

==== New commands ====
*<code>calculate structure ramachandran</code> is available for forcing the old method of calculation (it may have some advantages in certain contexts). It can be shortened to <code>calculate structure rama</code>.
* <code>calculate hBonds structure</code> does DSSP determining hydrogen bonds only.
* <code>set defaultStructureDSSP false</code> will change the defaults of file load and <code>calculate structure</code> to use the Ramachandran method; by default, it is ''true'' and uses the DSSP method.
* <code>set dsspCalculateHydrogenAlways false</code> does DSSP enforcing the use of backbone amide H atoms present in the file; by default, it is ''true'' and does standard DSSP, ignoring all backbone amide H atoms present in the file and using rough approximations instead.
* <code>show dssp</code> displays a report of the DSSP calculation.
* <code>set debug</code> displays a verbose DSSP calculation.
* <code>save structure s1</code> (for example) saves in memory the current structure assignment with the internal ID 's1' for later restoring.
* <code>restore structure s1</code> restores the previously saved assignment.

=== Selection and coloring ===
Atom expressions related to secondary structure:
<table>
<tr>
<th>by name</th>
<th></th>
<th>by index</th>
<th></th>
</tr>
<tr>
<td></td>
<td></td>
<td><code>select structure=-1</code></td>
<td>matches non-protein parts (ligands,...)</td>
</tr>
<tr>
<td></td>
<td></td>
<td><code>select structure=0</code></td>
<td>matches unstructured protein (random coil)</td>
</tr>
<tr>
<td><code>select turn</code></td>
<td>or</td>
<td><code>select structure=1</code></td>
<td>matches turns</td>
</tr>
<tr>
<td><code>select sheet</code></td>
<td>or</td>
<td><code>select structure=2</code></td>
<td>matches beta strands</td>
</tr>
<tr>
<td><code>select helix</code></td>
<td>or</td>
<td><code>select structure=3</code></td>
<td>includes all 3 types of helices</td>
</tr>
<tr>
<td><code>select dna</code></td>
<td>or </td>
<td><code>select structure=4</code></td>
<td>matches DNA</td>
</tr>
<tr>
<td><code>select rna</code></td>
<td>or </td>
<td><code>select structure=5</code></td>
<td>matches RNA</td>
</tr>
<tr>
<td><code>select carbohydrate</code></td>
<td>or</td>
<td><code>select structure=6</code></td>
<td>matches carbohydrates</td>
</tr>
<tr>
<td><code>select helix310</code></td>
<td>or</td>
<td><code>select substructure=7</code></td>
<td>matches 3<sub>10</sub> helices (DSSP code "G"). Example: 2JC9.pdb</td>
</tr>
<tr>
<td><code>select helixAlpha</code></td>
<td>or</td>
<td><code>select substructure=8</code></td>
<td>matches alpha helices (DSSP code "H")</td>
</tr>
<tr>
<td><code>select helixPi</code></td>
<td>or</td>
<td><code>select substructure=9</code></td>
<td>matches pi helices (DSSP code "I").<br>Example: 2JC9.pdb after <code>calculate structure</code></td>
</tr>
</table>
''Note'': it is safe to use <code>substructure</code> for all selection commands. For ex., <code>select substructure=2</code> selects beta strands, but <code>select structure=8</code> selects nothing.

* <code>color structure</code> now applies different shades of color to alpha, 3<sub>10</sub> and pi helices, in addition to the classic coloring of helices, beta strands, turns and unstructured portions (see colors below).

* <code>label %[substructure]</code> labels with the texts: helixalpha, helix310, helixpi, sheet, turn, none

* <code>structure helixAlpha</code> ....
* <code>structure helix310</code> ....
* <code>structure helixPi</code> ....

=== Reference section ===
* W. Kabsch, C. Sander (1983) Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features. ''Biopolymers'' '''22''': 2577-2637. doi:10.1002/bip.360221211
* DSSP Home Page: http://swift.cmbi.ru.nl/gv/dssp/

* M. Levitt, J. Greer (1977) Automatic identification of secondary structure in globular proteins. ''J. Mol. Biol.'' '''114''': 181-239. doi:10.1016/0022-2836(77)90207-8

* Codes for structure types:
{| cellpadding="4" cellspacing="0" border="1" style="margin-left:3ex;"
! !! PDB<br>code !! SPSS<br>code !! Jmol codes<br> substructure and name
|-
| right-handed alpha (default) || 1 || H || 8 helixalpha <span style="color:#FF0080;">&#x2588;</span>
|-
| right-handed omega || 2 || ||
|-
| right-handed pi || 3 || G || 7 helixpi <span style="color:#600080;">&#x2588;</span>
|-
| right-handed gamma || 4 || ||
|-
| right-handed 3<sub>10</sub> || 5 || I || 9 helix310 <span style="color:#A00080;">&#x2588;</span>
|-
| left-handed alpha || 6 || ||
|-
| left-handed omega || 7 || ||
|-
| left-handed gamma || 8 || ||
|-
| 2<sub>7</sub> ribbon/helix || 9 || ||
|-
| polyproline || 10 || ||
|-
| (beta) strand  || 0 first<br>1 parallel<br>-1 antiparallel || E || 2  sheet <span style="color:#FFC800;">&#x2588;</span>
|-
| turn || || T || 1 turn <span style="color:#6080FF;">&#x2588;</span>
|-
| isolated beta-bridge residue || || B || 
|-
| bend || || S || 
|-
| none (unstructured protein) || || - || 0 none <span style="color:#FFFFFF;border:1px solid black;">&#x2588;</span>
|-
| DNA || || || 4 dna <span style="color:#AE00FE;">&#x2588;</span>
|-
| RNA || || || 5 rna <span style="color:#FD0162;">&#x2588;</span>
|-
| carbohydrate || || || 6 carbohydrate <span style="color:#A6A6FA;">&#x2588;</span>
|-
| other*  || || || -1 <span style="color:#808080;">&#x2588;</span>
|}
: Helix types in PDB format: columns 39-40 of HELIX record, right-justified.
: Strand sense in PDB format: columns 39-40 of SHEET record, right-justified.
: (*) Other = not protein, not nucleic, not carbohydrate.

== Pretty cartoons ==
Apply this script to get cartoons with a better-than-default look (it will take some more time to rotate them around):
 set hermiteLevel -4; set ribbonAspectRatio 12;
or
 set hermiteLevel 4; set ribbonAspectRatio 12;
which accelerates rotation (the thickness disappears while manipulating with the mouse)