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

************************************************************