Re: [gmx-users] how to modify opls force field
fit_tone at 163.com
Mon Apr 18 04:04:48 CEST 2005
I have also modified opsl force field as the way which you mention.I think what you have done is OK! As for your problem, I think it is caused by the modified residue "MSO" is not included in the file "/share/top/aminoacids.dat",and the programe does not know it is a amino acid. you may add the name "MSO" to the file and try again.good luck!
> Dear Gromacs users,
> I tried to use opsl force field to simulation a pentapeptide solution.
> The residue in my pentapeptide is very similar to Methionine. The only
> difference is the atom between CG and CE is O, not S (in Methionine).
> Because this residue is not included in ffoplsaa.rtp, I made a new .rtp
> file by adding a new residue. I named it MSO. The atom type for O in MSO
> is opsl_999 (I treated it as opsl_180, O in the dialkyl ether). I also
> modified ffoplsaa.atp, ffoplsaabon.itp, ffoplsaa.hdb, ffoplsaa.itp and
> ffoplsaanb.itp by adding new parameters. In short, I made my own force
> field for my pentapeptide, and I gave it a new name. However, when I
> tried to use pdb2gmx to convert my pdb file, I got the error message as
> No N- or C-terminus found: this chain appears to contain no protein
> Fatal error: Atom -C not found in residue MSO1 while adding hydrogens
> Because I just added more parameters to opsl force field to get my own
> force field, I think my new force field should work for Methionine
> pentapeptide, so I changeed my pdb file to Methionine pentapeptide. When
> I tried pdb2gmx, it worked fine.
> Another test I have done is I modified ffoplsaa.rtp by changing the
> definition of MET, namely I changed SD to OD in [MET] section (and
> changed other things accordingly). I also changed SD to OD in MET
> section in the file ffoplsaa.hdb. Then I tried pdb2gmx, it worked as
> Could anyone tell me what's wrong with the modification by adding one
> more residue type? Any helps are highly appreciated. If you need more
> information, please let me know.
> Thanks in advance!
> Dongsheng Zhang
> gmx-users mailing list
> gmx-users at gromacs.org
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-users-request at gromacs.org.
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the gromacs.org_gmx-users