Hello Modellers,
I've been trying to run my first modelling, but cannot get any model so far. Instead, the log file tells me:
Read the alignment from file : hpabc.ali Total number of alignment positions: 144 # Code #_Res #_Segm PDB_code Name ------------------------------------------------------------------------------- 1 1G9L 144 1 1G9L PABC 2 wPABC 72 1 wPABC PABC
TOP_________> 121 107 CHECK_ALIGNMENT
check_a_343_> >> BEGINNING OF COMMAND openf5__224_> Open 11 OLD SEQUENTIAL ./1G9L.atm rdpdb___303E> No atoms were read from the specified input PDB file: 1) Possibly because an incorrect/non-existent PDB file is specified. 2) Possibly because the segment is specified incorrectly in the alignment file or by the TOP variable MODEL_SEGMENT. That is, the beginning residue number and/or chain id in MODEL_SEGMENT may not be found in the input PDB file; MODEL_SEGMENT:1 : To find out more, switch on maximal output by 'SET OUTPUT_CONTROL = 1 1 1 1 2' rdabrk__288W> Protein not accepted: 1 rdabrk__290E> Number of residues in the alignment and pdb files are different: 144 0 For alignment entry: 1 recover ____E> ERROR_STATUS >= STOP_ON_ERROR: 1 1
My alignment file is:
>P1;1G9L structureN:1G9L:1 : :144 : :PABC:homo sapiens: : GPLGSAAAATPAVRTVPQYKYAAGVRNPQQHLNAQPQVTMQQPAVHVQGQEPLTASMLAS APPQEQKQMLGERLFPLIQAMHPTLAGKITGMLLEIDNSELLHMLESPESLRSKVDEAVA VLQAHQAKEAAQKAVNSATGVPTV* >P1;wPABC sequence:wPABC:1 : :144 : :PABC:wheat: : --------------------------------------------------EPLTASMLAS APPQEQKQMLGERLFPLIQAMHPTLAGKITGMLLEIDNSELLHMLESPESLRSKVDEAVA VL----------------------*
Which I trust is the right way to do it. The top file has nothing that deviates from the example file. My pdb file was just renamed from 1G9L.pdb to 1G9L.atm, its content was not edited or changed, and it was located in the right directory.
So what am I missing there?
Any suggestion welcomed!
Thanks, Stephane