AMBER

From: chandran karunakaran <ckaru2000.yahoo.com>
Date: Sun, 20 Apr 2003 18:13:32 -0700 (PDT)

Dear AMBER users,

  I am a new user. I ran tleap for
my protein containing two metal atoms
in it.

  I got the following error message.
In my protein, there is no bond between
Zn-C, i dont know why it asks.
I would appreciate if anybody give me
some hint for tackling this.

Thanks in advance

Sincerely Dr.C.Karunakaran

 

Welcome to LEaP!
Sourcing leaprc:
/home/amkumar/soft/amber6/dat/leap/cmd/leaprc
Log file: ./leap.log
Loading parameters:
/home/amkumar/soft/amber6/dat/parm94.dat
Loading library:
/home/amkumar/soft/amber6/dat/leap/lib/all_nucleic94.lib
Loading library:
/home/amkumar/soft/amber6/dat/leap/lib/all_amino94.lib
Loading library:
/home/amkumar/soft/amber6/dat/leap/lib/all_aminoct94.lib
Loading library:
/home/amkumar/soft/amber6/dat/leap/lib/all_aminont94.lib
Loading library:
/home/amkumar/soft/amber6/dat/leap/lib/ions94.lib
Loading library:
/home/amkumar/soft/amber6/dat/leap/lib/water.lib
Sourcing: ./leap.in
Loading parameters: ./frcmod.pcy
Reading force field mod type file (frcmod)
Loading Prep file: ./hicu.in
Loading Prep file: ./asz.in
Loading PDB file: ./1spd.h.pdb
Splitting pdb file residue 48 into HIC and HIE
Splitting pdb file residue 63 into HIC and HIE
 2 residues had naming problems.
 Residues split, so res numbers will not correspond to
pdb.
  Added missing heavy atom: .R<HIC 48>.A<CB 5>
  Added missing heavy atom: .R<HIC 48>.A<C 17>
  Added missing heavy atom: .R<HIC 48>.A<CG 8>
  Added missing heavy atom: .R<HIC 48>.A<O 18>
  Added missing heavy atom: .R<HIC 48>.A<ND1 9>
  Added missing heavy atom: .R<HIC 48>.A<CD2 15>
  Added missing heavy atom: .R<HIC 48>.A<CE1 11>
  Added missing heavy atom: .R<HIC 48>.A<NE2 13>
  Added missing heavy atom: .R<HIE 49>.A<CA 3>
  Added missing heavy atom: .R<HIE 49>.A<N 1>
  Added missing heavy atom: .R<HIC 64>.A<CA 3>
  Added missing heavy atom: .R<HIC 64>.A<CB 5>
  Added missing heavy atom: .R<HIC 64>.A<C 17>
  Added missing heavy atom: .R<HIC 64>.A<CG 8>
  Added missing heavy atom: .R<HIC 64>.A<O 18>
  Added missing heavy atom: .R<HIC 64>.A<ND1 9>
  Added missing heavy atom: .R<HIC 64>.A<CD2 15>
  Added missing heavy atom: .R<HIC 64>.A<CE1 11>
  Added missing heavy atom: .R<HIC 64>.A<NE2 13>
  Added missing heavy atom: .R<HIE 65>.A<N 1>
  Added missing heavy atom: .R<HIC 156>.A<N 1>
  Added missing heavy atom: .R<HIC 156>.A<CA 3>
  Added missing heavy atom: .R<HIC 156>.A<CB 5>
  Added missing heavy atom: .R<HIC 156>.A<C 17>
  Added missing heavy atom: .R<HIC 156>.A<CG 8>
  Added missing heavy atom: .R<HIC 156>.A<O 18>
  Added missing heavy atom: .R<HIC 156>.A<ND1 9>
  Added missing heavy atom: .R<HIC 156>.A<CD2 15>
  Added missing heavy atom: .R<HIC 156>.A<CE1 11>
  Added missing heavy atom: .R<HIC 156>.A<NE2 13>
  Added missing heavy atom: .R<ASZ 157>.A<C 11>
  Added missing heavy atom: .R<ASZ 157>.A<CA 3>
  Added missing heavy atom: .R<ASZ 157>.A<O 12>
  Added missing heavy atom: .R<ASZ 157>.A<N 1>
  Added missing heavy atom: .R<ASZ 157>.A<CB 5>
  Added missing heavy atom: .R<ASZ 157>.A<CG 8>
  Added missing heavy atom: .R<ASZ 157>.A<OD1 9>
  Added missing heavy atom: .R<ASZ 157>.A<OD2 10>
  total atoms in file: 2190
  Leap added 67 missing atoms according to residue
