[gmx-users] hbonds(non covalent) contraint

[Martina Bertsch, PhD]

Dear Lakshmi,

First, define a custom group of atoms you wish to restrain using make_ndx:

make_ndx -f my_protein.gro -n index.ndx

For example, to define a group containing C-alpha atoms of residues 
100-125, at the prompt you would type:

 >r 100-125 & a CA
 >q

Next, use genpr to generate a position restraint itp file corresponding 
to that custom group:

genpr -f my_protein.gro -n index.ndx -o my_res.itp

and select the desired group at the prompt.

Please, see the Manual or the online Reference for details.

Regards,
Martina



Lakshmi Padmavathi wrote:

>Dear gromacs users,
>           I have a protein which contains 12 helices. while doing simulation 
>I would like to constrain or restraint the hbonds between the two turns of 
>helix to constraint the alpha- helix structure.
>
>how can i do this? to write the distance restrints, i have to write hundres of 
>atoms, is there any other was to maintain the helix structure trought the 
>simulation.
>
>very thankful for your suggestions.
>
>thank you very much.
>
>sincerely,
>
>         
>
>_______________________________________________
>gmx-users mailing list
>gmx-users at gromacs.org
>http://www.gromacs.org/mailman/listinfo/gmx-users
>Please don't post (un)subscribe requests to the list. Use the 
>www interface or send it to gmx-users-request at gromacs.org.
>
>