Re: [AMBER] Hydrogen bond analysis problem

From: Sumra Wajid Abbasi 30-FBAS/MSBI/F09 <sumra.msbi30.iiu.edu.pk>
Date: Thu, 19 Sep 2013 11:22:05 +0500

Dear Daniel

Even after making corrections in donor and acceptor mask definition same
error occurred

Regards


  \-/
  -/- PTRAJ: a utility for processing trajectory files
  /-\
  \-/ Version: "AMBER 12.0 integrated" (4/2012)
  -/- Executable is: "ptraj"
  /-\ Running on 1 processor(s)
  \-/ Residue labels:

 MET GLY SER SER GLU ASP GLN ALA TYR ARG
 LEU LEU ASN ASP TYR ALA ASN GLY PHE MET
 VAL SER GLN VAL LEU PHE ALA ALA CYS GLU
 LEU GLY VAL PHE ASP LEU LEU ALA GLU ALA
 PRO GLY PRO LEU ASP VAL ALA ALA VAL ALA
 ALA GLY VAL ARG ALA SER ALA HIE GLY THR
 GLU LEU LEU LEU ASP ILE CYS VAL SER LEU
 LYS LEU LEU LYS VAL GLU THR ARG GLY GLY
 LYS ALA PHE TYR ARG ASN THR GLU LEU SER
 SER ASP TYR LEU THR THR VAL SER PRO THR
 SER GLN CYS SER MET LEU LYS TYR MET GLY
 ARG THR SER TYR ARG CYS TRP GLY HIE LEU
 ALA ASP ALA VAL ARG GLU GLY ARG ASN GLN
 TYR LEU GLU THR PHE GLY VAL PRO ALA GLU
 GLU LEU PHE THR ALA ILE TYR ARG SER GLU
 GLY GLU ARG LEU GLN PHE MET GLN ALA LEU
 GLN GLU VAL TRP SER VAL ASN GLY ARG SER
 VAL LEU THR ALA PHE ASP LEU SER VAL PHE
 PRO LEU MET CYS ASP LEU GLY GLY THR TRP
 ILE LYS LEU GLU THR ILE ILE LEU SER LYS
 LEU SER GLN GLY GLN LYS THR LYS HIE ARG
 VAL PHE SER LEU ILE GLY GLY ALA GLY ALA
 LEU ALA LYS GLU CYS MET SER LEU TYR PRO
 GLY CYS LYS ILE THR VAL PHE ASP ILE PRO
 GLU VAL VAL TRP THR ALA LYS GLN HIE PHE
 SER PHE GLN GLU GLU GLU GLN ILE ASP PHE
 GLN GLU GLY ASP PHE PHE LYS ASP PRO LEU
 PRO GLU ALA ASP LEU TYR ILE LEU ALA ARG
 VAL LEU HID ASP TRP ALA ASP GLY LYS CYS
 SER HIE LEU LEU GLU ARG ILE TYR HID THR
 CYS LYS PRO GLY GLY GLY ILE LEU VAL ILE
 GLU SER LEU LEU ASP GLU ASP ARG ARG GLY
 PRO LEU LEU THR GLN LEU TYR SER LEU ASN
 MET LEU VAL GLN THR GLU GLY GLN GLU ARG
 THR PRO THR HIE TYR HIE MET LEU LEU SER
 SER ALA GLY PHE ARG ASP PHE GLN PHE LYS
 LYS THR GLY ALA ILE TYR ASP ALA ILE LEU
 ALA ARG LYS UNK WAT WAT WAT WAT WAT WAT
 WAT WAT WAT WAT WAT WAT WAT WAT WAT WAT
 ...
 WAT WAT WAT


PTRAJ: Processing input from "STDIN" ...

PTRAJ: trajin prod_1-8.binpos
  Checking coordinates: prod_1-8.binpos

PTRAJ: donor mask :198.O
Mask [:198.O] represents 1 atoms

PTRAJ: donor mask :198.N
Mask [:198.N] represents 1 atoms

PTRAJ: donor mask :196.O
Mask [:196.O] represents 1 atoms

PTRAJ: donor mask :196.N
Mask [:196.N] represents 1 atoms

PTRAJ: donor mask :197.N
Mask [:197.N] represents 1 atoms

PTRAJ: donor mask :249.O4
Mask [:249.O4] represents 0 atoms !!!NO ATOMS DETECTED!!!
WARNING in ptraj, donor: No atoms selected (:249.O4), ignoring...

PTRAJ: acceptor mask :374.N1 :196.H
Mask [:374.N1] represents 1 atoms
Mask [:196.H] represents 1 atoms

PTRAJ: acceptor mask :374.O1 :210.HH12
Mask [:374.O1] represents 1 atoms
Mask [:210.HH12] represents 1 atoms

PTRAJ: acceptor mask :374.O3 :210.HH12
Mask [:374.O3] represents 1 atoms
Mask [:210.HH12] represents 1 atoms

PTRAJ: acceptor mask :374.O3 :210.HH22
Mask [:374.O3] represents 1 atoms
Mask [:210.HH22] represents 1 atoms