templates:
       38 Heavy
       25 H / lone pairs
       4 unknown element
bond: Argument #1 is type String must be of type:
[atom]
usage: bond <atom1> <atom2> [order]
bond: Argument #1 is type String must be of type:
[atom]
usage: bond <atom1> <atom2> [order]
bond: Argument #1 is type String must be of type:
[atom]
usage: bond <atom1> <atom2> [order]
bond: Argument #1 is type String must be of type:
[atom]
usage: bond <atom1> <atom2> [order]
bond: Argument #1 is type String must be of type:
[atom]
usage: bond <atom1> <atom2> [order]
bond: Argument #1 is type String must be of type:
[atom]
usage: bond <atom1> <atom2> [order]
bond: Argument #1 is type String must be of type:
[atom]
usage: bond <atom1> <atom2> [order]
addIons: 1st Ion & target are the same charge:
     can't neutralize.
  Solute vdw bounding box: 48.390 37.633
44.449
  Total bounding box for atom centers: 54.390 43.633
50.449
  Solvent unit box: 18.774 18.774
18.774
  Total vdw box size: 57.635 46.888
53.455 angstroms.
  Volume: 144455.615 A^3
  Total mass 63469.948 amu, Density 0.730 g/cc
  Added 2596 residues.
Checking Unit.
WARNING: There is a bond of 0.353458 angstroms
between:
------- .R<HIC 48>.A<C 17> and .R<HIE 49>.A<N 1>
WARNING: There is a bond of 26.077951 angstroms
between:
------- .R<HIC 156>.A<C 17> and .R<ASZ 157>.A<N 1>
WARNING: There is a bond of 3.248675 angstroms
between:
------- .R<HIC 156>.A<NE2 13> and .R<HIC 156>.A<CD2
15>
WARNING: There is a bond of 28.702998 angstroms
between:
------- .R<HIC 156>.A<ND1 9> and .R<HIC 156>.A<CU 10>
WARNING: The unperturbed charge of the unit: 2.000000
is not zero.

 -- ignoring the warnings.

Building topology.
Building atom parameters.
Building bond parameters.
Could not find bond parameter for: C - ZN
Could not find bond parameter for: C - ZN
Building angle parameters.
Could not find angle parameter: ZN - C - N
Could not find angle parameter: O - C - ZN
Could not find angle parameter: CT - C - ZN
Could not find angle parameter: O - C - ZN
Could not find angle parameter: CT - C - ZN
Building proper torsion parameters.
Building improper torsion parameters.
old PREP-specified impropers:
 <HIC 46>: -M CA N H
 <HIC 46>: CA +M C O
 <HIC 46>: CE1 CD2 NE2 HE2
 <HIC 46>: CG NE2 CD2 HD2
 <HIC 46>: ND1 NE2 CE1 HE1
 <HIC 46>: ND1 CD2 CG CB
 <HIC 48>: -M CA N H
 <HIC 48>: CA +M C O
 <HIC 48>: CE1 CD2 NE2 HE2
 <HIC 48>: CG NE2 CD2 HD2
 <HIC 48>: ND1 NE2 CE1 HE1
 <HIC 48>: ND1 CD2 CG CB
 <HIC 64>: -M CA N H
 <HIC 64>: CA +M C O
 <HIC 64>: CE1 CD2 NE2 HE2
 <HIC 64>: CG NE2 CD2 HD2
 <HIC 64>: ND1 NE2 CE1 HE1
 <HIC 64>: ND1 CD2 CG CB
 <ASZ 85>: -M CA N H
 <ASZ 85>: CA +M C O
 <ASZ 85>: CB OD1 CG OD2
 <HIC 156>: -M CA N H
 <HIC 156>: CA +M C O
 <HIC 156>: CE1 CD2 NE2 HE2
 <HIC 156>: CG NE2 CD2 HD2
 <HIC 156>: ND1 NE2 CE1 HE1
 <HIC 156>: ND1 CD2 CG CB
 <ASZ 157>: -M CA N H
 <ASZ 157>: CA +M C O
 <ASZ 157>: CB OD1 CG OD2
 total 448 improper torsions applied
 30 improper torsions in old prep form
Building H-Bond parameters.
Parameter file was not saved.


__________________________________________________
Do you Yahoo!?
The New Yahoo! Search - Faster. Easier. Bingo
http://search.yahoo.com
Received on Mon Apr 21 2003 - 02:53:01 PDT
Custom Search