Jianghi I can answer some of your questions. See below. Jianghai Zhu wrote:
Hi,
Got some questions in using phenix.refine.
1) Can phenix.refine use the ligand library generated by refmac5?
Yes, it should.
2) How do I use elbow to generate the ligand library when I have more than one kind of ligands? When I used --residue=xxx more than once, elbow will only generated library for the last ligand. If I use elbow multiple times to generate multiple libraries, how do I read the libraries in? Should I combine them to one file?
If you have a more recent version of eLBOW you can create a single restraints file for a single PDB file thus elbow.builder protein_and_ligands.pdb --do-all --opt In older versions, you should be able to elbow.builder protein_and_ligands.pdb --residue=LG1,LG2,LG3 --opt to generate the residues you desire. Use elbow.join_cif_files to combine single CIF files into one file. To read more than one CIF file into phenix.refine simply add the file names to the command line phenix.refine filename.pdb filename.mtz ligand1.cif ligand2.cif
3) When I ran phenix.refine, I found the following message in the log.
Number of residues, atoms: 32, 418 Unusual residues: {'NAG': 22, 'MAN': 10} Classifications: {'undetermined': 32} Link IDs: {None: 31} Unresolved non-hydrogen bonds: 32 Unresolved non-hydrogen angles: 64 Unresolved non-hydrogen dihedrals: 22 Unresolved non-hydrogen chiralities: 32 Number of residues, atoms: 18, 28 Unusual residues: {'GOL': 2, ' MG': 2, ' CA': 14} Classifications: {'undetermined': 18} Link IDs: {None: 17}
I already gave the cif library for NAG. Looks like nothing happened. How come phenix.refine doesn't recognize Glycerol, Mg, and Ca?
The most recent version of phenix.refine will recognise the single metal atoms but currently you need to generate the CIF for glycerol with eLBOW. I have attached a CIF file for GOL using the standard atomic naming from the PDB MSDChem database.
4) The refinement (TLS + ML + B individual) went through, I got reasonable R, Rfree, rmsdBOND, rmsdANGLE. But the B factors are pretty low. The B factor of the backbone is much lower than the side chain, some have numbers like 4. Some metal atoms also have B factors around 4. What did I do wrong?
Thanks.
Jianghai _______________________________________________ phenixbb mailing list [email protected] http://www.phenix-online.org/mailman/listinfo/phenixbb
-- Nigel W. Moriarty Building 64R0246B, Physical Biosciences Division Lawrence Berkeley National Laboratory Berkeley, CA 94720-8235 Phone : 510-486-5709 Fax : 510-486-5909 Email : [email protected] Web : CCI.LBL.gov # electronic Ligand Builder and Optimisation Workbench (eLBOW) # - a module of PHENIX version 1.25b (Mon Sep 17 16:26:00 2006) # - file written: Thu Dec 14 08:55:16 2006 # # Quantum optimisation: True # Method: AM1 # SMILES string: OCC(O)CO # Template file: /net/cci/filer1/share/rdb_resources/hicup/g/gol_clean.pdb # 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 LIG LIG 'Unknown ' ligand 14 6 . # data_comp_LIG # 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 _chem_comp_atom.x _chem_comp_atom.y _chem_comp_atom.z LIG O1 O OH1 . 0.0208 0.0709 0.0355 LIG C1 C CH2 . 1.4325 -0.0011 0.0109 LIG C2 C CH1 . 1.9132 0.0138 1.4647 LIG O2 O OH1 . 1.3682 -1.1074 2.1354 LIG C3 C CH2 . 3.4402 -0.1109 1.4474 LIG O3 O OH1 . 3.8600 -0.1930 2.7956 LIG 1H01 H HOH1 . -0.2646 0.0606 -0.8845 LIG 1H02 H HCH2 . 1.8594 0.8904 -0.5225 LIG 2H02 H HCH2 . 1.7739 -0.9421 -0.4960 LIG 1H03 H HCH1 . 1.6018 0.9730 1.9593 LIG 1H04 H HOH1 . 0.9761 -0.7765 2.9506 LIG 1H05 H HCH2 . 3.8935 0.7878 0.9507 LIG 2H05 H HCH2 . 3.7316 -1.0360 0.8822 LIG 1H06 H HOH1 . 4.8200 -0.2706 2.7757 # 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 LIG C1 O1 single 1.41370 0.02 LIG C2 C1 single 1.53131 0.02 LIG O2 C2 single 1.41559 0.02 LIG C3 C2 single 1.53210 0.02 LIG O3 C3 single 1.41441 0.02 LIG 1H01 O1 single 0.96335 0.02 LIG 1H02 C1 single 1.12320 0.02 LIG 2H02 C1 single 1.12196 0.02 LIG 1H03 C2 single 1.12330 0.02 LIG 1H04 O2 single 0.96319 0.02 LIG 1H05 C3 single 1.12248 0.02 LIG 2H05 C3 single 1.12263 0.02 LIG 1H06 O3 single 0.96326 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 LIG C2 C1 O1 107.24269 3.0 LIG 1H02 C1 O1 110.32761 3.0 LIG 2H02 C1 O1 110.75942 3.0 LIG 1H01 O1 C1 106.18066 3.0 LIG O2 C2 C1 108.73907 3.0 LIG C3 C2 C1 107.54137 3.0 LIG 1H03 C2 C1 109.83557 3.0 LIG 1H02 C1 C2 108.89201 3.0 LIG 2H02 C1 C2 109.97070 3.0 LIG 1H04 O2 C2 106.59896 3.0 LIG O3 C3 C2 106.84813 3.0 LIG 1H05 C3 C2 110.01935 3.0 LIG 2H05 C3 C2 109.35677 3.0 LIG C3 C2 O2 108.94178 3.0 LIG 1H03 C2 O2 111.15686 3.0 LIG 1H03 C2 C3 110.53511 3.0 LIG 1H06 O3 C3 106.30838 3.0 LIG 1H05 C3 O3 110.39062 3.0 LIG 2H05 C3 O3 110.79298 3.0 LIG 2H02 C1 1H02 109.60233 3.0 LIG 2H05 C3 1H05 109.39863 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 LIG Var_01 O2 C2 C1 O1 60.33884 10.0 1 LIG Var_02 C3 C2 C1 O1 178.16082 10.0 1 LIG Var_03 O3 C3 C2 C1 -175.86192 10.0 1 LIG Var_04 O3 C3 C2 O2 -58.17137 10.0 1