- phenixbb - phenix-online.org

refinement / alpha_beta
by Kay Diederichs 05 Oct '10

05 Oct '10

Dear Phenix developers, I'd like to refine (using phenix.refine) a _very_ incomplete model (obtained by molecular replacement with one subunit of the complete structure) against the native amplitudes, with the goal to obtain meaningful phases and figures_of_merit which should serve to find heavy atom sites in derivative datasets. If I "just refine" in phenix.refine, then the R-factors are in the 50s and the program estimates the coordinate error to be very large. Thus I would like to tell phenix.refine that the coordinates are actually quite accurate, but that a lot of the model is simply missing. I found alpha_beta free_reflections_per_bin= 140 number_of_macromolecule_atoms_absent= 225 n_atoms_included= 0 bf_atoms_absent= 15.0 final_error= 0.0 absent_atom_type= "O" method= *est calc estimation_algorithm= *analytical iterative verbose= -1 interpolation= True fix_scale_for_calc_option= None number_of_waters_absent= 613 sigmaa_estimator kernel_width_free_reflections= 100 kernel_on_chebyshev_nodes= True number_of_sampling_points= 20 number_of_chebyshev_terms= 10 use_sampling_sum_weights= True in the documentation but I have trouble understanding its meanings and defaults - but I do think this should be used. Maybe I should switch to method=calc and specify number_of_macromolecule_atoms_absent= as well as bf_atoms_absent= ? Can anyone advise me, please? thank you, Kay

3 2

best strategy using two datasets?
by Manoj saxena 05 Oct '10

05 Oct '10

Hi , I have data 2.4 A for a protein ( Se-Met derivative) unit cell dimensions. 53.098,60.207,142.382. I have used this data to solve the partial structure with phenix.autosol and phenix.autobuild. The model from this data was then used as MR model in phaser with another data set from a crystal at 2.2 A (group P212121 ) .unit cell dimensions are 53.86,61.47,142.49. Currently my R free is 23%. My question is 1) Is there a a better way to use the two datasets. can I merge the reflection files from two datasets and does using a merged reflection file have an advantage over the strategy I used. 2) Although not directly related to phenix I was wondering can we use the merged reflection file for PDB depositions? Thanks . I will greatly appreciate any suggestions and guidance in this regard Manoj

2 1

Re: [phenixbb] Help on multi-start-simulated-annealing sigmaA-weighted 2fofc map - WORKING EXAMPLE
by Pavel Afonine 03 Oct '10

03 Oct '10

Hi Fengyun, no problems, here is the updated script and inputs files to illustrate how it works: http://cci.lbl.gov/~afonine/mssa2/ All it needs is a bunch of PDB files and corresponding map coefficients. Let me know if you have any questions. Pavel. On 9/9/10 7:48 PM, fn1(a)rice.edu wrote: > Hi Pavel, > > The script works fine. I hope you could slightly revise the script for > me. > > Since right now, I could run the multi-sa on different nodes at the > same time. So I could get the refined models in the name of sa-#.pdb > within a short time. > > So I could begin with the different models and combine them to > "all.pdb" and calculate the average map after that. But I don't know > how to revise the run.py to do that. > > Wait for your help. > Thanks! > Fengyun > > 引用 Pavel Afonine <pafonine(a)lbl.gov>： > >> Hi Fengyun, >> >> my previous email actually doens't answer your question - I realized >> this after I pushed Sent button -:) >> >> To really answer your question I just wrote a Python script that does >> the following: >> - run multi-start Simulated Annealing, >> - combine all refined models into one multi-model PDB file (models >> separated by MODEL-ENDMDL records), and >> - compute averaged 2mFo-DFc map. >> >> The complete working example is here: >> >> http://cci.lbl.gov/~afonine/mssa/ >> >> You have two options at this point: >> >> - take this example, slightly change inputs by editing the run.py >> file (change input data and model file names, number of SA runs, etc) >> and run it as following: >> >> phenix.python run.py >> >> (since I spent 15 minutes on writing this script it's obviously not >> thoroughly tested or parameter-optimized, although I believe it >> should do the job right) >> >> or >> >> - wait a few days for one of the next PHENIX nightly builds where the >> above script will be available as a user-friendly either >> phenix.multi_start_sa command or an option of phenix.maps (I haven't >> decided yet what is better). >> >> Let me know if you have any questions. >> >> Pavel. >> >>> I read on the paper that they use CNS (model_map.inp) to calculate >>> the avarage sigmaA-weighted 2fofc map. Right now I have set up the >>> simulated annealing with different random seeds. I wonder whether >>> phenix could calculate the average map with the multiple coordinates >>> from different starts? >> >> >> > > >

