- phenixbb - phenix-online.org

ADP: fix_wxu and wxu_scale
by fulvio saccoccia 25 Jan '11

25 Jan '11

Dear Phenix users and developers, I'm trying to refine a 3 A resolution structure. In order to better understand how can I use wxu_scale and fix_wxu options in phenix.refine, may you explain me their use, since I have a lot of problems with ADPs? In fact, I have big variations between atoms of the same residue. So, is it better to use wxu_scale or to impose fix_wxu to refine the ADPs? Thanks to all Fulvio Saccoccia

2 1

Problem with loop fitting
by Kendall Nettles 25 Jan '11

25 Jan '11

Dear PHENIX developers, I am trying fit_loops for the first time and getting error message: "No suitable gap found matching residues" I'm not sure what the problem is. Best regards, Kendall Nettles knettles$ phenix.fit_loops pdb_in=Liganded_overall_best-coot-1.pdb mtz_in=cycle_best_refine_map_coeffs_1.mtz seq_file=seq.dat start=331 end=336 chain_id=A # fit_loops # # Fit loops in a PDB file with gaps # Type phenix.doc for help Input sequence file is seq.dat Reading sequence from seq.dat Fit_loops: Fit a missing loop in a model. Map from: TEMP_cycle_best_refine_map_coeffs_1.mtz using labin FP=2FOFCWT PHIB=PH2FOFCWT Model from: Liganded_overall_best-coot-1.pdb LABIN LINE: FP=2FOFCWT PHIB=PH2FOFCWT No suitable gap found matching residues 331 - 336 chain: A No final model obtained

2 3

Ligandfit bug report
by Kendall Nettles 25 Jan '11

25 Jan '11

We are generating pdb files with chemdraw3D, which turn out to not be "real" pdb files. PHENIX.ELBOW fixes them, but leaves some of the atom ID numbers blank, such as Oxygen #2, and Carbon #18: HETATM 17 C16 LIG 1 -83.510 -0.856 0.582 1.00 20.00 C HETATM 18 O01 LIG 1 -86.640 -2.401 -0.357 1.00 20.00 O HETATM 19 O LIG 1 -75.199 2.925 2.145 1.00 20.00 O HETATM 20 C17 LIG 1 -83.058 -0.188 1.872 1.00 20.00 C HETATM 21 C LIG 1 -76.676 0.780 -1.835 1.00 20.00 C After editing it to add in the missing numbers, I re-ran elbow, and found that the resulting optimized ligand had imploded, producing altered atom connectivity. If I ran it with --final_geometry it came out OK. If you would like files to reproduce the error, I am happy to send them. Best, Kendall Nettles

2 1

autobuild and buccaneer reference
by Bryan Lepore 21 Jan '11

21 Jan '11

1.7-629 it appears that autobuild does not have buccaneer implemented such that _buccaneer's_ reference structure method is used - N.B. i do not mean the "reference structure" as used by resolve/autobuild. any comment appreciated, thanks -Bryan

2 1

Stereochemistry statistics
by miket＠chem.ucla.edu 20 Jan '11

20 Jan '11

Hi All, I have a very naive question, but nobody I have asked has been able to give me a confident answer. When phenix calculates the stereochemistry statistics for a model, two of the statistics given are RMS(chirality) and RMS(planarity). What specific measurements are being calculated to produce these statistics, and what are the units for the numbers given? For example, is RMS(planarity) telling me the average distance out of plane for an atom that should be in plane? If this is the case I imagine the units are angstroms? Thanks for the explanation, Mike Thompson -- Michael C. Thompson Graduate Student Biochemistry & Molecular Biology Division Department of Chemistry & Biochemistry University of California, Los Angeles miket(a)chem.ucla.edu

3 2

Ten days left to apply, and many seats available, for RapiData 2011, our course on Data Collection and Structure Solving at the NSLS.
by Robert Sweet 20 Jan '11

20 Jan '11

This may be the last reminder that we are offering RapiData 2011, the thirteenth offering of our popular course: Rapid Data Collection and Structure Solving at the NSLS: A Practical Course in Macromolecular X-Ray Diffraction Measurement The course will be held 3-8 April 2011. Students could be at any level from advanced undergraduate to full professor. The course should accommodate 48 students total. All students are encouraged to bring their own specimens for data collection, and to bring old data for the data-reduction and structure-solving tutorials. Please read the Course Announcement at http://www.bnl.gov/RapiData/. You'll see that many experts in the field will be available for lectures and tutorials. You'll find the application materials on the Course Application tab at this site. For the ninth time we will hold a short lecture course on the fundamentals of crystallography for roughly five hours on Sunday 3 April. The body of the RapiData course really requires that students have a healthy knowledge of crystallography. For potential students who have some experience but are shaky about fundamentals, this course will help. There will be a small additional fee for the fundamentals course, to pay for Saturday night accomodations and food on Sunday morning and noon. Latin American Scientists: Several scholarships are available, from the International Union of Crystallography, to pay partial travel and subsistence costs for Latin-American students and junior faculty (under 40 yrs). Please apply for the course, and then contact R. Sweet (sweet(a)bnl.gov) if you are interested in applying for a scholarship. In accordance with the standards of the International Union of Crystallography, we observe the basic policy of non-discrimination, affirming the right and freedom of scientists to associate in international scientific activity without regard to such factors as citizenship, religion, creed, political stance, ethnic origin, race, colour, language, age, or gender, in accordance with the Statutes of the International Council for Science. At this course no barriers will exist beyond the application procedure that would prevent the participation of bona fide scientists. Please apply or send your students to our course, Bob Sweet, Sal Sclafani, and Alex Soares Course Announcement at http://www.bnl.gov/RapiData/ ========================================================================= Robert M. Sweet E-Dress: sweet(a)bnl.gov Group Leader, PXRR: Macromolecular ^ (that's L Crystallography Research Resource at NSLS not 1) http://px.nsls.bnl.gov/ Biology Dept Brookhaven Nat'l Lab. Phones: Upton, NY 11973 631 344 3401 (Office) U.S.A. 631 344 2741 (Facsimile) =========================================================================

1 0

Submitting autobuild to cluster
by Kendall Nettles 20 Jan '11

20 Jan '11

Dear PHENIX developers, I'm preparing a bunch autobuild runs. If I submit 20 autobuild jobs to our 100 node cluster, what would you recommend for nproc and nbatch? How about for 10 jobs at once? Thanks! Kendall Nettles

2 1

Defining an Intermolecular link
by Ryan Timothy Vanderlinden 19 Jan '11

19 Jan '11

Hi, I am working on a complex of three molecules in which, I believe, there is a covalent bond between the sulfur atom of a cysteine on molecule "A" and the c-terminal carbon of molecule "C". There are other structures deposited that have a similar linkage. I have looked at these structures and included a replica description of the supposed link in my input pdb file for refinement. I have redefined this c-terminal residue as a ligand and used phenix elbow do create a .cif file for it. Even when I use the modified pdb file that describes the linkage and the cif file that describes the ligand, the intermolecular link falls apart and the ligand (c-terminal amino acid) separates from the rest of the chain. The resulting mtz file still shows density (2Fo-Fc) supporting the link but the ligand is not in the density. The fact that the ligand moves away from where I would expect it is evidence that modeling the intermolecular linkage is really incorrect or I have not included all the proper files to define this linkage and hold it together during refinement. What are all the necessary files required hold an intermolecular linkage of this nature together during refinement? Are there certain commands that I might be leaving out from my refine script that are necessary to recognize such linkages? Thanks in advance, Ryan VanderLinden University of Utah

3 2

labelit documentation
by Bryan Lepore 19 Jan '11

19 Jan '11

in 1.7-629 (not the latest ), I cannot find any labelit documentation (however, i found some on the labelit site OK). could anyone comment on that - e.g. is there a timeframe, etc. thanks, -Bryan p.s: the interesting CCN article compelled me to look.

2 1

Autobuild: problems with modified amino acids
by Kendall Nettles 19 Jan '11

19 Jan '11

I'm running autobuild from a molecular replacement solution and having two problems: 1. methionines are converted to MSE 2. I'd like to include a .cif file for b-mercapto-ethanol adducted cysteines: CME. The .cif file is one that works for phenix.refine (listed below). I suspect I am not writing the .eff file correctly: refinement.pdb_interpretation.apply_cif_modification { data_mod = cme residue_selection = chain A and resid 530 residue_selection = chain B and resid 381 } Thanks! Kendall Nettles # electronic Ligand Builder and Optimisation Workbench (eLBOW) # - a module of PHENIX version 1.6-289 (Mon Jan 21 12:43:00 2010) # - file written: Sat May 22 10:41:08 2010 # # Input file: /AutoBuild_run_1_/TEMP0/LIGANDS.pdb # Residue: CME # data_comp_list loop_ _chem_comp.id _chem_comp.three_letter_code _chem_comp.name _chem_comp.group _chem_comp.number_atoms_all _chem_comp.number_atoms_nh _chem_comp.desc_level CME CME 'Unknown ' ligand 21 10 . # data_comp_CME # loop_ _chem_comp_atom.comp_id _chem_comp_atom.atom_id _chem_comp_atom.type_symbol _chem_comp_atom.type_energy _chem_comp_atom.partial_charge CME N N NH2 . CME C C C1 . CME O O O . CME CA C CH1 . CME CB C CH2 . CME SG S S2 . CME SD S S2 . CME CE C CH2 . CME CZ C CH2 . CME OH O OH1 . CME H H HNH2 . CME H2 H HNH2 . CME HC1 H H . CME HA H HCH1 . CME HB2 H HCH2 . CME HB3 H HCH2 . CME HE2 H HCH2 . CME HE3 H HCH2 . CME HZ2 H HCH2 . CME HZ3 H HCH2 . CME HH1 H HOH1 . # loop_ _chem_comp_bond.comp_id _chem_comp_bond.atom_id_1 _chem_comp_bond.atom_id_2 _chem_comp_bond.type _chem_comp_bond.value_dist _chem_comp_bond.value_dist_esd CME N CA single 1.450 0.02 CME C O double 1.224 0.02 CME C CA single 1.492 0.02 CME CA CB single 1.497 0.02 CME CB SG single 1.780 0.02 CME SG SD single 2.018 0.02 CME SD CE single 1.783 0.02 CME CE CZ single 1.485 0.02 CME CZ OH single 1.374 0.02 CME H N single 0.913 0.02 CME H2 N single 0.914 0.02 CME HC1 C single 0.965 0.02 CME HA CA single 0.975 0.02 CME HB2 CB single 0.966 0.02 CME HB3 CB single 0.973 0.02 CME HE2 CE single 0.972 0.02 CME HE3 CE single 0.973 0.02 CME HZ2 CZ single 0.976 0.02 CME HZ3 CZ single 0.968 0.02 CME HH1 OH single 0.880 0.02 # loop_ _chem_comp_angle.comp_id _chem_comp_angle.atom_id_1 _chem_comp_angle.atom_id_2 _chem_comp_angle.atom_id_3 _chem_comp_angle.value_angle _chem_comp_angle.value_angle_esd CME C CA N 109.05 3.0 CME CB CA N 110.14 3.0 CME CB CA C 110.73 3.0 CME CA C O 119.98 3.0 CME SG CB CA 117.67 3.0 CME SD SG CB 109.52 3.0 CME CE SD SG 109.45 3.0 CME CZ CE SD 117.64 3.0 CME OH CZ CE 117.63 3.0 CME HA CA N 109.56 3.0 CME HA CA C 109.15 3.0 CME HC1 C O 120.00 3.0 CME H N CA 109.48 3.0 CME H2 N CA 109.47 3.0 CME HC1 C CA 120.00 3.0 CME HB2 CB CA 107.69 3.0 CME HB3 CB CA 107.72 3.0 CME HA CA CB 108.19 3.0 CME HB2 CB SG 107.73 3.0 CME HB3 CB SG 107.72 3.0 CME HE2 CE SD 107.78 3.0 CME HE3 CE SD 107.74 3.0 CME HZ2 CZ CE 107.71 3.0 CME HZ3 CZ CE 107.80 3.0 CME HE2 CE CZ 107.71 3.0 CME HE3 CE CZ 107.70 3.0 CME HH1 OH CZ 109.49 3.0 CME HZ2 CZ OH 107.74 3.0 CME HZ3 CZ OH 107.69 3.0 CME H2 N H 109.47 3.0 CME HB3 CB HB2 107.97 3.0 CME HE3 CE HE2 107.93 3.0 CME HZ3 CZ HZ2 107.94 3.0 # loop_ _chem_comp_tor.comp_id _chem_comp_tor.id _chem_comp_tor.atom_id_1 _chem_comp_tor.atom_id_2 _chem_comp_tor.atom_id_3 _chem_comp_tor.atom_id_4 _chem_comp_tor.value_angle _chem_comp_tor.value_angle_esd _chem_comp_tor.period CME Var_01 O C CA N 133.07 30.0 2 CME Var_02 SG CB CA N 179.49 30.0 3 CME Var_03 SG CB CA C 58.80 30.0 3 CME Var_04 CB CA C O -105.59 30.0 3 CME Var_05 SD SG CB CA 161.56 30.0 3 CME Var_06 CE SD SG CB -172.66 30.0 3 CME Var_07 CZ CE SD SG 177.33 30.0 3 CME Var_08 OH CZ CE SD 45.37 30.0 3 CME Var_09 HC1 C CA N -45.35 30.0 2 CME Var_10 HB2 CB CA N -58.63 30.0 3 CME Var_11 HB3 CB CA N 57.60 30.0 3 CME Var_12 H N CA C -165.67 30.0 3 CME Var_13 H2 N CA C -45.66 30.0 3 CME Var_14 HB2 CB CA C -179.32 30.0 3 CME Var_15 HB3 CB CA C -63.09 30.0 3 CME Var_16 HA CA C O 13.41 30.0 3 CME Var_17 H N CA CB 72.63 30.0 3 CME Var_18 H2 N CA CB -167.36 30.0 3 CME Var_19 HC1 C CA CB 75.99 30.0 2 CME Var_20 HA CA CB SG -60.78 30.0 2 CME Var_21 HE2 CE SD SG 55.43 30.0 3 CME Var_22 HE3 CE SD SG -60.81 30.0 3 CME Var_23 HB2 CB SG SD 39.69 30.0 1 CME Var_24 HB3 CB SG SD -76.55 30.0 1 CME Var_25 HZ2 CZ CE SD 167.25 30.0 3 CME Var_26 HZ3 CZ CE SD -76.51 30.0 3 CME Var_27 HH1 OH CZ CE 173.81 30.0 3 CME Var_28 HE2 CE CZ OH 167.31 30.0 3 CME Var_29 HE3 CE CZ OH -76.51 30.0 3 CME Var_30 HA CA N H -46.26 30.0 3 CME Var_31 HA CA N H2 73.74 30.0 3 CME Var_32 HA CA C HC1 -165.01 30.0 3 CME Var_33 HB2 CB CA HA 61.10 30.0 3 CME Var_34 HB3 CB CA HA 177.33 30.0 3 CME Var_35 HZ2 CZ CE HE2 -70.81 30.0 3 CME Var_36 HZ3 CZ CE HE2 45.43 30.0 3 CME Var_37 HZ2 CZ CE HE3 45.37 30.0 3 CME Var_38 HZ3 CZ CE HE3 161.61 30.0 3 CME Var_39 HH1 OH CZ HZ2 51.95 30.0 3 CME Var_40 HH1 OH CZ HZ3 -64.25 30.0 3 # loop_ _chem_comp_chir.comp_id _chem_comp_chir.id _chem_comp_chir.atom_id_centre _chem_comp_chir.atom_id_1 _chem_comp_chir.atom_id_2 _chem_comp_chir.atom_id_3 _chem_comp_chir.volume_sign CME chir_01 CA N C CB both #

4 5