Re: [AMBER] donor and acceptor masks

From: Daniel Roe <daniel.r.roe.gmail.com>
Date: Sat, 17 Oct 2015 11:43:58 -0600

On Sat, Oct 17, 2015 at 10:56 AM, Hadeer ELHabashy
<hadeer.elhabashi.gmail.com> wrote:
> That was great, I figured out a problem was the atoms numbering
> what if I want to use residue names instead of numbers
> is this correct ?!
> Acceptor MET O Donorh ALA H donor ALA N is this equilvelant to this

Have a look at the Amber mask syntax that cpptraj uses (Amber 15
manual 29.1.6). The 'hbond' command actually does a reasonable job of
detecting donors and acceptors on its own. For example, If you just
wanted to look at hydrogen bonds between residues named GLU and
residues named LYS you could do something like:

hbond :GLU,LYS out nhb.dat avgout avghb.dat

Hope this helps,

-Dan

>
>
>
> On Sat, Oct 17, 2015 at 6:41 PM, Thomas Cheatham <tec3.utah.edu> wrote:
>
>>
>> Start simple, and try to figure out the syntax for one command... Look at
>> all the warnings and messages.
>>
>> A good step was also looking at the visualization with VMD, but what does
>> it mean that VMD "detected" hydrogen bondings? By eye?
>>
>> I would also try calculation a simple distance and plotting that vs. time
>> (and also the distribution) to see if 3.5 angstroms is the appropriate
>> cutoff for these interactions.
>>
>> distance d1 out d1.dat :1.N :6.OXT
>>
>> > Actually when I visualize the system with VMD , I can visualize many
>> > hydrogen bondings
>> > but AMBER did not detected it
>> >
>> > > finally I ran this successfully, but I got all Zeros and that is
>> strange!
>> > > I ran with NVT with Langevin thermostat and all interaction calculated
>> > >
>> > > hbond out N_H1_OXT.out donormask :1.N donorhmask :1.H1 acceptormask
>> > > :6.OXT dist 3.5 series h1
>> > > hbond out N_H2_OXT.out donormask :1.N donorhmask :1.H2 acceptormask
>> > > :6.OXT dist 3.5 series h2
>> > > hbond out N_H3_OXT.out donormask :1.N donorhmask :1.H3 acceptormask
>> > > :6.OXT dist 3.5 series h3
>> > >
>> > > what if I want to get the life time of the bond ?
>> > > should just write life time in the command ?!
>>
>> _______________________________________________
>> AMBER mailing list
>> AMBER.ambermd.org
>> http://lists.ambermd.org/mailman/listinfo/amber
>>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber



-- 
-------------------------
Daniel R. Roe, PhD
Department of Medicinal Chemistry
University of Utah
30 South 2000 East, Room 307
Salt Lake City, UT 84112-5820
http://home.chpc.utah.edu/~cheatham/
(801) 587-9652
(801) 585-6208 (Fax)
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sat Oct 17 2015 - 11:00:03 PDT
Custom Search