Hello all, I had a question about unrestrained refinement. I have several data sets (right near 1 A resolution) that are nearly fully refined, but to finish I would like to run one round of unrestrained refinement. At this point, I want to focus on unrestrained refinement of just the protein, and not any waters. How would the syntax for refining just the chains of the protein (and not water) look in the command line? I also want to include anisotropic refinement of B-factors in this last refinement run. Should it look something like: strategy=individual_sites+occupanices+individual_adp \ adp.individual.anisotropic="not water and not element H" wc=0 wu=0 Then what should I add after this to include just my protein chains and omit waters from the unrestrained refinement? Also, if I wanted to not refine the waters at all in this round of refinement, what should I add to the command line to make that happen? Thank you in advance as always. Regards, Louis
Hi Louis,
I had a question about unrestrained refinement. I have several data sets (right near 1 A resolution) that are nearly fully refined, but to finish I would like to run one round of unrestrained refinement.
at might be a valid thing to do at subatomic resolution.
At this point, I want to focus on unrestrained refinement of just the protein, and not any waters. How would the syntax for refining just the chains of the protein (and not water) look in the command line?
I also want to include anisotropic refinement of B-factors in this last refinement run.
At this resolution you should always refine anisotropic ADP for all (including solvent) but H atoms.
Should it look something like: strategy=individual_sites+occupanices+individual_adp \ adp.individual.anisotropic="not water and not element H" wc=0 wu=0
It depends how you run phenix.refine: using the GUI or typing commands from the command line or using parameter files. From the command line something like phenix.refine model.pdb data.mtz sites.individual="(chain A and resseq 1:10) or (chain C and resname ALA)" will refine coordinates of atoms in residues from 1 to 10 in chain A, and coordinates of atoms in ALA residues of chain C. The coordinates of the other atoms will not be refined at all. You can do the same for ADPs. wc=0 will turn off geometry restraints (globally) wu=0 will turn off ADP restraints (globally). Probably, this is easier to do using a parameter file. For an example, see page 58 here: http://www.phenix-online.org/presentations/latest/pavel_phenix_refine.pdf Please let me know if you have any questions or problems with this. Pavel
participants (2)
-
Louis Lazar
-
Pavel Afonine