[gmx-users] Abnormal popc struture after inserting protein

minnale minnale_gnos at rediffmail.com
Wed Jul 2 06:22:38 CEST 2008


  


Thanks to Justin for comment about my problem
I will tell you clearly
1. Initially POPC bilayer taken from Tieleman website and run the simulation for 5ns.
2. Protein also has taken and simulated for 9ns actually this contain 10 helices among that 10 helices two helices have taken, because i want only 2 and generated .top file for 2 helices by using *pdb2gmx command. 
3. By using VMD orient 2 helices properly in POPC bilayer
4. BY using genbox command these helices embedded into POPC
5. Minimisation has done, it ran fine without error but it didnt reach FMAX<1000
6. Ran restarin with fc1000 showed the abnrmal POPC and water molecules stucture.
  
   Earlier I have run similar system with ffgmx FF with similar time 50 ps restrain that structure pretty fine.The idea of simulating system with short time is if it runs fine with short time I will give long step. 
 
Toplogy file is 
 
 ; Include forcefield parameters
#include "ffoplsaa.itp"
#include "popc.itp"

[ moleculetype ]
; Name            nrexcl
Protein             3

[ atoms ]
;   nr       type  resnr residue  atom   cgnr     charge       mass  typeB    chargeB      massB
     1   opls_287      1    ALA      N      1       -0.3    14.0067   ; qtot -0.3
     2   opls_290      1    ALA     H1      1       0.33      1.008   ; qtot 0.03
     3   opls_290      1    ALA     H2      1       0.33      1.008   ; qtot 0.36
     4   opls_290      1    ALA     H3      1       0.33      1.008   ; qtot 0.69
     5  opls_293B      1    ALA     CA      1       0.25     12.011   ; qtot 0.94
     6   opls_140      1    ALA     HA      1       0.06      1.008   ; qtot 1
     7   opls_135      1    ALA     CB      2      -0.18     12.011   ; qtot 0.82
     8   opls_140      1    ALA    HB1      2       0.06      1.008   ; qtot 0.88
     9   opls_140      1    ALA    HB2      2       0.06      1.008   ; qtot 0.94
    10   opls_140      1    ALA    HB3      2       0.06      1.008   ; qtot 1
    11   opls_235      1    ALA      C      3        0.5     12.011   ; qtot 1.5
    12   opls_236      1    ALA      O      3       -0.5    15.9994   ; qtot 1
    13   opls_238      2    VAL      N      4       -0.5    14.0067   ; qtot 0.5
    14   opls_241      2    VAL      H      4        0.3      1.008   ; qtot 0.8
    15  opls_224B      2    VAL     CA      4       0.14     12.011   ; qtot 0.94
    16   opls_140      2    VAL     HA      4       0.06      1.008   ; qtot 1
   .
   .
   .
   .
   .
   .
   .
   .
   .
   .
 ; Include Position restrain protein
#ifdef POSRES_PROTEIN
#include "posre.itp"
#endif

;Include Position restrain lipid
#ifdef POSRES_LIPID
#include "lip_posre.itp"
#endif

; Include water topology
#include "spc.itp"

#ifdef POSRES_WATER
; Position restraint for each water oxygen
[ position_restraints ]
;  i funct       fcx        fcy        fcz
   1    1       1000       1000       1000
#endif

; Include generic topology for ions
#include "ions.itp"

[ system ]
; Name
GROtesk MACabre and Sinister in water

[ molecules ]
; Compound        #mols
Protein             1
POPC               94
SOL              2434
CL-                 3

Any comments wiill be appreciated.

Thanks in advance

   




On Tue, 01 Jul 2008 Justin A.Lemkul wrote :
>
>
>minnale wrote:
>>
>>
>>   So I dont boughter about those sentences because I confirmed in list archives that I can proceed further steps and moreover em.gro file is      fine, so I went for restrain, here I am getting abnormal POPC structure
>>means POPC tails are tilting, water molecules structure also disturbed.
>>
>>Could you please tell me where iam doing mistake?
>>
>
>Not really.  You haven't told us much.  What parameters are you using (for the protein and the lipids)?  How did you insert your protein?
>
>Also realize that 50 ps is *extremely* short in the realm of membrane protein MD.  Even after several hundred picoseconds things may look a bit strange.  Prevailing wisdom is that between 25-30 ns of MD are required for the lipids to reach equilibrium (although they should look reasonable after just a few ns, but parameters like area per headgroup may not level off for a while).
>
>-Justin
>
>>my pr.mdp filetitle              =  protein in popc restrained
>>define              =  -DPOSRES_protein
>>constraints        =  all-bonds
>>integrator          =  md
>>dt                  =  0.002    ; ps !
>>nsteps              =  10000    ; total 50 ps.
>>nstcomm            =  1
>>nstxout            =  50
>>nstvout            =  1000
>>nstfout            =  0
>>nstlog              =  10
>>nstenergy          =  10
>>nstlist            =  10
>>ns_type            =  grid
>>rlist              =  0.9
>>coulombtype        =  PME
>>rcoulomb            =  0.9
>>rvdw                =  1.4
>>pbc                =  xyz
>>; Berendsen temperature coupling is on in three groups
>>Tcoupl              =  Berendsen
>>tc-grps            =  POPC    Protein  SOL_CL-
>>tau_t              =  0.1      0.1      0.1
>>ref_t              =  310      310      310
>>; Anisotropic pressure coupling is now on
>>Pcoupl              =  berendsen
>>pcoupltype          =  anisotropic
>>tau_p              =  1.0      1.0    1.0    0        0      0
>>compressibility    =  4.5e-5  4.5e-5  4.5e-5  0        0      0
>>ref_p              =  1.0      1.0    1.0    0        0      0
>>; Energy monitoring
>>energygrps          =  POPC  Protein  SOL_CL-
>>; Generate velocites is on at 300 K.
>>gen_vel            =  yes
>>gen_temp            =  310.0
>>gen_seed            =  173529
>>
>>Thanks in advance.                       
>>
>>
>>------------------------------------------------------------------------
>>
>>_______________________________________________
>>gmx-users mailing list    gmx-users at 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-request at gromacs.org.
>>Can't post? Read http://www.gromacs.org/mailing_lists/users.php
>
>-- ========================================
>
>Justin A. Lemkul
>Graduate Research Assistant
>Department of Biochemistry
>Virginia Tech
>Blacksburg, VA
>jalemkul[at]vt.edu | (540) 231-9080
>http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin
>
>========================================
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20080702/65faab62/attachment.html>


More information about the gromacs.org_gmx-users mailing list