AMBER Archive (2005)

Subject: AMBER: hbond detailed analysis over time?

From: pascal.baillod_at_epfl.ch
Date: Mon May 23 2005 - 13:04:10 CDT


Dear users and developers,

I am using ptraj to analyse hbonds in an rna system. I get the mean hbond
occupancy, distance and angle for the donors and acceptors that I specify, as
well as the schematic time evolution. I would be interested in getting a
detailed time evolution printed to a separate file, which must somehow be
possible.. Here below are examples of an input I use and the respective output I
obtained. Thank you very much in advance!

Pascal Baillod

-----INPUT-----------------------------------------------------------------

trajin ../traj.crd
acceptor mask :12_at_N3 :12_at_H3
donor mask :1_at_O6
hbond series hb time 0.75 distance 3.0 out 12Un3-1Go6.out

-----OUTPUT-----------------------------------------------------------------
  

HBOND SUMMARY:
      Data was saved to series hb, output to file 12Un3-1Go6.out, time interval
is 0.75
 
      data was sorted, intra-residue interactions are NOT included,
      Distance cutoff is 3.00 angstroms, angle cutoff is 120.00 degrees
      Hydrogen bond information dumped for occupancies > 0.00
 
  Dumping schematic of time series after each h-bond, key follows:
   | . - o x * @ |
      0-5% 5-20% 20-40% 40-60% 60-80% 80-95% 95-100% occupancy
 
        DONOR ACCEPTORH ACCEPTOR
  atom# :res_at_atom atom# :res_at_atom atom# :res_at_atom %occupied distance
angle lifetime maxocc
| 17 :1_at_O6 | 379 :12_at_H3 378 :12_at_N3 | 14.27 2.874 ( 0.08)
27.14 (13.24) 1.8 ( 2.4) 19 |--..-... |

*******************************************************************************
Pascal Baillod (PhD student)
*******************************************************************************
Swiss Federal Institute of Technology EPFL Tel: +41-(0)21-693-0322
Institute of Chemical Sciences and Engineering , Fax: +41-(0)21-693-0320
Laboratory of Computational Chemistry and Biochemistry pascal.baillod_at_epfl.ch
Room BCH 4121, Avenue Forel, http://lcbcpc21.epfl.ch
CH-1015 Lausanne
*******************************************************************************
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu