-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathD0.5.inp
170 lines (151 loc) · 8.92 KB
/
D0.5.inp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
structure No_Nter_Gnrh_1_SC.psf
coordinates No_Nter_Gnrh_1_SC.pdb
set temp 300;
set outputname D0.5;
source step5_input.str
set inputname D0.4;
outputname $outputname;
binCoordinates $inputname.coor;
binVelocities $inputname.vel;
extendedSystem $inputname.xsc;
#firsttimestep 260000;
restartfreq 1000; # 1000 steps = every 2ps
dcdfreq 5000;
dcdUnitCell yes; # the file will contain unit cell info in the style of
# charmm dcd files. if yes, the dcd files will contain
# unit cell information in the style of charmm DCD files.
xstFreq 5000; # XSTFreq: control how often the extended systen configuration
# will be appended to the XST file
outputEnergies 125; # 125 steps = every 0.25ps
# The number of timesteps between each energy output of NAMD
outputTiming 1000; # The number of timesteps between each timing output shows
# time per step and time to completion
# Force-Field Parameters
paraTypeCharmm on; # We're using charmm type parameter file(s)
# multiple definitions may be used but only one file per definition
parameters toppar/par_all36m_prot.prm
parameters toppar/par_all36_na.prm
parameters toppar/par_all36_carb.prm
parameters toppar/par_all36_lipid.prm
parameters toppar/par_all36_cgenff.prm
parameters toppar/par_interface.prm
parameters toppar/toppar_all36_moreions.str
parameters toppar/toppar_all36_nano_lig.str
parameters toppar/toppar_all36_nano_lig_patch.str
parameters toppar/toppar_all36_synthetic_polymer.str
parameters toppar/toppar_all36_synthetic_polymer_patch.str
parameters toppar/toppar_all36_polymer_solvent.str
parameters toppar/toppar_water_ions.str
parameters toppar/toppar_dum_noble_gases.str
parameters toppar/toppar_ions_won.str
parameters toppar/cam.str
parameters toppar/toppar_all36_prot_arg0.str
parameters toppar/toppar_all36_prot_c36m_d_aminoacids.str
parameters toppar/toppar_all36_prot_fluoro_alkanes.str
parameters toppar/toppar_all36_prot_heme.str
parameters toppar/toppar_all36_prot_na_combined.str
parameters toppar/toppar_all36_prot_retinol.str
parameters toppar/toppar_all36_prot_model.str
parameters toppar/toppar_all36_prot_modify_res.str
parameters toppar/toppar_all36_na_nad_ppi.str
parameters toppar/toppar_all36_na_rna_modified.str
parameters toppar/toppar_all36_lipid_sphingo.str
parameters toppar/toppar_all36_lipid_archaeal.str
parameters toppar/toppar_all36_lipid_bacterial.str
parameters toppar/toppar_all36_lipid_cardiolipin.str
parameters toppar/toppar_all36_lipid_cholesterol.str
parameters toppar/toppar_all36_lipid_dag.str
parameters toppar/toppar_all36_lipid_inositol.str
parameters toppar/toppar_all36_lipid_lnp.str
parameters toppar/toppar_all36_lipid_lps.str
parameters toppar/toppar_all36_lipid_mycobacterial.str
parameters toppar/toppar_all36_lipid_miscellaneous.str
parameters toppar/toppar_all36_lipid_model.str
parameters toppar/toppar_all36_lipid_prot.str
parameters toppar/toppar_all36_lipid_tag.str
parameters toppar/toppar_all36_lipid_yeast.str
parameters toppar/toppar_all36_lipid_hmmm.str
parameters toppar/toppar_all36_lipid_detergent.str
parameters toppar/toppar_all36_lipid_ether.str
parameters toppar/toppar_all36_lipid_oxidized.str
parameters toppar/toppar_all36_carb_glycolipid.str
parameters toppar/toppar_all36_carb_glycopeptide.str
parameters toppar/toppar_all36_carb_imlab.str
parameters toppar/toppar_all36_label_spin.str
parameters toppar/toppar_all36_label_fluorophore.str
# Nonbonded Parameters
exclude scaled1-4 # non-bonded exclusion policy to use "none,1-2,1-3,1-4,or scaled1-4"
# 1-2: all atoms pairs that are bonded are going to be ignored
# 1-3: 3 consecutively bonded are excluded
# scaled1-4: include all the 1-3, and modified 1-4 interactions
# electrostatic scaled by 1-4scaling factor 1.0
# vdW special 1-4 parameters in charmm parameter file.
1-4scaling 1.0
switching on
vdwForceSwitching on; # New option for force-based switching of vdW
# if both switching and vdwForceSwitching are on CHARMM force
# switching is used for vdW forces.
# You have some freedom choosing the cutoff
cutoff 12.0; # may use smaller, maybe 10., with PME
switchdist 8.0; # cutoff - 2.
# switchdist - where you start to switch
# cutoff - where you stop accounting for nonbond interactions.
# correspondence in charmm:
# (cutnb,ctofnb,ctonnb = pairlistdist,cutoff,switchdist)
pairlistdist 14.0; # stores the all the pairs with in the distance it should be larger
# than cutoff( + 2.)
stepspercycle 20; # 20 redo pairlists every ten steps
pairlistsPerCycle 2; # 2 is the default
# cycle represents the number of steps between atom reassignments
# this means every 20/2=10 steps the pairlist will be updated
# Integrator Parameters
timestep 1.0; # fs/step
rigidBonds water; # Bound constraint all bonds involving H are fixed in length
nonbondedFreq 1; # nonbonded forces every step
fullElectFrequency 1; # PME every step
# Constant Temperature Control ONLY DURING EQUILB
reassignFreq 500; # reassignFreq: use this to reassign velocity every 500 steps
reassignTemp $temp;
if { $boxtype == "hexa" } {
set wrapnearst on
} else {
set wrapnearst off
}
wrapWater on; # wrap water to central cell
wrapAll on; # wrap other molecules too
wrapNearest $wrapnearst; # use for non-rectangular cells (wrap to the nearest image)
# PME (for full-system periodic electrostatics)
PME yes;
PMEInterpOrder 6; # interpolation order (spline order 6 in charmm)
PMEGridSpacing 1.0; # maximum PME grid space / used to calculate grid size
# Pressure and volume control
useGroupPressure yes; # use a hydrogen-group based pseudo-molecular viral to calcualte pressure and
# has less fluctuation, is needed for rigid bonds (rigidBonds/SHAKE)
useFlexibleCell yes; # yes for anisotropic system like membrane
useConstantRatio yes; # keeps the ratio of the unit cell in the x-y plane constant A=B
langevin on
langevinDamping 5.0
langevinTemp $temp
langevinHydrogen off
# constant pressure
langevinPiston on
langevinPistonTarget 1.01325
langevinPistonPeriod 50.0
langevinPistonDecay 25.0
langevinPistonTemp $temp
constraints on
consexp 2
consref restraints/prot_posres.ref
conskfile restraints/prot_posres.ref
conskcol B
constraintScaling 2.5
# planar restraint
exec sed -e "s/Constant \$fc/Constant 2/g" step5_input.colvar.str > restraints/$outputname.col
colvars on
colvarsConfig restraints/$outputname.col
# dihedral restraint
exec sed -e "s/\$FC/100/g" restraints/dihe.txt > restraints/$outputname.dihe
extraBonds yes
extraBondsFile restraints/$outputname.dihe
numsteps 90000000
run 125000