2 3

Job Opening Announcement
by Lance Westerhoff 01 Oct '10

01 Oct '10

Hello All- We have a number of computational chemist and computational crystallographer post doc positions available. One has a significant focus on x-ray refinement method development, and the other one will explore computational chemistry method research and development. We are excited to fill these openings with talented individuals, so if you think you would be a good fit, we invite you to take a look at our careers website below and submit your CV (along with cover letter) to careers(a)quantumbioinc.com http://www.quantumbioinc.com/careers/ Thank you for your time! -Lance ____________________ Lance M. Westerhoff, Ph.D. General Manager QuantumBio Inc. WWW: http://www.quantumbioinc.com Email: lance(a)quantumbioinc.com Phone: 814-235-6908 Fax: 814-235-6909

1 0

using custom_param file from the gui
by Sebastiano Pasqualato 01 Oct '10

01 Oct '10

Hi all, I have a custom parameter file that I used to use running phenix.refine from the command line. Is there a way to use it from the gui? ciao, s -- Sebastiano Pasqualato, PhD IFOM-IEO Campus Dipartimento di Oncologia Sperimentale Istituto Europeo di Oncologia via Adamello, 16 20139 - Milano Italy tel +39 02 9437 5094 fax +39 02 9437 5990

1 1

question on the weight for adp
by fn1＠rice.edu 30 Sep '10

30 Sep '10

Hi everyone, I have some questions the weight for adp. My case is that some part of the model will be refined to large B values (about 100) if the default value of wxu_scale (1.0) is used. The other part of the model will have average B value of about 40, which the wilson B is about 35. The corresponding 2FOFC map has very poor density for the region with large B values. If the value of wxu_scale set to be 0.1, that part with large B will be refined to about 80 for B. The corresponding 2FOFC map has some better (but not good enough) density for this region. My question is that can I keep the wxu_scale at the small value for the refinement? Or finally do I need to loose this value? Thanks! Fengyun

4 5

Seeking Postdoctoral Researcher
by Francois Berenger 29 Sep '10

29 Sep '10

Laboratory: Zhang Initiative Research Unit, Advanced Science Institute, RIKEN, Japan. (Unit Leader: Dr. Kam Zhang) Job title and Job description: Postdoctoral Researchers Job description: We are seeking a highly motivated and experienced computational structural biologist to join the Zhang Initiative Research Unit at RIKEN. The successful candidate will use computational tools to tackle problems in the area of protein folding, structure prediction, and structure-based drug design. Our ability to predict the structure of proteins from their primary sequences will greatly facilitate our understanding of the important biological functions that proteins play in their host systems. The 3D structure information will also facilitate the rationale design of drugs for the treatment of various diseases with unmet medical needs. Qualifications: This position requires a PhD in areas related to biophysics, biochemistry, computational chemistry, bioinformatics, or structural biology with experience in employing computational tools to solve biological or chemical problems. Experiences in protein folding, protein structure prediction, protein structure analysis, protein crystallography, or structure-based drugs are highly desired. Proficiency in object-oriented programming such as C++, Java, scripting languages such as Python, Perl are also desired. Experiences in high performance computing, network computing, or grid computing would be a plus. The ideal candidate should possess a track record of accomplishments demonstrating technical proficiency, independent thinking, and scientific creativity. The research environment will be international and hence communication in English is encouraged. The complete job ad, as well as how to apply can be viewed here: http://www.riken.jp/engn/r-world/info/recruit/k100817_s_asi.html Best regards, Francois.

1 0

Ligand restraints editing with REEL -summary
by Joseph Brock 28 Sep '10

28 Sep '10

Hi all, I just wanted to post a quick summary of how my problem was resolved. I later realised that the 'native glutahione' molecule I was transferring ligand restraints from was actually produced by eLBOW released under the Phenix package v1.4, and was therefore defined differently to those restraints created by v1.6.4-486. I was advised by Nigel (author of eLBOW and REEL) to download the latest nightly build (dev-541) and create a clean subdirectory. Copying my model.pdb and data.mtz file into this, running phenix.ready_set on these files and running phenix.refine with the resulting .pdb and .cif files, resulted in a ligand that was correctly restrained to refine into the density extremely well. Hope this is of help to anyone with a similar problem. Cheers, Joe. > From: phenixbb-request(a)phenix-online.org > Subject: phenixbb Digest, Vol 58, Issue 18 > To: phenixbb(a)phenix-online.org > Date: Sat, 18 Sep 2010 12:00:01 -0700 > > Send phenixbb mailing list submissions to > phenixbb(a)phenix-online.org > > To subscribe or unsubscribe via the World Wide Web, visit > http://phenix-online.org/mailman/listinfo/phenixbb > or, via email, send a message with subject or body 'help' to > phenixbb-request(a)phenix-online.org > > You can reach the person managing the list at > phenixbb-owner(a)phenix-online.org > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of phenixbb digest..." > > > Today's Topics: > > 1. Ligand restraints editing with REEL (Joseph Brock) > 2. Re: Ligand restraints editing with REEL (Ralf W. Grosse-Kunstleve) > 3. Resolve NCS masks (Florian Schmitzberger) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sat, 18 Sep 2010 08:09:23 +0000 > From: Joseph Brock <joeylives2ride(a)hotmail.com> > To: <phenixbb(a)phenix-online.org> > Subject: [phenixbb] Ligand restraints editing with REEL > Message-ID: <SNT120-W580E16EA4DC524ED59A3209B7C0(a)phx.gbl> > Content-Type: text/plain; charset="iso-8859-1" > > > Hi all, > > I'm working on a 1.5A structure that contains a novel ligand within the active site, which is an Arsenic derivative of Glutathione, using the Phenix GUI, version 1.6.4-486. > > I have been using the 'Molinspiration' webpage (http://www.molinspiration.com/docu/webme/) to generate a smiles string, which I then feed into ELBOW. Although the structure/chemical restraints of Glutathione should be unchanged by the covalent linkage to arsenic through the cysteinyl sulfur, the resulting restraints file consistently restricts one of the dihedrals to a value that is incompatible with the (quite clear) electron density. This occurs regardless of wether i use AM1 optimization, and is unresponsive to using native Glutathione (which refines into the corresponding density extremely well) as an input geometry file with any of the available 'Use geometry for' options (final geometry, initial geometry, etc). > > I then tried opening the output file in REEL and manually changing the corresponding dihedrals to those of glutathione, however after running simple optimization, the resulting restraints file has the same problem, and upon choosing the 'transer geometry and restraints' option (or any option) midway through running the 'eLBOW optimization' or 'AM1 optimization' options, Phenix consistently crashes. > > Any advice on how to transfer the origional geometry of native Glutathione to my new ligand effectively would be greatly appreciated. I'm a little confused as to the abilities of REEL, and how to affect them. > > Thanks so much in advance. > > Cheers, > > -Joe. > > >

1 0

mlhl refinement
by Florian Schmitzberger 28 Sep '10

28 Sep '10

