Hi Maia,

you can try to build a new test set of reflections using "thin shells" algorithm. PHENIX can't create such test set, so you need to use an outside program that can generate it. Make sure you generate the same number of test reflections.

However, and this is a *big* however, you will have to remove the "memory" from newly generated test set. How? Remove all waters, reset all B-factors to average one, shake coordinates to achieve rmsd~1A and then run 10-20 macro-cycles of refinement in phenix.refine (or until the R-factors stop changing), that will include:
- simulated-annealing;
- automatic water picking;
- refinement of individual coordinates and B-factors;
- after you have done this run, run the second run of simple xyz and b-factors refinement using "optimize_wxc=true optimize_wxu=true".

Repeat the above steps with and w/o NCS. What do you get?

I would also do the above extensive model shaking/re-refining using the original test set, just to make sure that this will not be enough to make the Rfree-Rwork, Rwork and Rfree change significantly without resorting to a different test set.

Pavel.

PS> I also use 10% of test reflections...

On 5/5/09 2:14 PM, Maia Cherney wrote:

Thank you Nat, Pavel,

Actually, the NCS approaches a 3-fold crystallographic symmetry axis 
with 3 types of protein molecules plus 2 restraint groups for DNA molecules.

The result is the same with "optimize_wxc=true optimize_wxu=true" 
options, Rfree is going up without NCS by ~1.5%. Rwork is practically 
the same.

Should I increase the test set (now it's 5%)? or should I leave NCS=true 
for the final pdb?

Maia


Nathaniel Echols wrote:

On May 4, 2009, at 6:48 PM, Maia Cherney wrote:

The resolution is 2.15 A. The NCS was always on during the refinement
until we got low R factors (19.2% and 21.2%). Then the NCS was turned
off for the final refinement and the R factors increased, which is
strange as they should be going down when you apply less restraints.

This isn't necessarily true - regardless of which program you use to  
refine, R-work isn't actually the refinement target, and minimization  
algorithms aren't foolproof either.  When you remove the restraints  
you are effectively decreasing the observation:parameter ratio, which  
increases the risk of overfitting.  If nothing else I would expect the  
gap between Rwork and Rfree to increase, but I've also seen both  
increase when NCS restraints were removed.

It may be possible that the NCS restraints are biasing your test set,  
leading to the suspiciously low Rfree-Rwork gap, but I thought this  
only applied when the NCS approached crystallographic symmetry, and  
there's no such thing as fivefold crystallographic symmetry.   
Hopefully someone else on the list can clarify this.

-Nat

-------------------
Nathaniel Echols
Lawrence Berkeley Lab
510-486-5136
[email protected]



_______________________________________________
phenixbb mailing list
[email protected]
http://www.phenix-online.org/mailman/listinfo/phenixbb


_______________________________________________
phenixbb mailing list
[email protected]
http://www.phenix-online.org/mailman/listinfo/phenixbb