| | 1 | == Surface Exposed Residues == |
| | 2 | |
| | 3 | Here is how to select surface exposed residues that are hydrophobic. |
| | 4 | |
| | 5 | Use [http://rbvi.ucsf.edu/chimerax/docs/user/commands/measure.html#sasa measure sasa] to define an "area" attribute of residues, then assign residue hydrophobicities from a file with the [http://rbvi.ucsf.edu/chimerax/docs/user/commands/defattr.html defattr] command, and then use a single [http://rbvi.ucsf.edu/chimerax/docs/user/commands/atomspec.html#attributes select] command to intersect these two attributes. |
| | 6 | |
| | 7 | In the [http://rbvi.ucsf.edu/chimerax/docs/user/commands/defattr.html#examples defattr example files] I already happened to include a file kdHydrophobicity.txt to assign Kyte-Doolittle hydrophobity, as well as for several other hydrophobicity scales. Just download whichever one(s) you want. |
| | 8 | |
| | 9 | The commands could be something like: |
| | 10 | |
| | 11 | {{{ |
| | 12 | open 6g2j |
| | 13 | measure sasa protein |
| | 14 | defattr ~/Desktop/kdHydrophobicity.txt |
| | 15 | select ::area>25 & ::kdHydrophobicity>0.0 |
| | 16 | surface protein enclose protein |
| | 17 | color sel magenta |
| | 18 | }}} |
| | 19 | |
| | 20 | ... which looks like this ... |
| | 21 | |
| | 22 | [[Image(6gj2_kd.png)]] |
| | 23 | |
| | 24 | ... in pretty good agreement with MLP coloring, e.g. |
| | 25 | |
| | 26 | |
| | 27 | {{{ |
| | 28 | mlp surfaces #1.3 |
| | 29 | }}} |
| | 30 | |
| | 31 | [[Image(6gj2_mlp.png)]] |
| | 32 | |
| | 33 | In the select command you may want to experiment with the surface-area and hydrophobicity cutoffs, as well as which hydrophobicity scale you use. I didn't try a lot of values for surface-area cutoff, but here is with all residues with values >25 colored red, reasonable at least as a first approximation: |
| | 34 | |
| | 35 | [[Image(6gj2_exposed.png)]] |
| | 36 | |
| | 37 | Elaine Meng |