| | 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 | |