[phenixbb] Phenix.refine bonds/angles rmsd

[guanrj at foxmail.com]

Hi Pavel

Thanks for your suggestion. I have another question:

I am refining another 3.3 A resolution structure. In coot I manually adjusted the model
to have >95% residues in preferred region in ramanchandran plot, but after refinement  
in Phenix only <80% residues in the preferred region and a lot more in the outlier region.
This was wen I used the default refinement parameters "phenix.refine mtz_file pdb_file".
R and R_free are 0.22/0.27 in this case. 

Then when I used "phenix.refine mtz_file pdb_file wxc_scale=0.1 wxu_scale=0.3", I can 
get better results in terms of ramanchandran plot (say, ~92% in preferred region, and a
little more in the outlier). R and R_free factors are 0.25/0.30.

I guess at low resolution (like in my case, 3.3 A) I have to give more weights to geometry 
restrains. Did I do it in the right way, or did I missing anything?

Thanks

Rongjin Guan 




= = = = = = = = On 2010-04-13 16:42:21 You wrote = = = = = = = = 
Hi Rongjin,
  
when I use phenix.refine to refine my structure, the rmsd's of bonds and angles are 
pretty large. My data is 2.3 A resolution.  

REMARK
REMARK ******************** REFINEMENT SUMMARY: QUICK FACTS *******************
REMARK Start: r_work = 0.2838 r_free = 0.2789 bonds = 0.029 angles = 2.002
REMARK Final: r_work = 0.2407 r_free = 0.2466 bonds = 0.029 angles = 2.002
REMARK ************************************************************************
REMARK


really (just wondering)? 
What I would be really worry is the difference between Rfree and Rwork, which apparently is terribly small, indicating that there is something wrong with either free-R flags or refinement strategy.


I tried many options but none worked so far. Here are the options that I tried:

phenix.refine mtz_file pdb_file strategy=tls tls_group_selections.params wxc_scale=0.3 wxu_scale=0.8

phenix.refine mtz_file pdb_file strategy=tls tls_group_selections.params fix_wxc=0.3 fix_wxu=0.8 

phenix.refine mtz_file pdb_file strategy=tls tls_group_selections.params fix_wxc=0.3 fix_wxu=0.8 optimize_wxc=false optimize_wxu=false

These gave the same refinement results. Any suggestions?

It's not surprising, because in all refinements above you did not refine the coordinates at all. The keyword "strategy=tls" means you refined B-factors only (more precisely, you did constrained anisotropic refinement, called TLS refinement). And therefore using wxc_scale or fix_wxc or  optimize_wxc would not make any difference. Also, optimize_wxc is set to False by default, so you don't need to explicitly say it.

I would run it like this:

phenix.refine mtz_file pdb_file strategy=tls+individual_adp+individual_sites tls_group_selections.params optimize_wxc=true

and hope that this fixes the "problem".

Let me know if you have any questions!
Pavel.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://phenix-online.org/pipermail/phenixbb/attachments/20100414/e9e18b60/attachment-0003.htm>