PLUMED-NEST

Project ID: plumID:21.036
Source: DHH1N_sampling/BEMD/CV5/plumed.4.dat
Originally used with PLUMED version: 2.5.2
Stable: zipped raw stdout - zipped raw stderr - stderr
Master: zipped raw stdout - zipped raw stderr - stderr

Click on the labels of the actions for more information on what each action computes

#RESTART

MOLINFO
This command is used to provide information on the molecules that are present in your system. More details
 
STRUCTURE
a file in pdb format containing a reference structure
=Protein.pdb 
MOLTYPE
 what kind of molecule is contained in the pdb file - usually not needed since protein/RNA/DNA are compatible
=protein

The MOLINFO action with label  calculates something
WHOLEMOLECULES
This action is used to rebuild molecules that can become split by the periodic boundary conditions. More details
 
ENTITY0
the atoms that make up a molecule that you wish to align
=1-709

RANDOM_EXCHANGES
Set random pattern for exchanges. More details

cv5The ALPHARMSD action with label cv5 calculates the following quantities:
 Quantity   
 Type   
 Description  
cv5
scalar
if LESS_THAN is present the RMSD distance between each residue and the ideal alpha helix.  If LESS_THAN is not present the number of residue segments where the structure is similar to an alpha helix: 
ALPHARMSD
Probe the alpha helical content of a protein structure. This action is a shortcut and it has hidden defaults. More details
 
RESIDUES
this command is used to specify the set of residues that could conceivably form part of the secondary structure
=1-45 
TYPE
 the manner in which RMSD alignment is performed
=DRMSD 
R_0
The r_0 parameter of the switching function
=0.08 
NN
 The n parameter of the switching function
=2 
MM
 The m parameter of the switching function
=4 
NOPBC
 ignore the periodic boundary conditions
cv5: 
ALPHARMSD
Probe the alpha helical content of a protein structure. This action is a shortcut and uses the defaults shown here. More details
 
RESIDUES
this command is used to specify the set of residues that could conceivably form part of the secondary structure
=1-45 
TYPE
 the manner in which RMSD alignment is performed
=DRMSD 
R_0
The r_0 parameter of the switching function
=0.08 
NN
 The n parameter of the switching function
=2 
MM
 The m parameter of the switching function
=4 
NOPBC
 ignore the periodic boundary conditions
  
D_0
 The d_0 parameter of the switching function
=0.0
# PLUMED interprets the command:
# cv5: ALPHARMSD RESIDUES=1-45 TYPE=DRMSD R_0=0.08 NN=2 MM=4 NOPBC
# as follows (Click the red comment above to revert to the short version of the input):
cv5_rmsdThe SECONDARY_STRUCTURE_RMSD action with label cv5_rmsd calculates the following quantities:
 Quantity   
 Type   
 Description  
cv5_rmsd
vector
a vector containing the rmsd distance between each of the residue segments and the reference structure: 
SECONDARY_STRUCTURE_RMSD
Calclulate the distance between segments of a protein and a reference structure of interest More details
 
BONDLENGTH
the length to use for bonds
=0.17 
ATOMS
this is the full list of atoms that we are investigating
=1,5,7,8,9,10,12,14,19,20,21,23,25,38,39,40,42,44,52,53,54,56,58,66,67,68,70,72,80,81,82,84,86,100,101,102,104,106,114,115,116,118,120,128,129,130,132,134,142,143,144,146,148,156,157,158,160,162,170,171,172,174,176,181,182,183,185,187,195,196,197,199,201,209,210,211,213,215,221,222,223,225,227,240,241,242,244,246,252,253,254,256,258,276,277,278,280,282,288,289,290,292,294,312,313,314,316,318,334,335,336,338,340,348,349,350,352,354,358,359,360,362,364,377,378,379,381,383,391,392,393,395,397,410,411,412,416,418,424,425,426,428,430,446,447,448,450,452,468,469,470,472,474,480,481,482,484,486,494,495,496,498,500,518,519,520,524,526,532,533,534,536,538,549,550,551,553,555,563,564,565,567,569,575,576,577,579,581,587,588,589,591,593,603,604,605,607,609,622,623,624,626,628,636,637,638,640,642,650,651,652,654,656,672,673,674,676,678,679,680,681,683,685,693,694 
TYPE
 the manner in which RMSD alignment is performed
=DRMSD 
NOPBC
 ignore the periodic boundary conditions
 
SEGMENT1
this is the lists of atoms in the segment that are being considered
=0,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 
SEGMENT2
this is the lists of atoms in the segment that are being considered
=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 
SEGMENT3
this is the lists of atoms in the segment that are being considered
=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 
SEGMENT4
this is the lists of atoms in the segment that are being considered
=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 
SEGMENT5
this is the lists of atoms in the segment that are being considered
=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 
STRUCTURE1
the reference structure
=0.733,0.519,5.298,1.763,0.81,4.301,3.166,0.543,4.881,1.527,-0.045,3.053,1.646,0.436,1.928,1.18,-1.312,3.254,0.924,-2.203,2.126,0.65,-3.626,2.626,-0.239,-1.711,1.261,-0.19,-1.815,0.032,-1.28,-1.172,1.891,-2.416,-0.661,1.127,-3.548,-0.217,2.056,-1.964,0.529,0.276,-2.364,0.659,-0.88,-1.13,1.391,0.856,-0.62,2.565,0.148,0.228,3.439,1.077,0.231,2.129,-1.032,0.179,2.733,-2.099,1.028,1.084,-0.833,1.872,0.593,-1.919,2.85,-0.462,-1.397,1.02,0.02,-3.049,1.317,0.227,-4.224,-0.051,-0.684,-2.696,-0.927,-1.261,-3.713,-1.933,-2.219,-3.074,-1.663,-0.171,-4.475,-1.916,-0.296,-5.673     # Action input conctinues with 35 further SEGMENTn keywords, 
cv5_ltThe LESS_THAN action with label cv5_lt calculates the following quantities:
 Quantity   
 Type   
 Description  
cv5_lt
vector
the vector obtained by doing an element-wise application of a function that is one if the input is less than a threshold to the input vectors: 
LESS_THAN
Use a switching function to determine how many of the input variables are less than a certain cutoff. More details
 
ARG
the values input to this function
=cv5_rmsd 
SWITCH
This keyword is used if you want to employ an alternative to the continuous swiching function defined above
={RATIONAL R_0=0.08 D_0=0.0 NN=2 MM=4}
cv5The SUM action with label cv5 calculates the following quantities:
 Quantity   
 Type   
 Description  
cv5
scalar
the sum of all the elements in the input vector: 
SUM
Calculate the sum of the arguments More details
 
ARG
the values input to this function
=cv5_lt 
PERIODIC
if the output of your function is periodic then you should specify the periodicity of the function
=NO
# --- End of included input --- 
lwallThe LOWER_WALLS action with label lwall calculates the following quantities:
 Quantity   
 Type   
 Description  
lwall.bias
scalar
the instantaneous value of the bias potential
lwall.force2
scalar
the instantaneous value of the squared force due to this bias potential: 
LOWER_WALLS
Defines a wall for the value of one or more collective variables, More details
 
ARG
the arguments on which the bias is acting
=cv5 
AT
the positions of the wall
=2.0 
KAPPA
the force constant for the wall
=30.0
uwallThe UPPER_WALLS action with label uwall calculates the following quantities:
 Quantity   
 Type   
 Description  
uwall.bias
scalar
the instantaneous value of the bias potential
uwall.force2
scalar
the instantaneous value of the squared force due to this bias potential: 
UPPER_WALLS
Defines a wall for the value of one or more collective variables, More details
 
ARG
the arguments on which the bias is acting
=cv5 
AT
the positions of the wall
=15.0 
KAPPA
the force constant for the wall
=30.0

metadThe METAD action with label metad calculates the following quantities:
 Quantity   
 Type   
 Description  
metad.bias
scalar
the instantaneous value of the bias potential: 
METAD
Used to performed metadynamics on one or more collective variables. More details
 
ARG
the labels of the scalars on which the bias will act
=cv5 
PACE
the frequency for hill addition
=2500 
HEIGHT
the heights of the Gaussian hills
=0.3 
SIGMA
the widths of the Gaussian hills
=0.1 
FILE
 a file in which the list of added hills is stored
=HILLS 
GRID_MIN
the lower bounds for the grid
=0.0 
GRID_MAX
the upper bounds for the grid
=30.0 
GRID_SPACING
the approximate grid spacing (to be used as an alternative or together with GRID_BIN)
=0.025

PRINT
Print quantities to a file. More details
 
STRIDE
 the frequency with which the quantities of interest should be output
=10 
ARG
the labels of the values that you would like to print to the file
=cv5,metad.bias,uwall.bias,lwall.bias 
FILE
the name of the file on which to output these quantities
=COLVAR
#PRINT STRIDE=10 ARG=cv5,metad.bias FILE=COLVAR