AMBER Archive (2005)

Subject: AMBER: hbond analysis in ptaj/carnal

From: Jana Precechtelova (janap_at_ncbr.chemi.muni.cz)
Date: Thu May 26 2005 - 11:25:20 CDT


Dear Amber users,

OBJECTIVE:
==========
I would like to carry out a hydrogen boding
analysis for a simple DNA d(CGCGAATTCGCG)2
in order to find out how many hydrogen bonds
are created between water and phosphate
oxygens O1P/O2P and in order to get a statistics
like: e.g.

over 10% of the trajectory, there are 3 H-bonds
to O1P and 2 H-bonds to O2P (the 3/2 configuration)

I would like to use the O1P/O2P...H(wat) distance
and O1P/O2P...H-O(wat) angle to restrict the
search for H-bonds.

PROBLEMS:
=========
1. The distance cutoff in Amber is defined with
the help of the heavy atoms distance, not
with the H-bond distance.

2. The angle definition in Amber doesn't correspond
to the standard H-bond angle definition.
In ptraj/carnal, the angle is defined as H-DONOR-ACCEPTOR. However, the
standard definition
is DONOR-H-ACCEPTOR.

In other words, neither the distance nor the
angle definition in Amber allows me to employ
the geometry restrictions for H-bonds that I want.

QUESTIONS:
==========
1. Is there any way in Amber how to get
a hydrogen bonding statistics using the
the H-bond distance and the DONOR-H-ACCEPTOR
angle as geometry criteria for hydrogen bonding?

2. If there is no way haw to treat the
problems 1 and 2 at the same time, is it at least
possible to overcome one of the two problems?

3. I am not sure about the definitions od the
DONOR and ACCEPTOR in ptraj and carnal.
In ptraj, DONOR seemes to be defined as a lone pair
donor, while in carnal DONOR seems to be
defined as a donor of the hydrogen atom.
It would mean that a DONOR in ptraj is
actually an ACCEPTOR i carnal and vice versa,
am I right?

Many thanks for your suggestions

Regards,
Jana

=============================================================
         Jana Precechtelova (Ph.D. student)
-------------------------------------------------------------
         Institut fuer Anorganische Chemie
         Bayerische Julius-Maximilians Universitaet Wuerzburg
         Am Hubland, 97074 Wuerzburg, Germany
                    &
         NMR Laboratory
         National Centre for Biomolecular Research (NCBR)
         Faculty of Science, Masaryk University
         Kotlarska 2, 611 37 Brno, Czech Republic
=============================================================
         email: janap_at_ncbr.chemi.muni.cz
=============================================================

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu