Property Selectors: Difference between revisions
Line 160: | Line 160: | ||
==User Notes== | ==User Notes== | ||
===Select by PDB Assigned Atom Number=== | |||
To select atoms by their PDB ATOM number, use the '''id''' selector: | |||
<source lang="python"> | <source lang="python"> | ||
select Nt, id 1-30 | select Nt, id 1-30 | ||
Line 167: | Line 169: | ||
*To select all residues with alternate locations (alt, alt loc), simply do: | ===Alternate Locations=== | ||
====Select All Alternate Locations==== | |||
*To select '''all''' residues with alternate locations (alt, alt loc), simply do: | |||
<source lang="python"> | <source lang="python"> | ||
select aa, not alt "" | select aa, not alt "" | ||
Line 173: | Line 177: | ||
Try it on 1CBN. | Try it on 1CBN. | ||
* To select and remove all atoms with alternate locations (alt loc, | ====Select All Alternate Locations Except the A-Locations==== | ||
* To select and remove all atoms with alternate locations (ALT, alt loc) that aren't the '''A''' record try the following. This removes all alt locations that are not the original '''A''' location. For example given the following pdb record, | |||
<source lang="python" line="1"> | |||
ATOM 949 N ARG A 124 6.039 8.484 9.313 1.00 9.82 N | |||
ATOM 950 CA ARG A 124 5.839 7.071 8.980 1.00 10.05 C | |||
ATOM 951 C ARG A 124 6.009 6.809 7.493 1.00 9.18 C | |||
ATOM 952 O ARG A 124 5.803 7.715 6.648 1.00 9.06 O | |||
ATOM 953 CB ARG A 124 4.437 6.660 9.425 1.00 10.61 C | |||
ATOM 954 CG AARG A 124 4.381 6.539 10.955 0.50 12.54 C | |||
ATOM 955 CG BARG A 124 4.170 6.743 10.930 0.50 14.58 C | |||
ATOM 956 CD AARG A 124 2.996 6.418 11.561 0.50 11.86 C | |||
ATOM 957 CD BARG A 124 2.700 7.091 11.242 0.50 17.41 C | |||
ATOM 958 NE AARG A 124 2.110 7.506 11.144 0.50 13.01 N | |||
ATOM 959 NE BARG A 124 2.229 6.746 12.588 0.50 19.94 N | |||
ATOM 960 CZ AARG A 124 0.912 7.729 11.654 0.50 12.21 C | |||
ATOM 961 CZ BARG A 124 2.123 5.510 13.067 0.50 21.42 C | |||
ATOM 962 NH1AARG A 124 0.433 6.958 12.627 0.50 12.93 N | |||
ATOM 963 NH1BARG A 124 2.470 4.470 12.319 0.50 24.35 N | |||
ATOM 964 NH2AARG A 124 0.185 8.735 11.188 0.50 11.89 N | |||
ATOM 965 NH2BARG A 124 1.658 5.310 14.313 0.50 18.28 N | |||
</source> | |||
the command below will remove the (not A so) B-locations (lines 7, 9, 11, 13, 15, 17). | |||
<source lang="python"> | <source lang="python"> | ||
# select & remove all non A altlocs | |||
remove not (alt ''+A) | remove not (alt ''+A) | ||
# reset the PDB information | |||
alter all, alt='' | alter all, alt='' | ||
</source> | </source> |
Revision as of 09:02, 9 September 2008
Property Selectors
PyMOL reads data files written in PDB, MOL/SDF, Macromodel, ChemPy Model, and Tinker XYZ formats. Some of the data fields in these formats allow PyMOL to assign properties to atoms. You can group and select atoms according to these properties using property selectors and identifiers: the selectors correspond to the fields in the data files, and the identifiers correspond to the target words to match, or the target numbers to compare.
General
The items in a list of identifiers are separated by plus signs (+) only. Do not add spaces within a list of identifiers. The selector resi takes (+)-separated lists of identifiers, as in
PyMOL> select nterm, resi 1+2+3
or, alternatively, it may take a range given with a dash
PyMOL> select nterm, resi 1-3
However, you will get an error message if you try to combine a list and a range in an identifier to a resi as in
PyMOL> select mistake, resi 1-3+6 <--This is WRONG
The identifier for a blank field in an input file is and empty pair of quotes
PyMOL> select unstruct, ss "" # A named selection is created # to contain all atoms that are not assigned # a secondary structure.
Property Selector Table
Most property selectors select matches to their identifiers
Matching Property Selector | Short Form Selector | Identifier and Example |
---|---|---|
symbol | e. |
chemical-symbol-list
PyMOL> select polar, symbol o+n |
name | n. |
atom-name-list
PyMOL> select carbons, name ca+cb+cg+cd |
resn | r. |
residue-name-list
PyMOL> select aas, resn asp+glu+asn+gln or list of up to 2-letter codes for nucleic acids
<B>PyMOL></B> select bases, resn a+g |
resi | i. |
residue-identifier-list
PyMOL> select mults10, resi 1+10+100+1000 residue-identifier-range
<B>PyMOL></B> select nterm, resi 1-10 |
alt | alt |
alternate-conformation-identifier-list
PyMOL> select altconf, alt a+"" |
chain | c. |
chain-identifier-list
PyMOL> select firstch, chain a |
segi | s. |
segment-identifier-list
PyMOL> select ligand, segi lig |
flag | f. |
flag-number
PyMOL> select f1, flag 0 |
numeric_type | nt. |
type-number
PyMOL> select type1, nt. 5 |
text_type | tt. |
type-string
PyMOL> select subset, text_type HA+HC |
id | id |
external-index-number
PyMOL> select idno, id 23 |
index | idx. |
internal-index-number
<PyMOL> select intid, index 11 |
ss | ss |
secondary-structure-type
PyMOL> select allstrs, ss h+s+l+"" |
Numeric Selector Table
Other property selectors select by comparison to numeric identifiers
Numeric Selector | Short Form | Argument and Example |
---|---|---|
b | b |
comparison-operator b-factor-value
PyMOL> select fuzzy, b > 10 |
q | q |
comparison-operator occupancy-value
PyMOL> select lowcharges, q <0.50 |
formal_charge | fc. |
comparison-operator formal charge-value
PyMOL> select doubles, fc. = -1 |
partial_charge | pc. |
comparison-operator partial charge-value
PyMOL> select hicharges, pc. > 1 |
User Notes
Select by PDB Assigned Atom Number
To select atoms by their PDB ATOM number, use the id selector:
select Nt, id 1-30
makes a new selection called Nt and puts in it the first 30 atoms (assuming your PDB starts numbering at 1).
Alternate Locations
Select All Alternate Locations
- To select all residues with alternate locations (alt, alt loc), simply do:
select aa, not alt ""
Try it on 1CBN.
Select All Alternate Locations Except the A-Locations
- To select and remove all atoms with alternate locations (ALT, alt loc) that aren't the A record try the following. This removes all alt locations that are not the original A location. For example given the following pdb record,
ATOM 949 N ARG A 124 6.039 8.484 9.313 1.00 9.82 N
ATOM 950 CA ARG A 124 5.839 7.071 8.980 1.00 10.05 C
ATOM 951 C ARG A 124 6.009 6.809 7.493 1.00 9.18 C
ATOM 952 O ARG A 124 5.803 7.715 6.648 1.00 9.06 O
ATOM 953 CB ARG A 124 4.437 6.660 9.425 1.00 10.61 C
ATOM 954 CG AARG A 124 4.381 6.539 10.955 0.50 12.54 C
ATOM 955 CG BARG A 124 4.170 6.743 10.930 0.50 14.58 C
ATOM 956 CD AARG A 124 2.996 6.418 11.561 0.50 11.86 C
ATOM 957 CD BARG A 124 2.700 7.091 11.242 0.50 17.41 C
ATOM 958 NE AARG A 124 2.110 7.506 11.144 0.50 13.01 N
ATOM 959 NE BARG A 124 2.229 6.746 12.588 0.50 19.94 N
ATOM 960 CZ AARG A 124 0.912 7.729 11.654 0.50 12.21 C
ATOM 961 CZ BARG A 124 2.123 5.510 13.067 0.50 21.42 C
ATOM 962 NH1AARG A 124 0.433 6.958 12.627 0.50 12.93 N
ATOM 963 NH1BARG A 124 2.470 4.470 12.319 0.50 24.35 N
ATOM 964 NH2AARG A 124 0.185 8.735 11.188 0.50 11.89 N
ATOM 965 NH2BARG A 124 1.658 5.310 14.313 0.50 18.28 N
the command below will remove the (not A so) B-locations (lines 7, 9, 11, 13, 15, 17).
# select & remove all non A altlocs
remove not (alt ''+A)
# reset the PDB information
alter all, alt=''