[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

dihedral restraints



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