I am performing energy minimization on my molecule.
I obtained the parameters of the molecule from CGenFF server.
Since the molecule is protonated, I added chloride ion to make the system neutral.
Following is the command which I have used: gmx genion -s SOLVATED.tpr -nname CL -nn 1 -neutral -p topol-test.top -o NEUTRALIZED.gro
I am getting same problem. The complex.gro file is fine, but after energy minimizations there are some different bonds forming same as yours. The no of atoms are same
This setting is likely inappropriate. It is only for infinite species like graphene, CNT, etc.
Your cutoffs are also wrong for the CHARMM (or CGenFF) force field, but this is not the cause of distortion.
Presumably you have converted the coordinates and topology to GROMACS format with the conversion script we supply via Alex MacKerell’s website? Do the initial coordinates look correct, prior to minimization? Did grompp issue any warnings about mismatched atom names, which implies an incompatibility between coordinates and topology?
WARNING 1 [file topol-test.top, line 388]:
You are using Ewald electrostatics in a system with net charge. This can
lead to severe artifacts, such as ions moving into regions with low
dielectric, due to the uniform background charge. We suggest to
neutralize your system with counter ions, possibly in combination with a
physiological salt concentration.
The initial coordinates were fine before the energy minimization.
The collapse of the ring and the appearance of an H atom in the physical location where a heavy atom should be still suggests an issue with the topology. The partial net charge is a little troubling but could be simply due to rounding error and may not be a concern. It’s certainly not the cause of such distortion in the connectivity.
It would be easier to provide helpful advice if either (or both) of you would upload all input files for grompp (coordinates, all relevant topologies, .mdp file, and any necessary index file) for others to investigate.
Prof. Lemkul,
I am pasting a link containing my files. I have used the same steps as provided in your tutorial, except I have used SwissParam for ligand topology generation. I have pasted my ligand.mol2 as well as itp files.If you require any other files for debugging kindly let me know.
I think Justin probably already answered this correctly a few posts back. It seems likely that your topology either has atoms mixed up intrinsically (hopefully an unlikely scenario), or the order in which atoms appear in the topology is not the same as in your coordinates. If you gave the CGenFF server a molecule and got parameters, did you also get a structure? Is it possible that CGenFF gave you back atoms in a different order? If so, suggest to use the CGenFF output structure instead of the CGenFF input structure in your EM input.
The mol2 file of ligand was generated, Swissparam was used for generating the pdb and itp file.
The topol top files was created like this
;
; File ‘topol.top’ was generated
; By user: unknown (1000)
; On host: LAPTOP-FS0RKSQQ
;
;
; This is a standalone topology file
;
; Created by:
; :-) GROMACS - gmx pdb2gmx, 2020.1-Ubuntu-2020.1-1 (-:
;
; Executable: /usr/bin/gmx
; Data prefix: /usr
; Working dir: /mnt/c/Users/Savio Cardoza/Desktop/GROMACS systems/7lif_1_sp_b
; Command line:
; gmx pdb2gmx -f protein.pdb -o protein_processed.gro -ignh
; Force field was read from current directory or a relative path - path added.
;
; Include forcefield parameters #include “./charmm36-mar2019.ff/forcefield.itp”
; Include ligand topology #include “lig.itp”
; Include water topology #include “./charmm36-mar2019.ff/spc.itp”
#ifdef POSRES_WATER
; Position restraint for each water oxygen
[ position_restraints ]
; i funct fcx fcy fcz
1 1 1000 1000 1000 #endif
; Include topology for ions #include “./charmm36-mar2019.ff/ions.itp”
[ system ]
; Name
lig in water
[ molecules ]
; Compound #mols
lig 1
SOL 797
gmx pdb2gmx was used to create a .gro file for ligand.
solvent and ions were added (however no ions were added by gromacs)
gmx grompp -f em.mdp -c solv_ions.gro -p topol.top -o em.tpr AND
gmx mdrun -v -deffnm em WERE RUN
I got same distortion as earlier.
7.CHANGED emtol = 5000.0, in the em.mdp file,. I again ran energy minimization and although the energy is not converged but after 4 steps i saw the em.gro file and structure did not distort.
Is the steep descent algorithm incompatible for the molecule? and can I restrain the molecule? How can we detect issues in the topology of our ligand as Prof. Lemkul pointed out?
@abhinav.sri89 the issue is your small-molecule topology. It defines bonds and pairs but no angles or dihedrals. The entire PYP molecule collapses in on itself because the topology is simply wrong. It appears it was manually created; don’t do this. Use the cgenff_charmm2gmx.py script from Alex MacKerell’s website (as in my protein-ligand tutorial) and it will generate an appropriate topology. You also have #include statements within the [system] directive, which I have no idea how it works, but you should generate a more sensible and syntactically correct topology. This one is essentially unusable and you should start over.
@SavioC you did not provide all the necessary files.
-------------------------------------------------------
Program: gmx grompp, version 2021.3
Source file: src/gromacs/gmxpreprocess/gmxcpp.cpp (line 298)
Fatal error:
Topology include file "topol_Protein_chain_A.itp" not found
For more information and tips for troubleshooting, please check the GROMACS
website at http://www.gromacs.org/Documentation/Errors
-------------------------------------------------------
