<DIV>
<P class="MsoNormal" style="MARGIN: 0cm 0cm 0pt"><SPAN lang="EN-US"><FONT face="Times New Roman">Hi All,</FONT></SPAN></P>
<P class="MsoNormal" style="MARGIN: 0cm 0cm 0pt; TEXT-INDENT: 10.5pt; mso-char-indent-count: 1.0"><SPAN lang="EN-US"><FONT face="Times New Roman">I met a problem when I try to unfold a protein using Gromacs, It seemed the protein cannot be totally unfolded!</FONT></SPAN></P>
<P class="MsoNormal" style="MARGIN: 0cm 0cm 0pt; TEXT-INDENT: 10.5pt; mso-char-indent-count: 1.0"><SPAN lang="EN-US"><FONT face="Times New Roman">The simulated system has one Engrailed Homeodomain (En) protein (a three helix bundle protein with 54 residues, 629 atoms), total 4848 spce waters, and 7 Cl- used to neutralize the system in a 5.752(nm)^3 water BOX. I use the NTP ensemble with T=498K and P=26atm. The system has 1nm thick water in each side of the En protein, and the density of the system has been adjusted to <st1:chmetcnv w:st="on" UnitName="g" SourceValue=".829" HasSpace="True" Negative="False" NumberType="1" TCSC="0">0.829 g</st1:chmetcnv>/cm3.</FONT></SPAN></P>
<P class="MsoNormal" style="MARGIN: 0cm 0cm 0pt; TEXT-INDENT: 10.5pt; mso-char-indent-count: 1.0"><SPAN lang="EN-US"><FONT face="Times New Roman">The simulation lasted 24ns. The helixes disappeared at about 4ns. And after that some beta sheet formed in the N terminal of the protein. However, the protein kept in a compact state during the 24ns simulation. The radius of gyration of the protein fluctuated around 1.1nm during the simulation.</FONT></SPAN></P>
<P class="MsoNormal" style="MARGIN: 0cm 0cm 0pt; TEXT-INDENT: 10.5pt; mso-char-indent-count: 1.0"><SPAN lang="EN-US"><FONT face="Times New Roman">I've also noticed similar simulation done by others. For example, David Becka and Valerie Daggett reviewed their En protein unfolding in paper "Methods for molecular dynamics simulations of protein folding/unfolding in solution. Methods 34 (2004) 112¨C120". The simulations were performed with ENCAD and ilmm, and used an 0.8nm cutoff range. And the ensemble they used is NVE as I know. A stretched unfolding state occurred at about 5ns in their 60ns simulation in 498K, with little helix structure.</FONT></SPAN></P>
<P class="MsoNormal" style="MARGIN: 0cm 0cm 0pt; TEXT-INDENT: 10.5pt; mso-char-indent-count: 1.0"><SPAN lang="EN-US"><FONT face="Times New Roman">I was wondering whether the difference is caused by using different MD software and force field, or by some wrong parameters in my .mdp file. I've also conducted another 18ns simulation, and the result is almost the same. I list he mdp file below. Any comment is appreciated!</FONT></SPAN></P>
<P class="MsoNormal" style="MARGIN: 0cm 0cm 0pt; TEXT-INDENT: 10.5pt; mso-char-indent-count: 1.0"><SPAN lang="EN-US"><FONT face="Times New Roman"></FONT></SPAN> </P>
<P class="MsoNormal" style="MARGIN: 0cm 0cm 0pt; TEXT-INDENT: 10.5pt; mso-char-indent-count: 1.0"><SPAN lang="EN-US">;---------------------------------------------------------------------<BR>integrator = md<BR>dt = 0.002<BR>nsteps = 3000000<BR>nstxout = 500<BR>nstvout = 500<BR>nstlog = 500<BR>nstenergy = 500<BR>nstxtcout = 1000</SPAN></P>
<P class="MsoNormal" style="MARGIN: 0cm 0cm 0pt; TEXT-INDENT: 10.5pt; mso-char-indent-count: 1.0"><SPAN lang="EN-US">energygrps = protein non-protein</SPAN></P>
<P class="MsoNormal" style="MARGIN: 0cm 0cm 0pt; TEXT-INDENT: 10.5pt; mso-char-indent-count: 1.0"><SPAN lang="EN-US">nstlist = 10<BR>ns_type = grid<BR>pbc = xyz<BR>rlist = 1.0</SPAN></P>
<P class="MsoNormal" style="MARGIN: 0cm 0cm 0pt; TEXT-INDENT: 10.5pt; mso-char-indent-count: 1.0"><SPAN lang="EN-US">coulombtype = PME<BR>rcoulomb = 1.0<BR>vdwtype = cut-off<BR>rvdw = 1.4<BR>fourierspacing = 0.12<BR>fourier_nx = 0<BR>fourier_ny = 0<BR>fourier_nz = 0<BR>pme_order = 4<BR>ewald_rtol = 1e-5<BR>optimize_fft = yes</SPAN></P>
<P class="MsoNormal" style="MARGIN: 0cm 0cm 0pt; TEXT-INDENT: 10.5pt; mso-char-indent-count: 1.0"><SPAN lang="EN-US">tcoupl = v-rescale<BR>tc_grps = protein non-protein<BR>tau_t = 0.1 0.1<BR>ref_t = 498 498<BR>Pcoupl = Parrinello-Rahman<BR>pcoupltype = isotropic<BR>tau_p = 0.5<BR>compressibility = 4.5e-5<BR>ref_p = 26.0</SPAN></P>
<P class="MsoNormal" style="MARGIN: 0cm 0cm 0pt; TEXT-INDENT: 10.5pt; mso-char-indent-count: 1.0"><SPAN lang="EN-US">gen_vel = yes<BR>gen_temp = 498<BR>gen_seed = 173529</SPAN></P>
<P class="MsoNormal" style="MARGIN: 0cm 0cm 0pt; TEXT-INDENT: 10.5pt; mso-char-indent-count: 1.0"><SPAN lang="EN-US">constraints = hbonds<BR>lincs_order = 10<BR>;---------------------------------------------------------------------</SPAN></P></DIV><br><br><span title="neteasefooter"><span id="netease_mail_footer"></span></span>