Dear All, I am running Phenix version 1.6.4.-486 on imac. I just wanted to refine (phenix.refine) a model with MLHL target and experimental phase information. Strange thing is when I load the mtz file, as written out from Phenix Autosol (overall_best_denmod_map_coeffs.mtz) into the GUI, the Phenix refine GUI does not seem to any phase information. Only picks up FP and SIGFP labels; but not the Hendrickson-Lattman coefficients, which are present. The following message comes up when I try to chose mlhl refinement: " You have selected 'mlhl' as the refinement target, but no experimental phases were supplied. If you wish to use phase restraints in refinement, please add or select experimental phases, or choose another target function". I guess it is not such a major problem as I can save the .eff file and manually edit it to include the HL coefficients and then run Phenix (as long as Phenix then recognizes them from the mtz file). I would also like to refine f' and f'', as the data are from a selenomethionine derivatized protein. From what I can see the default in Phenix.Autosol seems to be not to keep the F+/- and Dano labels for the output mtz file? So I suppose I do need to provide the initial mtz file or is there a way to change the output options for the coefficients in Phenix Autosol? Cheers, Florian ----------------------------------------------------------- Florian Schmitzberger Biological Chemistry and Molecular Pharmacology Harvard Medical School 250 Longwood Avenue, SGM 130 Boston, MA 02115, US Tel: 001 617 432 5602

4 4

Does Phenix have the library or program to maintain bond angle and bond length to create an extended chain?
by Rojan Shrestha 27 Sep '10

27 Sep '10

Hello: I am wondering whether Phenix has the library or program to maintain the ideal bond length and bond angle. The torsion angles of residue and xyz of each atoms are already known. I want to juxtapose residue or residues to create an extended chain. Thanks Regards, Rojan

4 3

CCP4 Study Weekend 2011: "Model Building and Refinement and Validation"
by ronan.keegan＠stfc.ac.uk 24 Sep '10

24 Sep '10

Dear All, We're pleased to announce that registration for this coming January's CCP4 Study Weekend meeting is now open. The meeting will take place from January 5th to 7th, 2011 at University of Warwick in the UK. The topic of the meeting will be "Model building and Refinement and Validation". This year's organisers, Roberto Steiner and Bernhard Rupp, have arranged the usual mix of interesting and stimulating lectures from some of the foremost speakers in the field along with the opportunity to see demonstrations of the latest software developments from CCP4 during the lunch-time bytes sessions. For more information on the program and registration details please go to the Study Weekend website: http://www.cse.scitech.ac.uk/events/CCP4_2011/ We look forward to seeing you in Warwick. Best wishes, The CCP4 team -- Scanned by iCritical.

1 0

separate amino acid chains
by Florian Schmitzberger 23 Sep '10

23 Sep '10

Dear All, Quick question regarding the handling of different amino acid chains in the Phenix Autobuild GUI. From the manual I understand that different amino acid chains are separated by > ? I pasted the my amino acid chains in the GUI separated with > ; I don't think in my case Phenix recognizes the different chains as such, but as one contiguous protein. Am I missing anything obvious? Regards, Florian ----------------------------------------------------------- Florian Schmitzberger Biological Chemistry and Molecular Pharmacology Harvard Medical School 250 Longwood Avenue, SGM 130 Boston, MA 02115, US Tel: 001 617 432 5602

2 1

autobuild : use of ncs
by Lepore, Bryan 22 Sep '10

22 Sep '10

i read and see that autobuild finds ncs automatically and uses it. i have one .pdb with one copy of the model for segid's A and B, and chainid's A and B. however, i notice sometimes in density modification that ncs has "NO" flags as to its use, and the output .pdb's have a single chain in them. refinement has sensible .eff's and such for ncs, there's find_ncs_* output, but some .log's suggest it is not specified or used. so just generally, how is ncs used - e.g. is it strict, which would explain a single chain the output .pdb? -bryan

2 3

PDB Tools and CNS formatted structure factors
by Joseph Noel 22 Sep '10

22 Sep '10

