AMBER Archive (2003)

Subject: AMBER

From: chandran karunakaran (ckaru2000_at_yahoo.com)
Date: Tue Apr 15 2003 - 18:55:16 CDT


Dear AMBER users,

  I am interested in molecular dynamics of
mutant copper proteins. Before I do that,
I tried with plastocyanin. Thanks for
correcting me to include Cu in HIC
residue. Afterwards I ran and got this
output. Could any one suggest me what
should I do? Sorry for disturbance.
 
  With thanks

  Sincerely yours
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: ./leapp.in
Loading parameters: ./frcmodp.pcy
Reading force field mod type file (frcmod)
Loading Prep file: ./mem.in
Loading Prep file: ./hicu.in
Loading PDB file: ./1plc.H.pdb
  Added missing heavy atom: .R<HIC 100>.A<N 1>
  Added missing heavy atom: .R<HIC 100>.A<CA 3>
  Added missing heavy atom: .R<HIC 100>.A<CB 5>
  Added missing heavy atom: .R<HIC 100>.A<C 17>
  Added missing heavy atom: .R<HIC 100>.A<CG 8>
  Added missing heavy atom: .R<HIC 100>.A<O 18>
  Added missing heavy atom: .R<HIC 100>.A<ND1 9>
  Added missing heavy atom: .R<HIC 100>.A<CD2 15>
  Added missing heavy atom: .R<HIC 100>.A<CE1 11>
  Added missing heavy atom: .R<HIC 100>.A<NE2 13>
  total atoms in file: 1442
  Leap added 18 missing atoms according to residue
templates:
       10 Heavy
       7 H / lone pairs
       1 unknown element
Checking Unit.
WARNING: There is a bond of 3.223311 angstroms
between:
------- .R<HIC 37>.A<CU 10> and .R<MET 92>.A<SD 11>
WARNING: There is a bond of 31.183065 angstroms
between:
------- .R<HIC 100>.A<ND1 9> and .R<HIC 100>.A<CU 10>
WARNING: The unperturbed charge of the unit: -6.000000
is not zero.

 -- ignoring the warnings.

Building topology.
Building atom parameters.
Building bond parameters.
Building angle parameters.
Could not find angle parameter: S - CU - S
Building proper torsion parameters.
Building improper torsion parameters.
old PREP-specified impropers:
 <HIC 37>: -M CA N H
 <HIC 37>: CA +M C O
 <HIC 37>: CE1 CD2 NE2 HE2
 <HIC 37>: CG NE2 CD2 HD2
 <HIC 37>: ND1 NE2 CE1 HE1
 <HIC 37>: ND1 CD2 CG CB
 <HIC 100>: -M CA N H
 <HIC 100>: CA +M C O
 <HIC 100>: CE1 CD2 NE2 HE2
 <HIC 100>: CG NE2 CD2 HD2
 <HIC 100>: ND1 NE2 CE1 HE1
 <HIC 100>: ND1 CD2 CG CB
 total 293 improper torsions applied
 12 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