GROMACS version:
GROMACS modification: Yes/No
How to generate topology file from pdb.
In my case I am using an organic molecule and base to study the simulation. But when I am attempting to generate topology file I am getting error like residue ligand not found in the residue topology database. For example i am trying to simulate Naphthalene i am getting error like this.Can anyone explain how to resolve this issue and advice on how to modify my PDB file to ensure compatibility with GROMACS for the topology generation process.
Have you tried converting your pdb into mol2, modifying bond order and then submitting this mol2 structure in the CGenFF web browser?
No. But i converted pdb to mol2 using swissparam. Still it doesn’t work
Have you edited the header line in the mol2 file?
Like, if we open the MOL2 file using any text editor, the molecule header will look like this, “****”, we have to replace this with our ligand name and make it uniform throughout the MOL2 file.
Have you done that?
No what i am conveying is i am not even able to get a mol2 file. With the mol2 file only i can get itp file. But i can’t able to convert pdb to mol2 file.
If you have openbabel installed in your system, or any other tools like PyRx, Avogadro, you an easily convert PDB to MOL2.
Thank you.I did converted it using obabel. But can you tell me from where i can get the itp and other files because i am going to write topol file manually.
Also you mentioned CGenFF web browser. Is that a liscenced one?
That is a free website.
You can register yourself, and generate the ligand topology using the mol2 file you generated.
But you need to prepare the mol2 file before submitting it to CGenFF.
You can follow this tutorial here ([link].
Thanks i simulated using CGenFF but for some molecule the ligand topology is not generating in CGenFF . I can only able to download str files that shows error like
! “penalty” is the highest penalty score of the associated parameters.
! Penalties lower than 10 indicate the analogy is fair; penalties between 10
! and 50 mean some basic validation is recommended; penalties higher than
! 50 indicate poor analogy and mandate extensive validation/optimization.
how should i correct this error.
readpar debug: read 683 bonds, 2501 angles, 7460 dihedrals, 203 impropers
and 375437 characters worth of comments.
readmol2 warning: renaming non-unique atom H to H1 .
readmol2 warning: renaming non-unique atom H to H2 .
Now processing molecule water.pd …
paranalogy debug: finding analogy for OGTIP3 HGTIP3
paranalogy debug: finding analogy for HGTIP3 OGTIP3 HGTIP3
paranalogy debug: “smart algorithm” found no matching parameters;
switching to “exhaustive algorithm”
paranalogy debug: finding analogy for OGTIP3 HGTIP3
paranalogy debug: finding analogy for OGTIP3 HGTIP3
paranalogy debug: finding analogy for HGTIP3 OGTIP3 HGTIP3
paranalogy debug: “smart algorithm” found no matching parameters;
switching to “exhaustive algorithm”
this is the error file i got
That doesn’t seem like an error to me.
Check your penalty score. Most likely, if you penalty score is between 10-50 it will still work. If it’s very high like beyond 100, try optimizing your ligand once and generate the topology again. That might helps.
I am manually building and optimizing my molecule using Gaussian software. I tried to reoptimize my molecule but still my penalty score is going up to 182. I want to reduce it below 10 so that I will write the topology file using str file. Can you suggest me how should I reduce the score?
Citing CGenFF tools:
it should be noted that the CGenFF penalty scores are based on analogy to known parameters.
And also very important is that
Accordingly, a parameter with a high penalty may yield acceptable conformational energies, while a parameter already in CGenFF (i.e., penalty = 0) may not yield an acceptable energy surface due to, for example, the terminal rings about the associated bond creating a significantly different context than that in which the dihedral parameter was initially optimized.
As such, probably getting better charges with QM methods won’t lower your score, as the comparison will be again done against the same database and it’s not a measure of how good the parameters are in general, but of how similar are they to something already known (and ideally good). So what you can do is i) live with high penalty scores or ii) refine the molecule with the tools I linked or other QM/MM methods. Whatever is your choice, you will then need to benchmark the parametrized molecule to check whether the parameters reproduce the main properties that you need to investigate your scientific problem.
Thank you for your clarification. But it is showing error like this
python2 cgenff_charmm2gmx_py2.py CLA clath_fix.mol2 clath_fix.str charmm36-jul2022.ff
NOTE 1: Code tested with python 2.7.12. Your version: 2.7.18 (default, Dec 9 2024, 18:47:23)
[GCC 11.4.0]
NOTE 2: Please be sure to use the same version of CGenFF in your simulations that was used during parameter generation:
–Version of CGenFF detected in clath_fix.str : 4.6
–Version of CGenFF detected in charmm36-jul2022.ff/forcefield.doc : 4.6
NOTE 3: To avoid duplicated parameters, do NOT select the ‘Include parameters that are already in CGenFF’ option when uploading a molecule into CGenFF.
what FF do i have to use?
The force field i have installed in the latest version i cant able to find any new versions in this?
I am having a problem generating an all-atom topology for a protein with an amino acid residue of 6425 on Charmgui. the error that comes is “Error: float() argument must be a string or a number, not 'NoneType”
It looks like CHARMM-GUI ran into an issue because it’s expecting a number but found something missing instead. This usually happens when there’s a problem with the residue numbering (like having an unusually high number such as 6425), missing atom coordinates a non-standard amino acid or a formatting issue in the PDB file. To fix this you can check if the residue numbering is correct and make sure the PDB file is complete and properly formatted and confirm that the residue exists in the CHARMM force field.
[quote=“MDGRO, post:15, topic:11341, full:true”]
Thank you for your clarification. But it is showing error like this
python2 cgenff_charmm2gmx_py2.py CLA clath_fix.mol2 clath_fix.str charmm36-jul2022.ff
NOTE 1: Code tested with python 2.7.12. Your version: 2.7.18 (default, Dec 9 2024, 18:47:23)
[GCC 11.4.0]
NOTE 2: Please be sure to use the same version of CGenFF in your simulations that was used during parameter generation:
–Version of CGenFF detected in clath_fix.str : 4.6
–Version of CGenFF detected in charmm36-jul2022.ff/forcefield.doc : 4.6
NOTE 3: To avoid duplicated parameters, do NOT select the ‘Include parameters that are already in CGenFF’ option when uploading a molecule into CGenFF.
what FF do i have to use?
The force field i have installed in the latest version i cant able to find any new versions in this?
Any updates in this question?@obZehn@nilanjana88
cla.top (7.6 KB)
This is my topology file but i am facing the above error in this.
Basically what i am trying to do is MD simulation of my molecule which interacts with water at low temperature . For that i am manually editing the topol file(following protein-ligand simulation module in gromacs tutorial). In the place of protein i am taking my manually build optimized molecule and in the place of ligand i am taking water. So basically i am trying to study the interaction between two molecules in the certain temperature.
I followed all the steps but in gmx grompp -f ions.mdp -c solv.gro -p topol.top -o ions.tpr i am facing the above mentioned error. I have attached top file kindly check it.