VMD-L Mailing List

From: Ajasja Ljubetič (ajasja.ljubetic_at_gmail.com)
Date: Mon May 30 2016 - 03:11:25 CDT

Hi!

Thanks Chitrak and Norman! I tried both suggestions but still did not
manage to wrap/join everything in one unit.
It's funny, but just rerunning the problematic part of simulation with
"wrapAll off" will probably be the quickest solution (in terms of my
personal time). Ahh the joy of computer power:)

Best regards,
Ajasja

On 30 May 2016 at 07:51, Norman Geist <norman.geist_at_uni-greifswald.de>
wrote:

> Hey,
>
>
>
> If your solute is split by the periodic boundaries and are not connected
> by bonds, neither wrap –center com –centersel “your solute” nor pbc join
> can correctly work. But if there’s a single atom, which can represent the
> center of your solute so it will completely fit into the periodic cell, you
> can use that atom for –centersel
>
>
>
> pbc wrap –all –compound res –center com –centersel „index 1234“
>
>
>
>
>
> Norman Geist
>
>
>
> *Von:* owner-vmd-l_at_ks.uiuc.edu [mailto:owner-vmd-l_at_ks.uiuc.edu] *Im
> Auftrag von *Chitrak Gupta
> *Gesendet:* Freitag, 27. Mai 2016 19:48
> *An:* Ajasja Ljubetič <ajasja.ljubetic_at_gmail.com>
> *Cc:* vmd-l_at_ks.uiuc.edu
> *Betreff:* Re: vmd-l: Wrapping problem
>
>
>
> Hi Ajasja,
>
>
>
> Have you tried something like this?
>
>
>
> pbc wrap -all -sel "water or ions" -compound res -center com -centersel
> "protein"
>
>
>
>
>
> This one works if:
>
> 1. All your chains are proteins
>
> 2. Only other things you have in your system apart from these chains come
> under the category "water or ions".
>
>
>
> You could modify this based on your needs, but you get the basic point.
> Also, you could also change "-compound res" to "-compound fragment".
>
>
>
> Hope this helps,
>
> Chitrak.
>
>
>
> On Fri, May 27, 2016 at 12:10 PM, Ajasja Ljubetič <
> ajasja.ljubetic_at_gmail.com> wrote:
>
> Dear All!
>
> I have run a simulation in NAMD using "wrapAll on". I'm simulating a
> complex of three chains.
> In one simulation (out of 14) I get large jumps, such as these
> <https://gfycat.com/CapitalHollowAmmonite> (link to movie)
> <https://gfycat.com/CapitalHollowAmmonite>:
>
> I have tried "pbc unwrap", "pbc wrap" and "pbc join", but Ihave yet to find
> the right combination of these command to correctly unwrap all three chains
>
>
>