Hi,
we are trying to generate 3D-models for a short peptide, on the basis of two template structures. We would like to restrain the mainchain dihedral angles in the central region of our peptide during the modelling procedure.
To do this, we tried to modify the following input file which allowed us to use the 'model' routine and some special restraints at the same time. ======================================================================= #Homology modelling by the MODELLER TOP routine 'model'.
INCLUDE # Include the predefined TOP routines SET OUTPUT_CONTROL = 2 2 2 2 2 SET ALNFILE = 'ali1.pir' # alignment filename SET KNOWNS = 'sero_c' 'sero_o' # codes of the templates SET SEQUENCE = 'agen7' #code of the target SET STARTING_MODEL= 1 # index of the first model SET ENDING_MODEL = 1 # index of the last model CALL ROUTINE = 'model' # do homology modelling STOP SUBROUTINE ROUTINE = 'special_restraints' SET ADD_RESTRAINTS = on MAKE_RESTRAINTS RESTRAINT_TYPE = 'strand', RESIDUE_IDS = '3' '6' RETURN END_SUBROUTINE
==========================================================================
We thought it would be possible to modify the 'special_restraints' routine on the basis of the example file enforcing a cis conformation for a proline residue:
=========================================================================# Example for: ADD_RESTRAINT, DELETE_RESTRAINT
# This will enforce cis conformation for Pro-56.
# Make a model and stereochemical restraints: DEFINE_STRING VARIABLES = ATOM_IDS1 ATOM_IDS2 READ_TOPOLOGY FILE = '$(LIB)/top_heav.lib' READ_PARAMETERS FILE = '$(LIB)/par.lib' READ_MODEL FILE = '1fas' SEQUENCE_TO_ALI ATOM_FILES = '1fas', ALIGN_CODES = '1fas' SEQUENCE_TO_ALI ADD_SEQUENCE = on, ATOM_FILES = ATOM_FILES '1fas.ini', ; ALIGN_CODES = ALIGN_CODES '1fas-ini' GENERATE_TOPOLOGY SEQUENCE = '1fas-ini' TRANSFER_XYZ BUILD_MODEL INITIALIZE_XYZ = off MAKE_RESTRAINTS RESTRAINT_TYPE = 'stereo'
# Change the Pro-56 restraint from trans to cis: CALL ROUTINE = 'cispeptide', ATOM_IDS1 = 'O:56' 'C:56' 'N:57' 'CA:57', ; ATOM_IDS2 = 'CA:56' 'C:56' 'N:57' 'CA:57' WRITE_RESTRAINTS FILE = '1fas.rsr' ENERGY
SUBROUTINE ROUTINE = 'cispeptide' # Delete the old restraint on the same atoms: DELETE_RESTRAINT ATOM_IDS = ATOM_IDS1 # Add the new restraint: ADD_RESTRAINT RESTRAINT_PARAMETERS = 3 1 3 3 4 2 0 3.141593 0.087
DELETE_RESTRAINT ATOM_IDS = ATOM_IDS2 ADD_RESTRAINT RESTRAINT_PARAMETERS = 3 1 3 3 4 2 0 0.0 0.087
RETURN END_SUBROUTINE
=========================================================================
Doing this, we succeeded in writing the correct restraints file, but the 'model' routine stopped and no new model was generated.
We thought maybe it might be necessary to proceed in two steps: 1. writing the restraints file 2. reading this restraint file prior to run the 'model' routine.
Any help or suggestion would be appreciated.
Thank you
Isabelle