July 2009 - phenixbb - phenix-online.org

Re: [phenixbb] phenix.automr ensemble/component question
by Thomas C. Terwilliger 24 Jul '10

24 Jul '10

Hi Scott, Your first command looks good to me. It is saying that the composition of the crystal is 2 copies of each chain: component_1.component_copies=2 component_2.component_copies=2 and you want to search for 2 copies of each of your search models: ensemble_1.copies_to_find=2 ensemble_2.copies_to_find=2 All the best, Tom T > Hi Phenixers, > > I'm a bit confused about the proper syntax for my phenix.automr job. > This is probably a simple concept, but I'm running on 3 hours of sleep > as my kid was up all night, and my synapses aren't working so good. > > I expect 2 copies of my protein per AU > > My protein is composed of two domains > > I have a single PDB search model for each domain > > is either of these commands correct? > > phenix.automr data=mydata.sca \ > ensemble_1.coords=domain1.pdb \ > ensemble_1.identity=95 \ > ensemble_1.copies_to_find=2 \ > component_1.mass=40579 \ > component_1.component_copies=2 \ > ensemble_2.coords=domain2.pdb \ > ensemble_2.identity=95 \ > ensemble_2.copies_to_find=2 \ > component_2.mass=13410 \ > component_2.component_copies=2 > > > or > > > phenix.automr data=mydata.sca \ > ensemble_1.coords=domain1.pdb \ > ensemble_1.identity=95 \ > ensemble_1.copies_to_find=1 \ > component_1.mass=40579 \ > component_1.component_copies=2 \ > ensemble_2.coords=domain2.pdb \ > ensemble_2.identity=95 \ > ensemble_2.copies_to_find=1 \ > component_2.mass=13410 \ > component_2.component_copies=2 > > > essentially I'm not really using an ensemble of search models, I only > have a single PDB for each domain. So should ensemble_1.copies_to_find > and ensemble_2.copies_to_find be set to 1 or 2? > > > Thanks, > Scott > > ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ > Scott Classen, Ph.D. > SIBYLS Beamline 12.3.1 > http://bl1231.als.lbl.gov > Advanced Light Source > Lawrence Berkeley National Laboratory > 1 Cyclotron Rd > MS6R2100 > Berkeley, CA 94720 > C) 510.206.4418 > O) 510.495.2697 > Beamline) 510.495.2134 > ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ > > > > > _______________________________________________ > phenixbb mailing list > phenixbb(a)phenix-online.org > http://www.phenix-online.org/mailman/listinfo/phenixbb >

2 2

Re: [phenixbb] sigma cutoff
by Kay Diederichs 01 Dec '09

01 Dec '09

> > > Hi, > > According to the Phenix manual, phenix.refine uses no sigma cutoff by > default. Since I haven't changed this manually, why do I see in my pdb > file MIN(FOBS/SIGMA_FOBS) 1.340? > > Thanks for help, > Karolina > > This one I can answer because I asked myself the same question recently. The answer is that Phenix does not impose any sigma cutoff, but rather your data are such that the minimum value of Fobs/Sigma_Fobs is 1.34 . You can verify this for yourself by going through your reflection file, line by line, and calculating Fobs/Sigma_Fobs for each reflection. You will find that all your reflections have a Fobs/Sigma_Fobs larger or equal to 1.34 . So it has nothing to with phenix; phenix just reports a property of your data. An easier way than doing it "manually" is using CCP4 sftools for calculating a new column "ratio", writing out the resulting MTZ file, and looking at this with mtzdmp. HTH, Kay -- Kay Diederichs http://strucbio.biologie.uni-konstanz.de email: Kay.Diederichs(a)uni-konstanz.de Tel +49 7531 88 4049 Fax 3183 Fachbereich Biologie, Universitaet Konstanz, Box M647, D-78457 Konstanz

2 1

Autosol error
by Francis E Reyes 01 Aug '09

01 Aug '09

