|
Hi, I have a problem creating a model for a small protein with two zinc atoms. I use the ADD_RESTRAINT commands to create distance restrictions between the zincs and the residues attached to them. Looking at the log file, it looks that the restrictions are read without problems but at the end some ini file is read and I get an Atom index is out of bounds message. I am attaching the top, alignment and log file. I would appreciate any help.
Thanks Alfredo Cardenas
Top file:
# Homology modelling by the MODELLER TOP routine 'model'.
INCLUDE # Include the predefined TOP routines
SET OUTPUT_CONTROL = 1 1 1 1 1 # uncomment to produce a large log file SET ALNFILE = 'alignment1.txt' # alignment filename SET KNOWNS = '1PTQ' # codes of the templates SET SEQUENCE = 'pkci' # code of the target #SET TOPOLOGY_MODEL = 1 SET HETATM_IO = on SET TOPLIB = '$(LIB)/top.lib' SET PARLIB = '$(LIB)/par.lib'
#SET ATOM_FILES_DIRECTORY = './:../atom_files' # directories for input atom files SET STARTING_MODEL= 1 # index of the first model SET ENDING_MODEL = 1 # index of the last model # (determines how many models to calculate) CALL ROUTINE = 'model' # do homology modelling
SUBROUTINE ROUTINE = 'special_restraints' ADD_RESTRAINT ATOM_IDS = 'SG:14' 'ZN:55', ; RESTRAINT_PARAMETERS = 3 1 1 27 2 2 0 3.0 0.1 ADD_RESTRAINT ATOM_IDS = 'SG:17' 'ZN:55', ; RESTRAINT_PARAMETERS = 3 1 1 27 2 2 0 3.0 0.1 ADD_RESTRAINT ATOM_IDS = 'NE2:39' 'ZN:55', ; RESTRAINT_PARAMETERS = 3 1 1 27 2 2 0 3.0 0.1 ADD_RESTRAINT ATOM_IDS = 'SG:42' 'ZN:55', ; RESTRAINT_PARAMETERS = 3 1 1 27 2 2 0 3.0 0.1 ADD_RESTRAINT ATOM_IDS = 'SG:31' 'ZN:54', ; RESTRAINT_PARAMETERS = 3 1 1 27 2 2 0 3.0 0.1 ADD_RESTRAINT ATOM_IDS = 'SG:34' 'ZN:54', ; RESTRAINT_PARAMETERS = 3 1 1 27 2 2 0 3.0 0.1 ADD_RESTRAINT ATOM_IDS = 'SG:50' 'ZN:54', ; RESTRAINT_PARAMETERS = 3 1 1 27 2 2 0 3.0 0.1 ADD_RESTRAINT ATOM_IDS = 'NE2:53' 'ZN:54', ; RESTRAINT_PARAMETERS = 3 1 1 27 2 2 0 3.0 0.1 RETURN END_SUBROUTINE
Align file:
C; alignment of catalytic domain ok PKCiota and PKCtheta >P1;1PTQ structureX:1PTQ:::::::: HRFKVYNYMSPTFCDHCGSLLWGLVKQGLKCEDCGMNVHHKCREKVANLC---/--*
>P1;pkci sequence:pkci:::::::: HTFQAKRFNRRAHCAICTDRIWGLGRQGYKCINCKLLVHKKCHKLVTIECGRH/zz*
Log file:
MODELLER 7v7, Sep 12, 2004 09:15pm
PROTEIN STRUCTURE MODELLING BY SATISFACTION OF SPATIAL RESTRAINTS
Copyright(c) 1989-2004 Andrej Sali All Rights Reserved
Written by A. Sali with help from B. Webb, M.S. Madhusudhan, M-Y. Shen, M.A. Marti-Renom, N. Eswar, F. Alber, B. Oliva, A. Fiser, R. Sanchez, B. Yerkovich, A. Badretdinov, F. Melo, J.P. Overington, E. Feyfant University of California, San Francisco, USA Rockefeller University, New York, USA Harvard University, Cambridge, USA Imperial Cancer Research Fund, London, UK Birkbeck College, University of London, London, UK
Kind, OS, HostName, Kernel, Processor: 4, Win2000 build 2195 Service Pack 4, C000826, uni, x86 Family 15 Model 2 Stepping 7 Date and time of compilation : 09/13/2004 17:17:22 Job starting time (YY/MM/DD HH:MM:SS): 2005/01/28 14:19:34.867
getprog_531W> ROUTINE redefined: special_restraints TOP_________> 108 746 SET ALNFILE = 'alignment1.txt'
TOP_________> 109 747 SET KNOWNS = '1PTQ'
TOP_________> 110 748 SET SEQUENCE = 'pkci'
TOP_________> 111 749 SET HETATM_IO = ON
TOP_________> 112 750 SET TOPLIB = '$(LIB)/top.lib'
TOP_________> 113 751 SET PARLIB = '$(LIB)/par.lib'
TOP_________> 114 752 SET STARTING_MODEL = 1
TOP_________> 115 753 SET ENDING_MODEL = 1
TOP_________> 116 754 CALL ROUTINE = 'model'
TOP_________> 117 419 CALL ROUTINE = 'getnames'
TOP_________> 118 531 STRING_IF STRING_ARGUMENTS = MODEL 'undefined', OPERATION; = 'EQ', THEN = 'STRING_OPERATE OPERATION = CONCATENA; TE, STRING_ARGUMENTS = SEQUENCE .ini, RESULT = MODEL'
TOP_________> 119 532 STRING_IF STRING_ARGUMENTS = CSRFILE 'undefined', OPERATI; ON = 'EQ', THEN = 'STRING_OPERATE OPERATION = CONCATE; NATE, STRING_ARGUMENTS = SEQUENCE .rsr, RESULT = CSRFILE; '
TOP_________> 120 533 STRING_OPERATE OPERATION = 'CONCATENATE', ; STRING_ARGUMENTS = SEQUENCE '.sch', RESULT = SCHFILE
TOP_________> 121 534 STRING_OPERATE OPERATION = 'CONCATENATE', ; STRING_ARGUMENTS = SEQUENCE '.mat', RESULT = MATRIX_FI; LE
TOP_________> 122 535 SET ROOT_NAME = SEQUENCE
TOP_________> 123 536 RETURN
TOP_________> 124 420 CALL ROUTINE = 'homcsr'
TOP_________> 125 112 READ_ALIGNMENT FILE = ALNFILE, ALIGN_CODES = KNOWNS SEQUE; NCE
Dynamically allocated memory at amaxseq [B,kB,MB]: 84067 82.097 0.080 openf5__224_> Open 20 OLD SEQUENTIAL alignment1.txt openf5__224_> Open 13 OLD SEQUENTIAL ./\1PTQ.atm
Dynamically allocated memory at amaxbnd [B,kB,MB]: 1516107 1480.573 1.446 openf5__224_> Open 20 OLD SEQUENTIAL alignment1.txt openf5__224_> Open 13 OLD SEQUENTIAL ./\1PTQ.atm read_al_374_> Non-standard residue type,position,sequence: z 54 2 read_al_374_> Non-standard residue type,position,sequence: z 55 2
Read the alignment from file : alignment1.txt Total number of alignment positions: 55
# Code #_Res #_Segm PDB_code Name ------------------------------------------------------------------------------- 1 1PTQ 50 1 1PTQ 2 pkci 55 2 pkci TOP_________> 126 113 STRING_IF STRING_ARGUMENTS = TOP_VERSION 'accelrys', OPER; ATION = 'ne', THEN = 'GO_TO __ACCELRYS1'
TOP_________> 127 117 CHECK_ALIGNMENT
check_a_343_> >> BEGINNING OF COMMAND openf5__224_> Open 11 OLD SEQUENTIAL ./\1PTQ.atm
check_ali___> Checking the sequence-structure alignment.
Implied target CA(i)-CA(i+1) distances longer than 8.0 angstroms:
ALN_POS TMPL RID1 RID2 NAM1 NAM2 DIST ---------------------------------------------- END OF TABLE check_a_344_> << END OF COMMAND TOP_________> 128 118 CALL ROUTINE = GENERATE_METHOD
TOP_________> 129 83 READ_ALIGNMENT FILE = ALNFILE, ALIGN_CODES = KNOWNS
openf5__224_> Open 20 OLD SEQUENTIAL alignment1.txt openf5__224_> Open 13 OLD SEQUENTIAL ./\1PTQ.atm openf5__224_> Open 20 OLD SEQUENTIAL alignment1.txt openf5__224_> Open 13 OLD SEQUENTIAL ./\1PTQ.atm
Read the alignment from file : alignment1.txt Total number of alignment positions: 50
# Code #_Res #_Segm PDB_code Name ------------------------------------------------------------------------------- 1 1PTQ 50 1 1PTQ TOP_________> 130 84 READ_ALIGNMENT FILE = ALNFILE, ALIGN_CODES = KNOWNS
openf5__224_> Open 20 OLD SEQUENTIAL alignment1.txt openf5__224_> Open 13 OLD SEQUENTIAL ./\1PTQ.atm openf5__224_> Open 20 OLD SEQUENTIAL alignment1.txt openf5__224_> Open 13 OLD SEQUENTIAL ./\1PTQ.atm
Read the alignment from file : alignment1.txt Total number of alignment positions: 50
# Code #_Res #_Segm PDB_code Name ------------------------------------------------------------------------------- 1 1PTQ 50 1 1PTQ TOP_________> 131 85 IF ARGUMENTS = INITIAL_MALIGN3D 0, OPERATION = 'EQ', THEN; = 'GO_TO NO_INITIAL_MALIGN3D'
TOP_________> 132 88 READ_ALIGNMENT FILE = ALNFILE, ALIGN_CODES = KNOWNS SEQUE; NCE
openf5__224_> Open 20 OLD SEQUENTIAL alignment1.txt openf5__224_> Open 13 OLD SEQUENTIAL ./\1PTQ.atm openf5__224_> Open 20 OLD SEQUENTIAL alignment1.txt openf5__224_> Open 13 OLD SEQUENTIAL ./\1PTQ.atm read_al_374_> Non-standard residue type,position,sequence: z 54 2 read_al_374_> Non-standard residue type,position,sequence: z 55 2
Read the alignment from file : alignment1.txt Total number of alignment positions: 55
# Code #_Res #_Segm PDB_code Name ------------------------------------------------------------------------------- 1 1PTQ 50 1 1PTQ 2 pkci 55 2 pkci TOP_________> 133 89 READ_TOPOLOGY FILE = TOPLIB
openf5__224_> Open 11 OLD SEQUENTIAL ${MODINSTALL7v7}/modlib\/top.lib read_to_238_> Reading CHARMM residue topology file version: 22 1 openf5__224_> Open 11 UNKNOWN SEQUENTIAL ${MODINSTALL7v7}/modlib/models.lib TOP_________> 134 90 READ_PARAMETERS FILE = PARLIB
openf5__224_> Open 11 OLD SEQUENTIAL ${MODINSTALL7v7}/modlib\/par.lib
Dynamically allocated memory at amattacns [B,kB,MB]: 1516657 1481.110 1.446 openf5__224_> Open 11 OLD SEQUENTIAL ${MODINSTALL7v7}/modlib\/par.lib rdparf__232_> parameters BONDS ANGLS DIHEDS IMPROPS MRFP MODE 227 561 661 112 0 0 TOP_________> 135 91 CALL ROUTINE = 'create_topology'
TOP_________> 136 106 GENERATE_TOPOLOGY ADD_SEQUENCE = OFF
getf_______W> RTF restraint not found in the atoms list: residue type, indices: 7 53 atom names : C +N atom indices : 884 0
getf_______W> RTF restraint not found in the atoms list: residue type, indices: 7 1 atom names : N -C CA H atom indices : 1 0 3 2
getf_______W> RTF restraint not found in the atoms list: residue type, indices: 7 53 atom names : C CA +N O atom indices : 884 871 0 885
mkilst______> segment topology constructed from sequence and RTF: segments residues atoms bonds angles dihedrals impropers: 1 53 885 894 0 0 284
patch_______> segment topology patched using RTF: 1 ; HIS ; NTER segments residues atoms bonds angles dihedrals impropers: 1 53 887 896 1609 2350 284
patch_______> segment topology patched using RTF: 53 ; HIS ; CTER segments residues atoms bonds angles dihedrals impropers: 1 53 888 897 1611 2353 285
genseg______> segment topology constructed from sequence and RTF: segments residues atoms bonds angles dihedrals impropers: 1 53 888 897 1611 2353 285
mkilst______> segment topology constructed from sequence and RTF: segments residues atoms bonds angles dihedrals impropers: 2 55 890 897 1611 2353 285
genseg______> segment topology constructed from sequence and RTF: segments residues atoms bonds angles dihedrals impropers: 2 55 890 897 1611 2353 285 TOP_________> 137 107 CALL ROUTINE = 'default_patches'
TOP_________> 138 526 READ_ALIGNMENT FILE = ALNFILE, ALIGN_CODES = KNOWNS SEQUE; NCE
openf5__224_> Open 20 OLD SEQUENTIAL alignment1.txt openf5__224_> Open 13 OLD SEQUENTIAL ./\1PTQ.atm openf5__224_> Open 20 OLD SEQUENTIAL alignment1.txt openf5__224_> Open 13 OLD SEQUENTIAL ./\1PTQ.atm read_al_374_> Non-standard residue type,position,sequence: z 54 2 read_al_374_> Non-standard residue type,position,sequence: z 55 2
Read the alignment from file : alignment1.txt Total number of alignment positions: 55
# Code #_Res #_Segm PDB_code Name ------------------------------------------------------------------------------- 1 1PTQ 50 1 1PTQ 2 pkci 55 2 pkci TOP_________> 139 527 PATCH_SS_TEMPLATES
TOP_________> 140 528 RETURN
TOP_________> 141 108 CALL ROUTINE = 'special_patches'
TOP_________> 142 523 RETURN
TOP_________> 143 109 RETURN
TOP_________> 144 92 TRANSFER_XYZ CLUSTER_CUT = -1.0
transfe_506_> MODEL is an average of all templates. transfe_511_> Number of templates for coordinate transfer: 1 After transfering coordinates of the equivalent template atoms, there are defined, undefined atoms in MODEL: 328 562 TOP_________> 145 93 BUILD_MODEL INITIALIZE_XYZ = OFF
build___467W> All coordinates in MODEL have been assigned. build___468W> Some coordinates in MODEL are still undefined. build___464W> Inventing the MODEL coordinates. TOP_________> 146 94 WRITE_MODEL FILE = MODEL
openf5__224_> Open 14 UNKNOWN SEQUENTIAL pkci.ini wrpdb2__568_> Residues, atoms, selected atoms: 55 890 890 TOP_________> 147 95 RETURN
TOP_________> 148 119 IF ARGUMENTS = EXIT_STAGE 2, OPERATION = 'EQ', THEN = 'RE; TURN'
TOP_________> 149 120 IF ARGUMENTS = CREATE_RESTRAINTS 0, OPERATION = 'EQ', THE; N ='GO_TO __SKIP_RSRS'
TOP_________> 150 121 CALL ROUTINE = 'mkhomcsr'
TOP_________> 151 126 MAKE_RESTRAINTS RESTRAINT_TYPE = 'stereo', ADD_RESTRAINTS; = OFF
Dynamically allocated memory at amprmcns [B,kB,MB]: 5191489 5069.813 4.951 make_re_417_> Restraint type to be calculated: stereo r_stere_607W> Cannot find params in params file: CHARMM atoms : NY CA CPT H IUPAC atoms : NE1 CD1 CE2 HE1 Atom indices : 367 365 369 368 Residues : TRP TRP TRP TRP addprm__440W> Adding params (mean,force,period): 1.0268 60.0000 0 For atoms: NY CA CPT H r_stere_606_> Stereochemical restraints were constructed from RTF & PRMF. Added bond,angle,dihedral,improper restraints : 897 1611 1936 285 Total number of restraints before, now : 0 4729 make_re_422_> Number of previous, current restraints : 0 4729 make_re_423_> Number of previous, current selected restraints: 0 4729 TOP_________> 152 127 READ_ALIGNMENT FILE = ALNFILE, ALIGN_CODES = KNOWNS SEQUE; NCE
openf5__224_> Open 20 OLD SEQUENTIAL alignment1.txt openf5__224_> Open 13 OLD SEQUENTIAL ./\1PTQ.atm openf5__224_> Open 20 OLD SEQUENTIAL alignment1.txt openf5__224_> Open 13 OLD SEQUENTIAL ./\1PTQ.atm read_al_374_> Non-standard residue type,position,sequence: z 54 2 read_al_374_> Non-standard residue type,position,sequence: z 55 2
Read the alignment from file : alignment1.txt Total number of alignment positions: 55
# Code #_Res #_Segm PDB_code Name ------------------------------------------------------------------------------- 1 1PTQ 50 1 1PTQ 2 pkci 55 2 pkci TOP_________> 153 128 MAKE_RESTRAINTS RESTRAINT_TYPE = 'phi-psi_binormal', ADD_; RESTRAINTS = ON
Dynamically allocated memory at amprmcns [B,kB,MB]: 5191489 5069.813 4.951 make_re_417_> Restraint type to be calculated: phi-psi_binormal openf5__224_> Open 11 OLD SEQUENTIAL ${MODINSTALL7v7}/modlib/mdt.ini openf5__224_> Open 11 OLD SEQUENTIAL ${MODINSTALL7v7}/modlib/mnch1.bin initmdt_400_> Distance function type: 1 openf5__224_> Open 11 OLD SEQUENTIAL ${MODINSTALL7v7}/modlib/mnch1.mdt irddata_401_> USER SYMMETRY, SYMMETRY: 1 T SYMMETRIC = .T. ==> NALN*NALN/2 SYMMETRIC = .F. ==> NALN*NALN All protein pairs always generated. getdata_643_> Protein accepted: 1PTQ getdata_289_> Proteins (all/accepted): 1 1 make_re_422_> Number of previous, current restraints : 4729 4780 make_re_423_> Number of previous, current selected restraints: 4729 4780 TOP_________> 154 129 SET SPLINE_RANGE = 4.0, SPLINE_DX = 0.3, SPLINE_MIN_POINT; S = 5
TOP_________> 155 130 MAKE_RESTRAINTS RESTRAINT_TYPE = 'omega_dihedral', ADD_RE; STRAINTS = ON
Dynamically allocated memory at amprmcns [B,kB,MB]: 5191489 5069.813 4.951 make_re_417_> Restraint type to be calculated: omega_dihedral openf5__224_> Open 11 OLD SEQUENTIAL ${MODINSTALL7v7}/modlib/mdt.ini openf5__224_> Open 11 OLD SEQUENTIAL ${MODINSTALL7v7}/modlib/omega.bin initmdt_400_> Distance function type: 1 openf5__224_> Open 11 OLD SEQUENTIAL ${MODINSTALL7v7}/modlib/omega.mdt irddata_401_> USER SYMMETRY, SYMMETRY: 1 T SYMMETRIC = .T. ==> NALN*NALN/2 SYMMETRIC = .F. ==> NALN*NALN All protein pairs always generated. getdata_643_> Protein accepted: 1PTQ getdata_289_> Proteins (all/accepted): 1 1 delete__442E> One or more atoms absent from MODEL: O: 53: C: 53: N: 54: CA: 54: delete__442E> One or more atoms absent from MODEL: O: 54: C: 54: N: 55: CA: 55: omgdel__425W> Unselected all O C +N +CA dihedrals: 52 make_re_422_> Number of previous, current restraints : 4780 4832 make_re_423_> Number of previous, current selected restraints: 4780 4780 TOP_________> 156 131 MAKE_RESTRAINTS RESTRAINT_TYPE = 'chi1_dihedral', ADD_RES; TRAINTS = ON
Dynamically allocated memory at amprmcns [B,kB,MB]: 5191489 5069.813 4.951 make_re_417_> Restraint type to be calculated: chi1_dihedral openf5__224_> Open 11 OLD SEQUENTIAL ${MODINSTALL7v7}/modlib/mdt.ini openf5__224_> Open 11 OLD SEQUENTIAL ${MODINSTALL7v7}/modlib/chi1234.bin initmdt_400_> Distance function type: 1 openf5__224_> Open 11 OLD SEQUENTIAL ${MODINSTALL7v7}/modlib/chi1.mdt irddata_401_> USER SYMMETRY, SYMMETRY: 1 T SYMMETRIC = .T. ==> NALN*NALN/2 SYMMETRIC = .F. ==> NALN*NALN All protein pairs always generated. getdata_643_> Protein accepted: 1PTQ getdata_289_> Proteins (all/accepted): 1 1 make_re_422_> Number of previous, current restraints : 4832 4878 make_re_423_> Number of previous, current selected restraints: 4780 4826 TOP_________> 157 132 MAKE_RESTRAINTS RESTRAINT_TYPE = 'chi2_dihedral', ADD_RES; TRAINTS = ON
Dynamically allocated memory at amprmcns [B,kB,MB]: 5191489 5069.813 4.951 make_re_417_> Restraint type to be calculated: chi2_dihedral openf5__224_> Open 11 OLD SEQUENTIAL ${MODINSTALL7v7}/modlib/mdt.ini openf5__224_> Open 11 OLD SEQUENTIAL ${MODINSTALL7v7}/modlib/chi1234.bin initmdt_400_> Distance function type: 1 openf5__224_> Open 11 OLD SEQUENTIAL ${MODINSTALL7v7}/modlib/chi2.mdt irddata_401_> USER SYMMETRY, SYMMETRY: 1 T SYMMETRIC = .T. ==> NALN*NALN/2 SYMMETRIC = .F. ==> NALN*NALN All protein pairs always generated. getdata_643_> Protein accepted: 1PTQ getdata_289_> Proteins (all/accepted): 1 1 make_re_422_> Number of previous, current restraints : 4878 4913 make_re_423_> Number of previous, current selected restraints: 4826 4861 TOP_________> 158 133 MAKE_RESTRAINTS RESTRAINT_TYPE = 'chi3_dihedral', ADD_RES; TRAINTS = ON
Dynamically allocated memory at amprmcns [B,kB,MB]: 5191489 5069.813 4.951 make_re_417_> Restraint type to be calculated: chi3_dihedral openf5__224_> Open 11 OLD SEQUENTIAL ${MODINSTALL7v7}/modlib/mdt.ini openf5__224_> Open 11 OLD SEQUENTIAL ${MODINSTALL7v7}/modlib/chi1234.bin initmdt_400_> Distance function type: 1 openf5__224_> Open 11 OLD SEQUENTIAL ${MODINSTALL7v7}/modlib/chi3.mdt irddata_401_> USER SYMMETRY, SYMMETRY: 1 T SYMMETRIC = .T. ==> NALN*NALN/2 SYMMETRIC = .F. ==> NALN*NALN All protein pairs always generated. getdata_643_> Protein accepted: 1PTQ getdata_289_> Proteins (all/accepted): 1 1 make_re_422_> Number of previous, current restraints : 4913 4928 make_re_423_> Number of previous, current selected restraints: 4861 4876 TOP_________> 159 134 MAKE_RESTRAINTS RESTRAINT_TYPE = 'chi4_dihedral', ADD_RES; TRAINTS = ON
Dynamically allocated memory at amprmcns [B,kB,MB]: 5191489 5069.813 4.951 make_re_417_> Restraint type to be calculated: chi4_dihedral openf5__224_> Open 11 OLD SEQUENTIAL ${MODINSTALL7v7}/modlib/mdt.ini openf5__224_> Open 11 OLD SEQUENTIAL ${MODINSTALL7v7}/modlib/chi1234.bin initmdt_400_> Distance function type: 1 openf5__224_> Open 11 OLD SEQUENTIAL ${MODINSTALL7v7}/modlib/chi4.mdt irddata_401_> USER SYMMETRY, SYMMETRY: 1 T SYMMETRIC = .T. ==> NALN*NALN/2 SYMMETRIC = .F. ==> NALN*NALN All protein pairs always generated. mdtrsr__446W> A potential that relies on one protein is used, yet you have at least one known structure available. MDT, not library, potential is used. getdata_643_> Protein accepted: 1PTQ getdata_289_> Proteins (all/accepted): 1 1 make_re_422_> Number of previous, current restraints : 4928 4940 make_re_423_> Number of previous, current selected restraints: 4876 4888 TOP_________> 160 135 SET SPLINE_RANGE = 4.0, SPLINE_DX = 0.7, SPLINE_MIN_POINT; S = 5
TOP_________> 161 136 SET RES_TYPES = 'STD'
TOP_________> 162 137 SET DISTANCE_RSR_MODEL = 5, MAXIMAL_DISTANCE = MAX_CA-CA_; DISTANCE
TOP_________> 163 138 SET RESIDUE_SPAN_RANGE = 2 99999, RESIDUE_SPAN_SIGN = ON
TOP_________> 164 139 SET RESTRAINT_GROUP = 9
TOP_________> 165 140 PICK_ATOMS PICK_ATOMS_SET = 2, ATOM_TYPES = 'CA'
Number of atoms to choose from, total : 890 890 Atom types to be searched for (ATOM_TYPES) : CA Residue types to be searched for (RES_TYPES) : STD Selection mode (SELECTION_MODE) : ATOM What to do with atoms & sets (SELECTION_STATUS): INITIALIZE SEGMENT search; residue range (2i5,2a5) : 1 55 1: 55: selatm__462_> Number of selected atoms : 53
List of segments of contiguous residues with at least one selected atom:
SEGMENT RESNUM AA --- RESNUM AA LEN 1 1 : HIS --- 53 : HIS 53
TOP_________> 166 141 PICK_ATOMS PICK_ATOMS_SET = 3, ATOM_TYPES = 'CA'
Number of atoms to choose from, total : 890 890 Atom types to be searched for (ATOM_TYPES) : CA Residue types to be searched for (RES_TYPES) : STD Selection mode (SELECTION_MODE) : ATOM What to do with atoms & sets (SELECTION_STATUS): INITIALIZE SEGMENT search; residue range (2i5,2a5) : 1 55 1: 55: selatm__462_> Number of selected atoms : 53
List of segments of contiguous residues with at least one selected atom:
SEGMENT RESNUM AA --- RESNUM AA LEN 1 1 : HIS --- 53 : HIS 53
TOP_________> 167 142 MAKE_RESTRAINTS RESTRAINT_TYPE = 'distance', ADD_RESTRAIN; TS = 'ON'
Dynamically allocated memory at amprmcns [B,kB,MB]: 5191489 5069.813 4.951 make_re_417_> Restraint type to be calculated: distance make_re_422_> Number of previous, current restraints : 4940 5639 make_re_423_> Number of previous, current selected restraints: 4888 5587 TOP_________> 168 143 SET DISTANCE_RSR_MODEL = 6, MAXIMAL_DISTANCE = MAX_N-O_DI; STANCE
TOP_________> 169 144 SET RESIDUE_SPAN_RANGE = 2 99999, RESIDUE_SPAN_SIGN = OFF
TOP_________> 170 145 SET RESTRAINT_GROUP = 10
TOP_________> 171 146 PICK_ATOMS PICK_ATOMS_SET = 2, ATOM_TYPES = 'N'
Number of atoms to choose from, total : 890 890 Atom types to be searched for (ATOM_TYPES) : N Residue types to be searched for (RES_TYPES) : STD Selection mode (SELECTION_MODE) : ATOM What to do with atoms & sets (SELECTION_STATUS): INITIALIZE SEGMENT search; residue range (2i5,2a5) : 1 55 1: 55: selatm__462_> Number of selected atoms : 53
List of segments of contiguous residues with at least one selected atom:
SEGMENT RESNUM AA --- RESNUM AA LEN 1 1 : HIS --- 53 : HIS 53
TOP_________> 172 147 PICK_ATOMS PICK_ATOMS_SET = 3, ATOM_TYPES = 'O'
Number of atoms to choose from, total : 890 890 Atom types to be searched for (ATOM_TYPES) : O Residue types to be searched for (RES_TYPES) : STD Selection mode (SELECTION_MODE) : ATOM What to do with atoms & sets (SELECTION_STATUS): INITIALIZE SEGMENT search; residue range (2i5,2a5) : 1 55 1: 55: selatm__462_> Number of selected atoms : 53
List of segments of contiguous residues with at least one selected atom:
SEGMENT RESNUM AA --- RESNUM AA LEN 1 1 : HIS --- 53 : HIS 53
TOP_________> 173 148 MAKE_RESTRAINTS RESTRAINT_TYPE = 'distance', ADD_RESTRAIN; TS = 'ON'
Dynamically allocated memory at amprmcns [B,kB,MB]: 5191489 5069.813 4.951 make_re_417_> Restraint type to be calculated: distance make_re_422_> Number of previous, current restraints : 5639 6463 make_re_423_> Number of previous, current selected restraints: 5587 6411 TOP_________> 174 149 SET DISTANCE_RSR_MODEL = 6, MAXIMAL_DISTANCE = MAX_SC-MC_; DISTANCE
TOP_________> 175 150 SET RESIDUE_SPAN_RANGE = 1 2, RESIDUE_SPAN_SIGN = OFF
TOP_________> 176 151 SET RESTRAINT_GROUP = 23, RESTRAINT_STDEV = 0.5 1.5
TOP_________> 177 152 PICK_ATOMS PICK_ATOMS_SET = 2, ATOM_TYPES = 'SDCH'
Number of atoms to choose from, total : 890 890 Atom types to be searched for (ATOM_TYPES) : SDCH Residue types to be searched for (RES_TYPES) : STD Selection mode (SELECTION_MODE) : ATOM What to do with atoms & sets (SELECTION_STATUS): INITIALIZE SEGMENT search; residue range (2i5,2a5) : 1 55 1: 55: selatm__462_> Number of selected atoms : 676
List of segments of contiguous residues with at least one selected atom:
SEGMENT RESNUM AA --- RESNUM AA LEN 1 1 : HIS --- 53 : HIS 53
TOP_________> 178 153 PICK_ATOMS PICK_ATOMS_SET = 3, ATOM_TYPES = 'MNCH'
Number of atoms to choose from, total : 890 890 Atom types to be searched for (ATOM_TYPES) : MNCH Residue types to be searched for (RES_TYPES) : STD Selection mode (SELECTION_MODE) : ATOM What to do with atoms & sets (SELECTION_STATUS): INITIALIZE SEGMENT search; residue range (2i5,2a5) : 1 55 1: 55: selatm__462_> Number of selected atoms : 212
List of segments of contiguous residues with at least one selected atom:
SEGMENT RESNUM AA --- RESNUM AA LEN 1 1 : HIS --- 53 : HIS 53
TOP_________> 179 154 MAKE_RESTRAINTS RESTRAINT_TYPE = 'distance', ADD_RESTRAIN; TS = 'ON'
Dynamically allocated memory at amprmcns [B,kB,MB]: 5191489 5069.813 4.951 make_re_417_> Restraint type to be calculated: distance make_re_422_> Number of previous, current restraints : 6463 7008 make_re_423_> Number of previous, current selected restraints: 6411 6956 TOP_________> 180 155 SET DISTANCE_RSR_MODEL = 6, MAXIMAL_DISTANCE = MAX_SC-SC_; DISTANCE
TOP_________> 181 156 SET RESIDUE_SPAN_RANGE = 2 99999, RESIDUE_SPAN_SIGN = ON
TOP_________> 182 157 SET RESTRAINT_GROUP = 26, RESTRAINT_STDEV = 0.5 2.0
TOP_________> 183 158 PICK_ATOMS PICK_ATOMS_SET = 2, ATOM_TYPES = 'SDCH'
Number of atoms to choose from, total : 890 890 Atom types to be searched for (ATOM_TYPES) : SDCH Residue types to be searched for (RES_TYPES) : STD Selection mode (SELECTION_MODE) : ATOM What to do with atoms & sets (SELECTION_STATUS): INITIALIZE SEGMENT search; residue range (2i5,2a5) : 1 55 1: 55: selatm__462_> Number of selected atoms : 676
List of segments of contiguous residues with at least one selected atom:
SEGMENT RESNUM AA --- RESNUM AA LEN 1 1 : HIS --- 53 : HIS 53
TOP_________> 184 159 PICK_ATOMS PICK_ATOMS_SET = 3, ATOM_TYPES = 'SDCH'
Number of atoms to choose from, total : 890 890 Atom types to be searched for (ATOM_TYPES) : SDCH Residue types to be searched for (RES_TYPES) : STD Selection mode (SELECTION_MODE) : ATOM What to do with atoms & sets (SELECTION_STATUS): INITIALIZE SEGMENT search; residue range (2i5,2a5) : 1 55 1: 55: selatm__462_> Number of selected atoms : 676
List of segments of contiguous residues with at least one selected atom:
SEGMENT RESNUM AA --- RESNUM AA LEN 1 1 : HIS --- 53 : HIS 53
TOP_________> 185 160 MAKE_RESTRAINTS RESTRAINT_TYPE = 'distance', ADD_RESTRAIN; TS = 'ON'
Dynamically allocated memory at amprmcns [B,kB,MB]: 5191489 5069.813 4.951 make_re_417_> Restraint type to be calculated: distance make_re_422_> Number of previous, current restraints : 7008 7226 make_re_423_> Number of previous, current selected restraints: 6956 7174 TOP_________> 186 161 CALL ROUTINE = 'hetatm_restraints'
TOP_________> 187 170 SET RESTRAINT_TYPE = 'distance'
TOP_________> 188 171 SET DISTANCE_RSR_MODEL = 7
TOP_________> 189 172 SET MAXIMAL_DISTANCE = 7.0
TOP_________> 190 173 SET ADD_RESTRAINTS = ON
TOP_________> 191 174 SET RESTRAINT_GROUP = 27
TOP_________> 192 175 SET RESTRAINT_STDEV = 0.2 0.0
TOP_________> 193 176 SET RESIDUE_SPAN_RANGE = 0 99999, RESIDUE_SPAN_SIGN = OFF
TOP_________> 194 177 PICK_ATOMS PICK_ATOMS_SET = 2, ATOM_TYPES = 'ALL', RES_TY; PES = 'ALL'
Number of atoms to choose from, total : 890 890 Atom types to be searched for (ATOM_TYPES) : ALL Residue types to be searched for (RES_TYPES) : ALL Selection mode (SELECTION_MODE) : ATOM What to do with atoms & sets (SELECTION_STATUS): INITIALIZE SEGMENT search; residue range (2i5,2a5) : 1 55 1: 55: selatm__462_> Number of selected atoms : 890
List of segments of contiguous residues with at least one selected atom:
SEGMENT RESNUM AA --- RESNUM AA LEN 1 1 : HIS --- 55 : ZN2 55
TOP_________> 195 178 PICK_ATOMS PICK_ATOMS_SET = 3, ATOM_TYPES = 'ALL', RES_TY; PES = 'HET'
Number of atoms to choose from, total : 890 890 Atom types to be searched for (ATOM_TYPES) : ALL Residue types to be searched for (RES_TYPES) : HET Selection mode (SELECTION_MODE) : ATOM What to do with atoms & sets (SELECTION_STATUS): INITIALIZE SEGMENT search; residue range (2i5,2a5) : 1 55 1: 55: selatm__462_> Number of selected atoms : 2
List of segments of contiguous residues with at least one selected atom:
SEGMENT RESNUM AA --- RESNUM AA LEN 1 54 : ZN2 --- 55 : ZN2 2
TOP_________> 196 179 MAKE_RESTRAINTS
Dynamically allocated memory at amprmcns [B,kB,MB]: 5191489 5069.813 4.951 make_re_417_> Restraint type to be calculated: distance make_re_422_> Number of previous, current restraints : 7226 7226 make_re_423_> Number of previous, current selected restraints: 7174 7174 TOP_________> 197 180 RETURN
TOP_________> 198 162 CALL ROUTINE = 'blk_restraints'
TOP_________> 199 183 SET RESTRAINT_TYPE = 'distance'
TOP_________> 200 184 SET DISTANCE_RSR_MODEL = 7
TOP_________> 201 185 SET MAXIMAL_DISTANCE = 10.0
TOP_________> 202 186 SET ADD_RESTRAINTS = ON
TOP_________> 203 187 SET RESTRAINT_GROUP = 27
TOP_________> 204 188 SET RESTRAINT_STDEV = 0.05 0.0
TOP_________> 205 189 SET RESIDUE_SPAN_RANGE = 0 0, RESIDUE_SPAN_SIGN = ON
TOP_________> 206 190 PICK_ATOMS PICK_ATOMS_SET = 2, ATOM_TYPES = 'ALL', RES_T; YPES = 'BLK'
Number of atoms to choose from, total : 890 890 Atom types to be searched for (ATOM_TYPES) : ALL Residue types to be searched for (RES_TYPES) : BLK Selection mode (SELECTION_MODE) : ATOM What to do with atoms & sets (SELECTION_STATUS): INITIALIZE SEGMENT search; residue range (2i5,2a5) : 1 55 1: 55: selatm__462_> Number of selected atoms : 0
List of segments of contiguous residues with at least one selected atom:
SEGMENT RESNUM AA --- RESNUM AA LEN
TOP_________> 207 191 PICK_ATOMS PICK_ATOMS_SET = 3, ATOM_TYPES = 'ALL', RES_T; YPES = 'BLK'
Number of atoms to choose from, total : 890 890 Atom types to be searched for (ATOM_TYPES) : ALL Residue types to be searched for (RES_TYPES) : BLK Selection mode (SELECTION_MODE) : ATOM What to do with atoms & sets (SELECTION_STATUS): INITIALIZE SEGMENT search; residue range (2i5,2a5) : 1 55 1: 55: selatm__462_> Number of selected atoms : 0
List of segments of contiguous residues with at least one selected atom:
SEGMENT RESNUM AA --- RESNUM AA LEN
TOP_________> 208 192 MAKE_RESTRAINTS
Dynamically allocated memory at amprmcns [B,kB,MB]: 5191489 5069.813 4.951 make_re_417_> Restraint type to be calculated: distance make_re_422_> Number of previous, current restraints : 7226 7226 make_re_423_> Number of previous, current selected restraints: 7174 7174 TOP_________> 209 193 SET RESTRAINT_STDEV = 0.2 0.0
TOP_________> 210 194 SET RESIDUE_SPAN_RANGE = 1 99999, RESIDUE_SPAN_SIGN = OFF
TOP_________> 211 195 PICK_ATOMS PICK_ATOMS_SET = 2, ATOM_TYPES = 'CA', RES_TY; PES = 'ALL'
Number of atoms to choose from, total : 890 890 Atom types to be searched for (ATOM_TYPES) : CA Residue types to be searched for (RES_TYPES) : ALL Selection mode (SELECTION_MODE) : ATOM What to do with atoms & sets (SELECTION_STATUS): INITIALIZE SEGMENT search; residue range (2i5,2a5) : 1 55 1: 55: selatm__462_> Number of selected atoms : 53
List of segments of contiguous residues with at least one selected atom:
SEGMENT RESNUM AA --- RESNUM AA LEN 1 1 : HIS --- 53 : HIS 53
TOP_________> 212 196 PICK_ATOMS PICK_ATOMS_SET = 3, ATOM_TYPES = 'ALL', RES_TY; PES = 'BLK'
Number of atoms to choose from, total : 890 890 Atom types to be searched for (ATOM_TYPES) : ALL Residue types to be searched for (RES_TYPES) : BLK Selection mode (SELECTION_MODE) : ATOM What to do with atoms & sets (SELECTION_STATUS): INITIALIZE SEGMENT search; residue range (2i5,2a5) : 1 55 1: 55: selatm__462_> Number of selected atoms : 0
List of segments of contiguous residues with at least one selected atom:
SEGMENT RESNUM AA --- RESNUM AA LEN
TOP_________> 213 197 MAKE_RESTRAINTS
Dynamically allocated memory at amprmcns [B,kB,MB]: 5191489 5069.813 4.951 make_re_417_> Restraint type to be calculated: distance make_re_422_> Number of previous, current restraints : 7226 7226 make_re_423_> Number of previous, current selected restraints: 7174 7174 TOP_________> 214 198 DELETE_ALIGNMENT
TOP_________> 215 199 RETURN
Dynamically allocated memory at amaxseq [B,kB,MB]: 5097041 4977.579 4.861 TOP_________> 216 163 CALL ROUTINE = 'special_restraints'
TOP_________> 217 756 ADD_RESTRAINT ATOM_IDS = 'SG:14' 'ZN:55', RESTRAINT_PARAM; ETERS = 3 1 1 27 2 2 0 3.0 0.1
Dynamically allocated memory at amprmcns [B,kB,MB]: 5097191 4977.726 4.861 TOP_________> 218 757 ADD_RESTRAINT ATOM_IDS = 'SG:17' 'ZN:55', RESTRAINT_PARAM; ETERS = 3 1 1 27 2 2 0 3.0 0.1
Dynamically allocated memory at amprmcns [B,kB,MB]: 5097191 4977.726 4.861 TOP_________> 219 758 ADD_RESTRAINT ATOM_IDS = 'NE2:39' 'ZN:55', RESTRAINT_PARA; METERS = 3 1 1 27 2 2 0 3.0 0.1
Dynamically allocated memory at amprmcns [B,kB,MB]: 5097191 4977.726 4.861 TOP_________> 220 759 ADD_RESTRAINT ATOM_IDS = 'SG:42' 'ZN:55', RESTRAINT_PARAM; ETERS = 3 1 1 27 2 2 0 3.0 0.1
Dynamically allocated memory at amprmcns [B,kB,MB]: 5097191 4977.726 4.861 TOP_________> 221 760 ADD_RESTRAINT ATOM_IDS = 'SG:31' 'ZN:54', RESTRAINT_PARAM; ETERS = 3 1 1 27 2 2 0 3.0 0.1
Dynamically allocated memory at amprmcns [B,kB,MB]: 5097191 4977.726 4.861 TOP_________> 222 761 ADD_RESTRAINT ATOM_IDS = 'SG:34' 'ZN:54', RESTRAINT_PARAM; ETERS = 3 1 1 27 2 2 0 3.0 0.1
Dynamically allocated memory at amprmcns [B,kB,MB]: 5097191 4977.726 4.861 TOP_________> 223 762 ADD_RESTRAINT ATOM_IDS = 'SG:50' 'ZN:54', RESTRAINT_PARAM; ETERS = 3 1 1 27 2 2 0 3.0 0.1
Dynamically allocated memory at amprmcns [B,kB,MB]: 5097191 4977.726 4.861 TOP_________> 224 763 ADD_RESTRAINT ATOM_IDS = 'NE2:53' 'ZN:54', RESTRAINT_PARA; METERS = 3 1 1 27 2 2 0 3.0 0.1
Dynamically allocated memory at amprmcns [B,kB,MB]: 5097191 4977.726 4.861 TOP_________> 225 764 RETURN
TOP_________> 226 164 CONDENSE_RESTRAINTS
delete__443_> Restraints marked for deletion were removed. Total number of restraints before, now: 7234 6845 TOP_________> 227 165 WRITE_RESTRAINTS FILE = CSRFILE
openf5__224_> Open 14 UNKNOWN SEQUENTIAL pkci.rsr TOP_________> 228 166 SET RESIDUE_SPAN_RANGE = -999 -999, RESIDUE_SPAN_SIGN = O; N
top_mod_451W> Requested schedule step not available: 0 SCHEDULE_STEP set to N_SCHEDULE. TOP_________> 229 167 RETURN
TOP_________> 230 122 LABEL __SKIP_RSRS
TOP_________> 231 123 RETURN
TOP_________> 232 421 IF ARGUMENTS = EXIT_STAGE 1, OPERATION = 'GE', THEN = 'RE; TURN'
TOP_________> 233 422 CALL ROUTINE = 'rd_restraints'
TOP_________> 234 276 READ_RESTRAINTS FILE = CSRFILE, ADD_RESTRAINTS = 'off'
openf5__224_> Open 11 OLD SEQUENTIAL pkci.rsr openf5__224_> Open 11 OLD SEQUENTIAL pkci.rsr rdcsr2__307_> Number of restraints read : 0 Number of excluded pairs read: 0 Number of pseudo atoms read : 0
Dynamically allocated memory at amprmcns [B,kB,MB]: 5097191 4977.726 4.861 openf5__224_> Open 11 OLD SEQUENTIAL pkci.rsr openf5__224_> Open 11 OLD SEQUENTIAL pkci.rsr rdcsr2__307_> Number of restraints read : 0 Number of excluded pairs read: 0 Number of pseudo atoms read : 0 rdcsrs__304_> Restraints in memory, selected restraints: 6845 6845 Explicitly excluded atom pairs in memory : 0 Pseudo atoms in memory : 0 TOP_________> 235 277 RETURN
TOP_________> 236 423 CALL ROUTINE = 'multiple_models'
TOP_________> 237 202 DO ID2 = STARTING_MODEL ENDING_MODEL 1
TOP_________> 238 203 SET FINAL_MODEL = 'default'
TOP_________> 239 204 CALL ROUTINE = 'single_model'
TOP_________> 240 215 SET MAX_ITERATIONS = MAX_VAR_ITERATIONS
TOP_________> 241 216 SET ID1 = 0
TOP_________> 242 217 SWITCH_TRACE FILE = 'default', FILE_EXT = '', FILE_ID = '; .D'
openf5__224_> Open 18 UNKNOWN SEQUENTIAL pkci.D00000001 TOP_________> 243 218 READ_MODEL FILE = MODEL
openf5__224_> Open 11 OLD SEQUENTIAL pkci.ini openf5__224_> Open 11 OLD SEQUENTIAL pkci.ini rdatm___297_> Segments, residues, atoms: 2 55 438 rdatm___298_> Segment: 1 1 53 436 rdatm___298_> Segment: 2 54 55 2 TOP_________> 244 219 CALL ROUTINE = 'select_atoms'
TOP_________> 245 288 PICK_ATOMS SELECTION_SEGMENT ='@:@' 'X:X', SELECTION_SEAR; CH ='segment', PICK_ATOMS_SET =1, RES_TYPES =; 'all', ATOM_TYPES ='all', SELECTION_FROM ='a; ll', SELECTION_STATUS ='initialize'
Number of atoms to choose from, total : 438 438 Atom types to be searched for (ATOM_TYPES) : all Residue types to be searched for (RES_TYPES) : all Selection mode (SELECTION_MODE) : ATOM What to do with atoms & sets (SELECTION_STATUS): initialize SEGMENT search; residue range (2i5,2a5) : 1 55 1: 55: selatm__462_> Number of selected atoms : 438
List of segments of contiguous residues with at least one selected atom:
SEGMENT RESNUM AA --- RESNUM AA LEN 1 1 : HIS --- 55 : ZN2 55
TOP_________> 246 289 RETURN
TOP_________> 247 220 CALL ROUTINE = RAND_METHOD
TOP_________> 248 272 RANDOMIZE_XYZ
randomi_498_> Atoms,selected atoms,random_seed,amplitude: 438 438 1 4.0000 randomi_496_> Amplitude is > 0; randomization is done. TOP_________> 249 273 RETURN
TOP_________> 250 221 MAKE_SCHEDULE
openf5__224_> Open 11 OLD SEQUENTIAL ${MODINSTALL7v7}/modlib/sched.lib TOP_________> 251 222 WRITE_SCHEDULE FILE = SCHFILE
openf5__224_> Open 14 UNKNOWN SEQUENTIAL pkci.sch TOP_________> 252 223 IF ARGUMENTS = WRITE_INTERMEDIATES 1, OPERATION = 'EQ', ; THEN = 'WRITE_MODEL FILE = default, FILE_EXT = PDB_E; XT', FILE_ID = '.B'
TOP_________> 253 224 DO IREPEAT = 1 REPEAT_OPTIMIZATION 1
TOP_________> 254 225 CALL ROUTINE = 'single_model_pass'
TOP_________> 255 252 DO SCHEDULE_STEP = 1 N_SCHEDULE 1
TOP_________> 256 253 OPERATE ARGUMENTS = ID1 1, OPERATION = 'SUM', RESULT = ID; 1
TOP_________> 257 254 PICK_RESTRAINTS ADD_RESTRAINTS = OFF
csrrng__299E> Atom index is out of bounds: 441 438 recover____E> ERROR_STATUS >= STOP_ON_ERROR: 1 1
Dynamically allocated memory at finish [B,kB,MB]: 5094825 4975.415 4.859 Starting time : 2005/01/28 14:19:34.867 Closing time : 2005/01/28 14:19:42.305 Total CPU time [seconds] : 5.45
************************************************************ Alfredo Cardenas Assistant Professor Department of Chemistry University of South Florida 4202 E. Fowler Ave, SCA 400 Tampa, FL 33620-5250 Phone: (813) 974-1591 Fax: (813) 974-1733 ************************************************************
|