[gmx-users] High rmsf value
jalemkul at vt.edu
Sat Apr 21 19:53:20 CEST 2018
On 4/21/18 1:13 PM, SAYAN BHATTACHARJEE wrote:
> Dear Lemkul,
> My complex have two protein( Chain A & B) as dimer bound to DNA. So, I have
> to make index consisting only chain A or B, excluding DNA, then center and
> output options will both be chain A or B in trjconv command?
> Is it possible to select both chains A & B?
The selections for centering and output are entirely separate. Center on
chain, output whatever makes sense for your needs.
> On Sat, 21 Apr 2018, 22:00 Justin Lemkul, <jalemkul at vt.edu> wrote:
>> On 4/21/18 12:24 PM, SAYAN BHATTACHARJEE wrote:
>>> I bave simulated protein_dna complex for 500ns following gromacs manuals,
>>> during visualisation in vmd its doesn't maintain pbc, moves out from
>> box. I
>>> have tried those following commands, but failed to make my complex within
>>> box throughout trajectory marging protein and dna only in index.
>>> I have tried,
>>> 1. Trjconv -f 500ns.xtc -s md_0_1.tpr -pbc cluster -n protein_dna.ndx -o
>>> center.xtc -ur compact -center
>>> Auto select protein_dna from input index for both centering and output.
>>> I have tried -pbc mol option too but it didn't work here.
>>> Visualisation in vmd is now fine(using cluster option and loading
>>> center.xtc trajectory), complex is now within box, but, in the next step
>>> when I tried nojump option and fit with reference initial structure by
>>> uising the following commands, pbc condition didn't maintain and protein
>>> motion was horrible and the rmsf shows very high value(using fit.xtc/
>>> nojump.xtc/ centre.xtc for rmsf calculation).
>>> 2.Trjconv -f center.xtc -s md_0_1.tpr -pbc nojump -n protein_dna.ndx -o
>>> 3.Trjconv -f nojump.xtc -s md_0_1.tpr -fit rot+trans -n protein_dna.ndx
>>> Please give suggestions in this regard.
>> If you have a complex, don't try to center on the entire thing. Choose
>> one chain of the complex and center on that. If elements of the complex
>> are split on opposing "faces" of the box, their center is still
>> (geometrically speaking) coincident with the actual center of the box,
>> so trjconv works properly but it's not what you want.
>> Justin A. Lemkul, Ph.D.
>> Assistant Professor
>> Virginia Tech Department of Biochemistry
>> 303 Engel Hall
>> 340 West Campus Dr.
>> Blacksburg, VA 24061
>> jalemkul at vt.edu | (540) 231-3129
>> 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.
Justin A. Lemkul, Ph.D.
Virginia Tech Department of Biochemistry
303 Engel Hall
340 West Campus Dr.
Blacksburg, VA 24061
jalemkul at vt.edu | (540) 231-3129
More information about the gromacs.org_gmx-users