All When issuing phenix.autosol native.data=scala_072809B_3_001.mtz deriv.data=scala_072909A_3_001.mtz deriv.atom_type=Ir deriv.sites=5 deriv.inano=noinano native.labels='IMEAN SIGIMEAN' deriv.labels='IMEAN SIGIMEAN' I get the error **************************************** AutoSol Input failed For SIR datasets please enter two datafiles or datasets.(native +deriv) Please treat MIR as a series of SIR datasets ...to specify a particular set of data you can specify one of the following (command-line format shown, for a script take out the '=' and use the keyword 'group_labels_list' instead): native.labels='IMEAN SIGIMEAN' deriv.labels='IMEAN SIGIMEAN' native.labels='I(+) SIGI(+) I(-) SIGI(-)' deriv.labels='I(+) SIGI(+) I(-) SIGI(-)' **************************************** Any ideas? thanks FR --------------------------------------------- Francis Reyes M.Sc. 215 UCB University of Colorado at Boulder gpg --keyserver pgp.mit.edu --recv-keys 67BA8D5D 8AE2 F2F4 90F7 9640 28BC 686F 78FD 6669 67BA 8D5D

2 2

combining models from multi-start simulated annealing
by James Whittle 30 Jul '09

30 Jul '09

Is there a way to use AutoBuild to combine a set of models created by multi-start simulated annealing? I've used a random number generator to seed different SA runs at two temperatures, as suggested on this bulletin board before ... but can combine the best parts of these models automatically? --James

2 1

Integrating Phenix with Torque
by Jeremy Mann 29 Jul '09

29 Jul '09

A new user of mine wants to run Phenix on two of our clusters. We use Torque/Maui as our batch scheduler and would like to integrated their usage of Phenix. They want to run something like: phenix.autobuild data=data.hkl seq_file=seq.dat model=model.pdb multiple_models=True multiple_models_number=20 n_cycle_rebuild_max=5 nbatch=20 nproc=20 The documentation mentions briefly the ability to use batch schedulers but doesn't give examples of submission scripts. -- Jeremy Mann jeremy(a)biochem.uthscsa.edu University of Texas Health Science Center Bioinformatics Core Facility http://www.bioinformatics.uthscsa.edu Phone: (210) 567-2672

2 1

Phaser RMS - What happens when it is too high
by Peter Grey 29 Jul '09

29 Jul '09

Hi all, In Phaser RMS is a crucial parameter. I had several cases where underestimation of RMS resulted in negative LLG. The cure was simply increasing the RMS. But I wonder, are there negative effects to a too high RMS. What are the dangers of giving Phaser a too high RMS ? Are there any indications in the output for that. This is also important with EM maps as models since in that case estimation of RMS is difficult. Many thanks in advance for your advice and insight, Peter.

2 1

phenix.refine issue
by Dennis Wolan 28 Jul '09

28 Jul '09

I am having an issue with running a simple phenix.refine that I have never encountered before. I get a list of issues and the program quits with: Traceback (most recent call last): File "/usr/local/phenix-1.4-3/phenix/phenix/command_line/ refine.py", line 9, in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/usr/local/phenix-1.4-3/phenix/phenix/refinement/ command_line.py", line 68, in run log=log) File "/usr/local/phenix-1.4-3/phenix/phenix/refinement/ command_line.py", line 231, in __init__ reflection_file_server = self.reflection_file_server() File "/usr/local/phenix-1.4-3/phenix/phenix/refinement/ command_line.py", line 566, in reflection_file_server err=self.log) File "/usr/local/phenix-1.4-3/iotbx/reflection_file_utils.py", line 393, in __init__ force_symmetry=self.force_symmetry)) File "/usr/local/phenix-1.4-3/iotbx/reflection_file_reader.py", line 197, in as_miller_arrays base_array_info=base_array_info) File "/usr/local/phenix-1.4-3/iotbx/mtz/__init__.py", line 434, in as_miller_arrays and flex.min(column_group.sigmas()) > 0): RuntimeError: min() argument is an empty array I have no idea if the issue is with my pdb or mtz. I'm sure it's something very simple that I just am not identifying. Thanks in advance, Dennis

2 1

phaser MR
by sbiswas2＠ncsu.edu 27 Jul '09

27 Jul '09

Hi all, I was trying to solve a structure by molecular replacement. I scaled the data in P4222 space group (resolution 2.7A) with two molecules in the assymmetric unit (molecule A and B) I ran phaser with my model and got a Zscore of 5.1. When I look at the map that I got from phaser I could easily see good electron density for both molecules, However upon inspection of the electron density map there were considerable interaction or clashes with molecule B and a symmetry atom. Molecule A had no clashes however with the symmetry atoms. I was wondering if anyone knows how to resolve this. Could it be a problem of space group. The statistics are good for space group P4222 and the I/sigI was good till 2.7A. Any advice is appreciated, Shya

3 6

question about SAD+MR
by crystallogrphy 25 Jul '09

25 Jul '09

