phenix.get_cc_mtz_pdb but allow for atom selections.
It'd be wonderful to pass in an atom selection (say backbone, or ligand or what have you) into phenix.get_cc_mtz_pdb . Is this functionality available in some other phenix util? Thanks! F --------------------------------------------- Francis E. Reyes M.Sc. 215 UCB University of Colorado at Boulder
Hi Francis, phenix.real_space_correlation allows to compute correlations between nearly any user-defined maps, can use x-ray or neutron scattering dictionaries, accepts atoms selections, detail level (cc per residue, atom etc). I don't remember if it is available in the GUI (Nat to confirm), if not then it would be good to have. Pavel On 3/3/12 8:05 AM, Francis E Reyes wrote:
It'd be wonderful to pass in an atom selection (say backbone, or ligand or what have you) into phenix.get_cc_mtz_pdb . Is this functionality available in some other phenix util?
Thanks! F
--------------------------------------------- Francis E. Reyes M.Sc. 215 UCB University of Colorado at Boulder
_______________________________________________ phenixbb mailing list [email protected] http://phenix-online.org/mailman/listinfo/phenixbb
On Sat, Mar 3, 2012 at 8:27 AM, Pavel Afonine
phenix.real_space_correlation
allows to compute correlations between nearly any user-defined maps, can use x-ray or neutron scattering dictionaries, accepts atoms selections, detail level (cc per residue, atom etc).
This doesn't let you define your own level of detail, e.g. CC for a domain, or a loop. It would need to be refactored to pick out map grid points corresponding to a user-input atom selection instead of each residue or each atom. This could theoretically be done for get_cc_mtz_pdb, but I think that uses RESOLVE internally to handle the origin shift, which makes it non-trivial to change. In this case it is probably overkill.
I don't remember if it is available in the GUI (Nat to confirm), if not then it would be good to have.
It's not in the GUI because I found it far too confusing - looking at the code again I am still not sure what some of the parameters mean. I assume this is because it's written to be flexible enough to handle two models or two maps in addition to model-vs-map, but this makes it far too cumbersome for routine use. The validation GUI includes the same residue-or-atom enumeration anyway, which is what 99% of users would use phenix.real_space_correlation for anyway. What Francis is asking for is something like this: phenix.some_awesome_new_command model.pdb data.mtz selection="rename LIG" which phenix.real_space_correlation cannot do, especially since it does not seem to be able to figure out what to do with model and data files automatically. -Nat
Hi Nat, to compute map CC one needs to define two maps. The source of maps can be PDB files, actual maps (for example, in Xplor or CCP4 or other format), Fourier map coefficients (that need to be FTed to get actual maps), or any combination of the above. To compute maps (as opposed to reading them in from input file), the options for input information are: file(s) with map coefficients or model. In latter case you also need to know which scattering dictionary to use: x-ray or neutron or custom user-defined. Also one needs to define which grid points to use: atom level, residue level or any user-defined (the latter is missing in phenix.real_space_correlation and can be trivially added). phenix.real_space_correlation automatically decides whether to use atom or residue level, depending on resolution (obviously, there is no much sense to compute cc per atom at say 4A resolution, while at 1A resolution cc per atom is justified and more informative than cc computed in any larger scale). Also one needs to handle special cases properly. One of many is (of course) hydrogen atoms, which you need to treat differently given the resolution and data type (x-ray or neutron). Since we already have at least three phenix.*** tools to do more or less the same job (plus phenix.model_vs_data can output it too), adding another awesome one wouldn't probably be that awesome unless you somehow make clear what of the above it can and cannot do (or make it *really awesome* so it figures it all out automatically!). Pavel
What does the validation in the GUI use? What is the CC being reported? atom or residue level? F On Mar 3, 2012, at 7:59 PM, Pavel Afonine wrote:
Also one needs to define which grid points to use: atom level, residue level or any user-defined (the latter is missing in phenix.real_space_correlation and can be trivially added). phenix.real_space_correlation automatically decides whether to use atom or residue level, depending on resolution (obviously, there is no much sense to compute cc per atom at say 4A resolution, while at 1A resolution cc per atom is justified and more informative than cc computed in any larger scale).
--------------------------------------------- Francis E. Reyes M.Sc. 215 UCB University of Colorado at Boulder
On Sat, Mar 3, 2012 at 7:18 PM, Francis E Reyes
What does the validation in the GUI use?
Internally, it's mostly the same code as phenix.real_space_correlation doing the work, it just skips the configuration input and uses a simpler routine with the existing PDB and f_model objects. Unless, by the way, you have a twinned structure, in which case this won't work at all due to inconsistencies in the map-related APIs.
What is the CC being reported? atom or residue level?
If you have resolution greater than 2.5 it will show both atom and residue correlations (I think), otherwise just residues. -Nat
Hi Francis, I have added the ability to pass in an atom selection to get_cc_mtz_pdb. You can now (starting with tomorrow's phenix build) say: phenix.get_cc_mtz_pdb model.pdb map.mtz atom_selection="resid 20:30" and this will give you your model-map correlation for residues 20-30 and also for each individual residue in this range. ( It will also give as "overall" the correlation between the entire map and the map calculated from residues 20-30 only; this is not very useful of course.) Let me know of any problems with this. All the best, Tom T ________________________________________ From: [email protected] [[email protected]] on behalf of Francis E Reyes [[email protected]] Sent: Saturday, March 03, 2012 9:05 AM To: PHENIX user mailing list Subject: [phenixbb] phenix.get_cc_mtz_pdb but allow for atom selections. It'd be wonderful to pass in an atom selection (say backbone, or ligand or what have you) into phenix.get_cc_mtz_pdb . Is this functionality available in some other phenix util? Thanks! F --------------------------------------------- Francis E. Reyes M.Sc. 215 UCB University of Colorado at Boulder _______________________________________________ phenixbb mailing list [email protected] http://phenix-online.org/mailman/listinfo/phenixbb
HI all, I want to use phenix.autobuild to calculate SA-OMIT map. BY default, it will be seperated into 16 regions. How could I set a different number of omit regions? Is there any difference of the final omit map by setting different regions? The phenix.autobuild manual says: If you are calculating an omit map, then use nproc=5 * number of omit regions (i.e., up to 100 or more, depending on how many processors you have). Does the * mean to use regions number to multiply 5? If my computer have 16 processors and I want to use 48 regions, so my command should be this: nproc=240. Am I right? Thank you! The information in this e-mail is intended only for the person to whom it is addressed. If you believe this e-mail was sent to you in error and the e-mail contains patient information, please contact the Partners Compliance HelpLine at http://www.partners.org/complianceline . If the e-mail was sent to you in error but does not contain patient information, please contact the sender and properly dispose of the e-mail.
Dear All: I am using phenix-1.7.2 to refine my structure. Now i am facing a problem with .cif file. I could not generate cif files for a modified nucleic acid and MPD at same time with phenix.elbow or phenix.ready_set. If I use phenix.elbow to generate cif file for modefied nucleic acid only, phenix .refine complains about MPD cif not found, same for MPD. If I combine two cif files together, phenix.refine also complains. Is there better to generate cif file containing both modified nucleic acid and MPD? Many thanks for your help, Xue
Maybe I'm missing something here, but don't you just import both cif files
(for nucleic acid, and MPD individually) into phenix.refine? Did you try
that?
Kelly Daughtry
*******************************************************
Kelly Daughtry, Ph.D.
Post-Doctoral Fellow, Raetz Lab
Biochemistry Department
Duke University
Alex H. Sands, Jr. Building
303 Research Drive
RM 250
Durham, NC 27710
P: 919-684-5178
*******************************************************
On Mon, Mar 5, 2012 at 12:29 PM, Xue-Yuan Pei
Dear All: I am using phenix-1.7.2 to refine my structure. Now i am facing a problem with .cif file. I could not generate cif files for a modified nucleic acid and MPD at same time with phenix.elbow or phenix.ready_set. If I use phenix.elbow to generate cif file for modefied nucleic acid only, phenix .refine complains about MPD cif not found, same for MPD. If I combine two cif files together, phenix.refine also complains. Is there better to generate cif file containing both modified nucleic acid and MPD? Many thanks for your help, Xue
______________________________**_________________ phenixbb mailing list [email protected] http://phenix-online.org/**mailman/listinfo/phenixbbhttp://phenix-online.org/mailman/listinfo/phenixbb
Yes, I did try to import both cif files (for nucleic acid and MPD individually) into phenix.refine, but with no luck. any suggestion? Regards, Xue On Mon, 5 Mar 2012, Kelly Daughtry wrote:
Maybe I'm missing something here, but don't you just import both cif files (for nucleic acid, and MPD individually) into phenix.refine? Did you try that?Kelly Daughtry
******************************************************* Kelly Daughtry, Ph.D. Post-Doctoral Fellow, Raetz Lab Biochemistry Department Duke University Alex H. Sands, Jr. Building 303 Research Drive RM 250 Durham, NC 27710 P: 919-684-5178 *******************************************************
On Mon, Mar 5, 2012 at 12:29 PM, Xue-Yuan Pei
wrote: Dear All: I am using phenix-1.7.2 to refine my structure. Now i am facing a problem with .cif file. I could not generate cif files for a modified nucleic acid and MPD at same time with phenix.elbow or phenix.ready_set. If I use phenix.elbow to generate cif file for modefied nucleic acid only, phenix .refine complains about MPD cif not found, same for MPD. If I combine two cif files together, phenix.refine also complains. Is there better to generate cif file containing both modified nucleic acid and MPD? Many thanks for your help, Xue _______________________________________________ phenixbb mailing list [email protected] http://phenix-online.org/mailman/listinfo/phenixbb
On Tue, Mar 6, 2012 at 3:05 AM, Xue-Yuan Pei
Yes, I did try to import both cif files (for nucleic acid and MPD individually) into phenix.refine, but with no luck. any suggestion?
Can you please send us the .eff file and log file for a representative run? You can email them to [email protected] if you don't want the entire list to see. thanks, Nat
Regards, Xue
On Mon, 5 Mar 2012, Kelly Daughtry wrote:
Maybe I'm missing something here, but don't you just import both cif files (for nucleic acid, and MPD individually) into phenix.refine? Did you try that?Kelly Daughtry
******************************************************* Kelly Daughtry, Ph.D. Post-Doctoral Fellow, Raetz Lab Biochemistry Department Duke University Alex H. Sands, Jr. Building 303 Research Drive RM 250 Durham, NC 27710 P: 919-684-5178 *******************************************************
On Mon, Mar 5, 2012 at 12:29 PM, Xue-Yuan Pei
wrote: Dear All: I am using phenix-1.7.2 to refine my structure. Now i am facing a problem with .cif file. I could not generate cif files for a modified nucleic acid and MPD at same time with phenix.elbow or phenix.ready_set. If I use phenix.elbow to generate cif file for modefied nucleic acid only, phenix .refine complains about MPD cif not found, same for MPD. If I combine two cif files together, phenix.refine also complains. Is there better to generate cif file containing both modified nucleic acid and MPD? Many thanks for your help, Xue _______________________________________________ phenixbb mailing list [email protected] http://phenix-online.org/mailman/listinfo/phenixbb
_______________________________________________ phenixbb mailing list [email protected] http://phenix-online.org/mailman/listinfo/phenixbb
Hi Chen, You can set the target number of omit regions with the keyword n_box_target: n_box_target= None You can tell the Wizard how many omit boxes to try and set up (but it will not necessarily choose your number because it has to be nicely divisible into boxes that fit your asymmetric unit). A suitable number is 24. The larger the number of boxes, the better the map will be, but the longer it will take to calculate the map. The number of processors is up to you. The suggestion of 5 times the number of boxes is the maximum number of processors that would be useful. All the best, Tom T ________________________________________ From: [email protected] [[email protected]] on behalf of Qiang Chen [[email protected]] Sent: Monday, March 05, 2012 10:03 AM To: PHENIX user mailing list Subject: [phenixbb] SA-OMIT map/omit regions HI all, I want to use phenix.autobuild to calculate SA-OMIT map. BY default, it will be seperated into 16 regions. How could I set a different number of omit regions? Is there any difference of the final omit map by setting different regions? The phenix.autobuild manual says: If you are calculating an omit map, then use nproc=5 * number of omit regions (i.e., up to 100 or more, depending on how many processors you have). Does the * mean to use regions number to multiply 5? If my computer have 16 processors and I want to use 48 regions, so my command should be this: nproc=240. Am I right? Thank you! The information in this e-mail is intended only for the person to whom it is addressed. If you believe this e-mail was sent to you in error and the e-mail contains patient information, please contact the Partners Compliance HelpLine at http://www.partners.org/complianceline . If the e-mail was sent to you in error but does not contain patient information, please contact the sender and properly dispose of the e-mail. _______________________________________________ phenixbb mailing list [email protected] http://phenix-online.org/mailman/listinfo/phenixbb
participants (7)
-
Francis E Reyes
-
Kelly Daughtry
-
Nathaniel Echols
-
Pavel Afonine
-
Qiang Chen
-
Terwilliger, Thomas C
-
Xue-Yuan Pei