|
Difference between revisions of "geninput.pl"
(10 intermediate revisions by 5 users not shown) | |||
Line 99: | Line 99: | ||
; -help : usage information | ; -help : usage information | ||
− | + | ; -d force res1:res2[=...] : residue-residue distance restraints | |
+ | ; -par name=value[,...] : SICHO force field parameters | ||
+ | ; -seq file : MONSSTER sequence file | ||
== Examples == | == Examples == |
Latest revision as of 02:22, 29 July 2009
Usage
usage: geninput.pl [options] options: [-d force res1:res2[=res1:res2 ...]] [-par random=val,ncycle=val,icycle=val, tsteps=val,resmin=val,resmax=val, temp=first[:last],softcore=val, central=val,stiff=val,pair=val, kdcore=val,hbond=val,short=val, burial=val,multibody=val,threebody=val] [-seq file]
Description
Generates an input file for MONSSTER.
Available parameters and their default values are shown
in the table below. The values can be overwritten by specifying
the parameter name and the new value with -par
Key | Explanation | Default |
random | random seed | random value |
ncycle | cycles/temp. step | 20 |
icycle | MC moves/cycle | 50 |
tsteps | temperature steps | 10 |
resmin | first fragment residue (fixed template) | 0 |
resmax | last fragment residue (fixed template) | 0 |
temp | temperature range | 2.5:1.0 |
softcore | soft core radius for large residues | 4.0 |
central | centrosymmetric potential scaling | 0.25 |
stiff | general stiffness potential scaling | 1.25 |
pair | pair potential scaling | 1.0 |
kdcore | hydrophobic core potential scaling | 0.35 |
hbond | H-bond cooperativity potential scaling | -0.3 |
short | short range potential scaling | 0.5 |
burial | burial potential scaling | 1.0 |
multibody | multi body potential scaling | 0.5 |
threebody | 3 body potential scaling | 0.25 |
Distance constraints between residues (e.g. if disulfide linkages are known) may be specified using -d with a force constant and a list of residue pairs as parameters. A MONSSTER sequence file may be given with -seq for translating the residue numbering accordingly.
Options
- -help
- usage information
- -d force res1:res2[=...]
- residue-residue distance restraints
- -par name=value[,...]
- SICHO force field parameters
- -seq file
- MONSSTER sequence file
Examples
geninput.pl
generates the default MONSSTER input file
520735 20 50 10 0 0 2.50 1.00 4.000 0.250 1.250 1.000 0.350 -0.300 0.500 1.000 0.500 0.250 0 0.0 0 0 0 0
geninput.pl -par ncycle=50,stiff=0.1 -d 5.0 10:34=19:26
generates a MONSSTER input file iwth 50 cycles per time step and a
modified scaling factor for the general stiffness potential term.
In addition sidechain restraints are specified between residues
10 and 34 as well as 19 and 26.
657610 50 50 10 0 0 2.50 1.00 4.000 0.250 0.100 1.000 0.350 -0.300 0.500 1.000 0.500 0.250 2 5.000000 10 34 19 26 0 0 ...