============================================================ | | | Crystallography & NMR System (CNS) | | CNSsolve | | | ============================================================ Version: 1.3 Status: General release ============================================================ Written by: A.T.Brunger, P.D.Adams, G.M.Clore, W.L.DeLano, P.Gros, R.W.Grosse-Kunstleve,J.-S.Jiang,J.M.Krahn, J.Kuszewski, M.Nilges, N.S.Pannu, R.J.Read, L.M.Rice, G.F.Schroeder, T.Simonson, G.L.Warren. Copyright (c) 1997-2010 Yale University ============================================================ Running on machine: master3.bw03.cabm.rutgers.edu (x86_64/Linux,64-bit) with 8 threads Program started by: gliu Program started at: 14:40:10 on 29-Apr-2011 ============================================================ FFT3C: Using FFTPACK4.1 CNSsolve>{+ file: generate.inp +} CNSsolve>{+ directory: general +} CNSsolve>{+ description: Generate molecular topology file and/or generate coordinates for unknown atoms +} CNSsolve>{+ comment: CNSsolve> For most cases the generate_easy.inp input file should be used. +} CNSsolve>{+ authors: Paul Adams and Axel Brunger +} CNSsolve>{+ copyright: Yale University +} CNSsolve> CNSsolve>{- Guidelines for using this file: CNSsolve> - all strings must be quoted by double-quotes CNSsolve> - logical variables (true/false) are not quoted CNSsolve> - do not remove any evaluate statements from the file -} CNSsolve> CNSsolve>{- Special patches will have to be entered manually at the relevant points CNSsolve> in the file - see comments throughout the file -} CNSsolve> CNSsolve>evaluate ($par_nonbonded = "_PARAM_") EVALUATE: symbol $PAR_NONBONDED set to "_PARAM_" (string) CNSsolve> CNSsolve>{- begin block parameter definition -} define( DEFINE> DEFINE>{============================== important =================================} DEFINE> DEFINE>{* Coordinates for molecules of the same type (eg. all protein, all DEFINE> nucleic acid etc) can be input in the same coordinate file if the DEFINE> different chains are separated by a TER card or each chain has DEFINE> a different segid or chainid. *} DEFINE> DEFINE>{* A break in a chain can be detected automatically or should be delimited DEFINE> by a BREAK card. In this case no patch (head, tail or link) will be DEFINE> applied between the residues that bound the chain break. *} DEFINE> DEFINE>{* If a segid is present in the coordinate file it will be read unless DEFINE> segid renaming is used below. If renaming is used then all chains in a DEFINE> coordinate file will be given the same segid. *} DEFINE> DEFINE>{* If a PDB chain identifier is present in the coordinate file then this DEFINE> can be used for the segid *} DEFINE> DEFINE>{* NB. All input PDB files must finish with an END statement *} DEFINE> DEFINE>{============ protein topology, linkage, and parameter files =============} DEFINE> DEFINE>{* topology files *} DEFINE>{===>} topology_infile_1="TOPOWAT:protein-allhdg5-4.top"; DEFINE>{===>} topology_infile_2="TOPOWAT:water-allhdg5-4.top"; DEFINE>{===>} topology_infile_3="TOPOWAT:ion.top"; DEFINE> DEFINE>{* linkage files for linear, continuous polymers (protein, DNA, RNA) *} DEFINE>{===>} link_infile_1="TOPOWAT:protein-allhdg5-4.link"; DEFINE>{===>} link_infile_2=""; DEFINE>{===>} link_infile_3=""; DEFINE> DEFINE>{* parameter files *} DEFINE>{===>} parameter_infile_1="TOPOWAT:protein-allhdg5-4.param"; DEFINE>{===>} parameter_infile_2="TOPOWAT:water-allhdg5-4.param"; DEFINE>{===>} parameter_infile_3="TOPOWAT:ion.param"; DEFINE> DEFINE>{====================== other linkages and modifications ==================} DEFINE> DEFINE>{* extra linkages and modifications by custom patches *} DEFINE>{===>} patch_infile=""; DEFINE> DEFINE>{============================ protein files ================================} DEFINE> DEFINE>{* Multiple coordinate files of the same type can be defined by duplicating DEFINE> all of the entries below and incrementing the file number *} DEFINE> DEFINE>{* protein coordinate file *} DEFINE>{===>} prot_coordinate_infile_1="template_or36.pdb"; DEFINE> DEFINE>{* rename segid *} DEFINE>{+ choice: true false +} DEFINE>{===>} prot_rename_1=false; DEFINE> DEFINE>{* new segid *} DEFINE>{===>} prot_segid_1=""; DEFINE> DEFINE>{* convert chainid to segid if chainid is non-blank *} DEFINE>{+ choice: true false +} DEFINE>{===>} prot_convert_1=true; DEFINE> DEFINE>{* separate chains by segid - a new segid starts a new chain *} DEFINE>{+ choice: true false +} DEFINE>{===>} prot_separate_1=true; DEFINE> DEFINE>{========================= nucleic acid files ==============================} DEFINE> DEFINE>{* Multiple coordinate files of the same type can be defined by duplicating DEFINE> all of the entries below and incrementing the file number *} DEFINE> DEFINE>{* nucleic acid coordinate file *} DEFINE>{===>} nucl_coordinate_infile_1=""; DEFINE> DEFINE>{* rename segid *} DEFINE>{+ choice: true false +} DEFINE>{===>} nucl_rename_1=false; DEFINE> DEFINE>{* new segid *} DEFINE>{===>} nucl_segid_1=""; DEFINE> DEFINE>{* convert chainid to segid if chainid is non-blank *} DEFINE>{+ choice: true false +} DEFINE>{===>} nucl_convert_1=true; DEFINE> DEFINE>{* separate chains by segid - a new segid starts a new chain *} DEFINE>{+ choice: true false +} DEFINE>{===>} nucl_separate_1=true; DEFINE> DEFINE>{============================= water files =================================} DEFINE> DEFINE>{* Multiple coordinate files of the same type can be defined by duplicating DEFINE> all of the entries below and incrementing the file number *} DEFINE> DEFINE>{* water coordinate file *} DEFINE>{===>} water_coordinate_infile_1=""; DEFINE> DEFINE>{* rename segid *} DEFINE>{+ choice: true false +} DEFINE>{===>} water_rename_1=false; DEFINE> DEFINE>{* new segid *} DEFINE>{===>} water_segid_1=""; DEFINE> DEFINE>{* convert chainid to segid if chainid is non-blank *} DEFINE>{+ choice: true false +} DEFINE>{===>} water_convert_1=true; DEFINE> DEFINE>{========================= carbohydrate files ==============================} DEFINE> DEFINE>{* Multiple coordinate files of the same type can be defined by duplicating DEFINE> all of the entries below and incrementing the file number *} DEFINE> DEFINE>{* carbohydrate coordinate file *} DEFINE>{===>} carbo_coordinate_infile_1=""; DEFINE> DEFINE>{* rename segid *} DEFINE>{+ choice: true false +} DEFINE>{===>} carbo_rename_1=false; DEFINE> DEFINE>{* new segid *} DEFINE>{===>} carbo_segid_1=""; DEFINE> DEFINE>{* convert chainid to segid if chainid is non-blank *} DEFINE>{+ choice: true false +} DEFINE>{===>} carbo_convert_1=true; DEFINE> DEFINE>{======================== prosthetic group files ===========================} DEFINE> DEFINE>{* Multiple coordinate files of the same type can be defined by duplicating DEFINE> all of the entries below and incrementing the file number *} DEFINE> DEFINE>{* prosthetic group coordinate file *} DEFINE>{===>} prost_coordinate_infile_1=""; DEFINE> DEFINE>{* rename segid *} DEFINE>{+ choice: true false +} DEFINE>{===>} prost_rename_1=false; DEFINE> DEFINE>{* new segid *} DEFINE>{===>} prost_segid_1=""; DEFINE> DEFINE>{* convert chainid to segid if chainid is non-blank *} DEFINE>{+ choice: true false +} DEFINE>{===>} prost_convert_1=true; DEFINE> DEFINE>{============================ ligand files =================================} DEFINE> DEFINE>{* Multiple coordinate files of the same type can be defined by duplicating DEFINE> all of the entries below and incrementing the file number *} DEFINE> DEFINE>{* ligand coordinate file *} DEFINE>{===>} lig_coordinate_infile_1=""; DEFINE> DEFINE>{* rename segid *} DEFINE>{+ choice: true false +} DEFINE>{===>} lig_rename_1=false; DEFINE> DEFINE>{* new segid *} DEFINE>{===>} lig_segid_1=""; DEFINE> DEFINE>{* convert chainid to segid if chainid is non-blank *} DEFINE>{+ choice: true false +} DEFINE>{===>} lig_convert_1=true; DEFINE> DEFINE>{============================== ions files =================================} DEFINE> DEFINE>{* Multiple coordinate files of the same type can be defined by duplicating DEFINE> all of the entries below and incrementing the file number *} DEFINE> DEFINE>{* ion coordinate file *} DEFINE>{===>} ion_coordinate_infile_1=""; DEFINE> DEFINE>{* rename segid *} DEFINE>{+ choice: true false +} DEFINE>{===>} ion_rename_1=false; DEFINE> DEFINE>{* new segid *} DEFINE>{===>} ion_segid_1=""; DEFINE> DEFINE>{* convert chainid to segid if chainid is non-blank *} DEFINE>{+ choice: true false +} DEFINE>{===>} ion_convert_1=true; DEFINE> DEFINE>{============================ renaming atoms ===============================} DEFINE> DEFINE>{* some atoms may need to be renamed in the topology database to conform DEFINE> to what is present in the coordinate file *} DEFINE> DEFINE>{* CNS uses O5 for NDG and A2G carbohydrates. PDB users O. DEFINE> What is it currently called in the coordinate file? *} DEFINE>{* this will not be changed if left blank *} DEFINE>{===>} O5_becomes="O"; DEFINE> DEFINE>{======================= automatic chain breaks ============================} DEFINE> DEFINE>{* automatically detect mainchain breaks in proteins/DNA/RNA based on distance *} DEFINE>{* the link at break points will be removed *} DEFINE>{+ choice: true false +} DEFINE>{===>} auto_break=true; DEFINE> DEFINE>{* cutoff distance in Angstroms for identification of breaks *} DEFINE>{* the default of 2.5A should be reasonable for most cases. If the input DEFINE> structure has bad geometry it may be necessary to increase this distance *} DEFINE>{===>} break_cutoff=2.5; DEFINE> DEFINE>{======================= automatic disulphide bonds ========================} DEFINE> DEFINE>{* automatically detect disulphide bonds based on distance *} DEFINE>{+ choice: true false +} DEFINE>{===>} auto_ss=true; DEFINE> DEFINE>{* cutoff distance in Angstroms for identification of disulphides *} DEFINE>{* the default of 3.0A should be reasonable for most cases. If the input DEFINE> structure has bad geometry it may be necessary to increase this distance *} DEFINE>{* Note: the atom positions for the sulfurs have to be known *} DEFINE>{===>} disulphide_dist=3.0; DEFINE> DEFINE>{========================= manual disulphide bonds =========================} DEFINE> DEFINE>{* Select pairs of cysteine residues that form disulphide bonds *} DEFINE>{* First 2 entries are the segid and resid of the first cysteine (CYS A). *} DEFINE>{* Second 2 entries are the segid and resid of the second cysteine (CYS B). *} DEFINE>{+ table: rows=12 numbered DEFINE> cols=5 "use" "segid CYS A" "resid CYS A" "segid CYS B" "resid CYS B" +} DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} ss_use_1=false; DEFINE>{===>} ss_i_segid_1=""; ss_i_resid_1=0; DEFINE>{===>} ss_j_segid_1=""; ss_j_resid_1=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} ss_use_2=false; DEFINE>{===>} ss_i_segid_2=""; ss_i_resid_2=0; DEFINE>{===>} ss_j_segid_2=""; ss_j_resid_2=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} ss_use_3=false; DEFINE>{===>} ss_i_segid_3=""; ss_i_resid_3=0; DEFINE>{===>} ss_j_segid_3=""; ss_j_resid_3=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} ss_use_4=false; DEFINE>{===>} ss_i_segid_4=""; ss_i_resid_4=0; DEFINE>{===>} ss_j_segid_4=""; ss_j_resid_4=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} ss_use_5=false; DEFINE>{===>} ss_i_segid_5=""; ss_i_resid_5=0; DEFINE>{===>} ss_j_segid_5=""; ss_j_resid_5=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} ss_use_6=false; DEFINE>{===>} ss_i_segid_6=""; ss_i_resid_6=0; DEFINE>{===>} ss_j_segid_6=""; ss_j_resid_6=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} ss_use_7=false; DEFINE>{===>} ss_i_segid_7=""; ss_i_resid_7=0; DEFINE>{===>} ss_j_segid_7=""; ss_j_resid_7=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} ss_use_8=false; DEFINE>{===>} ss_i_segid_8=""; ss_i_resid_8=0; DEFINE>{===>} ss_j_segid_8=""; ss_j_resid_8=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} ss_use_9=false; DEFINE>{===>} ss_i_segid_9=""; ss_i_resid_9=0; DEFINE>{===>} ss_j_segid_9=""; ss_j_resid_9=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} ss_use_10=false; DEFINE>{===>} ss_i_segid_10=""; ss_i_resid_10=0; DEFINE>{===>} ss_j_segid_10=""; ss_j_resid_10=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} ss_use_11=false; DEFINE>{===>} ss_i_segid_11=""; ss_i_resid_11=0; DEFINE>{===>} ss_j_segid_11=""; ss_j_resid_11=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} ss_use_12=false; DEFINE>{===>} ss_i_segid_12=""; ss_i_resid_12=0; DEFINE>{===>} ss_j_segid_12=""; ss_j_resid_12=0; DEFINE> DEFINE> DEFINE>{=========================== carbohydrate links ===========================} DEFINE> DEFINE>{* Select pairs of residues that are linked *} DEFINE>{* First entry is the name of the patch residue. *} DEFINE>{* Second and third entries are the resid and segid for the atoms DEFINE> referenced by "-" in the patch. *} DEFINE>{* Fourth and fifth entries are the resid and segid for the atoms DEFINE> referenced by "+" in the patch *} DEFINE>{+ table: rows=20 numbered DEFINE> cols=6 "use" "patch name" "segid -" "resid -" "segid +" "resid +" +} DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} carbo_use_1=false; DEFINE>{===>} carbo_patch_1="B1N"; DEFINE>{===>} carbo_i_segid_1="BBBB"; carbo_i_resid_1=401; DEFINE>{===>} carbo_j_segid_1="AAAA"; carbo_j_resid_1=56; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} carbo_use_2=false; DEFINE>{===>} carbo_patch_2="B1N"; DEFINE>{===>} carbo_i_segid_2="BBBB"; carbo_i_resid_2=402; DEFINE>{===>} carbo_j_segid_2="AAAA"; carbo_j_resid_2=182; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} carbo_use_3=false; DEFINE>{===>} carbo_patch_3=""; DEFINE>{===>} carbo_i_segid_3=""; carbo_i_resid_3=0; DEFINE>{===>} carbo_j_segid_3=""; carbo_j_resid_3=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} carbo_use_4=false; DEFINE>{===>} carbo_patch_4=""; DEFINE>{===>} carbo_i_segid_4=""; carbo_i_resid_4=0; DEFINE>{===>} carbo_j_segid_4=""; carbo_j_resid_4=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} carbo_use_5=false; DEFINE>{===>} carbo_patch_5=""; DEFINE>{===>} carbo_i_segid_5=""; carbo_i_resid_5=0; DEFINE>{===>} carbo_j_segid_5=""; carbo_j_resid_5=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} carbo_use_6=false; DEFINE>{===>} carbo_patch_6=""; DEFINE>{===>} carbo_i_segid_6=""; carbo_i_resid_6=0; DEFINE>{===>} carbo_j_segid_6=""; carbo_j_resid_6=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} carbo_use_7=false; DEFINE>{===>} carbo_patch_7=""; DEFINE>{===>} carbo_i_segid_7=""; carbo_i_resid_7=0; DEFINE>{===>} carbo_j_segid_7=""; carbo_j_resid_7=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} carbo_use_8=false; DEFINE>{===>} carbo_patch_8=""; DEFINE>{===>} carbo_i_segid_8=""; carbo_i_resid_8=0; DEFINE>{===>} carbo_j_segid_8=""; carbo_j_resid_8=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} carbo_use_9=false; DEFINE>{===>} carbo_patch_9=""; DEFINE>{===>} carbo_i_segid_9=""; carbo_i_resid_9=0; DEFINE>{===>} carbo_j_segid_9=""; carbo_j_resid_9=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} carbo_use_10=false; DEFINE>{===>} carbo_patch_10=""; DEFINE>{===>} carbo_i_segid_10=""; carbo_i_resid_10=0; DEFINE>{===>} carbo_j_segid_10=""; carbo_j_resid_10=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} carbo_use_11=false; DEFINE>{===>} carbo_patch_11=""; DEFINE>{===>} carbo_i_segid_11=""; carbo_i_resid_11=0; DEFINE>{===>} carbo_j_segid_11=""; carbo_j_resid_11=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} carbo_use_12=false; DEFINE>{===>} carbo_patch_12=""; DEFINE>{===>} carbo_i_segid_12=""; carbo_i_resid_12=0; DEFINE>{===>} carbo_j_segid_12=""; carbo_j_resid_12=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} carbo_use_13=false; DEFINE>{===>} carbo_patch_13=""; DEFINE>{===>} carbo_i_segid_13=""; carbo_i_resid_13=0; DEFINE>{===>} carbo_j_segid_13=""; carbo_j_resid_13=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} carbo_use_14=false; DEFINE>{===>} carbo_patch_14=""; DEFINE>{===>} carbo_i_segid_14=""; carbo_i_resid_14=0; DEFINE>{===>} carbo_j_segid_14=""; carbo_j_resid_14=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} carbo_use_15=false; DEFINE>{===>} carbo_patch_15=""; DEFINE>{===>} carbo_i_segid_15=""; carbo_i_resid_15=0; DEFINE>{===>} carbo_j_segid_15=""; carbo_j_resid_15=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} carbo_use_16=false; DEFINE>{===>} carbo_patch_16=""; DEFINE>{===>} carbo_i_segid_16=""; carbo_i_resid_16=0; DEFINE>{===>} carbo_j_segid_16=""; carbo_j_resid_16=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} carbo_use_17=false; DEFINE>{===>} carbo_patch_17=""; DEFINE>{===>} carbo_i_segid_17=""; carbo_i_resid_17=0; DEFINE>{===>} carbo_j_segid_17=""; carbo_j_resid_17=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} carbo_use_18=false; DEFINE>{===>} carbo_patch_18=""; DEFINE>{===>} carbo_i_segid_18=""; carbo_i_resid_18=0; DEFINE>{===>} carbo_j_segid_18=""; carbo_j_resid_18=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} carbo_use_19=false; DEFINE>{===>} carbo_patch_19=""; DEFINE>{===>} carbo_i_segid_19=""; carbo_i_resid_19=0; DEFINE>{===>} carbo_j_segid_19=""; carbo_j_resid_19=0; DEFINE> DEFINE>{+ choice: true false +} DEFINE>{===>} carbo_use_20=false; DEFINE>{===>} carbo_patch_20=""; DEFINE>{===>} carbo_i_segid_20=""; carbo_i_resid_20=0; DEFINE>{===>} carbo_j_segid_20=""; carbo_j_resid_20=0; DEFINE> DEFINE>{========================= generate parameters =============================} DEFINE> DEFINE>{* hydrogen flag - determines whether hydrogens will be output *} DEFINE>{* must be true for NMR, atomic resolution X-ray crystallography DEFINE> or modelling. Set to false for X-ray crystallographic DEFINE> applications in the initial stages of refinement. *} DEFINE>{+ choice: true false +} DEFINE>{===>} hydrogen_flag=true; DEFINE> DEFINE>{* selection of atoms to be deleted *} DEFINE>{* to delete all unknown atoms use: (not known) *} DEFINE>{===>} atom_delete=(none); DEFINE> DEFINE>{* selection of atoms for which coordinates will be (re-)generated *} DEFINE>{* to prevent generation of unknown atom use: (none) DEFINE> and set the "selection of atoms to be deleted" to (not known). *} DEFINE>{===>} atom_build=(not known); DEFINE> DEFINE>{* set bfactor flag *} DEFINE>{+ choice: true false +} DEFINE>{===>} set_bfactor=false; DEFINE> DEFINE>{* set bfactor value *} DEFINE>{===>} bfactor=15.0; DEFINE> DEFINE>{* set occupancy flag *} DEFINE>{+ choice: true false +} DEFINE>{===>} set_occupancy=false; DEFINE> DEFINE>{* set occupancy value *} DEFINE>{===>} occupancy=1.0; DEFINE> DEFINE>{============================= output files ================================} DEFINE> DEFINE>{* output structure file *} DEFINE>{===>} structure_outfile="cnsPDB/sa_cns_6.mtf"; DEFINE> DEFINE>{* output coordinate file *} DEFINE>{===>} coordinate_outfile="cnsPDB/sa_cns_6.pdb"; DEFINE> DEFINE>{===========================================================================} DEFINE>{ things below this line do not need to be changed unless } DEFINE>{ you need to apply patches - at the appropriate places marked } DEFINE>{===========================================================================} DEFINE> DEFINE> ) {- end block parameter definition -} CNSsolve> CNSsolve> checkversion 1.3 Program version= 1.3 File version= 1.3 CNSsolve> CNSsolve> evaluate ($log_level=quiet) Assuming literal string "QUIET" EVALUATE: symbol $LOG_LEVEL set to "QUIET" (string) CNSsolve> CNSsolve> {- read topology files -} CNSsolve> topology RTFRDR> evaluate ($counter=1) EVALUATE: symbol $COUNTER set to 1.00000 (real) RTFRDR> evaluate ($done=false) EVALUATE: symbol $DONE set to FALSE (logical) RTFRDR> while ( $done = false ) loop read NEXTCD: condition evaluated as true RTFRDR> if ( &exist_topology_infile_$counter = true ) then NEXTCD: condition evaluated as true RTFRDR> if ( &BLANK%topology_infile_$counter = false ) then NEXTCD: condition evaluated as true RTFRDR> @@&topology_infile_$counter ASSFIL: file protein-allhdg5-4.top opened. RTFRDR>remarks file toppar/protein-allhdg5-4.top RTFRDR>remarks RTFRDR>remarks file topallhdg.pro version 5.4 date 07-July-10 RTFRDR>remarks for file parallhdg.pro version 5.4 date 07-July-10 or later RTFRDR>remarks Geometric energy function parameters for distance geometry and RTFRDR>remarks simulated annealing. RTFRDR>remarks Author: Michael Nilges, EMBL Heidelberg; Institut Pasteur, Paris RTFRDR>remarks This file contains modifications from M. Williams, UCL London RTFRDR>remarks and multiple modifications for HADDOCK from A. Bonvin, Utrecht University RTFRDR>remarks Last modification 10-July-10 RTFRDR> RTFRDR>set message ? end eval ($old_message=$result) set echo ? end eval ($old_echo=$result) MESSage=NORM EVALUATE: symbol $OLD_MESSAGE set to "NORM" (string) ECHO=TRUE {ON} EVALUATE: symbol $OLD_ECHO set to TRUE (logical) RTFRDR>set message=off echo=off end RTFRDR> RTFRDR> end if RTFRDR> else RTFRDR> evaluate ($done=true) RTFRDR> end if RTFRDR> evaluate ($counter=$counter+1) EVALUATE: symbol $COUNTER set to 2.00000 (real) RTFRDR> end loop read RTFRDR> while ( $done = false ) loop read NEXTCD: condition evaluated as true RTFRDR> if ( &exist_topology_infile_$counter = true ) then NEXTCD: condition evaluated as true RTFRDR> if ( &BLANK%topology_infile_$counter = false ) then NEXTCD: condition evaluated as true RTFRDR> @@&topology_infile_$counter ASSFIL: file water-allhdg5-4.top opened. RTFRDR>remarks file toppar/water-allhdg5-4.top RTFRDR>remarks RTFRDR>remarks TOPH19.SOL RTFRDR>remarks ========== RTFRDR>remarks topology file for solvent molecules RTFRDR>remarks water models available: TIP3P model RTFRDR> RTFRDR>set message ? end eval ($old_message=$result) set echo ? end eval ($old_echo=$result) MESSage=NORM EVALUATE: symbol $OLD_MESSAGE set to "NORM" (string) ECHO=TRUE {ON} EVALUATE: symbol $OLD_ECHO set to TRUE (logical) RTFRDR>set message=off echo=off end RTFRDR> RTFRDR> end if RTFRDR> else RTFRDR> evaluate ($done=true) RTFRDR> end if RTFRDR> evaluate ($counter=$counter+1) EVALUATE: symbol $COUNTER set to 3.00000 (real) RTFRDR> end loop read RTFRDR> while ( $done = false ) loop read NEXTCD: condition evaluated as true RTFRDR> if ( &exist_topology_infile_$counter = true ) then NEXTCD: condition evaluated as true RTFRDR> if ( &BLANK%topology_infile_$counter = false ) then NEXTCD: condition evaluated as true RTFRDR> @@&topology_infile_$counter ASSFIL: file ion.top opened. RTFRDR>remarks file toppar/ion.top RTFRDR>remarks topology and masses for common ions RTFRDR>remarks Dingle atom ion residues are given the name of the element. RTFRDR>remarks By default the atom will be uncharged (eg. the residue MG will RTFRDR>remarks contain the atom called MG with zero charge). RTFRDR>remarks To use the charged species the charge state is appended to RTFRDR>remarks the atom name (eg to use MG2+ the residue name is MG2, and the RTFRDR>remarks atom name is MG+2 and has charge +2.0). RTFRDR>remarks NOTE: not all ionic species are represented RTFRDR>remarks PDA 02/09/99 RTFRDR> RTFRDR>set message ? end eval ($old_message=$result) set echo ? end eval ($old_echo=$result) MESSage=NORM EVALUATE: symbol $OLD_MESSAGE set to "NORM" (string) ECHO=TRUE {ON} EVALUATE: symbol $OLD_ECHO set to TRUE (logical) RTFRDR>set message=off echo=off end RTFRDR> RTFRDR> end if RTFRDR> else RTFRDR> evaluate ($done=true) RTFRDR> end if RTFRDR> evaluate ($counter=$counter+1) EVALUATE: symbol $COUNTER set to 4.00000 (real) RTFRDR> end loop read RTFRDR> while ( $done = false ) loop read NEXTCD: condition evaluated as true RTFRDR> if ( &exist_topology_infile_$counter = true ) then NEXTCD: condition evaluated as false RTFRDR> if ( &BLANK%topology_infile_$counter = false ) then RTFRDR> @@&topology_infile_$counter RTFRDR> end if RTFRDR> else RTFRDR> evaluate ($done=true) EVALUATE: symbol $DONE set to TRUE (logical) RTFRDR> end if RTFRDR> evaluate ($counter=$counter+1) EVALUATE: symbol $COUNTER set to 5.00000 (real) RTFRDR> end loop read RTFRDR> while ( $done = false ) loop read NEXTCD: condition evaluated as false RTFRDR> if ( &exist_topology_infile_$counter = true ) then RTFRDR> if ( &BLANK%topology_infile_$counter = false ) then RTFRDR> @@&topology_infile_$counter RTFRDR> end if RTFRDR> else RTFRDR> evaluate ($done=true) RTFRDR> end if RTFRDR> evaluate ($counter=$counter+1) RTFRDR> end loop read RTFRDR> end CNSsolve> CNSsolve> do (refy=0) (all) SELRPN: 0 atoms have been selected out of 0 CNSsolve> evaluate ($counter=1) EVALUATE: symbol $COUNTER set to 1.00000 (real) CNSsolve> evaluate ($done=false) EVALUATE: symbol $DONE set to FALSE (logical) CNSsolve> while ( $done = false ) loop prot NEXTCD: condition evaluated as true CNSsolve> if ( &exist_prot_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &BLANK%prot_coordinate_infile_$counter = false ) then NEXTCD: condition evaluated as true CNSsolve> do (refx=0) (all) SELRPN: 0 atoms have been selected out of 0 CNSsolve> segment SEGMENT> chain CHAIN> if ( &prot_convert_$counter = true ) then NEXTCD: condition evaluated as true CHAIN> convert=true CHAIN> end if CHAIN> if ( &prot_separate_$counter = true ) then NEXTCD: condition evaluated as true CHAIN> separate=true CHAIN> end if CHAIN> evaluate ($count=1) EVALUATE: symbol $COUNT set to 1.00000 (real) CHAIN> evaluate ($done=false) EVALUATE: symbol $DONE set to FALSE (logical) CHAIN> while ( $done = false ) loop read NEXTCD: condition evaluated as true CHAIN> if ( &exist_link_infile_$count = true ) then NEXTCD: condition evaluated as true CHAIN> if ( &BLANK%link_infile_$count = false ) then NEXTCD: condition evaluated as true CHAIN> @@&link_infile_$count ASSFIL: file protein-allhdg5-4.link opened. CHAIN>remarks file toppar/protein-allhdg5-4.link CHAIN>remarks CHAIN>remarks this is a macro to define standard protein peptide bonds CHAIN>remarks and termini to generate a protein sequence. CHAIN> CHAIN>set message ? end eval ($old_message=$result) set echo ? end eval ($old_echo=$result) MESSage=NORM EVALUATE: symbol $OLD_MESSAGE set to "NORM" (string) ECHO=TRUE {ON} EVALUATE: symbol $OLD_ECHO set to TRUE (logical) CHAIN>set message=off echo=off end CHAIN> CHAIN> end if CHAIN> else CHAIN> evaluate ($done=true) CHAIN> end if CHAIN> evaluate ($count=$count+1) EVALUATE: symbol $COUNT set to 2.00000 (real) CHAIN> end loop read CHAIN> while ( $done = false ) loop read NEXTCD: condition evaluated as true CHAIN> if ( &exist_link_infile_$count = true ) then NEXTCD: condition evaluated as true CHAIN> if ( &BLANK%link_infile_$count = false ) then NEXTCD: condition evaluated as false CHAIN> @@&link_infile_$count CHAIN> end if CHAIN> else CHAIN> evaluate ($done=true) CHAIN> end if CHAIN> evaluate ($count=$count+1) EVALUATE: symbol $COUNT set to 3.00000 (real) CHAIN> end loop read CHAIN> while ( $done = false ) loop read NEXTCD: condition evaluated as true CHAIN> if ( &exist_link_infile_$count = true ) then NEXTCD: condition evaluated as true CHAIN> if ( &BLANK%link_infile_$count = false ) then NEXTCD: condition evaluated as false CHAIN> @@&link_infile_$count CHAIN> end if CHAIN> else CHAIN> evaluate ($done=true) CHAIN> end if CHAIN> evaluate ($count=$count+1) EVALUATE: symbol $COUNT set to 4.00000 (real) CHAIN> end loop read CHAIN> while ( $done = false ) loop read NEXTCD: condition evaluated as true CHAIN> if ( &exist_link_infile_$count = true ) then NEXTCD: condition evaluated as false CHAIN> if ( &BLANK%link_infile_$count = false ) then CHAIN> @@&link_infile_$count CHAIN> end if CHAIN> else CHAIN> evaluate ($done=true) EVALUATE: symbol $DONE set to TRUE (logical) CHAIN> end if CHAIN> evaluate ($count=$count+1) EVALUATE: symbol $COUNT set to 5.00000 (real) CHAIN> end loop read CHAIN> while ( $done = false ) loop read NEXTCD: condition evaluated as false CHAIN> if ( &exist_link_infile_$count = true ) then CHAIN> if ( &BLANK%link_infile_$count = false ) then CHAIN> @@&link_infile_$count CHAIN> end if CHAIN> else CHAIN> evaluate ($done=true) CHAIN> end if CHAIN> evaluate ($count=$count+1) CHAIN> end loop read CHAIN> coordinates @@&prot_coordinate_infile_$counter ASSFIL: file template_or36.pdb opened. COOR>ATOM 1 N MET A 1 2.290 0.393 -1.117 1.00 1.21 N SEGMNT-info: chain termination due to END keyword. 134 residues were inserted into segid "A " CHAIN> end SEGMENT> end Status of internal molecular topology database: -> NATOM= 2282(MAXA= 1000000) NBOND= 2294(MAXB= 1000000) -> NTHETA= 4184(MAXT= 2000000) NGRP= 136(MAXGRP= 1000000) -> NPHI= 3549(MAXP= 2000000) NIMPHI= 1144(MAXIMP= 1000000) -> NNB= 1038(MAXNB= 1000000) CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) EVALUATE: symbol $COUNTER set to 2.00000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop prot CNSsolve> while ( $done = false ) loop prot NEXTCD: condition evaluated as false CNSsolve> if ( &exist_prot_coordinate_infile_$counter = true ) then CNSsolve> if ( &BLANK%prot_coordinate_infile_$counter = false ) then CNSsolve> do (refx=0) (all) CNSsolve> segment CNSsolve> chain CNSsolve> if ( &prot_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> if ( &prot_separate_$counter = true ) then CNSsolve> separate=true CNSsolve> end if CNSsolve> evaluate ($count=1) CNSsolve> evaluate ($done=false) CNSsolve> while ( $done = false ) loop read CNSsolve> if ( &exist_link_infile_$count = true ) then CNSsolve> if ( &BLANK%link_infile_$count = false ) then CNSsolve> @@&link_infile_$count CNSsolve> end if CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> evaluate ($count=$count+1) CNSsolve> end loop read CNSsolve> coordinates @@&prot_coordinate_infile_$counter CNSsolve> end CNSsolve> end CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop prot CNSsolve> CNSsolve> {- execute custom patches -} CNSsolve> if ( &BLANK%patch_infile = false ) then NEXTCD: condition evaluated as false CNSsolve> display CNSsolve> display executing external patch file CNSsolve> display CNSsolve> inline @&patch_infile CNSsolve> end if CNSsolve> CNSsolve> {* any special protein patches can be applied here *} {* doing it my way (RT) *} CNSsolve> evaluate ($HaveCis = "no") {* --- Apply possible CIS peptide patches --- *} EVALUATE: symbol $HAVECIS set to "no" (string) CNSsolve> if ( $HaveCis = "yes" ) then NEXTCD: condition evaluated as false CNSsolve> !CISpep info CNSsolve> end if CNSsolve> CNSsolve> evaluate ($HaveHisd = "no") {* --- Apply possible HISD peptide patches --- *} EVALUATE: symbol $HAVEHISD set to "no" (string) CNSsolve> if ( $HaveHisd = "yes" ) then NEXTCD: condition evaluated as false CNSsolve> !HISDpep info CNSsolve> end if CNSsolve> CNSsolve> evaluate ($HaveHise = "no") {* --- Apply possible HISE peptide patches --- *} EVALUATE: symbol $HAVEHISE set to "no" (string) CNSsolve> if ( $HaveHise = "yes" ) then NEXTCD: condition evaluated as false CNSsolve> !HISEpep info CNSsolve> end if CNSsolve> CNSsolve> evaluate ($HaveCysn = "no") {* --- Apply possible CYNH peptide patches --- *} EVALUATE: symbol $HAVECYSN set to "no" (string) CNSsolve> if ( $HaveCysn = "yes" ) then NEXTCD: condition evaluated as false CNSsolve> !CYSNpep info CNSsolve> end if CNSsolve> CNSsolve> evaluate ($HaveDisu = "no") {* --- Getting ready for S-S bridges --- *} EVALUATE: symbol $HAVEDISU set to "no" (string) CNSsolve> if ( $HaveDisu = "yes" ) then NEXTCD: condition evaluated as false CNSsolve> !SSBridge info CNSsolve> end if CNSsolve> CNSsolve> evaluate ($counter=1) EVALUATE: symbol $COUNTER set to 1.00000 (real) CNSsolve> evaluate ($done=false) EVALUATE: symbol $DONE set to FALSE (logical) CNSsolve> while ( $done = false ) loop prot NEXTCD: condition evaluated as true CNSsolve> if ( &exist_prot_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &BLANK%prot_coordinate_infile_$counter = false ) then NEXTCD: condition evaluated as true CNSsolve> coor COOR> if ( &prot_convert_$counter = true ) then NEXTCD: condition evaluated as true COOR> convert=true COOR> end if COOR> @@&prot_coordinate_infile_$counter ASSFIL: file template_or36.pdb opened. COOR>ATOM 1 N MET A 1 2.290 0.393 -1.117 1.00 1.21 N COOR>ATOM 2 CA MET A 1 3.114 -0.167 -2.182 1.00 32.31 C %READC-ERR: atom A 1 MET HN not found in molecular structure %READC-WRN: still 10 missing coordinates (in selected subset) CNSsolve> set echo=off end SELRPN: 1 atoms have been selected out of 2282 SHOW: sum over selected elements = 1.000000 NEXTCD: condition evaluated as false CNSsolve> if ( &prot_rename_$counter = true ) then NEXTCD: condition evaluated as false CNSsolve> do (segid=capitalize(&prot_segid_$counter)) (attr refy=$counter) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) EVALUATE: symbol $COUNTER set to 2.00000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop prot CNSsolve> while ( $done = false ) loop prot NEXTCD: condition evaluated as true CNSsolve> if ( &exist_prot_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as false CNSsolve> if ( &BLANK%prot_coordinate_infile_$counter = false ) then CNSsolve> coor CNSsolve> if ( &prot_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> @@&prot_coordinate_infile_$counter CNSsolve> set echo=off end CNSsolve> show sum(1) ( not(hydrogen) and not(known) ) CNSsolve> if ( $select = 0 ) then CNSsolve> display %INFO: There are no coordinates missing for non-hydrogen atoms CNSsolve> end if CNSsolve> set echo=on end CNSsolve> if ( &prot_rename_$counter = true ) then CNSsolve> do (segid=capitalize(&prot_segid_$counter)) (attr refy=$counter) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) EVALUATE: symbol $DONE set to TRUE (logical) CNSsolve> end if CNSsolve> end loop prot CNSsolve> while ( $done = false ) loop prot NEXTCD: condition evaluated as false CNSsolve> if ( &exist_prot_coordinate_infile_$counter = true ) then CNSsolve> if ( &BLANK%prot_coordinate_infile_$counter = false ) then CNSsolve> coor CNSsolve> if ( &prot_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> @@&prot_coordinate_infile_$counter CNSsolve> set echo=off end CNSsolve> show sum(1) ( not(hydrogen) and not(known) ) CNSsolve> if ( $select = 0 ) then CNSsolve> display %INFO: There are no coordinates missing for non-hydrogen atoms CNSsolve> end if CNSsolve> set echo=on end CNSsolve> if ( &prot_rename_$counter = true ) then CNSsolve> do (segid=capitalize(&prot_segid_$counter)) (attr refy=$counter) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop prot CNSsolve> CNSsolve> do (refy=0) (all) SELRPN: 2282 atoms have been selected out of 2282 CNSsolve> evaluate ($counter=1) EVALUATE: symbol $COUNTER set to 1.00000 (real) CNSsolve> evaluate ($done=false) EVALUATE: symbol $DONE set to FALSE (logical) CNSsolve> while ( $done = false ) loop nucl NEXTCD: condition evaluated as true CNSsolve> if ( &exist_nucl_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &BLANK%nucl_coordinate_infile_$counter = false ) then NEXTCD: condition evaluated as false CNSsolve> do (refx=0) (all) CNSsolve> segment CNSsolve> chain CNSsolve> if ( &nucl_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> if ( &nucl_separate_$counter = true ) then CNSsolve> separate=true CNSsolve> end if CNSsolve> evaluate ($count=1) CNSsolve> evaluate ($done=false) CNSsolve> while ( $done = false ) loop read CNSsolve> if ( &exist_link_infile_$count = true ) then CNSsolve> if ( &BLANK%link_infile_$count = false ) then CNSsolve> @@&link_infile_$count CNSsolve> end if CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> evaluate ($count=$count+1) CNSsolve> end loop read CNSsolve> coordinates @@&nucl_coordinate_infile_$counter CNSsolve> end CNSsolve> end CNSsolve> if ( &nucl_rename_$counter = true ) then CNSsolve> do (refy=$counter) (attr refx=9999) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) EVALUATE: symbol $COUNTER set to 2.00000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop nucl CNSsolve> while ( $done = false ) loop nucl NEXTCD: condition evaluated as true CNSsolve> if ( &exist_nucl_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as false CNSsolve> if ( &BLANK%nucl_coordinate_infile_$counter = false ) then CNSsolve> do (refx=0) (all) CNSsolve> segment CNSsolve> chain CNSsolve> if ( &nucl_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> if ( &nucl_separate_$counter = true ) then CNSsolve> separate=true CNSsolve> end if CNSsolve> evaluate ($count=1) CNSsolve> evaluate ($done=false) CNSsolve> while ( $done = false ) loop read CNSsolve> if ( &exist_link_infile_$count = true ) then CNSsolve> if ( &BLANK%link_infile_$count = false ) then CNSsolve> @@&link_infile_$count CNSsolve> end if CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> evaluate ($count=$count+1) CNSsolve> end loop read CNSsolve> coordinates @@&nucl_coordinate_infile_$counter CNSsolve> end CNSsolve> end CNSsolve> if ( &nucl_rename_$counter = true ) then CNSsolve> do (refy=$counter) (attr refx=9999) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) EVALUATE: symbol $DONE set to TRUE (logical) CNSsolve> end if CNSsolve> end loop nucl CNSsolve> while ( $done = false ) loop nucl NEXTCD: condition evaluated as false CNSsolve> if ( &exist_nucl_coordinate_infile_$counter = true ) then CNSsolve> if ( &BLANK%nucl_coordinate_infile_$counter = false ) then CNSsolve> do (refx=0) (all) CNSsolve> segment CNSsolve> chain CNSsolve> if ( &nucl_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> if ( &nucl_separate_$counter = true ) then CNSsolve> separate=true CNSsolve> end if CNSsolve> evaluate ($count=1) CNSsolve> evaluate ($done=false) CNSsolve> while ( $done = false ) loop read CNSsolve> if ( &exist_link_infile_$count = true ) then CNSsolve> if ( &BLANK%link_infile_$count = false ) then CNSsolve> @@&link_infile_$count CNSsolve> end if CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> evaluate ($count=$count+1) CNSsolve> end loop read CNSsolve> coordinates @@&nucl_coordinate_infile_$counter CNSsolve> end CNSsolve> end CNSsolve> if ( &nucl_rename_$counter = true ) then CNSsolve> do (refy=$counter) (attr refx=9999) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop nucl CNSsolve> CNSsolve> {* any special nucleic acid patches can be applied here *} CNSsolve> {===>} CNSsolve> {<===} CNSsolve> CNSsolve> evaluate ($counter=1) EVALUATE: symbol $COUNTER set to 1.00000 (real) CNSsolve> evaluate ($done=false) EVALUATE: symbol $DONE set to FALSE (logical) CNSsolve> while ( $done = false ) loop nucl NEXTCD: condition evaluated as true CNSsolve> if ( &exist_nucl_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &BLANK%nucl_coordinate_infile_$counter = false ) then NEXTCD: condition evaluated as false CNSsolve> coor CNSsolve> if ( &nucl_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> @@&nucl_coordinate_infile_$counter CNSsolve> set echo=off end CNSsolve> show sum(1) ( not(hydrogen) and not(known) ) CNSsolve> if ( $select = 0 ) then CNSsolve> display %INFO: There are no coordinates missing for non-hydrogen atoms CNSsolve> end if CNSsolve> set echo=on end CNSsolve> if ( &nucl_rename_$counter = true ) then CNSsolve> do (segid=capitalize(&nucl_segid_$counter)) (attr refy=$counter) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) EVALUATE: symbol $COUNTER set to 2.00000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop nucl CNSsolve> while ( $done = false ) loop nucl NEXTCD: condition evaluated as true CNSsolve> if ( &exist_nucl_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as false CNSsolve> if ( &BLANK%nucl_coordinate_infile_$counter = false ) then CNSsolve> coor CNSsolve> if ( &nucl_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> @@&nucl_coordinate_infile_$counter CNSsolve> set echo=off end CNSsolve> show sum(1) ( not(hydrogen) and not(known) ) CNSsolve> if ( $select = 0 ) then CNSsolve> display %INFO: There are no coordinates missing for non-hydrogen atoms CNSsolve> end if CNSsolve> set echo=on end CNSsolve> if ( &nucl_rename_$counter = true ) then CNSsolve> do (segid=capitalize(&nucl_segid_$counter)) (attr refy=$counter) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) EVALUATE: symbol $DONE set to TRUE (logical) CNSsolve> end if CNSsolve> end loop nucl CNSsolve> while ( $done = false ) loop nucl NEXTCD: condition evaluated as false CNSsolve> if ( &exist_nucl_coordinate_infile_$counter = true ) then CNSsolve> if ( &BLANK%nucl_coordinate_infile_$counter = false ) then CNSsolve> coor CNSsolve> if ( &nucl_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> @@&nucl_coordinate_infile_$counter CNSsolve> set echo=off end CNSsolve> show sum(1) ( not(hydrogen) and not(known) ) CNSsolve> if ( $select = 0 ) then CNSsolve> display %INFO: There are no coordinates missing for non-hydrogen atoms CNSsolve> end if CNSsolve> set echo=on end CNSsolve> if ( &nucl_rename_$counter = true ) then CNSsolve> do (segid=capitalize(&nucl_segid_$counter)) (attr refy=$counter) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop nucl CNSsolve> CNSsolve>!!! {- patching of RNA to DNA -} CNSsolve>!!! evaluate ($counter=0) CNSsolve>!!! for $id in id ( tag and (&dna_sele) ) loop dna CNSsolve>!!! evaluate ($counter=$counter+1) CNSsolve>!!! show (segid) (id $id) CNSsolve>!!! evaluate ($dna.segid.$counter=$result) CNSsolve>!!! show (resid) (id $id) CNSsolve>!!! evaluate ($dna.resid.$counter=$result) CNSsolve>!!! end loop dna CNSsolve>!!! evaluate ($dna.num=$counter) CNSsolve> CNSsolve>!!! evaluate ($counter=0) CNSsolve>!!! while ($counter < $dna.num) loop dnap CNSsolve>!!! evaluate ($counter=$counter+1) CNSsolve>!!! patch deox reference=nil=(segid $dna.segid.$counter and CNSsolve>!!!! resid $dna.resid.$counter) end CNSsolve>!!! end loop dnap CNSsolve> CNSsolve> if ( $log_level = verbose ) then NEXTCD: condition evaluated as false CNSsolve> set message=normal echo=on end CNSsolve> else CNSsolve> set message=off echo=off end CNSsolve> CNSsolve> CNSsolve> do (refy=0) (all) SELRPN: 2282 atoms have been selected out of 2282 CNSsolve> evaluate ($counter=1) EVALUATE: symbol $COUNTER set to 1.00000 (real) CNSsolve> evaluate ($done=false) EVALUATE: symbol $DONE set to FALSE (logical) CNSsolve> while ( $done = false ) loop water NEXTCD: condition evaluated as true CNSsolve> if ( &exist_water_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &BLANK%water_coordinate_infile_$counter = false ) then NEXTCD: condition evaluated as false CNSsolve> do (refx=0) (all) CNSsolve> segment CNSsolve> chain CNSsolve> if ( &water_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> coordinates @@&water_coordinate_infile_$counter CNSsolve> end CNSsolve> end CNSsolve> if ( &water_rename_$counter = true ) then CNSsolve> do (refy=$counter) (attr refx=9999) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) EVALUATE: symbol $COUNTER set to 2.00000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop water CNSsolve> while ( $done = false ) loop water NEXTCD: condition evaluated as true CNSsolve> if ( &exist_water_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as false CNSsolve> if ( &BLANK%water_coordinate_infile_$counter = false ) then CNSsolve> do (refx=0) (all) CNSsolve> segment CNSsolve> chain CNSsolve> if ( &water_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> coordinates @@&water_coordinate_infile_$counter CNSsolve> end CNSsolve> end CNSsolve> if ( &water_rename_$counter = true ) then CNSsolve> do (refy=$counter) (attr refx=9999) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) EVALUATE: symbol $DONE set to TRUE (logical) CNSsolve> end if CNSsolve> end loop water CNSsolve> while ( $done = false ) loop water NEXTCD: condition evaluated as false CNSsolve> if ( &exist_water_coordinate_infile_$counter = true ) then CNSsolve> if ( &BLANK%water_coordinate_infile_$counter = false ) then CNSsolve> do (refx=0) (all) CNSsolve> segment CNSsolve> chain CNSsolve> if ( &water_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> coordinates @@&water_coordinate_infile_$counter CNSsolve> end CNSsolve> end CNSsolve> if ( &water_rename_$counter = true ) then CNSsolve> do (refy=$counter) (attr refx=9999) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop water CNSsolve> CNSsolve> {* any special water patches can be applied here *} CNSsolve> {===>} CNSsolve> {<===} CNSsolve> CNSsolve> evaluate ($counter=1) EVALUATE: symbol $COUNTER set to 1.00000 (real) CNSsolve> evaluate ($done=false) EVALUATE: symbol $DONE set to FALSE (logical) CNSsolve> while ( $done = false ) loop water NEXTCD: condition evaluated as true CNSsolve> if ( &exist_water_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &BLANK%water_coordinate_infile_$counter = false ) then NEXTCD: condition evaluated as false CNSsolve> coor CNSsolve> if ( &water_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> @@&water_coordinate_infile_$counter CNSsolve> set echo=off end CNSsolve> show sum(1) ( not(hydrogen) and not(known) ) CNSsolve> if ( $select = 0 ) then CNSsolve> display %INFO: There are no coordinates missing for non-hydrogen atoms CNSsolve> end if CNSsolve> set echo=on end CNSsolve> if ( &water_rename_$counter = true ) then CNSsolve> do (segid=capitalize(&water_segid_$counter)) (attr refy=$counter) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) EVALUATE: symbol $COUNTER set to 2.00000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop water CNSsolve> while ( $done = false ) loop water NEXTCD: condition evaluated as true CNSsolve> if ( &exist_water_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as false CNSsolve> if ( &BLANK%water_coordinate_infile_$counter = false ) then CNSsolve> coor CNSsolve> if ( &water_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> @@&water_coordinate_infile_$counter CNSsolve> set echo=off end CNSsolve> show sum(1) ( not(hydrogen) and not(known) ) CNSsolve> if ( $select = 0 ) then CNSsolve> display %INFO: There are no coordinates missing for non-hydrogen atoms CNSsolve> end if CNSsolve> set echo=on end CNSsolve> if ( &water_rename_$counter = true ) then CNSsolve> do (segid=capitalize(&water_segid_$counter)) (attr refy=$counter) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) EVALUATE: symbol $DONE set to TRUE (logical) CNSsolve> end if CNSsolve> end loop water CNSsolve> while ( $done = false ) loop water NEXTCD: condition evaluated as false CNSsolve> if ( &exist_water_coordinate_infile_$counter = true ) then CNSsolve> if ( &BLANK%water_coordinate_infile_$counter = false ) then CNSsolve> coor CNSsolve> if ( &water_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> @@&water_coordinate_infile_$counter CNSsolve> set echo=off end CNSsolve> show sum(1) ( not(hydrogen) and not(known) ) CNSsolve> if ( $select = 0 ) then CNSsolve> display %INFO: There are no coordinates missing for non-hydrogen atoms CNSsolve> end if CNSsolve> set echo=on end CNSsolve> if ( &water_rename_$counter = true ) then CNSsolve> do (segid=capitalize(&water_segid_$counter)) (attr refy=$counter) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop water CNSsolve> CNSsolve> do (refy=0) (all) SELRPN: 2282 atoms have been selected out of 2282 CNSsolve> evaluate ($counter=1) EVALUATE: symbol $COUNTER set to 1.00000 (real) CNSsolve> evaluate ($done=false) EVALUATE: symbol $DONE set to FALSE (logical) CNSsolve> while ( $done = false ) loop carbo NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &BLANK%carbo_coordinate_infile_$counter = false ) then NEXTCD: condition evaluated as false CNSsolve> do (refx=0) (all) CNSsolve> segment CNSsolve> chain CNSsolve> if ( &carbo_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> coordinates @@&carbo_coordinate_infile_$counter CNSsolve> end CNSsolve> end CNSsolve> if ( &BLANK%O5_becomes = false ) then CNSsolve> do (name=&O5_becomes) (name O5 and ( resname NDG or resname A2G )) CNSsolve> end if CNSsolve> if ( &carbo_rename_$counter = true ) then CNSsolve> do (refy=$counter) (attr refx=9999) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) EVALUATE: symbol $COUNTER set to 2.00000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop carbo CNSsolve> while ( $done = false ) loop carbo NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as false CNSsolve> if ( &BLANK%carbo_coordinate_infile_$counter = false ) then CNSsolve> do (refx=0) (all) CNSsolve> segment CNSsolve> chain CNSsolve> if ( &carbo_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> coordinates @@&carbo_coordinate_infile_$counter CNSsolve> end CNSsolve> end CNSsolve> if ( &BLANK%O5_becomes = false ) then CNSsolve> do (name=&O5_becomes) (name O5 and ( resname NDG or resname A2G )) CNSsolve> end if CNSsolve> if ( &carbo_rename_$counter = true ) then CNSsolve> do (refy=$counter) (attr refx=9999) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) EVALUATE: symbol $DONE set to TRUE (logical) CNSsolve> end if CNSsolve> end loop carbo CNSsolve> while ( $done = false ) loop carbo NEXTCD: condition evaluated as false CNSsolve> if ( &exist_carbo_coordinate_infile_$counter = true ) then CNSsolve> if ( &BLANK%carbo_coordinate_infile_$counter = false ) then CNSsolve> do (refx=0) (all) CNSsolve> segment CNSsolve> chain CNSsolve> if ( &carbo_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> coordinates @@&carbo_coordinate_infile_$counter CNSsolve> end CNSsolve> end CNSsolve> if ( &BLANK%O5_becomes = false ) then CNSsolve> do (name=&O5_becomes) (name O5 and ( resname NDG or resname A2G )) CNSsolve> end if CNSsolve> if ( &carbo_rename_$counter = true ) then CNSsolve> do (refy=$counter) (attr refx=9999) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop carbo CNSsolve> CNSsolve> evaluate ($counter=1) EVALUATE: symbol $COUNTER set to 1.00000 (real) CNSsolve> evaluate ($done=false) EVALUATE: symbol $DONE set to FALSE (logical) CNSsolve> while ( $done = false ) loop carbo NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &BLANK%carbo_coordinate_infile_$counter = false ) then NEXTCD: condition evaluated as false CNSsolve> coor CNSsolve> if ( &carbo_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> @@&carbo_coordinate_infile_$counter CNSsolve> set echo=off end CNSsolve> show sum(1) ( not(hydrogen) and not(known) ) CNSsolve> if ( $select = 0 ) then CNSsolve> display %INFO: There are no coordinates missing for non-hydrogen atoms CNSsolve> end if CNSsolve> set echo=on end CNSsolve> if ( &carbo_rename_$counter = true ) then CNSsolve> do (segid=capitalize(&carbo_segid_$counter)) (attr refy=$counter) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) EVALUATE: symbol $COUNTER set to 2.00000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop carbo CNSsolve> while ( $done = false ) loop carbo NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as false CNSsolve> if ( &BLANK%carbo_coordinate_infile_$counter = false ) then CNSsolve> coor CNSsolve> if ( &carbo_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> @@&carbo_coordinate_infile_$counter CNSsolve> set echo=off end CNSsolve> show sum(1) ( not(hydrogen) and not(known) ) CNSsolve> if ( $select = 0 ) then CNSsolve> display %INFO: There are no coordinates missing for non-hydrogen atoms CNSsolve> end if CNSsolve> set echo=on end CNSsolve> if ( &carbo_rename_$counter = true ) then CNSsolve> do (segid=capitalize(&carbo_segid_$counter)) (attr refy=$counter) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) EVALUATE: symbol $DONE set to TRUE (logical) CNSsolve> end if CNSsolve> end loop carbo CNSsolve> while ( $done = false ) loop carbo NEXTCD: condition evaluated as false CNSsolve> if ( &exist_carbo_coordinate_infile_$counter = true ) then CNSsolve> if ( &BLANK%carbo_coordinate_infile_$counter = false ) then CNSsolve> coor CNSsolve> if ( &carbo_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> @@&carbo_coordinate_infile_$counter CNSsolve> set echo=off end CNSsolve> show sum(1) ( not(hydrogen) and not(known) ) CNSsolve> if ( $select = 0 ) then CNSsolve> display %INFO: There are no coordinates missing for non-hydrogen atoms CNSsolve> end if CNSsolve> set echo=on end CNSsolve> if ( &carbo_rename_$counter = true ) then CNSsolve> do (segid=capitalize(&carbo_segid_$counter)) (attr refy=$counter) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop carbo CNSsolve> CNSsolve> evaluate ($carc=1) EVALUATE: symbol $CARC set to 1.00000 (real) CNSsolve> evaluate ($done=false) EVALUATE: symbol $DONE set to FALSE (logical) CNSsolve> while ( $done = false ) loop cabr NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_use_$carc = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &carbo_use_$carc = true ) then NEXTCD: condition evaluated as false CNSsolve> evaluate ($segidtmp1=capitalize(&carbo_i_segid_$carc)) CNSsolve> evaluate ($segidtmp2=capitalize(&carbo_j_segid_$carc)) CNSsolve> patch &carbo_patch_$carc CNSsolve> reference=-=(segid $QUOTE%segidtmp1 and CNSsolve> resid &carbo_i_resid_$carc) CNSsolve> reference=+=(segid $QUOTE%segidtmp2 and CNSsolve> resid &carbo_j_resid_$carc) CNSsolve> end CNSsolve> end if CNSsolve> evaluate ($carc=$carc+1) EVALUATE: symbol $CARC set to 2.00000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop cabr CNSsolve> while ( $done = false ) loop cabr NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_use_$carc = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &carbo_use_$carc = true ) then NEXTCD: condition evaluated as false CNSsolve> evaluate ($segidtmp1=capitalize(&carbo_i_segid_$carc)) CNSsolve> evaluate ($segidtmp2=capitalize(&carbo_j_segid_$carc)) CNSsolve> patch &carbo_patch_$carc CNSsolve> reference=-=(segid $QUOTE%segidtmp1 and CNSsolve> resid &carbo_i_resid_$carc) CNSsolve> reference=+=(segid $QUOTE%segidtmp2 and CNSsolve> resid &carbo_j_resid_$carc) CNSsolve> end CNSsolve> end if CNSsolve> evaluate ($carc=$carc+1) EVALUATE: symbol $CARC set to 3.00000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop cabr CNSsolve> while ( $done = false ) loop cabr NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_use_$carc = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &carbo_use_$carc = true ) then NEXTCD: condition evaluated as false CNSsolve> evaluate ($segidtmp1=capitalize(&carbo_i_segid_$carc)) CNSsolve> evaluate ($segidtmp2=capitalize(&carbo_j_segid_$carc)) CNSsolve> patch &carbo_patch_$carc CNSsolve> reference=-=(segid $QUOTE%segidtmp1 and CNSsolve> resid &carbo_i_resid_$carc) CNSsolve> reference=+=(segid $QUOTE%segidtmp2 and CNSsolve> resid &carbo_j_resid_$carc) CNSsolve> end CNSsolve> end if CNSsolve> evaluate ($carc=$carc+1) EVALUATE: symbol $CARC set to 4.00000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop cabr CNSsolve> while ( $done = false ) loop cabr NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_use_$carc = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &carbo_use_$carc = true ) then NEXTCD: condition evaluated as false CNSsolve> evaluate ($segidtmp1=capitalize(&carbo_i_segid_$carc)) CNSsolve> evaluate ($segidtmp2=capitalize(&carbo_j_segid_$carc)) CNSsolve> patch &carbo_patch_$carc CNSsolve> reference=-=(segid $QUOTE%segidtmp1 and CNSsolve> resid &carbo_i_resid_$carc) CNSsolve> reference=+=(segid $QUOTE%segidtmp2 and CNSsolve> resid &carbo_j_resid_$carc) CNSsolve> end CNSsolve> end if CNSsolve> evaluate ($carc=$carc+1) EVALUATE: symbol $CARC set to 5.00000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop cabr CNSsolve> while ( $done = false ) loop cabr NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_use_$carc = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &carbo_use_$carc = true ) then NEXTCD: condition evaluated as false CNSsolve> evaluate ($segidtmp1=capitalize(&carbo_i_segid_$carc)) CNSsolve> evaluate ($segidtmp2=capitalize(&carbo_j_segid_$carc)) CNSsolve> patch &carbo_patch_$carc CNSsolve> reference=-=(segid $QUOTE%segidtmp1 and CNSsolve> resid &carbo_i_resid_$carc) CNSsolve> reference=+=(segid $QUOTE%segidtmp2 and CNSsolve> resid &carbo_j_resid_$carc) CNSsolve> end CNSsolve> end if CNSsolve> evaluate ($carc=$carc+1) EVALUATE: symbol $CARC set to 6.00000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop cabr CNSsolve> while ( $done = false ) loop cabr NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_use_$carc = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &carbo_use_$carc = true ) then NEXTCD: condition evaluated as false CNSsolve> evaluate ($segidtmp1=capitalize(&carbo_i_segid_$carc)) CNSsolve> evaluate ($segidtmp2=capitalize(&carbo_j_segid_$carc)) CNSsolve> patch &carbo_patch_$carc CNSsolve> reference=-=(segid $QUOTE%segidtmp1 and CNSsolve> resid &carbo_i_resid_$carc) CNSsolve> reference=+=(segid $QUOTE%segidtmp2 and CNSsolve> resid &carbo_j_resid_$carc) CNSsolve> end CNSsolve> end if CNSsolve> evaluate ($carc=$carc+1) EVALUATE: symbol $CARC set to 7.00000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop cabr CNSsolve> while ( $done = false ) loop cabr NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_use_$carc = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &carbo_use_$carc = true ) then NEXTCD: condition evaluated as false CNSsolve> evaluate ($segidtmp1=capitalize(&carbo_i_segid_$carc)) CNSsolve> evaluate ($segidtmp2=capitalize(&carbo_j_segid_$carc)) CNSsolve> patch &carbo_patch_$carc CNSsolve> reference=-=(segid $QUOTE%segidtmp1 and CNSsolve> resid &carbo_i_resid_$carc) CNSsolve> reference=+=(segid $QUOTE%segidtmp2 and CNSsolve> resid &carbo_j_resid_$carc) CNSsolve> end CNSsolve> end if CNSsolve> evaluate ($carc=$carc+1) EVALUATE: symbol $CARC set to 8.00000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop cabr CNSsolve> while ( $done = false ) loop cabr NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_use_$carc = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &carbo_use_$carc = true ) then NEXTCD: condition evaluated as false CNSsolve> evaluate ($segidtmp1=capitalize(&carbo_i_segid_$carc)) CNSsolve> evaluate ($segidtmp2=capitalize(&carbo_j_segid_$carc)) CNSsolve> patch &carbo_patch_$carc CNSsolve> reference=-=(segid $QUOTE%segidtmp1 and CNSsolve> resid &carbo_i_resid_$carc) CNSsolve> reference=+=(segid $QUOTE%segidtmp2 and CNSsolve> resid &carbo_j_resid_$carc) CNSsolve> end CNSsolve> end if CNSsolve> evaluate ($carc=$carc+1) EVALUATE: symbol $CARC set to 9.00000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop cabr CNSsolve> while ( $done = false ) loop cabr NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_use_$carc = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &carbo_use_$carc = true ) then NEXTCD: condition evaluated as false CNSsolve> evaluate ($segidtmp1=capitalize(&carbo_i_segid_$carc)) CNSsolve> evaluate ($segidtmp2=capitalize(&carbo_j_segid_$carc)) CNSsolve> patch &carbo_patch_$carc CNSsolve> reference=-=(segid $QUOTE%segidtmp1 and CNSsolve> resid &carbo_i_resid_$carc) CNSsolve> reference=+=(segid $QUOTE%segidtmp2 and CNSsolve> resid &carbo_j_resid_$carc) CNSsolve> end CNSsolve> end if CNSsolve> evaluate ($carc=$carc+1) EVALUATE: symbol $CARC set to 10.0000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop cabr CNSsolve> while ( $done = false ) loop cabr NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_use_$carc = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &carbo_use_$carc = true ) then NEXTCD: condition evaluated as false CNSsolve> evaluate ($segidtmp1=capitalize(&carbo_i_segid_$carc)) CNSsolve> evaluate ($segidtmp2=capitalize(&carbo_j_segid_$carc)) CNSsolve> patch &carbo_patch_$carc CNSsolve> reference=-=(segid $QUOTE%segidtmp1 and CNSsolve> resid &carbo_i_resid_$carc) CNSsolve> reference=+=(segid $QUOTE%segidtmp2 and CNSsolve> resid &carbo_j_resid_$carc) CNSsolve> end CNSsolve> end if CNSsolve> evaluate ($carc=$carc+1) EVALUATE: symbol $CARC set to 11.0000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop cabr CNSsolve> while ( $done = false ) loop cabr NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_use_$carc = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &carbo_use_$carc = true ) then NEXTCD: condition evaluated as false CNSsolve> evaluate ($segidtmp1=capitalize(&carbo_i_segid_$carc)) CNSsolve> evaluate ($segidtmp2=capitalize(&carbo_j_segid_$carc)) CNSsolve> patch &carbo_patch_$carc CNSsolve> reference=-=(segid $QUOTE%segidtmp1 and CNSsolve> resid &carbo_i_resid_$carc) CNSsolve> reference=+=(segid $QUOTE%segidtmp2 and CNSsolve> resid &carbo_j_resid_$carc) CNSsolve> end CNSsolve> end if CNSsolve> evaluate ($carc=$carc+1) EVALUATE: symbol $CARC set to 12.0000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop cabr CNSsolve> while ( $done = false ) loop cabr NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_use_$carc = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &carbo_use_$carc = true ) then NEXTCD: condition evaluated as false CNSsolve> evaluate ($segidtmp1=capitalize(&carbo_i_segid_$carc)) CNSsolve> evaluate ($segidtmp2=capitalize(&carbo_j_segid_$carc)) CNSsolve> patch &carbo_patch_$carc CNSsolve> reference=-=(segid $QUOTE%segidtmp1 and CNSsolve> resid &carbo_i_resid_$carc) CNSsolve> reference=+=(segid $QUOTE%segidtmp2 and CNSsolve> resid &carbo_j_resid_$carc) CNSsolve> end CNSsolve> end if CNSsolve> evaluate ($carc=$carc+1) EVALUATE: symbol $CARC set to 13.0000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop cabr CNSsolve> while ( $done = false ) loop cabr NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_use_$carc = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &carbo_use_$carc = true ) then NEXTCD: condition evaluated as false CNSsolve> evaluate ($segidtmp1=capitalize(&carbo_i_segid_$carc)) CNSsolve> evaluate ($segidtmp2=capitalize(&carbo_j_segid_$carc)) CNSsolve> patch &carbo_patch_$carc CNSsolve> reference=-=(segid $QUOTE%segidtmp1 and CNSsolve> resid &carbo_i_resid_$carc) CNSsolve> reference=+=(segid $QUOTE%segidtmp2 and CNSsolve> resid &carbo_j_resid_$carc) CNSsolve> end CNSsolve> end if CNSsolve> evaluate ($carc=$carc+1) EVALUATE: symbol $CARC set to 14.0000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop cabr CNSsolve> while ( $done = false ) loop cabr NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_use_$carc = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &carbo_use_$carc = true ) then NEXTCD: condition evaluated as false CNSsolve> evaluate ($segidtmp1=capitalize(&carbo_i_segid_$carc)) CNSsolve> evaluate ($segidtmp2=capitalize(&carbo_j_segid_$carc)) CNSsolve> patch &carbo_patch_$carc CNSsolve> reference=-=(segid $QUOTE%segidtmp1 and CNSsolve> resid &carbo_i_resid_$carc) CNSsolve> reference=+=(segid $QUOTE%segidtmp2 and CNSsolve> resid &carbo_j_resid_$carc) CNSsolve> end CNSsolve> end if CNSsolve> evaluate ($carc=$carc+1) EVALUATE: symbol $CARC set to 15.0000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop cabr CNSsolve> while ( $done = false ) loop cabr NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_use_$carc = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &carbo_use_$carc = true ) then NEXTCD: condition evaluated as false CNSsolve> evaluate ($segidtmp1=capitalize(&carbo_i_segid_$carc)) CNSsolve> evaluate ($segidtmp2=capitalize(&carbo_j_segid_$carc)) CNSsolve> patch &carbo_patch_$carc CNSsolve> reference=-=(segid $QUOTE%segidtmp1 and CNSsolve> resid &carbo_i_resid_$carc) CNSsolve> reference=+=(segid $QUOTE%segidtmp2 and CNSsolve> resid &carbo_j_resid_$carc) CNSsolve> end CNSsolve> end if CNSsolve> evaluate ($carc=$carc+1) EVALUATE: symbol $CARC set to 16.0000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop cabr CNSsolve> while ( $done = false ) loop cabr NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_use_$carc = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &carbo_use_$carc = true ) then NEXTCD: condition evaluated as false CNSsolve> evaluate ($segidtmp1=capitalize(&carbo_i_segid_$carc)) CNSsolve> evaluate ($segidtmp2=capitalize(&carbo_j_segid_$carc)) CNSsolve> patch &carbo_patch_$carc CNSsolve> reference=-=(segid $QUOTE%segidtmp1 and CNSsolve> resid &carbo_i_resid_$carc) CNSsolve> reference=+=(segid $QUOTE%segidtmp2 and CNSsolve> resid &carbo_j_resid_$carc) CNSsolve> end CNSsolve> end if CNSsolve> evaluate ($carc=$carc+1) EVALUATE: symbol $CARC set to 17.0000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop cabr CNSsolve> while ( $done = false ) loop cabr NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_use_$carc = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &carbo_use_$carc = true ) then NEXTCD: condition evaluated as false CNSsolve> evaluate ($segidtmp1=capitalize(&carbo_i_segid_$carc)) CNSsolve> evaluate ($segidtmp2=capitalize(&carbo_j_segid_$carc)) CNSsolve> patch &carbo_patch_$carc CNSsolve> reference=-=(segid $QUOTE%segidtmp1 and CNSsolve> resid &carbo_i_resid_$carc) CNSsolve> reference=+=(segid $QUOTE%segidtmp2 and CNSsolve> resid &carbo_j_resid_$carc) CNSsolve> end CNSsolve> end if CNSsolve> evaluate ($carc=$carc+1) EVALUATE: symbol $CARC set to 18.0000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop cabr CNSsolve> while ( $done = false ) loop cabr NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_use_$carc = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &carbo_use_$carc = true ) then NEXTCD: condition evaluated as false CNSsolve> evaluate ($segidtmp1=capitalize(&carbo_i_segid_$carc)) CNSsolve> evaluate ($segidtmp2=capitalize(&carbo_j_segid_$carc)) CNSsolve> patch &carbo_patch_$carc CNSsolve> reference=-=(segid $QUOTE%segidtmp1 and CNSsolve> resid &carbo_i_resid_$carc) CNSsolve> reference=+=(segid $QUOTE%segidtmp2 and CNSsolve> resid &carbo_j_resid_$carc) CNSsolve> end CNSsolve> end if CNSsolve> evaluate ($carc=$carc+1) EVALUATE: symbol $CARC set to 19.0000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop cabr CNSsolve> while ( $done = false ) loop cabr NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_use_$carc = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &carbo_use_$carc = true ) then NEXTCD: condition evaluated as false CNSsolve> evaluate ($segidtmp1=capitalize(&carbo_i_segid_$carc)) CNSsolve> evaluate ($segidtmp2=capitalize(&carbo_j_segid_$carc)) CNSsolve> patch &carbo_patch_$carc CNSsolve> reference=-=(segid $QUOTE%segidtmp1 and CNSsolve> resid &carbo_i_resid_$carc) CNSsolve> reference=+=(segid $QUOTE%segidtmp2 and CNSsolve> resid &carbo_j_resid_$carc) CNSsolve> end CNSsolve> end if CNSsolve> evaluate ($carc=$carc+1) EVALUATE: symbol $CARC set to 20.0000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop cabr CNSsolve> while ( $done = false ) loop cabr NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_use_$carc = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &carbo_use_$carc = true ) then NEXTCD: condition evaluated as false CNSsolve> evaluate ($segidtmp1=capitalize(&carbo_i_segid_$carc)) CNSsolve> evaluate ($segidtmp2=capitalize(&carbo_j_segid_$carc)) CNSsolve> patch &carbo_patch_$carc CNSsolve> reference=-=(segid $QUOTE%segidtmp1 and CNSsolve> resid &carbo_i_resid_$carc) CNSsolve> reference=+=(segid $QUOTE%segidtmp2 and CNSsolve> resid &carbo_j_resid_$carc) CNSsolve> end CNSsolve> end if CNSsolve> evaluate ($carc=$carc+1) EVALUATE: symbol $CARC set to 21.0000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop cabr CNSsolve> while ( $done = false ) loop cabr NEXTCD: condition evaluated as true CNSsolve> if ( &exist_carbo_use_$carc = true ) then NEXTCD: condition evaluated as false CNSsolve> if ( &carbo_use_$carc = true ) then CNSsolve> evaluate ($segidtmp1=capitalize(&carbo_i_segid_$carc)) CNSsolve> evaluate ($segidtmp2=capitalize(&carbo_j_segid_$carc)) CNSsolve> patch &carbo_patch_$carc CNSsolve> reference=-=(segid $QUOTE%segidtmp1 and CNSsolve> resid &carbo_i_resid_$carc) CNSsolve> reference=+=(segid $QUOTE%segidtmp2 and CNSsolve> resid &carbo_j_resid_$carc) CNSsolve> end CNSsolve> end if CNSsolve> evaluate ($carc=$carc+1) CNSsolve> else CNSsolve> evaluate ($done=true) EVALUATE: symbol $DONE set to TRUE (logical) CNSsolve> end if CNSsolve> end loop cabr CNSsolve> while ( $done = false ) loop cabr NEXTCD: condition evaluated as false CNSsolve> if ( &exist_carbo_use_$carc = true ) then CNSsolve> if ( &carbo_use_$carc = true ) then CNSsolve> evaluate ($segidtmp1=capitalize(&carbo_i_segid_$carc)) CNSsolve> evaluate ($segidtmp2=capitalize(&carbo_j_segid_$carc)) CNSsolve> patch &carbo_patch_$carc CNSsolve> reference=-=(segid $QUOTE%segidtmp1 and CNSsolve> resid &carbo_i_resid_$carc) CNSsolve> reference=+=(segid $QUOTE%segidtmp2 and CNSsolve> resid &carbo_j_resid_$carc) CNSsolve> end CNSsolve> end if CNSsolve> evaluate ($carc=$carc+1) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop cabr CNSsolve> CNSsolve> {* any special carbohydrate patches can be applied here *} CNSsolve> {===>} CNSsolve> {<===} CNSsolve> CNSsolve> do (refy=0) (all) SELRPN: 2282 atoms have been selected out of 2282 CNSsolve> evaluate ($counter=1) EVALUATE: symbol $COUNTER set to 1.00000 (real) CNSsolve> evaluate ($done=false) EVALUATE: symbol $DONE set to FALSE (logical) CNSsolve> while ( $done = false ) loop prost NEXTCD: condition evaluated as true CNSsolve> if ( &exist_prost_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &BLANK%prost_coordinate_infile_$counter = false ) then NEXTCD: condition evaluated as false CNSsolve> do (refx=0) (all) CNSsolve> segment CNSsolve> chain CNSsolve> if ( &prost_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> coordinates @@&prost_coordinate_infile_$counter CNSsolve> end CNSsolve> end CNSsolve> if ( &prost_rename_$counter = true ) then CNSsolve> do (refy=$counter) (attr refx=9999) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) EVALUATE: symbol $COUNTER set to 2.00000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop prost CNSsolve> while ( $done = false ) loop prost NEXTCD: condition evaluated as true CNSsolve> if ( &exist_prost_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as false CNSsolve> if ( &BLANK%prost_coordinate_infile_$counter = false ) then CNSsolve> do (refx=0) (all) CNSsolve> segment CNSsolve> chain CNSsolve> if ( &prost_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> coordinates @@&prost_coordinate_infile_$counter CNSsolve> end CNSsolve> end CNSsolve> if ( &prost_rename_$counter = true ) then CNSsolve> do (refy=$counter) (attr refx=9999) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) EVALUATE: symbol $DONE set to TRUE (logical) CNSsolve> end if CNSsolve> end loop prost CNSsolve> while ( $done = false ) loop prost NEXTCD: condition evaluated as false CNSsolve> if ( &exist_prost_coordinate_infile_$counter = true ) then CNSsolve> if ( &BLANK%prost_coordinate_infile_$counter = false ) then CNSsolve> do (refx=0) (all) CNSsolve> segment CNSsolve> chain CNSsolve> if ( &prost_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> coordinates @@&prost_coordinate_infile_$counter CNSsolve> end CNSsolve> end CNSsolve> if ( &prost_rename_$counter = true ) then CNSsolve> do (refy=$counter) (attr refx=9999) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop prost CNSsolve> CNSsolve> {* any special prosthetic group patches can be applied here *} CNSsolve> {===>} CNSsolve> {<===} CNSsolve> CNSsolve> evaluate ($counter=1) EVALUATE: symbol $COUNTER set to 1.00000 (real) CNSsolve> evaluate ($done=false) EVALUATE: symbol $DONE set to FALSE (logical) CNSsolve> while ( $done = false ) loop prost NEXTCD: condition evaluated as true CNSsolve> if ( &exist_prost_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &BLANK%prost_coordinate_infile_$counter = false ) then NEXTCD: condition evaluated as false CNSsolve> coor CNSsolve> if ( &prost_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> @@&prost_coordinate_infile_$counter CNSsolve> set echo=off end CNSsolve> show sum(1) ( not(hydrogen) and not(known) ) CNSsolve> if ( $select = 0 ) then CNSsolve> display %INFO: There are no coordinates missing for non-hydrogen atoms CNSsolve> end if CNSsolve> set echo=on end CNSsolve> if ( &prost_rename_$counter = true ) then CNSsolve> do (segid=capitalize(&prost_segid_$counter)) (attr refy=$counter) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) EVALUATE: symbol $COUNTER set to 2.00000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop prost CNSsolve> while ( $done = false ) loop prost NEXTCD: condition evaluated as true CNSsolve> if ( &exist_prost_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as false CNSsolve> if ( &BLANK%prost_coordinate_infile_$counter = false ) then CNSsolve> coor CNSsolve> if ( &prost_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> @@&prost_coordinate_infile_$counter CNSsolve> set echo=off end CNSsolve> show sum(1) ( not(hydrogen) and not(known) ) CNSsolve> if ( $select = 0 ) then CNSsolve> display %INFO: There are no coordinates missing for non-hydrogen atoms CNSsolve> end if CNSsolve> set echo=on end CNSsolve> if ( &prost_rename_$counter = true ) then CNSsolve> do (segid=capitalize(&prost_segid_$counter)) (attr refy=$counter) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) EVALUATE: symbol $DONE set to TRUE (logical) CNSsolve> end if CNSsolve> end loop prost CNSsolve> while ( $done = false ) loop prost NEXTCD: condition evaluated as false CNSsolve> if ( &exist_prost_coordinate_infile_$counter = true ) then CNSsolve> if ( &BLANK%prost_coordinate_infile_$counter = false ) then CNSsolve> coor CNSsolve> if ( &prost_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> @@&prost_coordinate_infile_$counter CNSsolve> set echo=off end CNSsolve> show sum(1) ( not(hydrogen) and not(known) ) CNSsolve> if ( $select = 0 ) then CNSsolve> display %INFO: There are no coordinates missing for non-hydrogen atoms CNSsolve> end if CNSsolve> set echo=on end CNSsolve> if ( &prost_rename_$counter = true ) then CNSsolve> do (segid=capitalize(&prost_segid_$counter)) (attr refy=$counter) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop prost CNSsolve> CNSsolve> do (refy=0) (all) SELRPN: 2282 atoms have been selected out of 2282 CNSsolve> evaluate ($counter=1) EVALUATE: symbol $COUNTER set to 1.00000 (real) CNSsolve> evaluate ($done=false) EVALUATE: symbol $DONE set to FALSE (logical) CNSsolve> while ( $done = false ) loop liga NEXTCD: condition evaluated as true CNSsolve> if ( &exist_lig_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &BLANK%lig_coordinate_infile_$counter = false ) then NEXTCD: condition evaluated as false CNSsolve> do (refx=0) (all) CNSsolve> segment CNSsolve> chain CNSsolve> if ( &lig_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> coordinates @@&lig_coordinate_infile_$counter CNSsolve> end CNSsolve> end CNSsolve> if ( &lig_rename_$counter = true ) then CNSsolve> do (refy=$counter) (attr refx=9999) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) EVALUATE: symbol $COUNTER set to 2.00000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop liga CNSsolve> while ( $done = false ) loop liga NEXTCD: condition evaluated as true CNSsolve> if ( &exist_lig_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as false CNSsolve> if ( &BLANK%lig_coordinate_infile_$counter = false ) then CNSsolve> do (refx=0) (all) CNSsolve> segment CNSsolve> chain CNSsolve> if ( &lig_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> coordinates @@&lig_coordinate_infile_$counter CNSsolve> end CNSsolve> end CNSsolve> if ( &lig_rename_$counter = true ) then CNSsolve> do (refy=$counter) (attr refx=9999) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) EVALUATE: symbol $DONE set to TRUE (logical) CNSsolve> end if CNSsolve> end loop liga CNSsolve> while ( $done = false ) loop liga NEXTCD: condition evaluated as false CNSsolve> if ( &exist_lig_coordinate_infile_$counter = true ) then CNSsolve> if ( &BLANK%lig_coordinate_infile_$counter = false ) then CNSsolve> do (refx=0) (all) CNSsolve> segment CNSsolve> chain CNSsolve> if ( &lig_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> coordinates @@&lig_coordinate_infile_$counter CNSsolve> end CNSsolve> end CNSsolve> if ( &lig_rename_$counter = true ) then CNSsolve> do (refy=$counter) (attr refx=9999) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop liga CNSsolve> CNSsolve> {* any special ligand patches can be applied here *} CNSsolve> {===>} CNSsolve> {<===} CNSsolve> CNSsolve> evaluate ($counter=1) EVALUATE: symbol $COUNTER set to 1.00000 (real) CNSsolve> evaluate ($done=false) EVALUATE: symbol $DONE set to FALSE (logical) CNSsolve> while ( $done = false ) loop liga NEXTCD: condition evaluated as true CNSsolve> if ( &exist_lig_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &BLANK%lig_coordinate_infile_$counter = false ) then NEXTCD: condition evaluated as false CNSsolve> coor CNSsolve> if ( &lig_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> @@&lig_coordinate_infile_$counter CNSsolve> set echo=off end CNSsolve> show sum(1) ( not(hydrogen) and not(known) ) CNSsolve> if ( $select = 0 ) then CNSsolve> display %INFO: There are no coordinates missing for non-hydrogen atoms CNSsolve> end if CNSsolve> set echo=on end CNSsolve> if ( &lig_rename_$counter = true ) then CNSsolve> do (segid=capitalize(&lig_segid_$counter)) (attr refy=$counter) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) EVALUATE: symbol $COUNTER set to 2.00000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop liga CNSsolve> while ( $done = false ) loop liga NEXTCD: condition evaluated as true CNSsolve> if ( &exist_lig_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as false CNSsolve> if ( &BLANK%lig_coordinate_infile_$counter = false ) then CNSsolve> coor CNSsolve> if ( &lig_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> @@&lig_coordinate_infile_$counter CNSsolve> set echo=off end CNSsolve> show sum(1) ( not(hydrogen) and not(known) ) CNSsolve> if ( $select = 0 ) then CNSsolve> display %INFO: There are no coordinates missing for non-hydrogen atoms CNSsolve> end if CNSsolve> set echo=on end CNSsolve> if ( &lig_rename_$counter = true ) then CNSsolve> do (segid=capitalize(&lig_segid_$counter)) (attr refy=$counter) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) EVALUATE: symbol $DONE set to TRUE (logical) CNSsolve> end if CNSsolve> end loop liga CNSsolve> while ( $done = false ) loop liga NEXTCD: condition evaluated as false CNSsolve> if ( &exist_lig_coordinate_infile_$counter = true ) then CNSsolve> if ( &BLANK%lig_coordinate_infile_$counter = false ) then CNSsolve> coor CNSsolve> if ( &lig_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> @@&lig_coordinate_infile_$counter CNSsolve> set echo=off end CNSsolve> show sum(1) ( not(hydrogen) and not(known) ) CNSsolve> if ( $select = 0 ) then CNSsolve> display %INFO: There are no coordinates missing for non-hydrogen atoms CNSsolve> end if CNSsolve> set echo=on end CNSsolve> if ( &lig_rename_$counter = true ) then CNSsolve> do (segid=capitalize(&lig_segid_$counter)) (attr refy=$counter) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop liga CNSsolve> CNSsolve> do (refy=0) (all) SELRPN: 2282 atoms have been selected out of 2282 CNSsolve> evaluate ($counter=1) EVALUATE: symbol $COUNTER set to 1.00000 (real) CNSsolve> evaluate ($done=false) EVALUATE: symbol $DONE set to FALSE (logical) CNSsolve> while ( $done = false ) loop ion NEXTCD: condition evaluated as true CNSsolve> if ( &exist_ion_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &BLANK%ion_coordinate_infile_$counter = false ) then NEXTCD: condition evaluated as false CNSsolve> do (refx=0) (all) CNSsolve> segment CNSsolve> chain CNSsolve> if ( &ion_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> coordinates @@&ion_coordinate_infile_$counter CNSsolve> end CNSsolve> end CNSsolve> if ( &ion_rename_$counter = true ) then CNSsolve> do (refy=$counter) (attr refx=9999) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) EVALUATE: symbol $COUNTER set to 2.00000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop ion CNSsolve> while ( $done = false ) loop ion NEXTCD: condition evaluated as true CNSsolve> if ( &exist_ion_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as false CNSsolve> if ( &BLANK%ion_coordinate_infile_$counter = false ) then CNSsolve> do (refx=0) (all) CNSsolve> segment CNSsolve> chain CNSsolve> if ( &ion_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> coordinates @@&ion_coordinate_infile_$counter CNSsolve> end CNSsolve> end CNSsolve> if ( &ion_rename_$counter = true ) then CNSsolve> do (refy=$counter) (attr refx=9999) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) EVALUATE: symbol $DONE set to TRUE (logical) CNSsolve> end if CNSsolve> end loop ion CNSsolve> while ( $done = false ) loop ion NEXTCD: condition evaluated as false CNSsolve> if ( &exist_ion_coordinate_infile_$counter = true ) then CNSsolve> if ( &BLANK%ion_coordinate_infile_$counter = false ) then CNSsolve> do (refx=0) (all) CNSsolve> segment CNSsolve> chain CNSsolve> if ( &ion_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> coordinates @@&ion_coordinate_infile_$counter CNSsolve> end CNSsolve> end CNSsolve> if ( &ion_rename_$counter = true ) then CNSsolve> do (refy=$counter) (attr refx=9999) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop ion CNSsolve> CNSsolve> {* any special ion patches can be applied here *} CNSsolve> {===>} CNSsolve> {<===} CNSsolve> CNSsolve> evaluate ($counter=1) EVALUATE: symbol $COUNTER set to 1.00000 (real) CNSsolve> evaluate ($done=false) EVALUATE: symbol $DONE set to FALSE (logical) CNSsolve> while ( $done = false ) loop ion NEXTCD: condition evaluated as true CNSsolve> if ( &exist_ion_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as true CNSsolve> if ( &BLANK%ion_coordinate_infile_$counter = false ) then NEXTCD: condition evaluated as false CNSsolve> coor CNSsolve> if ( &ion_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> @@&ion_coordinate_infile_$counter CNSsolve> set echo=off end CNSsolve> show sum(1) ( not(hydrogen) and not(known) ) CNSsolve> if ( $select = 0 ) then CNSsolve> display %INFO: There are no coordinates missing for non-hydrogen atoms CNSsolve> end if CNSsolve> set echo=on end CNSsolve> if ( &ion_rename_$counter = true ) then CNSsolve> do (segid=capitalize(&ion_segid_$counter)) (attr refy=$counter) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) EVALUATE: symbol $COUNTER set to 2.00000 (real) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop ion CNSsolve> while ( $done = false ) loop ion NEXTCD: condition evaluated as true CNSsolve> if ( &exist_ion_coordinate_infile_$counter = true ) then NEXTCD: condition evaluated as false CNSsolve> if ( &BLANK%ion_coordinate_infile_$counter = false ) then CNSsolve> coor CNSsolve> if ( &ion_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> @@&ion_coordinate_infile_$counter CNSsolve> set echo=off end CNSsolve> show sum(1) ( not(hydrogen) and not(known) ) CNSsolve> if ( $select = 0 ) then CNSsolve> display %INFO: There are no coordinates missing for non-hydrogen atoms CNSsolve> end if CNSsolve> set echo=on end CNSsolve> if ( &ion_rename_$counter = true ) then CNSsolve> do (segid=capitalize(&ion_segid_$counter)) (attr refy=$counter) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) EVALUATE: symbol $DONE set to TRUE (logical) CNSsolve> end if CNSsolve> end loop ion CNSsolve> while ( $done = false ) loop ion NEXTCD: condition evaluated as false CNSsolve> if ( &exist_ion_coordinate_infile_$counter = true ) then CNSsolve> if ( &BLANK%ion_coordinate_infile_$counter = false ) then CNSsolve> coor CNSsolve> if ( &ion_convert_$counter = true ) then CNSsolve> convert=true CNSsolve> end if CNSsolve> @@&ion_coordinate_infile_$counter CNSsolve> set echo=off end CNSsolve> show sum(1) ( not(hydrogen) and not(known) ) CNSsolve> if ( $select = 0 ) then CNSsolve> display %INFO: There are no coordinates missing for non-hydrogen atoms CNSsolve> end if CNSsolve> set echo=on end CNSsolve> if ( &ion_rename_$counter = true ) then CNSsolve> do (segid=capitalize(&ion_segid_$counter)) (attr refy=$counter) CNSsolve> end if CNSsolve> end if CNSsolve> evaluate ($counter=$counter+1) CNSsolve> else CNSsolve> evaluate ($done=true) CNSsolve> end if CNSsolve> end loop ion CNSsolve> CNSsolve> CNSsolve> {* any final patches can be applied here *} CNSsolve> {===>} CNSsolve> {<===} CNSsolve> CNSsolve> if (&hydrogen_flag=false) then NEXTCD: condition evaluated as false CNSsolve> delete selection=( hydrogen ) end CNSsolve> end if CNSsolve> CNSsolve> delete selection=( &atom_delete ) end SELRPN: 0 atoms have been selected out of 2282 SCRATC-warning: STORe selections erased. Status of internal molecular topology database: -> NATOM= 2282(MAXA= 1000000) NBOND= 2294(MAXB= 1000000) -> NTHETA= 4184(MAXT= 2000000) NGRP= 136(MAXGRP= 1000000) -> NPHI= 3549(MAXP= 2000000) NIMPHI= 1144(MAXIMP= 1000000) -> NNB= 1038(MAXNB= 1000000) CNSsolve> CNSsolve> {- read parameter files -} CNSsolve> {- here atom-based parameter files are possible as well since all atoms have been defined -} CNSsolve> parameter PARRDR> evaluate ($counter=1) EVALUATE: symbol $COUNTER set to 1.00000 (real) PARRDR> evaluate ($done=false) EVALUATE: symbol $DONE set to FALSE (logical) PARRDR> while ( $done = false ) loop read NEXTCD: condition evaluated as true PARRDR> if ( &exist_parameter_infile_$counter = true ) then NEXTCD: condition evaluated as true PARRDR> if ( &BLANK%parameter_infile_$counter = false ) then NEXTCD: condition evaluated as true PARRDR> @@¶meter_infile_$counter ASSFIL: file protein-allhdg5-4.param opened. PARRDR>remarks file toppar/protein-allhdg5-4.param PARRDR>remark PARRDR>remark for file protein-allhdg-ucl.top version UCL date 14-MAR-00 PARRDR>remark for file protein-allhdg-dih-ucl.top version UCL date 07-JUL-01 PARRDR>remark Geometric energy function parameters for distance geometry and PARRDR>remark simulated annealing. PARRDR>remark Original author: Michael Nilges, EMBL Heidelberg PARRDR>remark Modifications: Mark A. Williams, UCL London PARRDR>remark Several modifications for HADDOCK: Alexandre Bonvin, Utrecht Uni PARRDR> PARRDR>set message ? end eval ($old_message=$result) set echo ? end eval ($old_echo=$result) MESSage=NORM EVALUATE: symbol $OLD_MESSAGE set to "NORM" (string) ECHO=TRUE {ON} EVALUATE: symbol $OLD_ECHO set to TRUE (logical) PARRDR>set message=off echo=off end PARRDR> PARRDR> end if PARRDR> else PARRDR> evaluate ($done=true) PARRDR> end if PARRDR> evaluate ($counter=$counter+1) EVALUATE: symbol $COUNTER set to 2.00000 (real) PARRDR> end loop read PARRDR> while ( $done = false ) loop read NEXTCD: condition evaluated as true PARRDR> if ( &exist_parameter_infile_$counter = true ) then NEXTCD: condition evaluated as true PARRDR> if ( &BLANK%parameter_infile_$counter = false ) then NEXTCD: condition evaluated as true PARRDR> @@¶meter_infile_$counter ASSFIL: file water-allhdg5-4.param opened. PARRDR>remarks file toppar/water-allhdg5-4.param PARRDR>remarks PARRDR>remarks PARAM19.SOL (solvent parameters) PARRDR>remarks =========== PARRDR>remarks available: TIPS3P and DMSO model PARRDR> PARRDR>set message ? end eval ($old_message=$result) set echo ? end eval ($old_echo=$result) MESSage=NORM EVALUATE: symbol $OLD_MESSAGE set to "NORM" (string) ECHO=TRUE {ON} EVALUATE: symbol $OLD_ECHO set to TRUE (logical) PARRDR>set message=off echo=off end PARRDR> end if PARRDR> else PARRDR> evaluate ($done=true) PARRDR> end if PARRDR> evaluate ($counter=$counter+1) EVALUATE: symbol $COUNTER set to 3.00000 (real) PARRDR> end loop read PARRDR> while ( $done = false ) loop read NEXTCD: condition evaluated as true PARRDR> if ( &exist_parameter_infile_$counter = true ) then NEXTCD: condition evaluated as true PARRDR> if ( &BLANK%parameter_infile_$counter = false ) then NEXTCD: condition evaluated as true PARRDR> @@¶meter_infile_$counter ASSFIL: file ion.param opened. PARRDR>remarks file toppar/ion.param PARRDR>remarks nonbonded parameters for common ions PARRDR>remarks new parameters derived from literature for single atom species PARRDR>remarks PDA 02/09/99 PARRDR> PARRDR>set message ? end eval ($old_message=$result) set echo ? end eval ($old_echo=$result) MESSage=NORM EVALUATE: symbol $OLD_MESSAGE set to "NORM" (string) ECHO=TRUE {ON} EVALUATE: symbol $OLD_ECHO set to TRUE (logical) PARRDR>set message=off echo=off end PARRDR> PARRDR> end if PARRDR> else PARRDR> evaluate ($done=true) PARRDR> end if PARRDR> evaluate ($counter=$counter+1) EVALUATE: symbol $COUNTER set to 4.00000 (real) PARRDR> end loop read PARRDR> while ( $done = false ) loop read NEXTCD: condition evaluated as true PARRDR> if ( &exist_parameter_infile_$counter = true ) then NEXTCD: condition evaluated as false PARRDR> if ( &BLANK%parameter_infile_$counter = false ) then PARRDR> @@¶meter_infile_$counter PARRDR> end if PARRDR> else PARRDR> evaluate ($done=true) EVALUATE: symbol $DONE set to TRUE (logical) PARRDR> end if PARRDR> evaluate ($counter=$counter+1) EVALUATE: symbol $COUNTER set to 5.00000 (real) PARRDR> end loop read PARRDR> while ( $done = false ) loop read NEXTCD: condition evaluated as false PARRDR> if ( &exist_parameter_infile_$counter = true ) then PARRDR> if ( &BLANK%parameter_infile_$counter = false ) then PARRDR> @@¶meter_infile_$counter PARRDR> end if PARRDR> else PARRDR> evaluate ($done=true) PARRDR> end if PARRDR> evaluate ($counter=$counter+1) PARRDR> end loop read PARRDR> end %PARRDR-ERROR: duplication of bond OUF SUF %PARRDR-ERROR: duplication of angle OUF SUF OUF %PARRDR-ERROR: duplication of nonbonded entry SUF %PARRDR-ERROR: duplication of nonbonded entry OUF CNSsolve> CNSsolve> {- build atoms if requested -} CNSsolve> identity (store1) (&atom_build) SELRPN: 10 atoms have been selected out of 2282 CNSsolve> CNSsolve> show sum(1) (store1) SELRPN: 10 atoms have been selected out of 2282 SHOW: sum over selected elements = 10.000000 CNSsolve> evaluate ($tobuild=$result) EVALUATE: symbol $TOBUILD set to 10.0000 (real) CNSsolve> CNSsolve> if ( $tobuild > 0 ) then NEXTCD: condition evaluated as true CNSsolve> CNSsolve> fix selection=(not(store1)) end SELRPN: 2272 atoms have been selected out of 2282 CNSsolve> CNSsolve> for $id in id (tag and byres(store1)) loop avco SELRPN: 7 atoms have been selected out of 2282 FOR ID LOOP: symbol ID set to 1.00000 (real) CNSsolve> CNSsolve> show ave(x) (byres(id $id) and known) SELRPN: 16 atoms have been selected out of 2282 SHOW: average of selected elements = 4.953750 CNSsolve> evaluate ($ave_x=$result) EVALUATE: symbol $AVE_X set to 4.95375 (real) CNSsolve> show ave(y) (byres(id $id) and known) SELRPN: 16 atoms have been selected out of 2282 SHOW: average of selected elements = 0.056250 CNSsolve> evaluate ($ave_y=$result) EVALUATE: symbol $AVE_Y set to 0.562500E-01 (real) CNSsolve> show ave(z) (byres(id $id) and known) SELRPN: 16 atoms have been selected out of 2282 SHOW: average of selected elements = -2.956250 CNSsolve> evaluate ($ave_z=$result) EVALUATE: symbol $AVE_Z set to -2.95625 (real) CNSsolve> CNSsolve> do (x=$ave_x) (byres(id $id) and store1) SELRPN: 3 atoms have been selected out of 2282 CNSsolve> do (y=$ave_y) (byres(id $id) and store1) SELRPN: 3 atoms have been selected out of 2282 CNSsolve> do (z=$ave_z) (byres(id $id) and store1) SELRPN: 3 atoms have been selected out of 2282 CNSsolve> CNSsolve> end loop avco CNSsolve> for $id in id (tag and byres(store1)) loop avco FOR ID LOOP: symbol ID set to 2174.00 (real) CNSsolve> CNSsolve> show ave(x) (byres(id $id) and known) SELRPN: 17 atoms have been selected out of 2282 SHOW: average of selected elements = 14.636824 CNSsolve> evaluate ($ave_x=$result) EVALUATE: symbol $AVE_X set to 14.6368 (real) CNSsolve> show ave(y) (byres(id $id) and known) SELRPN: 17 atoms have been selected out of 2282 SHOW: average of selected elements = -3.139941 CNSsolve> evaluate ($ave_y=$result) EVALUATE: symbol $AVE_Y set to -3.13994 (real) CNSsolve> show ave(z) (byres(id $id) and known) SELRPN: 17 atoms have been selected out of 2282 SHOW: average of selected elements = 0.576294 CNSsolve> evaluate ($ave_z=$result) EVALUATE: symbol $AVE_Z set to 0.576294 (real) CNSsolve> CNSsolve> do (x=$ave_x) (byres(id $id) and store1) SELRPN: 1 atoms have been selected out of 2282 CNSsolve> do (y=$ave_y) (byres(id $id) and store1) SELRPN: 1 atoms have been selected out of 2282 CNSsolve> do (z=$ave_z) (byres(id $id) and store1) SELRPN: 1 atoms have been selected out of 2282 CNSsolve> CNSsolve> end loop avco CNSsolve> for $id in id (tag and byres(store1)) loop avco FOR ID LOOP: symbol ID set to 2192.00 (real) CNSsolve> CNSsolve> show ave(x) (byres(id $id) and known) SELRPN: 17 atoms have been selected out of 2282 SHOW: average of selected elements = 12.819941 CNSsolve> evaluate ($ave_x=$result) EVALUATE: symbol $AVE_X set to 12.8199 (real) CNSsolve> show ave(y) (byres(id $id) and known) SELRPN: 17 atoms have been selected out of 2282 SHOW: average of selected elements = 2.359059 CNSsolve> evaluate ($ave_y=$result) EVALUATE: symbol $AVE_Y set to 2.35906 (real) CNSsolve> show ave(z) (byres(id $id) and known) SELRPN: 17 atoms have been selected out of 2282 SHOW: average of selected elements = 1.277706 CNSsolve> evaluate ($ave_z=$result) EVALUATE: symbol $AVE_Z set to 1.27771 (real) CNSsolve> CNSsolve> do (x=$ave_x) (byres(id $id) and store1) SELRPN: 1 atoms have been selected out of 2282 CNSsolve> do (y=$ave_y) (byres(id $id) and store1) SELRPN: 1 atoms have been selected out of 2282 CNSsolve> do (z=$ave_z) (byres(id $id) and store1) SELRPN: 1 atoms have been selected out of 2282 CNSsolve> CNSsolve> end loop avco CNSsolve> for $id in id (tag and byres(store1)) loop avco FOR ID LOOP: symbol ID set to 2210.00 (real) CNSsolve> CNSsolve> show ave(x) (byres(id $id) and known) SELRPN: 17 atoms have been selected out of 2282 SHOW: average of selected elements = 17.885000 CNSsolve> evaluate ($ave_x=$result) EVALUATE: symbol $AVE_X set to 17.8850 (real) CNSsolve> show ave(y) (byres(id $id) and known) SELRPN: 17 atoms have been selected out of 2282 SHOW: average of selected elements = 1.996235 CNSsolve> evaluate ($ave_y=$result) EVALUATE: symbol $AVE_Y set to 1.99624 (real) CNSsolve> show ave(z) (byres(id $id) and known) SELRPN: 17 atoms have been selected out of 2282 SHOW: average of selected elements = -1.854647 CNSsolve> evaluate ($ave_z=$result) EVALUATE: symbol $AVE_Z set to -1.85465 (real) CNSsolve> CNSsolve> do (x=$ave_x) (byres(id $id) and store1) SELRPN: 1 atoms have been selected out of 2282 CNSsolve> do (y=$ave_y) (byres(id $id) and store1) SELRPN: 1 atoms have been selected out of 2282 CNSsolve> do (z=$ave_z) (byres(id $id) and store1) SELRPN: 1 atoms have been selected out of 2282 CNSsolve> CNSsolve> end loop avco CNSsolve> for $id in id (tag and byres(store1)) loop avco FOR ID LOOP: symbol ID set to 2228.00 (real) CNSsolve> CNSsolve> show ave(x) (byres(id $id) and known) SELRPN: 17 atoms have been selected out of 2282 SHOW: average of selected elements = 16.920647 CNSsolve> evaluate ($ave_x=$result) EVALUATE: symbol $AVE_X set to 16.9206 (real) CNSsolve> show ave(y) (byres(id $id) and known) SELRPN: 17 atoms have been selected out of 2282 SHOW: average of selected elements = 6.356235 CNSsolve> evaluate ($ave_y=$result) EVALUATE: symbol $AVE_Y set to 6.35624 (real) CNSsolve> show ave(z) (byres(id $id) and known) SELRPN: 17 atoms have been selected out of 2282 SHOW: average of selected elements = -4.749588 CNSsolve> evaluate ($ave_z=$result) EVALUATE: symbol $AVE_Z set to -4.74959 (real) CNSsolve> CNSsolve> do (x=$ave_x) (byres(id $id) and store1) SELRPN: 1 atoms have been selected out of 2282 CNSsolve> do (y=$ave_y) (byres(id $id) and store1) SELRPN: 1 atoms have been selected out of 2282 CNSsolve> do (z=$ave_z) (byres(id $id) and store1) SELRPN: 1 atoms have been selected out of 2282 CNSsolve> CNSsolve> end loop avco CNSsolve> for $id in id (tag and byres(store1)) loop avco FOR ID LOOP: symbol ID set to 2246.00 (real) CNSsolve> CNSsolve> show ave(x) (byres(id $id) and known) SELRPN: 17 atoms have been selected out of 2282 SHOW: average of selected elements = 17.942647 CNSsolve> evaluate ($ave_x=$result) EVALUATE: symbol $AVE_X set to 17.9426 (real) CNSsolve> show ave(y) (byres(id $id) and known) SELRPN: 17 atoms have been selected out of 2282 SHOW: average of selected elements = 2.124588 CNSsolve> evaluate ($ave_y=$result) EVALUATE: symbol $AVE_Y set to 2.12459 (real) CNSsolve> show ave(z) (byres(id $id) and known) SELRPN: 17 atoms have been selected out of 2282 SHOW: average of selected elements = -7.567235 CNSsolve> evaluate ($ave_z=$result) EVALUATE: symbol $AVE_Z set to -7.56724 (real) CNSsolve> CNSsolve> do (x=$ave_x) (byres(id $id) and store1) SELRPN: 1 atoms have been selected out of 2282 CNSsolve> do (y=$ave_y) (byres(id $id) and store1) SELRPN: 1 atoms have been selected out of 2282 CNSsolve> do (z=$ave_z) (byres(id $id) and store1) SELRPN: 1 atoms have been selected out of 2282 CNSsolve> CNSsolve> end loop avco CNSsolve> for $id in id (tag and byres(store1)) loop avco FOR ID LOOP: symbol ID set to 2264.00 (real) CNSsolve> CNSsolve> show ave(x) (byres(id $id) and known) SELRPN: 17 atoms have been selected out of 2282 SHOW: average of selected elements = 19.897471 CNSsolve> evaluate ($ave_x=$result) EVALUATE: symbol $AVE_X set to 19.8975 (real) CNSsolve> show ave(y) (byres(id $id) and known) SELRPN: 17 atoms have been selected out of 2282 SHOW: average of selected elements = 5.607000 CNSsolve> evaluate ($ave_y=$result) EVALUATE: symbol $AVE_Y set to 5.60700 (real) CNSsolve> show ave(z) (byres(id $id) and known) SELRPN: 17 atoms have been selected out of 2282 SHOW: average of selected elements = -10.797235 CNSsolve> evaluate ($ave_z=$result) EVALUATE: symbol $AVE_Z set to -10.7972 (real) CNSsolve> CNSsolve> do (x=$ave_x) (byres(id $id) and store1) SELRPN: 2 atoms have been selected out of 2282 CNSsolve> do (y=$ave_y) (byres(id $id) and store1) SELRPN: 2 atoms have been selected out of 2282 CNSsolve> do (z=$ave_z) (byres(id $id) and store1) SELRPN: 2 atoms have been selected out of 2282 CNSsolve> CNSsolve> end loop avco CNSsolve> for $id in id (tag and byres(store1)) loop avco CNSsolve> CNSsolve> show ave(x) (byres(id $id) and known) CNSsolve> evaluate ($ave_x=$result) CNSsolve> show ave(y) (byres(id $id) and known) CNSsolve> evaluate ($ave_y=$result) CNSsolve> show ave(z) (byres(id $id) and known) CNSsolve> evaluate ($ave_z=$result) CNSsolve> CNSsolve> do (x=$ave_x) (byres(id $id) and store1) CNSsolve> do (y=$ave_y) (byres(id $id) and store1) CNSsolve> do (z=$ave_z) (byres(id $id) and store1) CNSsolve> CNSsolve> end loop avco CNSsolve> CNSsolve> do (x=x+random(2.0)) (store1) SELRPN: 10 atoms have been selected out of 2282 CNSsolve> do (y=y+random(2.0)) (store1) SELRPN: 10 atoms have been selected out of 2282 CNSsolve> do (z=z+random(2.0)) (store1) SELRPN: 10 atoms have been selected out of 2282 CNSsolve> CNSsolve> {- start parameter for the side chain building -} CNSsolve> parameter PARRDR> nbonds NBDSET> rcon=20. nbxmod=-2 repel=0.9 wmin=0.1 tolerance=1. NBDSET> rexp=2 irexp=2 inhibit=0.25 NBDSET> end PARRDR> end %PARRDR-ERROR: duplication of bond OUF SUF %PARRDR-ERROR: duplication of angle OUF SUF OUF %PARRDR-ERROR: duplication of nonbonded entry SUF %PARRDR-ERROR: duplication of nonbonded entry OUF CNSsolve> CNSsolve> {- Friction coefficient, in 1/ps. -} CNSsolve> do (fbeta=100) (store1) SELRPN: 10 atoms have been selected out of 2282 CNSsolve> CNSsolve> evaluate ($bath=300.0) EVALUATE: symbol $BATH set to 300.000 (real) CNSsolve> evaluate ($nstep=500) EVALUATE: symbol $NSTEP set to 500.000 (real) CNSsolve> evaluate ($timestep=0.0005) EVALUATE: symbol $TIMESTEP set to 0.500000E-03 (real) CNSsolve> CNSsolve> do (refy=mass) (store1) SELRPN: 10 atoms have been selected out of 2282 CNSsolve> CNSsolve> do (mass=20) (store1) SELRPN: 10 atoms have been selected out of 2282 CNSsolve> CNSsolve> igroup interaction SELRPN> (store1) (store1 or known) SELRPN: 10 atoms have been selected out of 2282 SELRPN: 2282 atoms have been selected out of 2282 IGROup> end CNSsolve> CNSsolve> {- turn on initial energy terms -} CNSsolve> flags exclude * include bond angle vdw end CNSsolve> CNSsolve> minimize lbfgs nstep=50 nprint=10 end LBFGS: number of degrees of freedom= 30 NBONDS: generating intra-molecular exclusion list with mode=-2 MAKINB: mode -2 found 10 exclusions and 0 interactions(1-4) NBONDS: found 874 intra-atom interactions NBONDS: found 699 intra-atom interactions NBONDS: found 757 intra-atom interactions NBONDS: found 768 intra-atom interactions NBONDS: found 742 intra-atom interactions ------------------------------ cycle= 10 ----------------------------------- | Etotal =2965.183 grad(E)=37.667 E(BOND)=1235.078 E(ANGL)=51.459 | | E(VDW )=1678.646 | ------------------------------------------------------------------------------- NBONDS: found 751 intra-atom interactions ------------------------------ cycle= 20 ----------------------------------- | Etotal =2896.313 grad(E)=36.149 E(BOND)=1189.414 E(ANGL)=31.859 | | E(VDW )=1675.040 | ------------------------------------------------------------------------------- ------------------------------ cycle= 30 ----------------------------------- | Etotal =2893.550 grad(E)=36.083 E(BOND)=1174.237 E(ANGL)=32.898 | | E(VDW )=1686.414 | ------------------------------------------------------------------------------- ------------------------------ cycle= 40 ----------------------------------- | Etotal =2893.548 grad(E)=36.084 E(BOND)=1174.347 E(ANGL)=32.891 | | E(VDW )=1686.310 | ------------------------------------------------------------------------------- ------------------------------ cycle= 50 ----------------------------------- | Etotal =2893.548 grad(E)=36.084 E(BOND)=1174.346 E(ANGL)=32.892 | | E(VDW )=1686.310 | ------------------------------------------------------------------------------- LBFGS: normal termination - NSTEP limit reached CNSsolve> CNSsolve> do (vx=maxwell($bath)) (store1) SELRPN: 10 atoms have been selected out of 2282 CNSsolve> do (vy=maxwell($bath)) (store1) SELRPN: 10 atoms have been selected out of 2282 CNSsolve> do (vz=maxwell($bath)) (store1) SELRPN: 10 atoms have been selected out of 2282 CNSsolve> CNSsolve> flags exclude vdw include impr end CNSsolve> CNSsolve> dynamics cartesian Cartesian Dynamics> nstep=50 Cartesian Dynamics> timestep=$timestep Cartesian Dynamics> tcoupling=true temperature=$bath DCART: temperature coupling (TCOUpling) enabled Cartesian Dynamics> nprint=$nstep Cartesian Dynamics> cmremove=false Cartesian Dynamics> end -------------------------- Cartesian dynamics start --------------------------- | E(kin)+E(total)=1218.599 E(kin)=11.094 temperature=372.180 | | Etotal =1207.505 grad(E)=38.037 E(BOND)=1174.346 E(ANGL)=32.892 | | E(IMPR)=0.267 | ------------------------------------------------------------------------------- ----------------- final step= 50 at 0.02500 ps --------------------- | E(kin)+E(total)=335.598 E(kin)=309.999 temperature=10399.905 | | Etotal =25.599 grad(E)=5.329 E(BOND)=23.578 E(ANGL)=1.687 | | E(IMPR)=0.335 | ------------------------------------------------------------------------------- CENMAS: Information about center of free masses position [A] : 13.00647 2.60147 -3.60648 velocity [A/ps] : 0.74651 6.18435 4.49281 ang. mom. [amu A/ps] : -12501.68725 7174.13222 -6951.28192 kin. ener. [Kcal/mol] : 14.09867 CNSsolve> CNSsolve> flags include vdw end CNSsolve> CNSsolve> minimize lbfgs nstep=50 nprint=10 end LBFGS: number of degrees of freedom= 30 ------------------------------ cycle= 10 ----------------------------------- | Etotal =2893.832 grad(E)=36.145 E(BOND)=1177.228 E(ANGL)=33.488 | | E(IMPR)=0.157 E(VDW )=1682.959 | ------------------------------------------------------------------------------- ------------------------------ cycle= 20 ----------------------------------- | Etotal =2893.755 grad(E)=36.107 E(BOND)=1174.241 E(ANGL)=33.220 | | E(IMPR)=0.160 E(VDW )=1686.134 | ------------------------------------------------------------------------------- ------------------------------ cycle= 30 ----------------------------------- | Etotal =2893.755 grad(E)=36.107 E(BOND)=1174.159 E(ANGL)=33.219 | | E(IMPR)=0.160 E(VDW )=1686.216 | ------------------------------------------------------------------------------- ------------------------------ cycle= 40 ----------------------------------- | Etotal =2893.755 grad(E)=36.107 E(BOND)=1174.158 E(ANGL)=33.219 | | E(IMPR)=0.160 E(VDW )=1686.217 | ------------------------------------------------------------------------------- LBFGS: Line search terminated LBFGS: maximum number of function evaluations reached CNSsolve> CNSsolve> do (vx=maxwell($bath)) (store1) SELRPN: 10 atoms have been selected out of 2282 CNSsolve> do (vy=maxwell($bath)) (store1) SELRPN: 10 atoms have been selected out of 2282 CNSsolve> do (vz=maxwell($bath)) (store1) SELRPN: 10 atoms have been selected out of 2282 CNSsolve> CNSsolve> dynamics cartesian Cartesian Dynamics> nstep=50 Cartesian Dynamics> timestep=$timestep Cartesian Dynamics> tcoupling=true temperature=$bath DCART: temperature coupling (TCOUpling) enabled Cartesian Dynamics> nprint=$nstep Cartesian Dynamics> cmremove=false Cartesian Dynamics> end -------------------------- Cartesian dynamics start --------------------------- | E(kin)+E(total)=2903.620 E(kin)=9.866 temperature=330.970 | | Etotal =2893.755 grad(E)=36.107 E(BOND)=1174.158 E(ANGL)=33.219 | | E(IMPR)=0.160 E(VDW )=1686.217 | ------------------------------------------------------------------------------- ----------------- final step= 50 at 0.02500 ps --------------------- | E(kin)+E(total)=2910.748 E(kin)=5.144 temperature=172.580 | | Etotal =2905.604 grad(E)=36.301 E(BOND)=1140.178 E(ANGL)=41.231 | | E(IMPR)=1.004 E(VDW )=1723.192 | ------------------------------------------------------------------------------- CENMAS: Information about center of free masses position [A] : 13.03809 2.68127 -3.59412 velocity [A/ps] : -0.07758 -0.01786 0.28742 ang. mom. [amu A/ps] : 1136.59333 -472.14575 2557.02729 kin. ener. [Kcal/mol] : 0.02126 CNSsolve> CNSsolve> parameter PARRDR> nbonds NBDSET> rcon=2. nbxmod=-3 repel=0.75 NBDSET> end PARRDR> end %PARRDR-ERROR: duplication of bond OUF SUF %PARRDR-ERROR: duplication of angle OUF SUF OUF %PARRDR-ERROR: duplication of nonbonded entry SUF %PARRDR-ERROR: duplication of nonbonded entry OUF CNSsolve> CNSsolve> minimize lbfgs nstep=100 nprint=25 end LBFGS: number of degrees of freedom= 30 NBONDS: generating intra-molecular exclusion list with mode=-3 MAKINB: mode -3 found 30 exclusions and 0 interactions(1-4) NBONDS: found 736 intra-atom interactions ------------------------------ cycle= 25 ----------------------------------- | Etotal =0.422 grad(E)=0.146 E(BOND)=0.000 E(ANGL)=0.422 | | E(IMPR)=0.000 E(VDW )=0.000 | ------------------------------------------------------------------------------- ------------------------------ cycle= 50 ----------------------------------- | Etotal =0.403 grad(E)=0.119 E(BOND)=0.000 E(ANGL)=0.403 | | E(IMPR)=0.000 E(VDW )=0.000 | ------------------------------------------------------------------------------- ------------------------------ cycle= 75 ----------------------------------- | Etotal =0.403 grad(E)=0.119 E(BOND)=0.000 E(ANGL)=0.403 | | E(IMPR)=0.000 E(VDW )=0.000 | ------------------------------------------------------------------------------- ------------------------------ cycle= 100 ----------------------------------- | Etotal =0.403 grad(E)=0.119 E(BOND)=0.000 E(ANGL)=0.403 | | E(IMPR)=0.000 E(VDW )=0.000 | ------------------------------------------------------------------------------- LBFGS: normal termination - NSTEP limit reached CNSsolve> CNSsolve> do (vx=maxwell($bath)) (store1) SELRPN: 10 atoms have been selected out of 2282 CNSsolve> do (vy=maxwell($bath)) (store1) SELRPN: 10 atoms have been selected out of 2282 CNSsolve> do (vz=maxwell($bath)) (store1) SELRPN: 10 atoms have been selected out of 2282 CNSsolve> CNSsolve> dynamics cartesian Cartesian Dynamics> nstep=$nstep Cartesian Dynamics> timestep=$timestep Cartesian Dynamics> tcoupling=true temperature=$bath DCART: temperature coupling (TCOUpling) enabled Cartesian Dynamics> nprint=$nstep Cartesian Dynamics> cmremove=false Cartesian Dynamics> end -------------------------- Cartesian dynamics start --------------------------- | E(kin)+E(total)=14.194 E(kin)=13.791 temperature=462.666 | | Etotal =0.403 grad(E)=0.119 E(BOND)=0.000 E(ANGL)=0.403 | | E(IMPR)=0.000 E(VDW )=0.000 | ------------------------------------------------------------------------------- ----------------- final step= 500 at 0.25000 ps --------------------- | E(kin)+E(total)=14.086 E(kin)=8.763 temperature=293.994 | | Etotal =5.323 grad(E)=2.645 E(BOND)=0.657 E(ANGL)=3.897 | | E(IMPR)=0.769 E(VDW )=0.000 | ------------------------------------------------------------------------------- CENMAS: Information about center of free masses position [A] : 13.04425 2.57951 -3.63486 velocity [A/ps] : -0.00142 -0.73055 -0.02263 ang. mom. [amu A/ps] : -866.01056 -46.10955 -3118.68075 kin. ener. [Kcal/mol] : 0.12768 CNSsolve> CNSsolve> {- turn on all energy terms -} CNSsolve> flags include dihe ? end EFLAGS: the following energy flags are set EFLAGS: BOND ANGL DIHE IMPR VDW CNSsolve> CNSsolve> {- set repel to ~vdw radii -} CNSsolve> parameter PARRDR> nbonds NBDSET> repel=0.89 NBDSET> end PARRDR> end %PARRDR-ERROR: duplication of bond OUF SUF %PARRDR-ERROR: duplication of angle OUF SUF OUF %PARRDR-ERROR: duplication of nonbonded entry SUF %PARRDR-ERROR: duplication of nonbonded entry OUF CNSsolve> CNSsolve> minimize lbfgs nstep=500 nprint=50 end LBFGS: number of degrees of freedom= 30 NBONDS: generating intra-molecular exclusion list with mode=-3 MAKINB: mode -3 found 30 exclusions and 0 interactions(1-4) NBONDS: found 728 intra-atom interactions ------------------------------ cycle= 50 ----------------------------------- | Etotal =3.042 grad(E)=0.674 E(BOND)=0.030 E(ANGL)=0.580 | | E(DIHE)=0.002 E(IMPR)=0.001 E(VDW )=2.428 | ------------------------------------------------------------------------------- ------------------------------ cycle= 100 ----------------------------------- | Etotal =3.038 grad(E)=0.674 E(BOND)=0.031 E(ANGL)=0.577 | | E(DIHE)=0.002 E(IMPR)=0.000 E(VDW )=2.428 | ------------------------------------------------------------------------------- ------------------------------ cycle= 150 ----------------------------------- | Etotal =3.038 grad(E)=0.674 E(BOND)=0.031 E(ANGL)=0.577 | | E(DIHE)=0.002 E(IMPR)=0.000 E(VDW )=2.428 | ------------------------------------------------------------------------------- LBFGS: normal termination - gradient converged CNSsolve> CNSsolve> flags exclude * include bond angl impr dihe vdw end CNSsolve> CNSsolve> {- return masses to something sensible -} CNSsolve> do (mass=refy) (store1) SELRPN: 10 atoms have been selected out of 2282 CNSsolve> CNSsolve> do (vx=maxwell($bath)) (store1) SELRPN: 10 atoms have been selected out of 2282 CNSsolve> do (vy=maxwell($bath)) (store1) SELRPN: 10 atoms have been selected out of 2282 CNSsolve> do (vz=maxwell($bath)) (store1) SELRPN: 10 atoms have been selected out of 2282 CNSsolve> CNSsolve> dynamics cartesian Cartesian Dynamics> nstep=$nstep Cartesian Dynamics> timestep=$timestep Cartesian Dynamics> tcoupling=true temperature=$bath DCART: temperature coupling (TCOUpling) enabled Cartesian Dynamics> nprint=$nstep Cartesian Dynamics> cmremove=false Cartesian Dynamics> end -------------------------- Cartesian dynamics start --------------------------- | E(kin)+E(total)=9.507 E(kin)=6.469 temperature=217.012 | | Etotal =3.038 grad(E)=0.674 E(BOND)=0.031 E(ANGL)=0.577 | | E(DIHE)=0.002 E(IMPR)=0.000 E(VDW )=2.428 | ------------------------------------------------------------------------------- ----------------- final step= 500 at 0.25000 ps --------------------- | E(kin)+E(total)=19.245 E(kin)=9.557 temperature=320.612 | | Etotal =9.688 grad(E)=2.611 E(BOND)=1.606 E(ANGL)=4.561 | | E(DIHE)=0.049 E(IMPR)=0.055 E(VDW )=3.417 | ------------------------------------------------------------------------------- CENMAS: Information about center of free masses position [A] : 15.84564 6.16440 -7.88293 velocity [A/ps] : 1.49176 5.87862 1.19991 ang. mom. [amu A/ps] : 10.19680 112.83236 319.70412 kin. ener. [Kcal/mol] : 1.14519 CNSsolve> CNSsolve> {- some final minimisation -} CNSsolve> minimize lbfgs LBFGS> nstep=500 LBFGS> drop=40.0 LBFGS> nprint=50 LBFGS> end LBFGS: number of degrees of freedom= 30 ------------------------------ cycle= 50 ----------------------------------- | Etotal =3.027 grad(E)=0.673 E(BOND)=0.031 E(ANGL)=0.577 | | E(DIHE)=0.002 E(IMPR)=0.000 E(VDW )=2.416 | ------------------------------------------------------------------------------- ------------------------------ cycle= 100 ----------------------------------- | Etotal =3.026 grad(E)=0.673 E(BOND)=0.031 E(ANGL)=0.577 | | E(DIHE)=0.002 E(IMPR)=0.000 E(VDW )=2.417 | ------------------------------------------------------------------------------- LBFGS: normal termination - gradient converged CNSsolve> CNSsolve> print thres=0.02 bonds (atom-i |atom-j ) dist. equil. delta energy const. Number of violations greater 0.020: 0 RMS deviation= 0.002 CNSsolve> print thres=5. angles (atom-i |atom-j |atom-k ) angle equil. delta energy const. Number of violations greater 5.000: 0 RMS deviation= 0.435 CNSsolve> CNSsolve> fix selection=( none ) end SELRPN: 0 atoms have been selected out of 2282 CNSsolve> CNSsolve> end if CNSsolve> CNSsolve> set echo=false end SELRPN: 0 atoms have been selected out of 2282 SHOW: zero atoms selected NEXTCD: condition evaluated as true SELRPN: 0 atoms have been selected out of 2282 CNSsolve> CNSsolve> if (&set_bfactor=true) then NEXTCD: condition evaluated as false CNSsolve> do (b=&bfactor) ( all ) CNSsolve> else CNSsolve> show ave(b) (known and not(store1)) SELRPN: 2272 atoms have been selected out of 2282 SHOW: average of selected elements = 37.213116 CNSsolve> do (b=$result) (store1 and (attr b < 0.01)) SELRPN: 10 atoms have been selected out of 2282 CNSsolve> end if CNSsolve> CNSsolve> if (&set_occupancy=true) then NEXTCD: condition evaluated as false CNSsolve> do (q=&occupancy) ( all ) CNSsolve> end if CNSsolve> CNSsolve> set echo=false end SELRPN: 10 atoms have been selected out of 2282 SHOW: sum over selected elements = 10.000000 NEXTCD: condition evaluated as true SELRPN: 10 atoms have been selected out of 2282 FOR ID LOOP: symbol ID set to 16.0000 (real) SELRPN: 1 atoms have been selected out of 2282 ( A MET 1 HT1 ) A EVALUATE: symbol $SEGID set to "A" (string) SELRPN: 1 atoms have been selected out of 2282 ( A MET 1 HT1 ) MET EVALUATE: symbol $RESNAME set to "MET" (string) SELRPN: 1 atoms have been selected out of 2282 ( A MET 1 HT1 ) 1 EVALUATE: symbol $RESID set to "1" (string) SELRPN: 1 atoms have been selected out of 2282 ( A MET 1 HT1 ) HT1 EVALUATE: symbol $NAME set to "HT1" (string) FOR ID LOOP: symbol ID set to 17.0000 (real) SELRPN: 1 atoms have been selected out of 2282 ( A MET 1 HT2 ) A EVALUATE: symbol $SEGID set to "A" (string) SELRPN: 1 atoms have been selected out of 2282 ( A MET 1 HT2 ) MET EVALUATE: symbol $RESNAME set to "MET" (string) SELRPN: 1 atoms have been selected out of 2282 ( A MET 1 HT2 ) 1 EVALUATE: symbol $RESID set to "1" (string) SELRPN: 1 atoms have been selected out of 2282 ( A MET 1 HT2 ) HT2 EVALUATE: symbol $NAME set to "HT2" (string) FOR ID LOOP: symbol ID set to 18.0000 (real) SELRPN: 1 atoms have been selected out of 2282 ( A MET 1 HT3 ) A EVALUATE: symbol $SEGID set to "A" (string) SELRPN: 1 atoms have been selected out of 2282 ( A MET 1 HT3 ) MET EVALUATE: symbol $RESNAME set to "MET" (string) SELRPN: 1 atoms have been selected out of 2282 ( A MET 1 HT3 ) 1 EVALUATE: symbol $RESID set to "1" (string) SELRPN: 1 atoms have been selected out of 2282 ( A MET 1 HT3 ) HT3 EVALUATE: symbol $NAME set to "HT3" (string) FOR ID LOOP: symbol ID set to 2189.00 (real) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 129 HE2 ) A EVALUATE: symbol $SEGID set to "A" (string) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 129 HE2 ) HIS EVALUATE: symbol $RESNAME set to "HIS" (string) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 129 HE2 ) 129 EVALUATE: symbol $RESID set to "129" (string) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 129 HE2 ) HE2 EVALUATE: symbol $NAME set to "HE2" (string) FOR ID LOOP: symbol ID set to 2207.00 (real) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 130 HE2 ) A EVALUATE: symbol $SEGID set to "A" (string) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 130 HE2 ) HIS EVALUATE: symbol $RESNAME set to "HIS" (string) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 130 HE2 ) 130 EVALUATE: symbol $RESID set to "130" (string) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 130 HE2 ) HE2 EVALUATE: symbol $NAME set to "HE2" (string) FOR ID LOOP: symbol ID set to 2225.00 (real) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 131 HE2 ) A EVALUATE: symbol $SEGID set to "A" (string) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 131 HE2 ) HIS EVALUATE: symbol $RESNAME set to "HIS" (string) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 131 HE2 ) 131 EVALUATE: symbol $RESID set to "131" (string) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 131 HE2 ) HE2 EVALUATE: symbol $NAME set to "HE2" (string) FOR ID LOOP: symbol ID set to 2243.00 (real) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 132 HE2 ) A EVALUATE: symbol $SEGID set to "A" (string) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 132 HE2 ) HIS EVALUATE: symbol $RESNAME set to "HIS" (string) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 132 HE2 ) 132 EVALUATE: symbol $RESID set to "132" (string) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 132 HE2 ) HE2 EVALUATE: symbol $NAME set to "HE2" (string) FOR ID LOOP: symbol ID set to 2261.00 (real) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 133 HE2 ) A EVALUATE: symbol $SEGID set to "A" (string) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 133 HE2 ) HIS EVALUATE: symbol $RESNAME set to "HIS" (string) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 133 HE2 ) 133 EVALUATE: symbol $RESID set to "133" (string) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 133 HE2 ) HE2 EVALUATE: symbol $NAME set to "HE2" (string) FOR ID LOOP: symbol ID set to 2278.00 (real) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 134 HE2 ) A EVALUATE: symbol $SEGID set to "A" (string) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 134 HE2 ) HIS EVALUATE: symbol $RESNAME set to "HIS" (string) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 134 HE2 ) 134 EVALUATE: symbol $RESID set to "134" (string) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 134 HE2 ) HE2 EVALUATE: symbol $NAME set to "HE2" (string) FOR ID LOOP: symbol ID set to 2282.00 (real) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 134 OXT ) A EVALUATE: symbol $SEGID set to "A" (string) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 134 OXT ) HIS EVALUATE: symbol $RESNAME set to "HIS" (string) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 134 OXT ) 134 EVALUATE: symbol $RESID set to "134" (string) SELRPN: 1 atoms have been selected out of 2282 ( A HIS 134 OXT ) OXT EVALUATE: symbol $NAME set to "OXT" (string) CNSsolve> CNSsolve> set remarks=reset end CNSsolve> set remarks=accumulate end CNSsolve> CNSsolve> buffer message BUFFER> to=remarks BUFFER> dump BUFFER> end CNSsolve> CNSsolve> !write structure output=&structure_outfile end CNSsolve> CNSsolve> write coordinates format=PDBO output=&coordinate_outfile selection=( known ) end SELRPN: 2282 atoms have been selected out of 2282 ASSFIL: file sa_cns_6.pdb opened. CNSsolve> CNSsolve> stop ============================================================ Maximum dynamic memory allocation: 8841840 bytes Maximum dynamic memory overhead: 800 bytes Program started at: 14:40:10 on 29-Apr-2011 Program stopped at: 14:40:11 on 29-Apr-2011 CPU time used: 0.5839 seconds ============================================================