[gmx-users] Protein fragments during production simulation
sapnauser365 at gmail.com
Thu Apr 30 16:12:09 CEST 2015
I have some doubts regarding pbc being new to gromacs. Since my protein was
breaking into fragments, I introduced pbc options and found the trajectory
whole and centred in the water box, but this was done after the simulation
and just for taking the snapshots with the command prompt:
trjconv -s run.tpr -f run.xtc -pbc mol -center -skip 500 -sep -o step.pdb
Is this the correct way to do it or do I have to incorporate pbc in the
.mdp files and re-run the simulation.
Also, before that I also checked if equilibration was an issue. So I
extended equilibration of the system to 5ns and 10ns. The results(without
using pbc) after a 10ns production run for both 5ns and 10ns equilibration
run shows that the protein is whole, however it had started breaking in
fixed intervals at fixed amino acid residues and then becomes whole for
around 2 ns and again fragments. The number of fragments has reduced
comparatively to the 1ns equilibration run.
Is this also a pbc related issue, and can someone please put some light on
what exactly is the problem going on?
On Wed, Apr 29, 2015 at 10:46 AM, SAPNA BORAH <sapnauser365 at gmail.com>
> Hi Tsjerk,
> Thank your for your help. I could get over the issue. :)
> On Tue, Apr 28, 2015 at 11:44 AM, Tsjerk Wassenaar <tsjerkw at gmail.com>
>> Hi Sapna,
>> On Apr 28, 2015 6:26 AM, "SAPNA BORAH" <sapnauser365 at gmail.com> wrote:
>> > Dear all,
>> > I am trying to simulate a globular protein predicted by Ab-initio method
>> > using web servers, Robetta and Itasser. I am currently using Gromacs
>> > A problem has occured during the production md run, when my protein
>> > to break into fragments while it remains inside the water box.
>> > Energy minimisation is achieved in 15000 steps using steep gradient.
>> > Equilibration is run for 1ns.
>> > Production md is run for 100ns.
>> > My complete mdp file for production run in as follows:
>> > ------run.mdp------
>> > integrator = md
>> > nsteps = 50000000
>> > dt = 0.002
>> > nstlist = 10
>> > rlist = 1.0
>> > coulombtype = pme
>> > rcoulomb = 1.0
>> > cutoff-scheme = verlet
>> > vdw-type = cut-off
>> > rvdw = 1.0
>> > tcoupl = v-rescale
>> > tc-grps = protein non-protein
>> > tau-t = 0.1 0.1
>> > ref-t = 298 298
>> > nstxtcout = 1000
>> > nstenergy = 1000
>> > constraints = all-bonds
>> > ------------------
>> > Can you please help me with the issue?
>> > Thanking You.
>> > Regards,
>> > Sapna Mayuri Borah
>> > Research Student
>> > Tezpur University
>> > India
>> > --
>> > Gromacs Users mailing list
>> > * Please search the archive at
>> > http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
>> > posting!
>> > * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>> > * For (un)subscribe requests visit
>> > https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
>> > send a mail to gmx-users-request at gromacs.org.
>> Gromacs Users mailing list
>> * Please search the archive at
>> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
>> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>> * For (un)subscribe requests visit
>> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
>> send a mail to gmx-users-request at gromacs.org.
> Sapna Mayuri Borah
> Research student
> Tezpur University,
Sapna Mayuri Borah
More information about the gromacs.org_gmx-users