[gmx-users] Charge on drug
jalemkul at vt.edu
Mon Dec 19 14:01:19 CET 2016
On 12/19/16 5:31 AM, tasneem kausar wrote:
> Thanks for reply
> I have visualized solvate.gro and em.gro. Some bonds that are not present
> in box.gro is seen in em.gro.
> I am giving you the link of the file to visualize.
You have an amide group that distorts, bringing the carbonyl C closer to the
amide proton, which VMD then fictitiously assigns as a C-H bond. The C-N bond
is much too short. This suggests your topology is of insufficient quality for
doing simulations and requires refinement. There is, of course, no actual bond
there, as Mark has pointed out. Visualization software tries to guess where
bonds should be based on interatomic distances. The topology is definitive.
What you see on your screen is not.
> I addition to that I have done energy minimization multiple times with
> steepest decent and conjugate gradient method. Potential energy value is
> negative in every minimization. When I am trying position restrain core
> dump with multiple lincs warning occurs.
Your topology needs work. It will not be stable for any simulation.
> Please tell me how this problem would be resolved
> On Mon, Dec 19, 2016 at 1:52 PM, Mark Abraham <mark.j.abraham at gmail.com>
>> Your energy minimization literally cannot create bonds because those are
>> already defined by your topology. See
>> for what's probably going on.
>> On Mon, 19 Dec 2016 19:05 tasneem kausar <tasneemkausar12 at gmail.com>
>>> Dear All
>>> I am using Acpype server to generate topology file of drug molecule to
>>> perform protein drug MD simulation. Acpype have generated the itp file
>>> other necessary for MD simulations. I am doing MD on gromacs 5.1.4
>>> using amber99sb force field. When I have done energy minimization some
>>> unwanted bonds are formed in drug molecule. So the position restraint is
>>> not able to execute. The problem is due to the unusual bond in drug atom.
>>> When I used the itp file of PRODRG server position restraint goes
>>> successfully. I have used the gromos54a7 ff. But there is problem with
>>> charge present on the atom that comes from PRODRG.
>>> Please tell me how charge on the atom can be corrected in itp file of
>>> PRODRG. How the problem in energy minimization with the Acpype itp will
>>> Thanks in Advance
>>> With Regards
>>> Gromacs Users mailing list
>>> * Please search the archive at
>>> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
>>> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>>> * For (un)subscribe requests visit
>>> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
>>> send a mail to gmx-users-request at gromacs.org.
>> Gromacs Users mailing list
>> * Please search the archive at http://www.gromacs.org/
>> Support/Mailing_Lists/GMX-Users_List before posting!
>> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>> * For (un)subscribe requests visit
>> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
>> send a mail to gmx-users-request at gromacs.org.
Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow
Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 629
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201
jalemkul at outerbanks.umaryland.edu | (410) 706-7441
More information about the gromacs.org_gmx-users