nitu sharma wrote:
Hello Justin sir

sorry , The suggestion U given to me is not working. Actually I am a master student and just started to do simulation of membrane protein . Here no one did this work thats why i am facing more problem. sorry for putting question all the time in gmx-user list . But its only option for me to getting some information about my work.

OK, but you still need to be seen to be applying your own troubleshooting skills. If something doesn't work, try some alternatives. :-)

Also, learn the lesson that if you don't have someone locally who is useful to ask for basic help, don't take on the project!

Now, The the is still there . The main problem at this time is after inserting protein prorein into lipid bilayer by using inflategro script-

./inflategro.pl system.gro 4 DMP 14 inflated_system.gro 5 inflated_area.dat
this works and the structure I have got from here looks fine , in next step I am doing enegy minimisation of inflate_system.gro by using grompp- grompp -f em.mdp -c Inflated_system.gro -p topol.top -o inflated_system.tpr with mdrun- mdrun -v -s inflated_system.tpr -o inflated_system.trr -c inflated-em.pdb -e inflated-em.edr -g inflated-em.log This works well for 1000 nstep but structure I am getting from here the protein is totally distorted why?

Position restraints aren't magic. If the other forces on your protein are stronger then they'll overrule them. If they are very large, then you have a serious problem with your structure, or your model physics (i.e. the contents of your .itp files). To find out where the problem is coming from, you need to simplify your problem. Remove an element and see how things work differently. So try not using inflategro.pl, or changing its settings, or not using both your proteins, or not using any proteins, or not using position restraints.

I don't think you've ever mentioned where you got the contents of your .itp files from. If your structure after inflategro.pl is fine, and it's not after you apply some physics to it, then your physics is broken. We can't help more than that without a lot more information.

Mark

Even after applying strong position restrsain on protein.
The em.mdp file is-

pp                 =  /usr/bin/cpp
define              =  -DPOSRES
constraints         =  none
integrator          =  steep
nsteps              =  1000
;
;       Energy minimizing stuff
;
emtol               =  1000
emstep              =  0.0001

nstcomm             =  1
ns_type             =  grid
rlist               =  1
rcoulomb            =  1.0
rvdw                =  1.0
Tcoupl              =  no
Pcoupl              =  no
gen_vel             =  no
comm-mode           = Linear
pbc                 = XYZ

Is anything wrong in applying postion restrain in parameter file?

this is topology file-
; Include forcefield parameters
#include "ffG53a6_lipid.itp"

; Include chain topology
#include "topol_A.itp"
#include "topol_B.itp"

; Include DMPC chain topology
#include "dmpc.itp"

[ system ]
; Name
Protein in DMPC bilayer

[ molecules ]
; Compound        #mols
Protein_A           1
Protein_B           1
DMPC               120

Where I should make changes to make protein as rigid molecule I am unable to understand.

can u suggest me something ? I am asking again becoz the step of compression and minimisation is working but the protein structure become distorted , and I want to keep protein in its proper configuration.
Any suggestion specific to this will really helpful for me.

Nitu sharma
School of Life Sciences
Jawaherlal nehru University
New Delhi,India


------------------------------------------------------------------------

_______________________________________________
gmx-users mailing list    gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/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/mailing_lists/users.php
_______________________________________________
gmx-users mailing list    gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/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/mailing_lists/users.php

Reply via email to