<P>
&nbsp; <BR>
<BR>
<BR>
Thanks to Justin for comment about my problem<BR>
I will tell you clearly<BR>
1. Initially POPC bilayer taken from Tieleman website and run the simulation for 5ns.<BR>
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. <BR>
3. By using VMD orient 2 helices properly in POPC bilayer<BR>
4. BY using genbox command these helices embedded into POPC<BR>
5. Minimisation has done, it ran fine without error but it didnt reach FMAX&lt;1000<BR>
6. Ran restarin with fc1000 showed the abnrmal POPC and water molecules stucture.<BR>
&nbsp; <BR>
&nbsp;  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. <BR>
 <BR>
Toplogy file is <BR>
 <BR>
 ; Include forcefield parameters<BR>
#include &quot;ffoplsaa.itp&quot;<BR>
#include &quot;popc.itp&quot;<BR>
<BR>
[ moleculetype ]<BR>
; Name&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; nrexcl<BR>
Protein&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  3<BR>
<BR>
[ atoms ]<BR>
;&nbsp;  nr&nbsp; &nbsp; &nbsp;  type&nbsp; resnr residue&nbsp; atom&nbsp;  cgnr&nbsp; &nbsp;  charge&nbsp; &nbsp; &nbsp;  mass&nbsp; typeB&nbsp; &nbsp; chargeB&nbsp; &nbsp; &nbsp; massB<BR>
&nbsp; &nbsp;  1&nbsp;  opls_287&nbsp; &nbsp; &nbsp; 1&nbsp; &nbsp; ALA&nbsp; &nbsp; &nbsp; N&nbsp; &nbsp; &nbsp; 1&nbsp; &nbsp; &nbsp;  -0.3&nbsp; &nbsp; 14.0067&nbsp;  ; qtot -0.3<BR>
&nbsp; &nbsp;  2&nbsp;  opls_290&nbsp; &nbsp; &nbsp; 1&nbsp; &nbsp; ALA&nbsp; &nbsp;  H1&nbsp; &nbsp; &nbsp; 1&nbsp; &nbsp; &nbsp;  0.33&nbsp; &nbsp; &nbsp; 1.008&nbsp;  ; qtot 0.03<BR>
&nbsp; &nbsp;  3&nbsp;  opls_290&nbsp; &nbsp; &nbsp; 1&nbsp; &nbsp; ALA&nbsp; &nbsp;  H2&nbsp; &nbsp; &nbsp; 1&nbsp; &nbsp; &nbsp;  0.33&nbsp; &nbsp; &nbsp; 1.008&nbsp;  ; qtot 0.36<BR>
&nbsp; &nbsp;  4&nbsp;  opls_290&nbsp; &nbsp; &nbsp; 1&nbsp; &nbsp; ALA&nbsp; &nbsp;  H3&nbsp; &nbsp; &nbsp; 1&nbsp; &nbsp; &nbsp;  0.33&nbsp; &nbsp; &nbsp; 1.008&nbsp;  ; qtot 0.69<BR>
&nbsp; &nbsp;  5&nbsp; opls_293B&nbsp; &nbsp; &nbsp; 1&nbsp; &nbsp; ALA&nbsp; &nbsp;  CA&nbsp; &nbsp; &nbsp; 1&nbsp; &nbsp; &nbsp;  0.25&nbsp; &nbsp;  12.011&nbsp;  ; qtot 0.94<BR>
&nbsp; &nbsp;  6&nbsp;  opls_140&nbsp; &nbsp; &nbsp; 1&nbsp; &nbsp; ALA&nbsp; &nbsp;  HA&nbsp; &nbsp; &nbsp; 1&nbsp; &nbsp; &nbsp;  0.06&nbsp; &nbsp; &nbsp; 1.008&nbsp;  ; qtot 1<BR>
&nbsp; &nbsp;  7&nbsp;  opls_135&nbsp; &nbsp; &nbsp; 1&nbsp; &nbsp; ALA&nbsp; &nbsp;  CB&nbsp; &nbsp; &nbsp; 2&nbsp; &nbsp; &nbsp; -0.18&nbsp; &nbsp;  12.011&nbsp;  ; qtot 0.82<BR>
&nbsp; &nbsp;  8&nbsp;  opls_140&nbsp; &nbsp; &nbsp; 1&nbsp; &nbsp; ALA&nbsp; &nbsp; HB1&nbsp; &nbsp; &nbsp; 2&nbsp; &nbsp; &nbsp;  0.06&nbsp; &nbsp; &nbsp; 1.008&nbsp;  ; qtot 0.88<BR>
&nbsp; &nbsp;  9&nbsp;  opls_140&nbsp; &nbsp; &nbsp; 1&nbsp; &nbsp; ALA&nbsp; &nbsp; HB2&nbsp; &nbsp; &nbsp; 2&nbsp; &nbsp; &nbsp;  0.06&nbsp; &nbsp; &nbsp; 1.008&nbsp;  ; qtot 0.94<BR>
&nbsp; &nbsp; 10&nbsp;  opls_140&nbsp; &nbsp; &nbsp; 1&nbsp; &nbsp; ALA&nbsp; &nbsp; HB3&nbsp; &nbsp; &nbsp; 2&nbsp; &nbsp; &nbsp;  0.06&nbsp; &nbsp; &nbsp; 1.008&nbsp;  ; qtot 1<BR>
&nbsp; &nbsp; 11&nbsp;  opls_235&nbsp; &nbsp; &nbsp; 1&nbsp; &nbsp; ALA&nbsp; &nbsp; &nbsp; C&nbsp; &nbsp; &nbsp; 3&nbsp; &nbsp; &nbsp; &nbsp; 0.5&nbsp; &nbsp;  12.011&nbsp;  ; qtot 1.5<BR>
&nbsp; &nbsp; 12&nbsp;  opls_236&nbsp; &nbsp; &nbsp; 1&nbsp; &nbsp; ALA&nbsp; &nbsp; &nbsp; O&nbsp; &nbsp; &nbsp; 3&nbsp; &nbsp; &nbsp;  -0.5&nbsp; &nbsp; 15.9994&nbsp;  ; qtot 1<BR>
&nbsp; &nbsp; 13&nbsp;  opls_238&nbsp; &nbsp; &nbsp; 2&nbsp; &nbsp; VAL&nbsp; &nbsp; &nbsp; N&nbsp; &nbsp; &nbsp; 4&nbsp; &nbsp; &nbsp;  -0.5&nbsp; &nbsp; 14.0067&nbsp;  ; qtot 0.5<BR>
&nbsp; &nbsp; 14&nbsp;  opls_241&nbsp; &nbsp; &nbsp; 2&nbsp; &nbsp; VAL&nbsp; &nbsp; &nbsp; H&nbsp; &nbsp; &nbsp; 4&nbsp; &nbsp; &nbsp; &nbsp; 0.3&nbsp; &nbsp; &nbsp; 1.008&nbsp;  ; qtot 0.8<BR>
&nbsp; &nbsp; 15&nbsp; opls_224B&nbsp; &nbsp; &nbsp; 2&nbsp; &nbsp; VAL&nbsp; &nbsp;  CA&nbsp; &nbsp; &nbsp; 4&nbsp; &nbsp; &nbsp;  0.14&nbsp; &nbsp;  12.011&nbsp;  ; qtot 0.94<BR>
&nbsp; &nbsp; 16&nbsp;  opls_140&nbsp; &nbsp; &nbsp; 2&nbsp; &nbsp; VAL&nbsp; &nbsp;  HA&nbsp; &nbsp; &nbsp; 4&nbsp; &nbsp; &nbsp;  0.06&nbsp; &nbsp; &nbsp; 1.008&nbsp;  ; qtot 1<BR>
&nbsp;  .<BR>
&nbsp;  .<BR>
&nbsp;  .<BR>
&nbsp;  .<BR>
&nbsp;  .<BR>
&nbsp;  .<BR>
&nbsp;  .<BR>
&nbsp;  .<BR>
&nbsp;  .<BR>
&nbsp;  .<BR>
 ; Include Position restrain protein<BR>
#ifdef POSRES_PROTEIN<BR>
#include &quot;posre.itp&quot;<BR>
#endif<BR>
<BR>
;Include Position restrain lipid<BR>
#ifdef POSRES_LIPID<BR>
#include &quot;lip_posre.itp&quot;<BR>
#endif<BR>
<BR>
; Include water topology<BR>
#include &quot;spc.itp&quot;<BR>
<BR>
#ifdef POSRES_WATER<BR>
; Position restraint for each water oxygen<BR>
[ position_restraints ]<BR>
;&nbsp; i funct&nbsp; &nbsp; &nbsp;  fcx&nbsp; &nbsp; &nbsp; &nbsp; fcy&nbsp; &nbsp; &nbsp; &nbsp; fcz<BR>
&nbsp;  1&nbsp; &nbsp; 1&nbsp; &nbsp; &nbsp;  1000&nbsp; &nbsp; &nbsp;  1000&nbsp; &nbsp; &nbsp;  1000<BR>
#endif<BR>
<BR>
; Include generic topology for ions<BR>
#include &quot;ions.itp&quot;<BR>
<BR>
[ system ]<BR>
; Name<BR>
GROtesk MACabre and Sinister in water<BR>
<BR>
[ molecules ]<BR>
; Compound&nbsp; &nbsp; &nbsp; &nbsp; #mols<BR>
Protein&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  1<BR>
POPC&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  94<BR>
SOL&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 2434<BR>
CL-&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  3<BR>
<BR>
Any comments wiill be appreciated.<BR>
<BR>
Thanks in advance<BR>
<BR>
&nbsp;  <BR>
<BR>
<BR>
<BR>
<BR>
On Tue, 01 Jul 2008 Justin A.Lemkul wrote :<BR>
&gt;<BR>
&gt;<BR>
&gt;minnale wrote:<BR>
&gt;&gt;<BR>
&gt;&gt;<BR>
&gt;&gt;&nbsp;  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&nbsp; &nbsp; &nbsp; fine, so I went for restrain, here I am getting abnormal POPC structure<BR>
&gt;&gt;means POPC tails are tilting, water molecules structure also disturbed.<BR>
&gt;&gt;<BR>
&gt;&gt;Could you please tell me where iam doing mistake?<BR>
&gt;&gt;<BR>
&gt;<BR>
&gt;Not really.&nbsp; You haven't told us much.&nbsp; What parameters are you using (for the protein and the lipids)?&nbsp; How did you insert your protein?<BR>
&gt;<BR>
&gt;Also realize that 50 ps is *extremely* short in the realm of membrane protein MD.&nbsp; Even after several hundred picoseconds things may look a bit strange.&nbsp; 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).<BR>
&gt;<BR>
&gt;-Justin<BR>
&gt;<BR>
&gt;&gt;my pr.mdp filetitle&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; protein in popc restrained<BR>
&gt;&gt;define&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; -DPOSRES_protein<BR>
&gt;&gt;constraints&nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; all-bonds<BR>
&gt;&gt;integrator&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; md<BR>
&gt;&gt;dt&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; 0.002&nbsp; &nbsp; ; ps !<BR>
&gt;&gt;nsteps&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; 10000&nbsp; &nbsp; ; total 50 ps.<BR>
&gt;&gt;nstcomm&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; 1<BR>
&gt;&gt;nstxout&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; 50<BR>
&gt;&gt;nstvout&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; 1000<BR>
&gt;&gt;nstfout&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; 0<BR>
&gt;&gt;nstlog&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; 10<BR>
&gt;&gt;nstenergy&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; 10<BR>
&gt;&gt;nstlist&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; 10<BR>
&gt;&gt;ns_type&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; grid<BR>
&gt;&gt;rlist&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; 0.9<BR>
&gt;&gt;coulombtype&nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; PME<BR>
&gt;&gt;rcoulomb&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; 0.9<BR>
&gt;&gt;rvdw&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; 1.4<BR>
&gt;&gt;pbc&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; xyz<BR>
&gt;&gt;; Berendsen temperature coupling is on in three groups<BR>
&gt;&gt;Tcoupl&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; Berendsen<BR>
&gt;&gt;tc-grps&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; POPC&nbsp; &nbsp; Protein&nbsp; SOL_CL-<BR>
&gt;&gt;tau_t&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; 0.1&nbsp; &nbsp; &nbsp; 0.1&nbsp; &nbsp; &nbsp; 0.1<BR>
&gt;&gt;ref_t&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; 310&nbsp; &nbsp; &nbsp; 310&nbsp; &nbsp; &nbsp; 310<BR>
&gt;&gt;; Anisotropic pressure coupling is now on<BR>
&gt;&gt;Pcoupl&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; berendsen<BR>
&gt;&gt;pcoupltype&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; anisotropic<BR>
&gt;&gt;tau_p&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; 1.0&nbsp; &nbsp; &nbsp; 1.0&nbsp; &nbsp; 1.0&nbsp; &nbsp; 0&nbsp; &nbsp; &nbsp; &nbsp; 0&nbsp; &nbsp; &nbsp; 0<BR>
&gt;&gt;compressibility&nbsp; &nbsp; =&nbsp; 4.5e-5&nbsp; 4.5e-5&nbsp; 4.5e-5&nbsp; 0&nbsp; &nbsp; &nbsp; &nbsp; 0&nbsp; &nbsp; &nbsp; 0<BR>
&gt;&gt;ref_p&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; 1.0&nbsp; &nbsp; &nbsp; 1.0&nbsp; &nbsp; 1.0&nbsp; &nbsp; 0&nbsp; &nbsp; &nbsp; &nbsp; 0&nbsp; &nbsp; &nbsp; 0<BR>
&gt;&gt;; Energy monitoring<BR>
&gt;&gt;energygrps&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; POPC&nbsp; Protein&nbsp; SOL_CL-<BR>
&gt;&gt;; Generate velocites is on at 300 K.<BR>
&gt;&gt;gen_vel&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; yes<BR>
&gt;&gt;gen_temp&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; 310.0<BR>
&gt;&gt;gen_seed&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; 173529<BR>
&gt;&gt;<BR>
&gt;&gt;Thanks in advance.&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <BR>
&gt;&gt;<BR>
&gt;&gt;<BR>
&gt;&gt;------------------------------------------------------------------------<BR>
&gt;&gt;<BR>
&gt;&gt;_______________________________________________<BR>
&gt;&gt;gmx-users mailing list&nbsp; &nbsp; gmx-users@gromacs.org<BR>
&gt;&gt;http://www.gromacs.org/mailman/listinfo/gmx-users<BR>
&gt;&gt;Please search the archive at http://www.gromacs.org/search before posting!<BR>
&gt;&gt;Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-request@gromacs.org.<BR>
&gt;&gt;Can't post? Read http://www.gromacs.org/mailing_lists/users.php<BR>
&gt;<BR>
&gt;-- ========================================<BR>
&gt;<BR>
&gt;Justin A. Lemkul<BR>
&gt;Graduate Research Assistant<BR>
&gt;Department of Biochemistry<BR>
&gt;Virginia Tech<BR>
&gt;Blacksburg, VA<BR>
&gt;jalemkul[at]vt.edu | (540) 231-9080<BR>
&gt;http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin<BR>
&gt;<BR>
&gt;========================================<BR>
&gt;<BR>

</P>
<br><br>