-- Dan Thomas Major (at Dr. B. Fischer's lab) Bar-Ilan University Ramat-Gan, Israel
--- Begin Message ---I'm trying to perform ALIGN2D of a target sequence with the sequence of a protein with with X-ray structure, but I keep on getting errors. The protein has several subunits, and I want to align with subunit A only.
- To: modeller_usage <>
- Subject: ALIGN2D problem
- From: Dan Thomas Major <>
- Date: Mon, 12 Nov 2001 16:13:19 +0200
- Organization: Bar-Ilan University
The top file is:
READ_MODEL FILE = '1HZX.pdb', MODEL_SEGMENT '@:A' 'X:A', OUTPUT_CONTROL = 1 1 1 1
SEQUENCE_TO_ALI ALIGN_CODES = '1HZX'
READ_ALIGNMENT FILE = 'hp2y1.seq', ALIGN_CODES = ALIGN_CODES 'hp2y1', ;
ADD_SEQUENCE = on
READ_TOPOLOGY FILE = '/u/ch/majord/modeller/modeller6a/modlib/top_heav.lib'
ALIGN2D
WRITE_ALIGNMENT FILE = 'hp2y1-1HZX.ali', ALIGNMENT_FORMAT = 'PIR'
WRITE_ALIGNMENT FILE = 'hp2y1-1HZX.pap', ALIGNMENT_FORMAT = 'PAP', ;
ALIGNMENT_FEATURES = 'indices helix beta'The error reported is:
...
# Code #_Res #_Segm PDB_code Name
-------------------------------------------------------------------------------
1 1HZX 641 3 1HZX undefined
2 hp2y1 373 1 hp2y1
TOP_________> 4 4 READ_TOPOLOGY FILE = '/u/ch/majord/modeller/modeller6a/mo;
dlib/top_heav.lib'
openf5__224_> Open 11 OLD SEQUENTIAL ${MODINSTALL6a}/modlib//as1.sim.mat
rdrrwgh_268_> Number of residue types: 20
dispers_247E> Internal error: 640 641
recover____E> ERROR_STATUS >= STOP_ON_ERROR: 1 1Dynamically allocated memory at finish [B,kB,MB]: 23066567 22525.945 21.998
Starting time : 2001/11/12 13:53:59.313
Closing time : 2001/11/12 14:41:13.406
Total CPU time [seconds] : 0.00Thanks, Dan
-- Dan Thomas Major (at Dr. B. Fischer's lab) Bar-Ilan University Ramat-Gan, Israel
--- End Message ---