Is there anyway with the PDB Tools to renumber protein chains? I have a model that doesn't correspond to the actual sequence, off by a couple of residues and just looking for a quick way to change that. Also, which module in Phenix is used to read in and convert CNS formatted ascii structure factors to the mtz read by phenix other then using ccp4. ___________________________________________________________ Joseph P. Noel, Ph.D. Investigator, Howard Hughes Medical Institute Professor, The Jack H. Skirball Center for Chemical Biology and Proteomics The Salk Institute for Biological Studies 10010 North Torrey Pines Road La Jolla, CA 92037 USA Phone: (858) 453-4100 extension 1442 Cell: (858) 349-4700 Fax: (858) 597-0855 E-mail: noel(a)salk.edu Web Site (Salk): http://www.salk.edu/faculty/faculty_details.php?id=37 Web Site (HHMI): http://hhmi.org/research/investigators/noel.html ___________________________________________________________

4 3

FYI : autobuild : "refine_pdb_in.mtz is missing" etc.
by Lepore, Bryan 22 Sep '10

22 Sep '10

here's a google-able post of a *solved problem* for the archive : phenix 1.6.4-486 (latest official release as of 21sep2010) building ends (FINISHED) with error messages (see below). some associated job setup details : phenix.autobuild with CL flag refine_eff_file_list=my_refinement_params.eff cartesian or torsion annealing, or start_temperature work OK in building/refinement. two segid's in .pdb, one chain id etc. SOLUTION : use the latest build. i used phenix-dev-533 GNU/Linux 32-bit. and thanks especially to Tom Terwilliger for the solution. [key error messages]: Failed to carry out AutoBuild_multiple_models: Sorry a subprocess has failed... Failed to carry out AutoBuild_build_cycle: Sorry...the file AutoBuild_run_1_/TEMP0/refine_pdb_in.mtz is missing...

1 0

autobuild : annealing
by Lepore, Bryan 21 Sep '10

21 Sep '10

[ phenix 1.6.4-486 ] looks like autobuild annealing is Cartesian. are there any (simple) tricks to get torsion, b/c seems there are no flags for it. -bryan

2 1

deleting jobs in phenix
by Sebastiano Pasqualato 20 Sep '10

20 Sep '10

Hi all, sorry for the maybe silly question, but I don't seem to find a way to delete jobs from a project in phenix (I'm using ver. 1.6.2-473 on a MacBook Pro). The 'job history' tab on the main window does not allow that, so I was wondering if there's an easy way to do it. Thanks a lot for the great software suite, keep it up guys, ciao s. -- Sebastiano Pasqualato, PhD IFOM-IEO Campus Dipartimento di Oncologia Sperimentale Istituto Europeo di Oncologia via Adamello, 16 20139 - Milano Italy tel +39 02 9437 5094 fax +39 02 9437 5990

2 4

Phenix renewal support
by Paul Adams 20 Sep '10

20 Sep '10

Dear Phenix Users, I hope that Phenix has been helpful in your work. We are in the process of submitting a renewal for funding of the Phenix project and I would be most appreciative of your support for our renewal. If you would be willing to write me a short note, before September 18th, in the form of a signed letter (PDF) that indicates how much Phenix has helped your research, that would be incredibly helpful for our renewal. You can send your letter to me at PDAdams(a)lbl.gov Cheers, Paul -- 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 [ KLMoody(a)lbl.gov ] [ 1-510-495-2506 ] --

2 1

R factor
by zhangh1＠umbc.edu 19 Sep '10

19 Sep '10

Hi, I want to calculate Fc based on the PDB file and reproduce the R factor in it. I first tried phenix.model_vs_data. It optimized ksol bsol and b_cart, and the output R factor (0.2143) is close enough to that in PDB (0.215). However, it didn't generate mtz files containing the calculated Fc. Then, I take the optimized ksol bsol and b_cart as the input and tried phenix.fmodel. I wrote my own code to calcuate the R factor. However, I got a lower R factor (0.2078). I think the reason might be my R calculation formula (sigma(fo-fc)/sigma(fo)) is wrong; or model_vs_data did more than just did flat solvent correction and anisotropic scaling. So, I am wondering whether phenix can: (1) Do the similar thing as model_vs_data (hope can be somehow faster:), but generate Fc-containing mtz file; (2) Do a simple R factor calculate to see whether my code is wrong (maybe answered already:) Thanks! Best Regards, Hailiang

