[phenixbb] real_space_refine questions
Daniel Larsson
daniel.larsson at icm.uu.se
Mon Jun 12 01:49:10 PDT 2017
I'm have some questions about doing real space refinements using phenix.real_space_refine. I am refining a model for an intermediate-resolution (3.79 Å) cryo-EM map. As the starting point I use a model from a higher-resolution map (3.25 Å) of the same system in another conformation. My idea is to first do rigid body fitting of domains and then gradually relax constrains and restraints.
- Is it possible to restrain the position of certain parts of the model (i.e. "freeze" them so that they move at all) to do focused refinement of subdomains or loops?
- In the instructions it says that you can do both NCS constraints and restraints. How you do the latter? I can only do the former by default or by defining ncs_constraints=True and give a NCS definition file.
- When I do a run with run=rigid_body, some rigid body groups move in to each other and overlap. Is there no collision detection using this algorithm? If I do a run=minimization_global, things work a bit better in this regard, but is not really what I want (at this stage). Am I doing something wrong? Here is the run command and the rigid body definition file:
phenix.real_space_refine \
model.pdb \
map.ccp4 \
rigid.eff \
resolution=8.0 \
ncs_constraints=False \
rotamer_restraints=True \
ramachandran_restraints=True \
c_beta_restraints=True \
run=rigid_body \
macro_cycles=10
# rigid.eff:
refinement {
rigid_body {
group = chain A and resseq 101:266
group = chain A and resseq 273:388
group = chain B and resseq 101:266
group = chain B and resseq 273:388
group = chain C and resseq 104:266
group = chain C and resseq 273:388
}
}
Regards,
Daniel Larsson
More information about the phenixbb
mailing list