Re: [phenixbb] refinement of N-glycosylation
Dear Georg,
Maybe you need to use quotation marks? I also have n-glycosylation and I in
my .eff file there is:
}
apply_cif_link {
data_link = "NAG-ASN"
residue_selection_1 = "chain A and resname NAG and resid XXX"
residue_selection_2 = "chain A and resname ASN and resid XXX"
}
Hope this is useful. Please correct me if I am wrong.
Ling
On Wed, Feb 24, 2010 at 1:13 PM, Georg Mlynek
Dear All,
I´m quite new to crystallography and I´m refining a structure with ASN-NAG-NAG-BMA attached and already got help from Blaine who told me how to do it, but I´m making a small mistake somewhere and can´t figure out what I´m doing wrong. Must be just one paranthesis or so, I guess.
In command line just for testing I put:
phenix.refine refine_data.mtz refine_1-coot-1FrenumberedGlu.pdb cif_link.params BMA.cif strategy=individual_sites+individual_adp main.number_of_macro_cycles=2 --overwrite
where I put in cif_link.params as suggested in the mailing list http://www.phenix-online.org/pipermail/phenixbb/2009-November/002837.htmlwit... and without paranthesis at the beginning
}
apply_cif_link {
data_link = NAG-ASN
residue_selection_1 = chain A and resname NAG and resid 500
residue_selection_2 = chain A and resname ASN and resid 297
}
apply_cif_link {
data_link = BETA1-4
residue_selection_1 = chain A and resname NAG and resid 500
residue_selection_2 = chain A and resname NAG and resid 501
}
apply_cif_link {
data_link = BETA1-4
residue_selection_1 = chain A and resname NAG and resid 501
residue_selection_2 = chain A and resname BMA and resid 502
}
apply_cif_link {
data_link = ALPHA1-3
residue_selection_1 = chain A and resname BMA and resid 502
residue_selection_2 = chain A and resname MAN and resid 503
and suggested by Blaine with and without paranthesis at the beginning
refinement.pdb_interpretation.apply_cif_link {
data_link = NAG-ASN
residue_selection_1 = chain A and resname NAG and resid 500
residue_selection_2 = chain A and resname ASN and resid 297
}
refinement.pdb_interpretation.apply_cif_link {
data_link = BETA1-4
residue_selection_1 = chain A and resname NAG and resid 500
residue_selection_2 = chain A and resname NAG and resid 501
}
refinement.pdb_interpretation.apply_cif_link {
data_link = BETA1-4
residue_selection_1 = chain A and resname NAG and resid 501
residue_selection_2 = chain A and resname BMA and resid 502
}
refinement.pdb_interpretation.apply_cif_link {
data_link = ALPHA1-3
residue_selection_1 = chain A and resname BMA and resid 502
residue_selection_2 = chain A and resname MAN and resid 503
I see in the .geo file that the residues are not bonded, and also in the .eff file that these lines are not added there, which should be done I guess, I put the info also direct into the .def file but I get the same problem?
By the way it shouldn´t matter which name I give to this file as long as the extension is correct?
Thank you very much and all the best, Georg.
-----Ursprüngliche Nachricht----- Von: [email protected] [mailto: [email protected]] Im Auftrag von Ursula Schulze-Gahmen Gesendet: Dienstag, 23. Februar 2010 19:20 An: [email protected] Betreff: [phenixbb] weights optimization, add hydrogens
I am refining a 2.0 A structure in Phenix using mostly the graphical
interface. The maps are very good and the R factors are low
(0.208/0.173). But the geometry is not so great with bond deviations of
0.208 and Angle deviations of 1.85. What is the best thing to do to try
to improve the geometry. Should I try to optimize the weights?
The other question is about adding hydrogens to the model. I did add
them in a previous refinement cycle. After rebuilding the model in coot,
most of the model still has the hydrogens in the input file for the next
refinement cycle, but some residues and waters don't have hydrogens.
When I tried to add hydrogens again to the model, I got an error message
about atoms being to close. Does Phenix know where to add hydrogens if
the input model has already hydroegns on many residues?
Thanks
Ursula
--
Ursula Schulze-Gahmen, PhD.
QB3, Tjian Lab
MCB, 16 Barker Hall #3204
University of California Berkeley
Berkeley, CA 94720-3204
Phone: (510) 642 8258
_______________________________________________
phenixbb mailing list
http://phenix-online.org/mailman/listinfo/phenixbb
_______________________________________________ phenixbb mailing list [email protected] http://phenix-online.org/mailman/listinfo/phenixbb
Dear All, 1) thank you very much for the help with the glycosylations. Blaine was right yesterday that I just forgot to set the last brace. As this is also a common question, I thought maybe the manual I wrote now for this issue can be put also in the FAQs .(lots I wrote is just pasted from the manual. 1) check if you take the right sugars. BMA is not same as MAN for example. good papers are: pdb-care (PDB CArbohydrate REsidue check): a program to support annotation of complex carbohydrate structures in PDB files Thomas Lütteke and Claus-W von der Lieth Data mining the protein data bank: automatic detection and assignment of carbohydrate structures Thomas Lütteke, Martin Frank and Claus-W. von der Lieth 2) check which links you need. There are already lots of predefined links in phenix stored in the file mon_lib_list.cif The full path to this file can be obtained with the command: % phenix.where_mon_lib_list_cif All apply_cif_modification and apply_cif_link definitions will be included into the .def files. I.e. it is not necessary to specify the definitions again if further refinement runs are started with .def files. Note that all LINK, SSBOND, HYDBND, SLTBRG and CISPEP records in the input PDB files are ignored. 3) Create a file which you can call whatever you want even the extension doesn't matter. E.g.: cif_link.params Put in this file your links: (Which Atom of residue1 will be connected to which Atom in residue2 is already defined in mon_lib_list.cif) refinement.pdb_interpretation.apply_cif_link { data_link = NAG-ASN residue_selection_1 = chain A and resname NAG and resid 500 residue_selection_2 = chain A and resname ASN and resid 297 } refinement.pdb_interpretation.apply_cif_link { data_link = BETA1-4 residue_selection_1 = chain A and resname NAG and resid 500 residue_selection_2 = chain A and resname NAG and resid 501 } refinement.pdb_interpretation.apply_cif_link { data_link = BETA1-4 residue_selection_1 = chain A and resname NAG and resid 501 residue_selection_2 = chain A and resname BMA and resid 502 } refinement.pdb_interpretation.apply_cif_link { data_link = ALPHA1-3. residue_selection_1 = chain A and resname BMA and resid 502 residue_selection_2 = chain A and resname MAN and resid 503 } 4) run % phenix.refine model.pdb data.hkl cif_link.params 5)Trouble shooting a)run iotbx.phil cif_link.params (this checks if the syntax is right) b) check .eff file the links should be added there. You can also check the . geo file if you see anything strange. c)Having unknown to phenix.refine item in PDB file (novel ligand, etc... (BMA is also not there) ). phenix.refine uses the CCP4 Monomer Library as the source of stereochemical information for building geometry restraints and reporting statistics. If phenix.refine is unable to match an item in input PDB file against the Monomer Library it will stop with "Sorry" message explaining what to do and listing the problem atoms. If this happened, it is necessary to obtain a cif file (parameter file, describing unknown molecule) by either making it manually or having eLBOW program to generate it: phenix.elbow model.pdb --do-all --output=all_ligands this will ask eLBOW to inspect the model_new.pdb file, find all unknown items in it and create one cif file for them all_ligands.cif. Alternatively, one can specify a three-letters name for the unknown residue: phenix.elbow model.pdb --residue=MAN --output=man !Check the file if everything is ok! Once the cif file is created, the new run of phenix.refine will be: phenix.refine model.pdb data.pdb man.cif Consult eLBOW documentation for more details. Check the file 2) From the command line everything is perfect now, but as I also want to know how to do this with the gui (version 1.6-289), I added the files BMA.cif file and the apply_cif_link.params together with my coordinates and the mtz in the Input files. But I get following error message: Screenshot-PHENIX errorcif.png Is this a Bug? Thanks again and kind regards Georg.
participants (2)
-
Georg Mlynek
-
Ling Qin