2 7

generate Rfree flag for twinned crystals with 3 twin fractions
by Zheng Zhou 19 Sep '10

19 Sep '10

Dear all After reading the informative posts in mailing list, I started to deal with my first twinning dataset. The space group is P3. Xtriage found 3 twin laws. R vs R statistics indicated perfect twinning for two twin laws and one is twin with RPS. I realize that twin-related reflections should be in the same work or free group. Could phenix handle the Rfree flag correctly for me when there are more than one twin fraction identified? By applying one twin law, the Rfree dropped from 37% to 28% using phenix. When I use Refmac, the twin fractions are refined and the Rfree is down to 26%. I copied the Rfree flag from phenix to the ccp4 in this case. When I let ccp4 generate Rfree flag automatically, the Rfree is refined to 26% too. Should I label the Rfree at a higher space group? Thanks for your advice Best, Zheng (Joe) Zhou

2 1

Resolve NCS masks
by Florian Schmitzberger 18 Sep '10

18 Sep '10

Dear All, I would have a question regarding the use of masks in Resolve in Phenix for density modification: Is it possible to output the NCS and/or solvent masks Resolve is using during the Autosol run (similar as the options from the stand alone Resolve version) such as: mask_as_mtz ncs_mask_file mask.mtz protein_mask_file What are the best criteria (other than the final maps of course) within the output from Phenix.resolve to judge if the NCS (and/or solvent) mask was chosen adequately? (NCS_overlap?) I am mostly concerned if the masks chosen do include some of the more peripheral parts of my protein molecule and do not flatten these regions. Secondly, it is possible within Phenix Autosol wizard to supply a pdb file (with dummy atoms) for the mask for Resolve to construct the NCS mask? i.e. similar to set mask_pdb file= dummyatom.pdb Regards, Florian ----------------------------------------------------------- Florian Schmitzberger Biological Chemistry and Molecular Pharmacology Harvard Medical School 250 Longwood Avenue, SGM 130 Boston, MA 02115, US Tel: 001 617 432 5602

2 1

Ligand restraints editing with REEL
by Joseph Brock 18 Sep '10

18 Sep '10

Hi all, I'm working on a 1.5A structure that contains a novel ligand within the active site, which is an Arsenic derivative of Glutathione, using the Phenix GUI, version 1.6.4-486. I have been using the 'Molinspiration' webpage (http://www.molinspiration.com/docu/webme/) to generate a smiles string, which I then feed into ELBOW. Although the structure/chemical restraints of Glutathione should be unchanged by the covalent linkage to arsenic through the cysteinyl sulfur, the resulting restraints file consistently restricts one of the dihedrals to a value that is incompatible with the (quite clear) electron density. This occurs regardless of wether i use AM1 optimization, and is unresponsive to using native Glutathione (which refines into the corresponding density extremely well) as an input geometry file with any of the available 'Use geometry for' options (final geometry, initial geometry, etc). I then tried opening the output file in REEL and manually changing the corresponding dihedrals to those of glutathione, however after running simple optimization, the resulting restraints file has the same problem, and upon choosing the 'transer geometry and restraints' option (or any option) midway through running the 'eLBOW optimization' or 'AM1 optimization' options, Phenix consistently crashes. Any advice on how to transfer the origional geometry of native Glutathione to my new ligand effectively would be greatly appreciated. I'm a little confused as to the abilities of REEL, and how to affect them. Thanks so much in advance. Cheers, -Joe.

2 1

autobuild : MSE's keep coming back
by Lepore, Bryan 17 Sep '10

17 Sep '10

i used "semet=False" the .pdb has MET the .seq has M the output models, however, keep getting MSE's in them, complete with selenium. am i missing something? -bryan

5 6

Re: [phenixbb] Something that came up on the CCP4 list...
by Simon Kolstoe 16 Sep '10

16 Sep '10

