From: Mayne, Christopher G (
Date: Thu Sep 17 2015 - 16:21:12 CDT


I'm not quite sure what's going on here. It seems like it's looking for an atom type that is a concatenation of several atoms, and I'm not sure why it would do that since you're PSF seems OK. I can come back to this problem, but I see a couple of other issues here.

Based on the contents of your initial PAR file I can tell that you didn't include the par_all36_prot.prm file as an "associated parameter file". Many of the parameters that you have in your PAR file are defined by the standard CHARMM force field for proteins. Including this file when generating the initial parameter file will significantly reduce the number of missing parameters that you will have to work on.

You should check out the CGenFF topology file, specifically how residue GUA changes to MGUA. Comparing these two entries will be helpful in assessing the changes associated with methylating the nitrogens of guanidinium, e.g., how the partial atomic charges change, and it will give you proper LJ parameters for the terminal methyl. Many of the missing parameters should be defined, albeit in CGenFF terms (i.e., using CGenFF atom types rather than CHARMM protein atom types). It may be a little tedious, but you should be able to translate these.

Anything that remains missing won't required parameterizing the full molecule that you've constructed. At worst you can truncate the system to an N-methyl and N-propyl guanidinium.

Christopher Mayne

On Sep 17, 2015, at 2:39 PM, Amy Rice wrote:

Hi all,
I am attempting to parameterize a modified arginine residue following the ffTK screencasts, using VMD 1.9.2. I've been able to complete the first 4 steps without issue, but I'm running into some errors during the charge optimization, in both downhill mode and simulated annealing. The error message reads:
"can't read "ljPars(HA2 HA2 HA2 H)": no such element in array"

I noticed a similar problem was discussed previously (>), however as far as I can tell the atom types HA2 and H exist in both the psf and parm files. I have LJ parameters in my .par file for both of these atom types, so I'm not sure why I'm receiving an error. I've attached my preliminary psf, pdb, and par for reference. Note that I am also loading the charmm 36 protein prm as an input paramter file.

Thank you for the help,
- Amy

Amy Rice
Ph.D. Student
Physics Department
Illinois Institute of Technology