| 3 | == Examples for proposed surface command == |
| 4 | |
| 5 | Here is a suggestion for how the Chimera 2 surface command could work (TomG). |
| 6 | The basic idea is to surface by chain (no ligands, solvent, ions), but allow surfacing whole proteins as another option, and allow surfacing any atoms as a rare case. The rule to ignore ligands, solvent and ions unless only those are specified is tricky but I think necessary to handle the common case in a simple way. |
| 7 | |
| 8 | To surface each chain of #1 excluding ligands, solvent, ions: |
| 9 | |
| 10 | surface #1 |
| 11 | |
| 12 | To surface all atoms in one surface use the "enclose" option which automatically excludes ligands, solvent and ions: |
| 13 | |
| 14 | surface enclose #1 |
| 15 | |
| 16 | To surface just chain A of #1 which has chains A,B,C: |
| 17 | |
| 18 | surface #1/A |
| 19 | |
| 20 | To surface chain A but only show residues 10-20: |
| 21 | |
| 22 | surface #1/A:10-20 |
| 23 | |
| 24 | To surface the whole molecule excluding ligands, solvent, ions, but only show chain A residues 10-20 |
| 25 | |
| 26 | surface #1/A:10-20 enclose #1 |
| 27 | |
| 28 | To surface just the ligands. Above I said "surface #1" excludes ligands but "surface ligand" includes them -- if the specifier is empty when ligands, ions and solvent are thrown out but not empty when they are kept, then keep them. |
| 29 | |
| 30 | surface ligand |
| 31 | |
| 32 | To surface ligands only of molecule #1 when other molecules are open: |
| 33 | |
| 34 | surface #1 & ligand |
| 35 | |
| 36 | To surface #1 including main and ligands per chain use "ligands" keyword option: |
| 37 | |
| 38 | surface #1 ligands true |
| 39 | |
| 40 | To change which patch is being displayed of an already created surface. In this example the first command creates the chain A surface showing residue 10-20. The second command modifies the existing surface because it sees that the chain A surface already exists (exact same set of atoms surfaced): |
| 41 | |
| 42 | surface #1/a:10-20 |
| 43 | |
| 44 | surface #1/a:30-40 |
| 45 | |
| 46 | If you want two copies of a surface enclosing the exact same set of atoms |
| 47 | |
| 48 | surface #1 copy true |
| 49 | |
| 50 | Now how do you change the two surfaces to show different patches? |
| 51 | |
| 52 | surface #1/a:10-20 using #1.1 |
| 53 | |
| 54 | surface #1/a:30-40 using #1.2 |
| 55 | |
| 56 | What will the model numbers be for surfaces? |
| 57 | |
| 58 | surface #1 |
| 59 | |
| 60 | produces (the model name is "surfaces") |
| 61 | |
| 62 | #1.1 surfaces |
| 63 | |
| 64 | and if #1 has three chains these will be |
| 65 | |
| 66 | #1.1.1 A |
| 67 | |
| 68 | #1.1.2 B |
| 69 | |
| 70 | #1.1.3 C |
| 71 | |
| 72 | Probably model panel will not show that third level by default, you would need to expand the line "#1.1 surfaces" to see them. If you then make another surface |
| 73 | |
| 74 | surface enclose #1 |
| 75 | |
| 76 | you get another model (with no children since the enclose option produces a single surface) |
| 77 | |
| 78 | #1.2 surface |
| 79 | |
| 80 | and the #1.1 surfaces still exist and could be closed if desired with |
| 81 | |
| 82 | close #1.1 |
| 83 | |
| 84 | If you surface two or more molecules with one surface (rare situation) |
| 85 | |
| 86 | surface enclose #1 | #2 |
| 87 | |
| 88 | then the resulting surface model will be #3 (or the next available model number) instead of being a child of either model. |
| 89 | |
| 90 | |