Hi, I am using NCS restraints for 4 copies of a molecule at 3.05 Å. We are using the standard weights: coordinate_sigma = 0.05 b_factor_weight = 10 and the 4 molecules look identical. What are suggested values for looser restraints? One of my colleagues said that varying them did not make a difference- ncs was still very strict. Is there a way to make them less strict? Thanks, Timothy A. Springer, Ph.D.
Hi Tim, increasing the coordinate_sigma to 0.5 will make the restraint considerably looser. However, at 3.0A resolution I'd want to be convinced that this change made things better by looking at the R-free (i.e. I'd expect that tighter restraints were more appropriate at this resolution). If there are areas where there are detectable differences between the NCS copies you might want to exclude them from the restraint selection and keep the rest tight. Cheers, Paul On Dec 7, 2010, at 7:31 PM, Timothy Springer wrote:
Hi, I am using NCS restraints for 4 copies of a molecule at 3.05 Å. We are using the standard weights: coordinate_sigma = 0.05 b_factor_weight = 10 and the 4 molecules look identical. What are suggested values for looser restraints? One of my colleagues said that varying them did not make a difference- ncs was still very strict. Is there a way to make them less strict? Thanks, Timothy A. Springer, Ph.D.
_______________________________________________ phenixbb mailing list [email protected] http://phenix-online.org/mailman/listinfo/phenixbb
-- Paul Adams Deputy Division Director, Physical Biosciences Division, Lawrence Berkeley Lab Adjunct Professor, Department of Bioengineering, U.C. Berkeley Vice President for Technology, the Joint BioEnergy Institute Head, Berkeley Center for Structural Biology Building 64, Room 248 Tel: 1-510-486-4225, Fax: 1-510-486-5909 http://cci.lbl.gov/paul Lawrence Berkeley Laboratory 1 Cyclotron Road BLDG 64R0121 Berkeley, CA 94720, USA. Executive Assistant: Kathleen Moody [ [email protected] ][ 1-510-495-2506 ] --
Hi Timothy, the default values (coordinate_sigma = 0.05 , b_factor_weight = 10) should be good in most cases. It depends what you call "looser". Yes, you can make them looser by increasing coordinate_sigma (for refinement of individual coordinates), and decreasing b_factor_weight (for individual isotropic B-factor refinement). You can do it specifically for selected atoms, since you can define as many NCS related groups of atoms as you wish and each defined with their individual weights: http://phenix-online.org/documentation/refinement.htm#anch24 refinement.ncs.restraint_group { reference = chain A resid 1:4 selection = chain B and resid 1:3 selection = chain C coordinate_sigma = 0.05 b_factor_weight = 10 } refinement.ncs.restraint_group { reference = chain E selection = chain F coordinate_sigma = 0.025 b_factor_weight = 3 } refinement.ncs.restraint_group { reference = chain Z and backbone selection = chain F and backbone } refinement.ncs.restraint_group { reference = chain Z and sidechain selection = chain F and sidechain selection = (chain I or chain K) and sidechain } Let me know if you have any questions. Pavel. On 12/7/10 7:31 PM, Timothy Springer wrote:
Hi, I am using NCS restraints for 4 copies of a molecule at 3.05 Å. We are using the standard weights: coordinate_sigma = 0.05 b_factor_weight = 10 and the 4 molecules look identical. What are suggested values for looser restraints? One of my colleagues said that varying them did not make a difference- ncs was still very strict. Is there a way to make them less strict? Thanks, Timothy A. Springer, Ph.D.
Hi Tim, If you start from nearly identical copies (due to previous refinement), in the absence of a strong signal from the data, it is expected that the copies remain fairly similar even if you increase the sigma / lower the weight. We just added a new feature to phenix.refine that could be useful to introduce differences between the NCS copies, guided by the model density through a new real-space refinement protocol: strategy=individual_sites_real_space plus real_space_refinement.mode=lockit. At your resolution it is turned on by default. I.e. to try it out all you need to do is get dev-602 or higher and run phenix.refine again. Only trouble, those installers aren't generally visible since there were some minor problems, but watch out for a new one soon. Ralf
From: Timothy Springer
To: [email protected] Sent: Tue, December 7, 2010 7:31:44 PM Subject: [phenixbb] ncs restraints Hi, I am using NCS restraints for 4 copies of a molecule at 3.05 Å. We are using the standard weights: coordinate_sigma = 0.05 b_factor_weight = 10 and the 4 molecules look identical. What are suggested values for looser restraints? One of my colleagues said that varying them did not make a difference- ncs was still very strict. Is there a way to make them less strict? Thanks, Timothy A. Springer, Ph.D.
Thanks, I believe this is an important tip. Tim On Dec 8, 2010, at 12:25 AM, Ralf W. Grosse-Kunstleve wrote:
Hi Tim,
If you start from nearly identical copies (due to previous refinement), in the absence of a strong signal from the data, it is expected that the copies remain fairly similar even if you increase the sigma / lower the weight. We just added a new feature to phenix.refine that could be useful to introduce differences between the NCS copies, guided by the model density through a new real-space refinement protocol: strategy=individual_sites_real_space plus real_space_refinement.mode=lockit. At your resolution it is turned on by default. I.e. to try it out all you need to do is get dev-602 or higher and run phenix.refine again. Only trouble, those installers aren't generally visible since there were some minor problems, but watch out for a new one soon.
Ralf
From: Timothy Springer
To: [email protected] Sent: Tue, December 7, 2010 7:31:44 PM Subject: [phenixbb] ncs restraints Hi, I am using NCS restraints for 4 copies of a molecule at 3.05 Å. We are using the standard weights: coordinate_sigma = 0.05 b_factor_weight = 10 and the 4 molecules look identical. What are suggested values for looser restraints? One of my colleagues said that varying them did not make a difference- ncs was still very strict. Is there a way to make them less strict? Thanks, Timothy A. Springer, Ph.D.
_______________________________________________ phenixbb mailing list [email protected] http://phenix-online.org/mailman/listinfo/phenixbb
participants (4)
-
Paul Adams
-
Pavel Afonine
-
Ralf W. Grosse-Kunstleve
-
Timothy Springer