RE: Problems to parametrize molecules

From: Bennion, Brian (bennion1_at_llnl.gov)
Date: Wed Apr 15 2015 - 13:04:21 CDT

Hello,

can you state how you used the molefacture plugin. I need to know what steps you are applying to solve this problem, EXACTLY.

Molefacture will not automatically assign partial charges and in cases where the atom type is larger than 4 characters it will cut off the rest.

It APPEARS that you are reading in the pdb file and molefacture is using that atom type.

Please list the steps that you are following and the results you get from those steps and then what you THINK those steps should have done.

Brian

________________________________
From: Evandro Semighini [epsemighini_at_gmail.com]
Sent: Wednesday, April 15, 2015 9:00 AM
To: Bennion, Brian
Subject: Re: namd-l: Problems to parametrize molecules

I got the psf from charmm_gui, but is seems that it's not working with my molecules.

I tried now with another molecule and it doesn't recognize the bond parameters between two atoms that are bonded and listed at the BOND list of psf file.

There is another way to parametrize this molecules ?

Molefacture doesn't put the partial charges and misses the correct atom types (instead of CG311 it calls C3)

Thank you

2015-04-14 17:57 GMT-03:00 Bennion, Brian <bennion1_at_llnl.gov<mailto:bennion1_at_llnl.gov>>:
still guessing but the psf file is not correct and that is what counts. I am not sure how you got your psf file from the pdb or mol2 but molefacture might be a plugin to use just to see the difference in the psf file that it creates compared to the current psf file.

Brian

________________________________
From: Evandro Semighini [epsemighini_at_gmail.com<mailto:epsemighini_at_gmail.com>]
Sent: Tuesday, April 14, 2015 1:51 PM
To: namd-l_at_ks.uiuc.edu<mailto:namd-l_at_ks.uiuc.edu>; Bennion, Brian
Subject: Re: namd-l: Problems to parametrize molecules

Thank you Brian,

I thought that this as the problem too, but I don't know how to fix it, since the input files (pdb or mol2) have the correct bond orders.

I can change this with molefracture, before try the psfgen or the sites I tried ?

2015-04-14 17:45 GMT-03:00 Bennion, Brian <bennion1_at_llnl.gov<mailto:bennion1_at_llnl.gov>>:
I am just guessing here...

You need to pay attention to the what the error states. The errors suggest that the same atom is bonded to itself and it obviously cant find that in the parameter file.

________________________________
From: owner-namd-l_at_ks.uiuc.edu<mailto:owner-namd-l_at_ks.uiuc.edu> [owner-namd-l_at_ks.uiuc.edu<mailto:owner-namd-l_at_ks.uiuc.edu>] on behalf of Evandro Semighini [epsemighini_at_gmail.com<mailto:epsemighini_at_gmail.com>]
Sent: Tuesday, April 14, 2015 1:30 PM
To: namd-l_at_ks.uiuc.edu<mailto:namd-l_at_ks.uiuc.edu>
Subject: namd-l: Problems to parametrize molecules

Hello NAMD community

I am trying to run NAMD with the results of virtual screening assays I made, but I am not being able to parametrize the molecules properly.

I read everything I found about it but, unfortunately, it didn't worked.
I tried charmmgui and paramchem webservers to generate the psf files, but when I try to run NAMD, theres always some unknown atom type (this errors are from two different molecules)

FATAL ERROR: UNABLE TO FIND BOND PARAMETERS FOR CG301 CG301 (ATOMS 22 22)

FATAL ERROR: UNABLE TO FIND BOND PARAMETERS FOR NG2S1 NG2S1 (ATOMS 8 8)

This atom types are at the par_all36_cgenff.prm and I added all the parameter files from the package to run this simulations.
There were no high penalties for the atoms from the servers.

Can you please help me to find where I am missing or show me where I should go to parametrize this molecules for the simulations ?

Thank you,

Evandro Semighini
University of São Paulo - Brazil

This archive was generated by hypermail 2.1.6 : Thu Dec 31 2015 - 23:21:48 CST