|
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: |