RESP version 3.1
The novelties brought by version 3.1 are available from the readme.txt file.
RESP version 2.4
This documentation is based on the Amber8 & Amber9 manuals. Additional pieces of information have been added to update the former documentation version.
The novelties brought by version 2.4 are highlighted in yellow.
-7th area- intra-molecular charge constraint(s)
-8th area- inter-molecular charge constraint(s)
-9th area- inter-structure charge equivalencing
-2nd line up to n_atoms+1 line-
-n_atoms+2 line up to n_atoms+2+n_MEP_points line-
Description of some Fortran formats
Download the RESP program & patches
Usage:
resp [-O] -i input -e espot [-o output -p punch -q qin -t qout -w qwts -s esout -j espdb -y esqpotpdb -z esmpotpdb]
-O Overwrite output files if they exist.
The RESP (Restrained ElectroStatic Potential) program fits the quantum mechanically calculated molecular electrostatic potential (MEP) at molecular surfaces using an atom-centered point charge model. This method was developed primarily by Bayly. [1-3] A quantum mechanical program such as Gaussian, [4] GAMESS, [5] or Firefly [6] must be used to generate the MEP input(s) for RESP.
- Note that the Antechamber program automates most of this process: [7] use the -fo gcrt option to create a Gaussian input file; then run Gaussian; then use the -fi gout -c resp option to automatically create the RESP input file and run a two-stage fitting procedure. If you don't use Gaussian, you can still run respgen to automatically create the input files needed for RESP.
- Other alternatives are to use the RESP ESP charge Derive (R.E.D.) tools, [8,9], R.E.D. Server/R.E.D. Server Development [10] and PyRED [11] to create RESP inputs, generate force field libraries and to manage the calculations. Examples of RESP and ESP charge derivation and force field libraries are available in the RESP ESP charge DDataBase (R.E.DD.B.). [12] See q4md-forcefieldtools.org for information about these options.
file name | flag | fortran unit |
| purpose |
-------------------------------------------------------------------------------------------------------
input | -i | 5 | required | input options |
output | -o | 6 | always produced | output of results |
punch | -p | 7 | always produced | synopsis of results |
qin | -q | 3 | optional | replacement charges |
qout | -t | 19 | always produced | output of current charges |
espot | -e | 10 | required | input of MEP and coordinates |
|
|
|
| (note: these must be in atomic units) |
qwts | -w | 4 | optional | input of new weight factors |
esout | -s | 20 | optional | generated MEP values for new charges |
espdb | -j | 30 | always produced | PDB-like file with relative residual in the TempFactor field |
esqpotpdb | -y | 40 | always produced | PDB-like file with input MEP values in the TempFactor field (atomic units) |
esmpotpdb | -z | 50 | always produced | PDB-like file with MEP values for new charges in TempFactor field (atomic units) |
-------------------------------------------------------------------------------------------------------
Input included in the "-i" file (see examples at the end)
TITLE | a character string |
Begin with namelist " &cntrl"
inopt | = | 0 | normal run |
| = | 1 | cycle through a list of different "qwt" read from -w unit |
ioutopt | = | 0 | normal run |
| = | 1 | write restart info of new MEP to -s unit |
iqopt | = | 1 | reset all initial charges to zero (default) |
| = | 2 | read in new initial charges from -q unit |
| = | 3 | read in new initial charges from -q unit and perform averaging of those new initial charges according to "ivary" values (normally not used) |
nmep | = | the number of Molecular Electrostatic Potential(s) (MEP) in a multiple MEP fit (default 1) MEP(s) = orientation(s), conformation(s) or molecule(s) |
ihfree | = | 0 | all atoms are restrained |
| = | 1 | hydrogens not restrained (default) |
irstrnt | = | 0 | use of harmonic restraints (old style) |
| = | 1 | use of hyperbolic restraints (default) |
| = | 2 | only analysis of input charges; no charge fitting is carried out |
| = | 3 | use of quadratic restraints (Henchman & Essex[13]) |
qwt | = | normally use 0.0005 for Stage 1 - Amber force fields (default) use 0.01 for Stage 1 - Glycam force field |
| = | normally use 0.001 for Stage 2 - Amber force fields (default) no Stage 2 - Glycam force field |
End namelist " &cntrl" with " &end"
wtmol | relative weight for the structure if multiple structure fit (1.0 otherwise) |
subtitle for the structure (a character string)
charge | iuniq | total charge value & total number of atoms for the considered structure (in 2I5 format) |
|
one line for each atom:
Atomic number | ivary (in 2I5 format) |
ivary |
|
|
|
| = | 0 | current charge fitted independently of other centers |
| = | n | current charge fitted and equivalenced to that of center "n" |
| = | -1 | current charge frozen at "initial stage" value typically read in from -q unit |
*blank to end whatever the number of structure(s) "nmep" is
intra-molecular charge constraint(s) (in I5,F10.5 format) blank line if no constraint |
ngrp | = | the number of centers in the group of atoms associated with this constraint |
grpchg | = | charge value to which the associated group of atoms (given on the next card) is to be constrained |
-7.1th area-
imol | iatom | structure & atom indices (in 16I5 format) |
the list ("ngrp" long) of the atom indices of those atoms to be constrained to the charge specified on the previous line.
*blank to end
inter-molecular charge constraint(s)
same format as intra-molecular charge constraint(s) - see the 7th & 7.1th areas
*blank to end
multiple structure atom charge equivalencing
format is similar to 7th & 7.1 areas
ngrp (in I5 format)
and then, on separate lines:
imol | iatom (in 16I5 format) |
*blank to end
input of replacement charges if requested (in 8F10.6 format) (note same format as that produced by -t unit) |
input new weight factors if requested |
input | qwt | number of new weights to cycle thru (in I5 format) |
input | new weights ("nqwt" lines in F10.5 format) |
input of MEP and coordinates |
n_atoms n_MEP_points (total number of atoms & MEP points, in I5,I6 format)
-2nd line up to n_atoms+1 line-
atom coordinates X Y Z (in Bohrs, in 17X,3E16.7 format)
-n_atoms+2 line up to n_atoms+2+n_MEP_points line-
MEP & coordinates
qpot X Y Y (in a.u. & Bohrs, in 1X,4E16.7 format)
Description of some Fortran formats
2I5 format 2 blocks of 5 characters (I = integer)
16I5 format 16 blocks of 5 characters (I = integer)
F10.5 format A block of 10 characters long (F = float) with 5 digits after the decimal point
17X,3E16.7 format A block of 17 space characters long followed by 3 blocks of 16 characters long with 7 digits after the decimal point (using the E scientific notation)
1X,4E16.7 format A space character followed by 4 blocks of 16 characters long with 7 digits after the decimal point (using the scientific E notation)
Several examples of input and output files are available in $AMBERHOME/examples/resp_charge_fit; these should be consulted by those interested in running the program.
See the RESP ESP charge DDataBase as well. In each R.E.DD.B. project, the RESP inputs used in the reported charge fitting step are provided. See: http://q4md-forcefieldtools.org/REDDB/projects/REDDB-code/input$X.in (where "REDDB-code" and "$X" are the R.E.DD.B. code and the number of stages of the charge fitting step, respectively). Examples:
- Single orientation/single conformation/single molecule RESP fit (Amber force fields):
W-1 R.E.DD.B. project: input1.in & input2
- Multiple orientation/single conformation/multiple molecule RESP fit (Amber force fields):
W-46 R.E.DD.B. project: input1.in & input2.in
- Use of intra-molecular charge constraints in RESP charge fitting (Amber force fields):
F-1 R.E.DD.B. project: input1.in & input2.in
- Use of inter-molecular charge constraints & inter-molecular atom charge equivalencing in RESP charge fitting (Amber force fields):
F-60 R.E.DD.B. project: input1.in & input2.in
- Use of intra-molecular charge constraints, inter-molecular charge constraints & inter-molecular atom charge equivalencing in RESP charge fitting (Glycam force field):
F-71 R.E.DD.B. project: input1.in & input2.in
- Use of intra-molecular charge constraints, inter-molecular charge constraints & inter-molecular atom charge equivalencing in RESP charge fitting (Amber force fields):
F-74 R.E.DD.B. project: input1.in & input2.in
- Use of intra-molecular charge constraints, inter-molecular charge constraints & inter-molecular atom charge equivalencing in RESP charge fitting (Amber force fields):
F-85 R.E.DD.B. project: input1.in & input2.in
etc...
[1] C.I. Bayly, P. Cieplak, W.D. Cornell & P.A. Kollman. A Well-Behaved Electrostatic Potential Based Method Using Charge Restraints For Determining Atom-Centered Charges: The RESP Model. J. Phys. Chem. 1993, 97, 10269-10280.
[2] W.D. Cornell, P. Cieplak, C.I. Bayly & P.A. Kollman. Application of RESP charges to calculate conformational energies, hydrogen bond energies and free energies of solvation. J. Am. Chem. Soc. 1993, 115, 9620-9631.
[3] P. Cieplak, W.D. Cornell, C. Bayly & P.A. Kollman. Application of the multimolecule and multiconformational RESP methodology to biopolymers: Charge derivation for DNA, RNA and proteins. J. Comput. Chem. 1995, 16, 1357-1377.
[4] Gaussian.
[5] GAMESS.
[6] PC-GAMESS/Firefly.
[7] J. Wang, W. Wang, P. A. Kollman & D. A. Case, Automatic atom type and bond type perception in molecular mechanical calculations, J. Mol. Graph. Model. 2006, 25, 247-260, Antechamber.
[8] A. Pigache, P. Cieplak & F.-Y. Dupradeau, Automatic and highly reproducible RESP and ESP charge derivation: Application to the development of programs RED and X RED, 227th ACS National Meeting, Anaheim, CA, USA, March 28 - April 1, 2004, RED.
[9] F.-Y. Dupradeau, A. Pigache, T. Zaffran, C. Savineau, R. Lelong, N. Grivel, D. Lelong, W. Rosanski & P. Cieplak, The R.E.D. tools: Advances in RESP and ESP charge derivation and force field library building, Phys. Chem. Chem. Phys. 2010, 12, 7821-7839, [PMCID].
[10] E. Vanquelef, S. Simon, G. Marquant, E. Garcia, G. Klimerak, J. C. Delepine, P. Cieplak and F.-Y. Dupradeau, R.E.D. Server: a web service for deriving RESP and ESP charges and building force field libraries for new molecules and molecular fragment, Nucl. Acids Res. (Web server issue) 2011, W511-W517, RED Server & [PMCID].
[11] F. Wang, J.-P. Becker, P. Cieplak & F.-Y. Dupradeau, R.E.D. Python: Object oriented programming for Amber force fields, 247th ACS National Meeting, Dallas, TX, USA, March 18-20, 2014, RED Server Development.
[12] F.-Y. Dupradeau, C. Cezard, R. Lelong, E. Stanislawiak, J. Pecher, J. C. Delepine & P. Cieplak, R.E.DD.B.: A database for RESP and ESP atomic charges, and force field libraries, Nucl. Acids Res. (Database issue) 2008, D360-D367, REDDB & [PMCID].
[13] R.H. Henchman & J.W. Essex. Generation of OPLS-like Charges From Molecular Electrostatic Potential Using Restraints. J. Comp. Chem. 1999, 20, 483-498.
Download the RESP program & patches
. Download the standalone version of the RESP program version 2.1 from q4md-forfieldtools.org: the readme.txt file & package.
. Download a patch to increase the number of MEP points defined in the espot file (resp 2.1).
. Download the standalone version of the RESP program version 2.2 from q4md-forfieldtools.org: the readme.txt file & package.
. Download the standalone version of the RESP program version 2.4 from q4md-forfieldtools.org: the readme.txt file & package.
. Download the standalone version of the RESP program version 3.1 from q4md-forfieldtools.org: the readme.txt file & package.
Last update of this RESP documentation: November 30th, 2014