Just for the record I've had it make a big difference a number of times as well with both refmac & shelx. I was under the impression that phenix.refine was "cleverly" allowing me to use that trick by just starring the add riding hydrogens in the def file however it seems not - something worth knowing! Simon On 16 Sep 2010, at 20:01, Phil Jeffrey wrote: > I've had it make a big difference in certain circumstances. And in > other circumstances not so much. Definitely worth a try. I do not > use the technique routinely. > > Phil Jeffrey > Princeton > > > > On 9/16/10 11:49 AM, Simon Kolstoe wrote: >> So just to clarify further - if i am having problems with the final >> geometry of a model (eg clashscore as reported by molprobity), is >> it an >> appropriate trick to add hydrogens and then refine in phenix.refine >> using the riding hydrogens option, or does this not really change >> anything? >> >> Simon >> >> On 16 Sep 2010, at 14:26, Pavel Afonine wrote: >> >>> Hi Simon, >>> >>>> Could I just seek some clarification? I was under the impression >>>> that >>>> phenix.refine calculated riding hydrogens by default and took these >>>> into account as a restraint on the geometry? >>> >>> if you want hydrogen atoms to be used in refinement (as "riding >>> model" >>> or explicit) your need to make sure the PDB file that you give to >>> phenix.refine contains H atoms. That is you need to add them using >>> programs like phenix.reduce or phenix.ready_set or using >>> corresponding >>> functionality in the GUI. >>> >>> phenix.refine will not add them internally as part of refinement IF >>> hydrogen atoms are NOT present in input PDB file. >>> >>> Adding hydrogen atoms as 'riding model' does not add any refinable >>> parameters, but the scattering contribution from them is included >>> most >>> of the time. This is why we do not use H atoms silently but require >>> them to be in input PDB file - simply to make sure the reported >>> statistics (such as R-factors) is easily reproducible. >>> >>> All the best! >>> Pavel. >>> _______________________________________________ >>> phenixbb mailing list >>> phenixbb(a)phenix-online.org >>> http://phenix-online.org/mailman/listinfo/phenixbb >> >> _______________________________________________ >> phenixbb mailing list >> phenixbb(a)phenix-online.org >> http://phenix-online.org/mailman/listinfo/phenixbb >

2 1

LABIN labels for Autosol/resolve NCS scores
by Florian Schmitzberger 16 Sep '10

16 Sep '10

Dear All, I would have a number of Phenix related questions: Firstly, I used the anisotropy server ( to anisotropically truncate my mtz file with anomalous signal (SAD data). My question is if the labels (for F+/F- and SIGF+/SIGF- are correctly recognized by Phenix (Autosol). The columns labels are: F_ref for FP; QSIGF_ref for SIGFP; GF(+)pk for F(+); LSIGF(+)pk for SIGF(+), GF(-)pk for F(-); LSIGF(-)pk for SIGF(-). It seems QSIGF and GF(+)pk/LSIGF(+)pk etc. are , but I am uncertain if it recognizes F_ref, because for the LABIN INPUT LABELS it says "NONE" for"F". Phenix does not seem to detect the anistropy ("Not using aniso- corrected data files as the range of aniso b is only 19.52 and 'correct_aniso' is not set); so I suppose I should use the correct_aniso=True keyword. Secondly, what are good values for the evaluation of the NCS mask and NCS averaging Resolve is using in Phenix; NCS overlap values (e.g. in my case 1.1) and NCS_CORR (in my case 0.52) in the "resolve.scores" file. Does an NCS overlap value of greater than 1 make sense? I am working with a (anisotropic) dataset with 4.1/4.5 A resolution, 78 % solvent content, 2-fold NCS; Perhaps a naive question: It is not possible to use B-factor sharpened maps for autobuilding with Phenix Autobuild is it? Thank you for any comments! Regards, Florian ----------------------------------------------------------- Florian Schmitzberger Biological Chemistry and Molecular Pharmacology Harvard Medical School 250 Longwood Avenue, SGM 130 Boston, MA 02115, US Tel: 001 617 432 5602

2 3