---------------------------- Original Message ---------------------------- Subject: md with multiple ligands From: ansuman@localhost Date: Mon, February 18, 2013 8:28 pm To: gmx-users@gromacs.org --------------------------------------------------------------------------
Message: 8 >> Date: Wed, 13 Feb 2013 08:28:00 -0500 >> From: Justin Lemkul <[hidden email]> >> Subject: Re: [gmx-users] md with multiple ligands >> To: Discussion list for GROMACS users <[hidden email]> >> Message-ID: <[hidden email]> >> Content-Type: text/plain; charset=ISO-8859-1; format=flowed >> >> >> >> On 2/13/13 5:24 AM, Ansuman Biswas wrote: >>> >>> Dear gromacs users, >>> >>> I wish to run a protein ligand md in Gromacs 4.5.3 with charmm force >>> field. Its a dimeric protein with a MG ion and 2 ligands bound to two >>> active sites. I wish to run the dimer and there would be 4 ligands. I >>> generated the topology files (.itp) of ligands using SWISSPARAM. Before >>> running the energy minimization I changed the topol.top file accordingly >>> as mentioned in the SWISSPARAM website. I added all the itp files >>> together.But while running GROMPP I got an error message like, >>> >>> >>> WARNING 1 [file ADP.itp, line 24]: >>> Overriding atomtype NPYL >>> >>> >>> WARNING 2 [file ADP.itp, line 25]: >>> Overriding atomtype C5A >>> >>> >>> WARNING 3 [file ADP.itp, line 26]: >>> Overriding atomtype N5B >>> >>> >>> WARNING 4 [file ADP.itp, line 27]: >>> Overriding atomtype C5B >>> >>> >>> Generated 25425 of the 25425 non-bonded parameter combinations >>> Generating 1-4 interactions: fudge = 1 >>> Generated 22285 of the 25425 1-4 parameter combinations >>> >>> ------------------------------------------------------- >>> Program grompp, VERSION 4.5.3 >>> Source code file: toppush.c, line: 1071 >>> >>> Fatal error: >>> Atoms in the .top are not numbered consecutively from 1 (rather, atomnr >>> = >>> 1, while at->nr = 40) >>> >>> I know that the problem is related to the merging of 4 itp files. >>> Please let me how I should do that. >>> >> >> The ligands probably use common atom types and each .itp file likely has >> an >> [atomtypes] directive that introduces these new types. Probably the >> easiest >> thing to do is create a merged [atomtypes] directive that lists all of the >> necessary atom types so you don't list them separately in each file. One >> approach that I have used in the past is to create an .itp file that only >> has >> atom types in it, then call something like: >> >> #include "charmm27.ff/forcefield.itp" >> >> ; only has [atomtypes] >> #include "my_ligand_atomtypes.itp" >> >> ; these don't introduce [atomtypes] any more >> #include "ligand1.itp" >> #include "ligand2.itp" >> >> -Justin > But according to SWISSPARAM website I can only add one ligand.itp file. « [hide part of quote] Just a guess, but it's probably because each ligand.itp file introduces new atom types which give the conflicts above. > So, is it possible to add ligand1 and ligand2 itp files? > Sure, my method should work. Please try it. I have done the exact same thing with Amber topologies from acpype. > Will it be ok if I merge all the ligand pdbs and then get the combined itp > file using that from SWISSPARAM and follow the protocol for one ligand? > Merging the coordinate files will probably not work. If you have multiple molecules, you need multiple topologies. Justin Well, I have tried using the method you mentioned; i.e using 2 separate itp files and use one common atomtype at the beginning. But it seemed that gromacs was unable to read the next lig1 and lig2 itp files without atomtypes option. It gave error message. So I tried merging two itp files from swissparam so that only unique parameters remain in the combined.itp file. I could run Energy Minimization only with steepest descent but not with CG as CG could not reduce the maximum force beyond 5th order. After that when I tried running NVT, I found Lincs error and saw several bonds rotated beyond 30 deg, as well as domain decomposition error: constraint lengths are too long compared to the domain decomposition cell size *******More on that when I opened the em.pdb (em output) in pymol I found all the ligand atoms were jumbled up together. please help, Ansuman -- This message has been scanned for viruses and dangerous content by MailScanner, and is believed to be clean. -- gmx-users mailing list gmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting! * Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists