Dear Modellers,
I'm relatively new to Modeller, and I have a question which I hope has an easy answer. Whenever the pdb file for a solved structure has a gap in the residue numbers (so, for instance, it jumps from residue 300 to 315), I get an error saying that the number of residues in the alignment file and the pdb file are different. How can I fix this problem?
Thanks!
Amanda MacFarlane
Amanda MacFarlane Autoimmune Disease Group/Diabetes Ottawa Health Research Institute Lab N1 Box 221 501 Smyth Road Ottawa ON K1S 5B7
Hi,
The error message is unlikely to be due to numbering reasons. Modeller does not keep numbering unless specifically requested, even in that case there is no reason for error message like this. It just reads in sequentially all ATOM records and renumbers them: all resulting models will start with residue 1 or assign the same numbers as in PDB.
You probably have such an alignment file that does not correspond to the ATOM records. These - the residue sequence(s) and the corresponding coordinate entries - must match.
You can test what modeller 'sees' when 'reading' you structure in:
################# read_model FILE = 'MY.pdb' sequence_to_ali atom_files = 'MY.pdb', align_codes = 'MY.pdb' write_alignment file='outputmyseq.pir' #####################
this will output a pir formatted sequence entry for your protein coordinate file. Compare it with the seq. that you currently use in your alignment.
best wishes, Andras Fiser
On Mon, 2003-12-22 at 13:38, MacFarlane, Amanda J. wrote: > Dear Modellers, > > > > I'm relatively new to Modeller, and I have a question which I hope has > an easy answer. Whenever the pdb file for a solved structure has a > gap in the residue numbers (so, for instance, it jumps from residue > 300 to 315), I get an error saying that the number of residues in the > alignment file and the pdb file are different. How can I fix this > problem? > > > > Thanks! > > > > Amanda MacFarlane > > > > Amanda MacFarlane > > Autoimmune Disease Group/Diabetes > > Ottawa Health Research Institute > > Lab N1 Box 221 > > 501 Smyth Road > > Ottawa ON K1S 5B7 > >
Amanda,
A "/" in the alignment means a chain break for MODELLER. So two chain breaks (in both the template and target sequences) and the requisite number of gaps in the template sequence should work.
Eswar.
---
Eswar Narayanan, Ph.D
Mission Bay Genentech Hall 600 16th Street, Suite N474Q University of California, San Francisco San Francisco, CA 94143-2240 (CA 94158 for courier)
Tel +1 (415) 514-4233; Fax +1 (415) 514-4231
http://www.salilab.org/~eashwar
-----Original Message----- From: modeller_usage-bounces@salilab.org [mailto:modeller_usage-bounces@salilab.org] On Behalf Of MacFarlane, Amanda J. Sent: Monday, December 22, 2003 10:38 AM To: 'modeller_usage@salilab.org' Subject: [modeller_usage] Gaps in pdb file
Dear Modellers,
I'm relatively new to Modeller, and I have a question which I hope has an easy answer. Whenever the pdb file for a solved structure has a gap in the residue numbers (so, for instance, it jumps from residue 300 to 315), I get an error saying that the number of residues in the alignment file and the pdb file are different. How can I fix this problem?
Thanks!
Amanda MacFarlane
Amanda MacFarlane
Autoimmune Disease Group/Diabetes
Ottawa Health Research Institute
Lab N1 Box 221
501 Smyth Road
Ottawa ON K1S 5B7
participants (3)
-
Andras Fiser
-
Eswar Narayanan
-
MacFarlane, Amanda J.