PTRAJ: acceptor mask :374.O2 :210.HH11
Mask [:374.O2] represents 1 atoms
Mask [:210.HH11] represents 1 atoms

PTRAJ: acceptor mask :374.O3 :210.HH11
Mask [:374.O3] represents 1 atoms
Mask [:210.HH11] represents 1 atoms

PTRAJ: acceptor mask :374.O4 :98.H
Mask [:374.O4] represents 1 atoms
Mask [:98.H] represents 1 atoms

PTRAJ: acceptor mask :374.N2 :198.H
Mask [:374.N2] represents 1 atoms
Mask [:198.H] represents 1 atoms

PTRAJ: acceptor mask :374.N :198.H18
Mask [:374.N] represents 0 atoms !!!NO ATOMS DETECTED!!!
Mask [:198.H18] represents 0 atoms !!!NO ATOMS DETECTED!!!
WARNING in ptraj, acceptor: No heavy atom was selected (:374.N), ignoring...
WARNING in ptraj, acceptor: No hydrogen atom was selected (:198.H18),
ignoring...

PTRAJ: acceptor mask :374.N1 :196.H
Mask [:374.N1] represents 1 atoms
Mask [:196.H] represents 1 atoms

PTRAJ: acceptor mask :374.N1 :195.HG1
Mask [:374.N1] represents 1 atoms
Mask [:195.HG1] represents 1 atoms

PTRAJ: acceptor mask :374.O4 :195.HG1
Mask [:374.O4] represents 1 atoms
Mask [:195.HG1] represents 1 atoms

PTRAJ: acceptor mask :374.O2 :211.H
Mask [:374.O2] represents 1 atoms
Mask [:211.H] represents 1 atoms

PTRAJ: acceptor mask :374.O1 :198.H
Mask [:374.O1] represents 1 atoms
Mask [:198.H] represents 1 atoms

PTRAJ: acceptor mask :374.O4 :249.HE2
Mask [:374.O4] represents 1 atoms
Mask [:249.HE2] represents 1 atoms

PTRAJ: acceptor mask :374.N2 :196.H
Mask [:374.N2] represents 1 atoms
Mask [:196.H] represents 1 atoms

PTRAJ: hbond print .05 series hbt
[No output trajectory specified (trajout)]
  prod_1-8.binpos: 16004 frames.

PTRAJ: Successfully read the input file.
       Coordinate processing will occur on 16004 frames.
       Summary of I/O and actions follows:

INPUT COORDINATE FILES
  File (prod_1-8.binpos) is a BINPOS file with 16004 sets

NO OUTPUT COORDINATE FILE WAS SPECIFIED

ACTIONS
  1> HBOND saved to series hbt, inter-solvent hbonds will not be included,

      data will be sorted, intra-residue interactions will NOT be included,
      Distance cutoff is 3.00 angstroms, angle cutoff is 120.00 degrees
      Hydrogen bond information will be dumped for occupancies > 0.05
      Estimated memory usage for this hbond call: 13.74 MB
      donors: 5 acceptors: 15 frames: 16004



Processing BINPOS file prod_1-8.binpos

 1% ............ 25% ............ 50% ............ 75% ............ 100%


PTRAJ: Successfully read in 16004 sets and processed 16004 sets.

Dumping accumulated results (if any)

  HBOND SUMMARY:
      Data was saved to series hbt,
      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.05

  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.atom atom# :res.atom atom# :res.atom %occupied
distance angle lifetime maxocc
 ---------------- ---------------------------------
-------------------------------------


Sumra Wajid Abbasi
PhD Scholar
Computational Biology Lab, National Center for Bioinformatics
Quaid-e-Azam University, Islamabad-45320 Pakistan.


On Wed, Sep 18, 2013 at 8:16 PM, Daniel Roe <daniel.r.roe.gmail.com> wrote:

> Hi,
>
> On Wed, Sep 18, 2013 at 3:50 AM, Sumra Wajid Abbasi 30-FBAS/MSBI/F09
> <sumra.msbi30.iiu.edu.pk> wrote:
> > acceptor mask :196.H :374.N1
> > acceptor mask :210.HH12 :374.O1
> > acceptor mask :210.HH11 :374.O2
> > acceptor mask :98.H :374.O4
>
> The 'acceptor mask' expects the heavy atom first, so you should have:
>
> acceptor mask :374.N1 :196.H
>
> etc. The ordering will affect how the angle is calculated.
>
> > PTRAJ: donor mask :LEU198.O
> > Mask [:LEU198.O] represents 0 atoms !!!NO ATOMS DETECTED!!!
> > WARNING in ptraj, donor: No atoms selected (:LEU198.O), ignoring...
>
> This is not a valid mask. You probably just want ':198.O'.
>
> Hope this helps,
>
> -Dan
>
> --
> -------------------------
> Daniel R. Roe, PhD
> Department of Medicinal Chemistry
> University of Utah
> 30 South 2000 East, Room 201
> Salt Lake City, UT 84112-5820
> http://home.chpc.utah.edu/~cheatham/
> (801) 587-9652
> (801) 585-9119 (Fax)
>
> _______________________________________________
> 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
Received on Wed Sep 18 2013 - 23:30:02 PDT
Custom Search