-
The number of atoms (or united-atoms) in this molecule. Must be less than or equal to NUMAX.
-
The maximum number of atoms to regrow during a configurational-bias regrowth. It is safe to set this to
the same value as nunit, but some molecules are so large that almost none of the moves which regrow large portions are accepted.
In those cases setting a lower value of nmaxcbmc focuses on moves that are likely to be accepted.
.TRUE.enable the input of information about the pdb (protein data bank) atom name, amino acid number, amino acid name. This information is needed in order to use the cartoon feature of certain pdb viewers (such as Rasmol)..FALSE.disable the input of additional pdb information. This does not disable the output of pdb files, it just means thing like the Rasmol cartoon feature will not work properly.
-
The Towhee shorthand name for the forcefield that you want to use to build this molecule.
See Towhee Capabilities for a complete list of the forcefields included with the Towhee distribution.
- 'bond increment': the bond increment method is used to assign the charges on each atom and the parameter charge is not required. Note that this method is not functional for all forcefields.
- 'manual': the charges on each atom are manually specified in the towhee_input file as the parameter charge in the following section.
- 'none': a zero charge is assigned to each atom and the parameter charge is not required.
- '3-bond: single improper': atoms that are bonded to exactly three neighbors are assigned a single improper torsions automatically. The three neighbors are automatically assigned in increasing order of their nonbonded integer type number. Originally implemented for use with the MMFF94 force field. The improper torsion parameters are not required in the subsequent section.
- 'manual': the improper torsions each atom are manually specified in the subsequent section.
- 'none': no improper torsions are assigned for this molecule and the improper torsion parameters are not required in the subsequent section.
- 'standard' : determines matching interactions for a given forcefield by checking for exact type matches first, and then utilizing a procedure that checks for an increasing number of wildcards. This is the approach recommended for the majority of forcefields and is the default value used by most of the other molecule assembly options.
- 'stereoselective' : very similar to the 'standard' option, but also requires the input of a list of torsions that have a specific stereochemistry. This allows the process to distinguish between cis and trans dihedrals when combined with a forcefield that makes such a distinction. When using this option the special_torsion_count parameters are required after the usual input of the connectivity.
- 'MMFF checkdown' : determines matching interactions for a given forcefield using the procedure described in the MMFF94 forcefield papers. This uses an internal checkdown procedure based upon atom names that are implemented as a 10 character string that contains 5 2-character substrings.
This section is repeated for each atom in the molecule (determined by nunit)
-
Three parameters listed on a single line.
First, the number of the atom in order starting from atom number 1.
Then the character string that contains the atom type for this atom
(see the individual forcefield documentation summarized at Towhee Capabilities).
Finally, charge on this atom (only required if the charge_assignment is 'manual').
.TRUE.
-
Three parameters listed on a single line.
First, a four character string that is output in the pdb file. The precise spacing is important if you want most pdb viewers to recognize
the atoms as the pdb file is extremely specific.
Then, the number of each amino acid starting from the N-terminus (only relevant for polypeptides).
Finally, the three letter code for each amino acid, or other group (such as caps on the C or N termini).
- The first line under the vibration heading is the number of atoms that are bonded to the current atom. Must be a number between 0 and NNBOND.
- For each vibration a subsequent line contains the atom number (integer) of the bond partner and the
vibration order (character*10) of the bond. Commonly used vibration orders are listed here.
- 'any': matches any bond order listed in the force field file. This effectively disables the vibration order searching. This entry is only valid in the towhee_input file, and not in an ff*.F file.
- 'aromatic': aromatic bond
- 'single': single bond
- 'double': double bond
- 'triple': triple bond
- 'wild': matches with any type of bond, but only if there is not an exact match. This is only a valid entry in an ff*.F file, and not in the towhee_input file.
- The first line under the improper heading is the number of improper torsions (any form) that have the central atom located at the current atom. You must list all improper torsions that have the current atom at the central position, but you do not list improper torsions that only contain the current atom at one of the ends. Must be a number between 0 and MAXIMPROP.
- For each improper torsion specified above a line with the following information is required. A list of the three other atoms in the improper torsion followed by the appropriate improper type for the force field. The code currently only allows an improper torsion where the three other atoms are all bonded to the first atom. If you want to enable the automatic determination of the improper torsion force field parameters then set the improper force field type to 0.
- The first line under the special_torsion_count heading is the number of special torsion assignments that involve a torsion order other than 'wild'. If this value is not 0, then a subsequent line is required for each value of this parameter.
- For each value of the special_torsion_count specified above, a line with the following information is required.
A list of the integer values of the 4 atoms that make up the special torsion, followed by a character string that describes the special
torsion. The character string is checked against the torsion order specified in the various forcefield files in order to determine a match.
The torsions only need to be specified in a single order, so if there was a special torsion between atoms 1,2,3, and 4 it is only required to
list the 1 2 3 4 'special flag' direction and not the 4 3 2 1 'special flag' direction. The special torsion order character string depends
upon the force field, but typical values are listed here.
- 'wild' : default value if none is specified. This matches up with torsions that do not have a specified torsion order as the default output for the torsion order is 'wild'. There is no reason to go through the effort of listing this value since it is the default, but it would result in a valid match with the default torsions so it is allowed.
- 'cis' : matches up with a torsion that is in the cis conformation only. Typically occurs for double bonded carbons that are then bonded to dissimilar atoms.
- 'trans' : matches up with a torsion that is in the trans conformation only. Typically occurs for double bonded carbons that are then bonded to dissimilar atoms.
Return to the main towhee web page