[gmx-users] grompp is adding a bond in my model that should not be there

[Jackson Chief]

The missing residues correspond to a loop that connects two alpha helices.
 The C-terminal end of one helix, and N-terminal end of the other helix
contain a random coil.  After energy minimization, the coils are extended
in space and not connected.  I visualized this by looking at the .gro file
in VMD.  After 5ns of unrestrained MD, the two coils are clearly connected
by looking at the output .gro file in VMD.  Weird.



On Sun, Mar 25, 2012 at 9:48 PM, Mark Abraham <Mark.Abraham at anu.edu.au>wrote:

> On 26/03/2012 12:20 PM, Jackson Chief wrote:
>
>> I made a model of a receptor protein, bilayer, and solvent.  My protein
>> contains a 20 residue gap.  This gap corresponds to a region of the protein
>> that had been digested by trypsin before crystallization.  The trypsin
>> digestion has no affect on receptor activity experimentally. I performed
>> energy minimization without problem.  The protein looked like it should,
>> containing the gap.
>>
>
> You have to treat this gap somehow. Either you have to cap the peptide
> chains (see pdb2gmx -h), or model in the missing residues using some
> (non-GROMACS) software.
>
>
>
>> When I performed equilibration, the output had the C-terminus of one
>> protein fragment connected to the N-terminus of the other protein fragment.
>>  I don't know how this peptide bond was created, because it was not in the
>> input *.gro file to grompp.  Please help.
>>
>>
> .gro files have coordinates, never bonds. Since you haven't described how
> you are treating the gap, and haven't said how you've observed the
> "creation" of a peptide bond, it's hard to give specific guidance.
>
> Mark
> --
> gmx-users mailing list    gmx-users at gromacs.org
> http://lists.gromacs.org/**mailman/listinfo/gmx-users<http://lists.gromacs.org/mailman/listinfo/gmx-users>
> Please search the archive at http://www.gromacs.org/**
> Support/Mailing_Lists/Search<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-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/**Support/Mailing_Lists<http://www.gromacs.org/Support/Mailing_Lists>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20120325/540971f8/attachment.html>