Dear all, I have a Se-Met SAD data which diffracted to 2.8A (peak wavelength). The space group is P1. The map from SAD autosol and autobuild still gives crappy map, I can only see one continous alpha-helix density with all other density fragmented and difficult to identify. I suppose it is because of high mosaicity (1.5), low data redundency, and some radiation damage. Now I want to use mocular replace (sequence identify ~0.26) plus SAD data to give better map. 1. The result from MR is shown below and the density is poor as the sole SAD map: REMARK Log-Likelihood Gain: 121.374 REMARK RFZ=3.4 TFZ=100.0 PAK=0 LLG=29 RFZ=3.7 TFZ=5.4 PAK=0 LLG=69 RFZ=6.1 TFZ=4.0 PAK=0 LLG=97 RFZ=5.6 TFZ=5.0 PAK=0 LLG=118 LLG=121 2. I then tried to combine MR model and SAD data set phenix.autosol seq_file=zntb-sd.seq sites=14 atom_type=Se data=output-peak.sca sg="P1" cell="44.355 56.687 58.451 61.061 82.013 82.740" resolution=3.0 input_partpdb_file=AutoMR_run_1_/MR.1.pdb The program failed to run with such error message: Failed to carry out AutoSol_score_ha_solutions: [Errno 2] No such file or directory: 'AutoSol_run_3_/TEMP0/phaser_mtz.hand.mtz' I do not know where to get this phaser.mtz.hand.mtz file or how to generate it. Does anyone know? -- Qun Wan Department of Pharmacology, Case Western Reserve University 10900 Euclid Avenue, Cleveland, OH 44106 email: qunwan1(a)gmail.com lab phone: 216-368-3337

2 1

new gui on osx (10.5)
by Miguel Ortiz Lombardia 25 Jul '09

25 Jul '09

Dear all, We would like to run the new phenix GUIs from a remote computer. Both server and clients are OSX 10.5 (leopard) and they have the same versions of software (X11, python, phenix... all them pretty new versions) When trying to do so, we invariably get this error: This program needs access to the screen. Please run with 'pythonw', not 'python', and only when you are logged in on the main display of your Mac. I have searched the internet for similar problems and I found an old (one year) unresolved question on the same/similar issue: http://www.phenix-online.org/pipermail/phenixbb/2008-July/001129.html In our case, the /usr/local directory is not mounted over the NFS, it is local. As in that case, though, changing the LIBTBX_PYEXE_BASENAME to pythonw (for example in the phenix.autosol_gui wrapper) did not solve the issue. However, it is not clear to me that this variable has to be changed in a single place: it is found in many files. My question is: was this issue solved and I missed that? If not, anyone has a clue on how to solve the problem? Just to make it clear, the GUIs work perfectly locally in the server computer, that is, if not exported. Thanks a lot! Miguel -- Miguel Ortiz Lombardía Architecture et Fonction des Macromolécules Biologiques (UMR6098) CNRS, Universités d'Aix-Marseille I & II Case 932 163 Avenue de Luminy 13288 Marseille cedex 9 France Tel : +33(0) 491 82 55 93 Fax: +33(0) 491 26 67 20 e-mail: miguel.ortiz-lombardia(a)afmb.univ-mrs.fr Web: http://www.pangea.org/mol/spip.php?rubrique2

2 2

sigma cutoff
by dziuba＠amu.edu.pl 24 Jul '09

24 Jul '09

Hi, According to the Phenix manual, phenix.refine uses no sigma cutoff by default. Since I haven't changed this manually, why do I see in my pdb file MIN(FOBS/SIGMA_FOBS) 1.340? Thanks for help, Karolina

2 1

ENSEmble parameters when the model is an EM map
by Peter Grey 24 Jul '09

24 Jul '09

Dear benevolent experts, I am trying to use Phaser for molecular replacement with EM-map as a model. The EM map, originally in a small P1 cell, was put in large P1 using CCP4's MAPROT and a mask that was defined by a certain cut-off level. Could you please advise me how to derive the EXTEnt and CENTre parameters needed for ENSEmble building ? Does the simple way of giving the size of the original P1 as the EXTEnt and its centre as CENTre sound reasonable ? I would be grateful for any advice and comments, Peter.

2 3

anisotropic ADP factor refinement and optimize_wxc
by Joyce, Gordon M. (NIH/NIAID) [F] 22 Jul '09

22 Jul '09

Dear All, I have just recently started using phenix.refine and I hope I can get some advice. I am trying to refine my structure with anisotropic ADP factors while also optimizing the weighting parameters.When I give the following command, anisotropic refinement is carried out but optimization is not carried out as far as I can see (ie the rms bond and angles are not reduced, which is one of my concerns at the moment). phenix.refine protein1.pdb 072109checks.mtz ordered_solvent=true strategy=individual_adp adp.individual.anisotropic="not element H" optimize_wxc=true --overwrite pdb output values REMARK Start: r_work = 0.1952 r_free = 0.2388 bonds = 0.015 angles = 1.646 REMARK Final: r_work = 0.1413 r_free = 0.2448 bonds = 0.015 angles = 1.646 If I give the following command, optimization occurs and rms bond and angle values decrease. phenix.refine protein1.pdb 072109checks.mtz optimize_wxc=true --overwrite REMARK Start: r_work = 0.2048 r_free = 0.2388 bonds = 0.015 angles = 1.646 REMARK Final: r_work = 0.1917 r_free = 0.2352 bonds = 0.009 angles = 1.301 Is it possible to refine both anisotropic B factors and optimize ADP factors in the same run. Any advice would be great. Thanks, Gordon

2 2

Phenix.refine coordinate error estimate discrepancies
by Phil Jeffrey 22 Jul '09

22 Jul '09

There seems to be a large discrepancy between the likelihood-based coordinate error estimate depending on whether the intensities or structure factor magnitudes are used, which doesn't make any sense to me. If I do: % phenix.refine rebuild.pdb data_withI.mtz phenix.refine picks up Imean and SigImean from the MTZ and refines the structure xray_data { file_name = "data_withI.mtz" labels = "IMEAN,SIGIMEAN" r_free_flags { file_name = "data_withI.mtz" label = "FreeR_flag" test_flag_value = 0 } then I get: | maximum likelihood estimate for coordinate error: 1.39 A But if I create an MTZ file without the intensities by using CAD on the original file, and use the same coordinate file: % phenix.refine rebuild.pdb data_withoutI.mtz phenix.refine picks up the structure factors instead: xray_data { file_name = "data_withoutI.mtz" labels = "F,SIGF" r_free_flags { file_name = "data_withoutI.mtz" label = "FreeR_flag" test_flag_value = 0 } } then this seems to be the right ballpark for a 1.75 Angstrom structure.... | maximum likelihood estimate for coordinate error: 0.28 A Similar R-factors, but there is a factor of 10 difference in the scale reported here: I's: r_work= 0.2266 r_free= 0.2570 ksol= 0.41 Bsol= 50.00 scale= 0.097 F's: r_work= 0.2221 r_free= 0.2557 ksol= 0.41 Bsol= 50.00 scale= 0.965 Version: 1.4 Release tag: 3 cctbx tag: 2008_12_07_1353 Platform: mac-intel-osx osx-10.5.6 User: phil Cheers, Phil Jeffrey Princeton

2 1

Assertion Error when using kicked maps in solvent update
by Schubert, Carsten [PRDUS] 21 Jul '09

21 Jul '09

Hi In 1.4.106 I get the following Assertion Error when I use kicked maps (refinement.ordered_solvent.use_kick_maps=True) in the solvent update section: Traceback (most recent call last): File "/raidnet/xdata1/programs/phenix-1.4-106/phenix/phenix/command_line/refi ne.py", line 19, in <module> command_line.run(command_name="phenix.refine", args=sys.argv[1:]) File "/raidnet/xdata1/programs/phenix-1.4-106/phenix/phenix/refinement/comman d_line.py", line 89, in run call_back_handler=call_back_handler) File "/raidnet/xdata1/programs/phenix-1.4-106/phenix/phenix/refinement/driver .py", line 1125, in run call_back_handler = call_back_handler) File "/raidnet/xdata1/programs/phenix-1.4-106/phenix/phenix/refinement/strate gies.py", line 358, in refinement_machine log = log) File "/raidnet/xdata1/programs/phenix-1.4-106/cctbx_project/mmtbx/solvent/ord ered_solvent.py", line 205, in __init__ use_kick_map = self.params.use_kick_maps).peaks_mapped() File "/raidnet/xdata1/programs/phenix-1.4-106/cctbx_project/mmtbx/solvent/ord ered_solvent.py", line 394, in find_peaks log = self.log) File "/raidnet/xdata1/programs/phenix-1.4-106/cctbx_project/mmtbx/find_peaks. py", line 98, in __init__ crystal_gridding_tags = fft_map.tags() File "/raidnet/xdata1/programs/phenix-1.4-106/cctbx_project/cctbx/maptbx/__in it__.py", line 152, in tags return crystal_gridding_tags(self) File "/raidnet/xdata1/programs/phenix-1.4-106/cctbx_project/cctbx/maptbx/__in it__.py", line 158, in __init__ assert gridding.symmetry_flags() is not None AssertionError Running from the same .def file w/o kicked maps works fine. Let me know if you need more info to debug this. Cheers, Carsten

2 1

Duration of a phenix.elbow run
by Dirk 21 Jul '09

21 Jul '09

Hi, I'm trying to create restraints for a carbohydrate ligand using phenix.elbow. If I run it without the --opt command everything works fine, but I'm not completely happy with the geometry of my ligand after refinement. I tried to get an optimized cif by using --opt command, but the program just runs up to the following output: MoleculeClass : H:38 C:23 N: 1 O:19 (PDB format) 81 atoms 83 bonds 151 angles 238 dihedrals -------------------- Time Energy delta E RMS Max Step size and then nothing happens for several hours. Does the optimization needs such a long time and if I let it run overnight (or perhaps longer) I will get an output or is something going wrong? The input is a pdb (I tested several) and I tried both the stable phenix release (1.3) and the current version (1.4). Any help would be very appreciated Thanks Dirk

7 9

Unable to produce kicked maps
by Schubert, Carsten [PRDUS] 18 Jul '09

18 Jul '09

Hi , I am trying to make some kicked maps in 1.4.106. At the end of the run all I am getting is Sorry: Duplicate mtz_label_amplitude FOFCWT Any workaround, except not to use kicked maps? Here the diff_params output for easer debugging. Cheers Carsten <<diff.def>>

2 2

autosol
by sbiswas2＠ncsu.edu 17 Jul '09

17 Jul '09

Hi, I was trying to run autosol using gui however it crashed with the following message: Floating-point error (Python and libc call stacks above) This crash may be due to a problem in any imported Python module, including modules which are not part of the cctbx project. To disable the traps leading to this message, define these environment variables (e.g. assign the value 1): BOOST_ADAPTBX_FPE_DEFAULT BOOST_ADAPTBX_SIGNALS_DEFAULT Is this again a version problem? I am running version 1.4-6 Any help will be appreciated, Shya

2 1

rigid body refinement in phenix.refine
by Pavel Afonine 17 Jul '09

17 Jul '09

To PHENIX users who asked about this: I recall several people (unfortunately, I don't remember who exactly) asked me about the details of rigid-body refinement implementation in phenix.refine (MZ protocol) and I promised to let them know when the corresponding paper is published. So here it is: J. Appl. Cryst. (2009). 42, 607-615 Automatic multiple-zone rigid-body refinement with a large convergence radius P. V. Afonine, R. W. Grosse-Kunstleve, A. Urzhumtsev and P. D. Adams http://journals.iucr.org/j/issues/2009/04/00/dd5048/index.html Pavel.

1 0

error message
by zachcp 16 Jul '09

16 Jul '09

A Python error was detected. This is almost certainly a bug; please send an email containing this message to phenixbb(a)phenix-online.org so we can fix it. Thanks! SorryException : No array of observed xray data found. Traceback: Traceback (most recent call last): File "/usr/local/phenix-1.4-3/phenix/wxGUI2/Xtriage.py", line 158, in OnExcept raise event.data SorryException: No array of observed xray data found.

3 2

Average B-factor and ML error
by Leo Sazanov 16 Jul '09

16 Jul '09

Hi, phenix.refine outputs "WILSON B" value into PDB header after refinement, which is not the same as average B-factor of all the atoms in the structure, and not exactly the same as Wilson B-factor from SCALA scaling. How exactly is Wilson B-factor for phenix header calculated and would it be possible to output there also the average B-factor of all the atoms in the structure? Also, deposition into PDB requires only Luzzati coordinate errors, while phenix produces only ML based ones. Presumably ML estimate is better than Luzzati?, but it looks like often people just ignore PDB input for errors. It seems that it would it be useful in future to have ML option in the PDB submission process? Or should Phenix produce Luzzati errors as well? (I know that one can calculate av. B-fac in other programs or ask PDB to add some data to the deposition, but if the idea is to streamline the process...) Thanks, -- Dr. Leonid A. Sazanov Research group leader Medical Research Council Mitochondrial Biology Unit Wellcome Trust / MRC Building Hills Road Cambridge CB2 0XY WEB: www.mrc-mbu.cam.ac.uk Tel: +44-1223-252910 Fax: +44-1223-252915 -- Dr. Leonid A. Sazanov Research group leader Medical Research Council Mitochondrial Biology Unit Wellcome Trust / MRC Building Hills Road Cambridge CB2 0XY WEB: www.mrc-mbu.cam.ac.uk Tel: +44-1223-252910 Fax: +44-1223-252915

2 3

Anyone have a phenix.resolve script for after Phaser ?
by Phil Jeffrey 15 Jul '09

15 Jul '09

Although my usual m.o. is SHELX -> SHARP I'm exploring HYSS -> PHASER -> RESOLVE to see if the maps are any more interpretable. I think I've got good scripts for the first two steps but I'm not sure about my resolve script which stands at: phenix.resolve << EOD hklin myp3_auto.mtz LABIN FP=F SIGFP=SIGF PHIB=PHIB FOM=FOM HLA=HLA HLB=HLB HLC=HLC HLD=HLD hklout resolve.mtz solvent_content 0.45 seq_file myp3.pir EOD This runs, but I could not tell from delving into the documentation if this was ideal or not. Thanks Phil Jeffrey (a Luddite and not using the GUI) Princeton

2 2

Re: [phenixbb] refine occupancy of partially overlapping ligands
by Chris Ulens 15 Jul '09

15 Jul '09

Hi, I'm trying to refine the occupancy of 2 ligands that occupy the same binding pocket in 2 *partially* overlapping orientations. Let's assume I would define the ligands in 2 separate "chains" termed A and B. Could the experts point me to the proper syntax to define the correct group restraints to refine the occupancy? Thank you. Chris +++++++++++++++++++++++++++++++ Chris Ulens, Ph.D. Lab of Structural Neurobiology Division of Pharmacology Campus Gasthuisberg, ON1 Herestraat 49, PB 601 B-3000 Leuven Belgium Tel +32 16 345812 Fax +32 16 345699 Telepathy 2.69 GHz +++++++++++++++++++++++++++++++

2 1

can't find phaser_mtz.hand.mtz?
by Pei-Chun Lin 14 Jul '09

14 Jul '09

Hi, I am the beginner of crystallography and have a 3.5A map generated by MAD method using phenix GUI. I already traced and built 2/3 of my protein. I wish to use this partial model combined with SAD to generate a new map and see if the new map can help me to deal with the rest of the structure. Does it sound reasonable to you? However, after running AutoSol. I got this error message: Failed to carry out AutoSol_score_ha_solutions: [Errno 2] No such file or directory: 'AutoSol_run_1_/TEMP0/phaser_mtz.hand.mtz' What should I do now? Can't phenix automatically calculate it? PC _________________________________________________________________ Windows Live™: Keep your life in sync. Check it out! http://windowslive.com/explore?ocid=TXT_TAGLM_WL_t1_allup_explore_012009

2 1

tightening geometry
by Patrick Loll 11 Jul '09

11 Jul '09

Hi, I'm a little baffled. We have a 3 Å structure, and while the R values seem OK, I'd like to tighten up the rmsd values for the geometry. However, altering wxc_scale seems to have no effect. Specifics: In a clean directory, I run the following command phenix.refine X.mtz Y.pdb Z1.cif Z2.cif Z3.cif wxc_scale=0.??? strategy=group_adp (the cif files are libraries for different ligands) When I vary wxc_scale from 0.5 to 0.005, the geometry doesn't change at all, which makes no sense to me: For wxc_scale=0.5 r_work = 0.2672 r_free = 0.2900 bonds = 0.028 angles = 3.151 For wxc_scale=0.05 r_work = 0.2672 r_free = 0.2900 bonds = 0.028 angles = 3.151 For wxc_scale=0.005 r_work = 0.2672 r_free = 0.2900 bonds = 0.028 angles = 3.151 ...you get the idea. These were run with phenix-1.4-3. With phenix-1.4-56, we get subtly different r values, but the bonds and angles numbers are the same as shown above (to 3 decimal places!), and the numbers don't change upon changing wxc_scale. Am I missing something here? Is there a better way to tighten up the geometry? Thanks, Pat --------------------------------------------------------------------------------------- Patrick J. Loll, Ph. D. Professor of Biochemistry & Molecular Biology Director, Biochemistry Graduate Program Drexel University College of Medicine Room 10-102 New College Building 245 N. 15th St., Mailstop 497 Philadelphia, PA 19102-1192 USA (215) 762-7706 pat.loll(a)drexelmed.edu

3 2