log started: Mon Sep 23 19:11:54 2019 Log file: ./leap.log >> # >> # ----- leaprc for loading the ff14SB force field >> # ----- NOTE: this is designed for PDB format 3! >> # >> # load atom type hybridizations >> # >> addAtomTypes { >> { "H" "H" "sp3" } >> { "HO" "H" "sp3" } >> { "HS" "H" "sp3" } >> { "H1" "H" "sp3" } >> { "H2" "H" "sp3" } >> { "H3" "H" "sp3" } >> { "H4" "H" "sp3" } >> { "H5" "H" "sp3" } >> { "HW" "H" "sp3" } >> { "HC" "H" "sp3" } >> { "HA" "H" "sp3" } >> { "HP" "H" "sp3" } >> { "HZ" "H" "sp3" } >> { "OH" "O" "sp3" } >> { "OS" "O" "sp3" } >> { "O" "O" "sp2" } >> { "O2" "O" "sp2" } >> { "OP" "O" "sp2" } >> { "OW" "O" "sp3" } >> { "CT" "C" "sp3" } >> { "CX" "C" "sp3" } >> { "C8" "C" "sp3" } >> { "2C" "C" "sp3" } >> { "3C" "C" "sp3" } >> { "CH" "C" "sp3" } >> { "CS" "C" "sp2" } >> { "C" "C" "sp2" } >> { "CO" "C" "sp2" } >> { "C*" "C" "sp2" } >> { "CA" "C" "sp2" } >> { "CB" "C" "sp2" } >> { "CC" "C" "sp2" } >> { "CN" "C" "sp2" } >> { "CM" "C" "sp2" } >> { "CK" "C" "sp2" } >> { "CQ" "C" "sp2" } >> { "CD" "C" "sp2" } >> { "C5" "C" "sp2" } >> { "C4" "C" "sp2" } >> { "CP" "C" "sp2" } >> { "CI" "C" "sp3" } >> { "CJ" "C" "sp2" } >> { "CW" "C" "sp2" } >> { "CV" "C" "sp2" } >> { "CR" "C" "sp2" } >> { "CA" "C" "sp2" } >> { "CY" "C" "sp2" } >> { "C0" "Ca" "sp3" } >> { "MG" "Mg" "sp3" } >> { "N" "N" "sp2" } >> { "NA" "N" "sp2" } >> { "N2" "N" "sp2" } >> { "N*" "N" "sp2" } >> { "NP" "N" "sp2" } >> { "NQ" "N" "sp2" } >> { "NB" "N" "sp2" } >> { "NC" "N" "sp2" } >> { "NT" "N" "sp3" } >> { "NY" "N" "sp2" } >> { "N3" "N" "sp3" } >> { "S" "S" "sp3" } >> { "SH" "S" "sp3" } >> { "P" "P" "sp3" } >> { "LP" "" "sp3" } >> { "EP" "" "sp3" } >> { "F" "F" "sp3" } >> { "Cl" "Cl" "sp3" } >> { "Br" "Br" "sp3" } >> { "I" "I" "sp3" } >> } >> # >> # Load the main parameter set. >> # >> parm10 = loadamberparams parm10.dat Loading parameters: /home/suchetana/amber18/dat/leap/parm/parm10.dat Reading title: PARM99 + frcmod.ff99SB + frcmod.parmbsc0 + OL3 for RNA (UNKNOWN ATOM TYPE: Zn) (UNKNOWN ATOM TYPE: EP) >> frcmod14SB = loadamberparams frcmod.ff14SB Loading parameters: /home/suchetana/amber18/dat/leap/parm/frcmod.ff14SB Reading force field modification type file (frcmod) Reading title: ff14SB protein backbone and sidechain parameters >> # >> # Load main chain and terminating amino acid libraries >> # >> loadOff amino12.lib Loading library: /home/suchetana/amber18/dat/leap/lib/amino12.lib Loading: ALA Loading: ARG Loading: ASH Loading: ASN Loading: ASP Loading: CYM Loading: CYS Loading: CYX Loading: GLH Loading: GLN Loading: GLU Loading: GLY Loading: HID Loading: HIE Loading: HIP Loading: HYP Loading: ILE Loading: LEU Loading: LYN Loading: LYS Loading: MET Loading: PHE Loading: PRO Loading: SER Loading: THR Loading: TRP Loading: TYR Loading: VAL >> loadOff aminoct12.lib Loading library: /home/suchetana/amber18/dat/leap/lib/aminoct12.lib Loading: CALA Loading: CARG Loading: CASN Loading: CASP Loading: CCYS Loading: CCYX Loading: CGLN Loading: CGLU Loading: CGLY Loading: CHID Loading: CHIE Loading: CHIP Loading: CHYP Loading: CILE Loading: CLEU Loading: CLYS Loading: CMET Loading: CPHE Loading: CPRO Loading: CSER Loading: CTHR Loading: CTRP Loading: CTYR Loading: CVAL Loading: NHE Loading: NME >> loadOff aminont12.lib Loading library: /home/suchetana/amber18/dat/leap/lib/aminont12.lib Loading: ACE Loading: NALA Loading: NARG Loading: NASN Loading: NASP Loading: NCYS Loading: NCYX Loading: NGLN Loading: NGLU Loading: NGLY Loading: NHID Loading: NHIE Loading: NHIP Loading: NILE Loading: NLEU Loading: NLYS Loading: NMET Loading: NPHE Loading: NPRO Loading: NSER Loading: NTHR Loading: NTRP Loading: NTYR Loading: NVAL >> >> # >> # Define the PDB name map for the amino acids >> # >> addPdbResMap { >> { 0 "HYP" "HYP" } { 1 "HYP" "CHYP" } >> { 0 "ALA" "NALA" } { 1 "ALA" "CALA" } >> { 0 "ARG" "NARG" } { 1 "ARG" "CARG" } >> { 0 "ASN" "NASN" } { 1 "ASN" "CASN" } >> { 0 "ASP" "NASP" } { 1 "ASP" "CASP" } >> { 0 "CYS" "NCYS" } { 1 "CYS" "CCYS" } >> { 0 "CYX" "NCYX" } { 1 "CYX" "CCYX" } >> { 0 "GLN" "NGLN" } { 1 "GLN" "CGLN" } >> { 0 "GLU" "NGLU" } { 1 "GLU" "CGLU" } >> { 0 "GLY" "NGLY" } { 1 "GLY" "CGLY" } >> { 0 "HID" "NHID" } { 1 "HID" "CHID" } >> { 0 "HIE" "NHIE" } { 1 "HIE" "CHIE" } >> { 0 "HIP" "NHIP" } { 1 "HIP" "CHIP" } >> { 0 "ILE" "NILE" } { 1 "ILE" "CILE" } >> { 0 "LEU" "NLEU" } { 1 "LEU" "CLEU" } >> { 0 "LYS" "NLYS" } { 1 "LYS" "CLYS" } >> { 0 "MET" "NMET" } { 1 "MET" "CMET" } >> { 0 "PHE" "NPHE" } { 1 "PHE" "CPHE" } >> { 0 "PRO" "NPRO" } { 1 "PRO" "CPRO" } >> { 0 "SER" "NSER" } { 1 "SER" "CSER" } >> { 0 "THR" "NTHR" } { 1 "THR" "CTHR" } >> { 0 "TRP" "NTRP" } { 1 "TRP" "CTRP" } >> { 0 "TYR" "NTYR" } { 1 "TYR" "CTYR" } >> { 0 "VAL" "NVAL" } { 1 "VAL" "CVAL" } >> { 0 "HIS" "NHIS" } { 1 "HIS" "CHIS" } >> } >> >> # >> # assume that most often proteins use HIE >> # >> NHIS = NHIE >> HIS = HIE >> CHIS = CHIE > > source leaprc.lipid17 ----- Source: /home/suchetana/amber18/dat/leap/cmd/leaprc.lipid17 ----- Source of /home/suchetana/amber18/dat/leap/cmd/leaprc.lipid17 done >> # >> # Set log file >> # >> logFile leap.log log started: Mon Sep 23 19:11:54 2019 Log file: ./leap.log >> # >> # leaprc for loading the Lipid17 v1.1 force field >> # >> # Load atom type hybridizations >> # >> addAtomTypes { >> { "cA" "C" "sp3" } >> { "cB" "C" "sp2" } >> { "cC" "C" "sp2" } >> { "cD" "C" "sp3" } >> { "cP" "C" "sp3" } >> { "cR" "C" "sp3" } >> { "hA" "H" "sp3" } >> { "hB" "H" "sp3" } >> { "hE" "H" "sp3" } >> { "hL" "H" "sp3" } >> { "hN" "H" "sp3" } >> { "hO" "H" "sp3" } >> { "hR" "H" "sp3" } >> { "hS" "H" "sp3" } >> { "hX" "H" "sp3" } >> { "nA" "N" "sp3" } >> { "nN" "N" "sp2" } >> { "oC" "O" "sp2" } >> { "oH" "O" "sp3" } >> { "oO" "O" "sp2" } >> { "oP" "O" "sp2" } >> { "oR" "O" "sp3" } >> { "oS" "O" "sp3" } >> { "oT" "O" "sp3" } >> { "pA" "P" "sp3" } >> } >> # >> # Load the Lipid17 parameter set. >> # >> lipid17 = loadamberparams lipid17.dat Loading parameters: /home/suchetana/amber18/dat/leap/parm/lipid17.dat Reading title: AMBER Lipid17 v1.1 Force Field, A. Skjevik, C. Dickson, B. Madej, I.R. Gould, R.C. Walker* >> # >> # Load the Lipid17 master lib file. >> # >> loadoff lipid17.lib Loading library: /home/suchetana/amber18/dat/leap/lib/lipid17.lib Loading: AR Loading: CHL Loading: DHA Loading: LAL Loading: MY Loading: OL Loading: PA Loading: PC Loading: PE Loading: PGR Loading: PH- Loading: PS Loading: ST >> > > source leaprc.gaff ----- Source: /home/suchetana/amber18/dat/leap/cmd/leaprc.gaff ----- Source of /home/suchetana/amber18/dat/leap/cmd/leaprc.gaff done >> logFile leap.log log started: Mon Sep 23 19:11:54 2019 Log file: ./leap.log >> # >> # ----- leaprc for loading the general Amber Force field. >> # This file is mostly for use with Antechamber >> # >> # load atom type hybridizations >> # >> addAtomTypes { >> { "h1" "H" "sp3" } >> { "h2" "H" "sp3" } >> { "h3" "H" "sp3" } >> { "h4" "H" "sp3" } >> { "h5" "H" "sp3" } >> { "ha" "H" "sp3" } >> { "hc" "H" "sp3" } >> { "hn" "H" "sp3" } >> { "ho" "H" "sp3" } >> { "hp" "H" "sp3" } >> { "hs" "H" "sp3" } >> { "hw" "H" "sp3" } >> { "hx" "H" "sp3" } >> { "o" "O" "sp2" } >> { "o2" "O" "sp2" } >> { "oh" "O" "sp3" } >> { "os" "O" "sp3" } >> { "ow" "O" "sp3" } >> { "c" "C" "sp2" } >> { "c1" "C" "sp2" } >> { "c2" "C" "sp2" } >> { "c3" "C" "sp3" } >> { "ca" "C" "sp2" } >> { "cc" "C" "sp2" } >> { "cd" "C" "sp2" } >> { "ce" "C" "sp2" } >> { "cf" "C" "sp2" } >> { "cg" "C" "sp2" } >> { "ch" "C" "sp2" } >> { "cp" "C" "sp2" } >> { "cq" "C" "sp2" } >> { "cu" "C" "sp2" } >> { "cv" "C" "sp2" } >> { "cx" "C" "sp2" } >> { "cy" "C" "sp2" } >> { "cz" "C" "sp2" } >> { "n" "N" "sp2" } >> { "n1" "N" "sp2" } >> { "n2" "N" "sp2" } >> { "n3" "N" "sp3" } >> { "n4" "N" "sp3" } >> { "na" "N" "sp2" } >> { "nb" "N" "sp2" } >> { "nc" "N" "sp2" } >> { "nd" "N" "sp2" } >> { "ne" "N" "sp2" } >> { "nf" "N" "sp2" } >> { "nh" "N" "sp2" } >> { "no" "N" "sp2" } >> { "s" "S" "sp2" } >> { "s2" "S" "sp2" } >> { "s3" "S" "sp3" } >> { "s4" "S" "sp3" } >> { "s6" "S" "sp3" } >> { "sh" "S" "sp3" } >> { "ss" "S" "sp3" } >> { "sx" "S" "sp3" } >> { "sy" "S" "sp3" } >> { "p2" "P" "sp2" } >> { "p3" "P" "sp3" } >> { "p4" "P" "sp3" } >> { "p5" "P" "sp3" } >> { "pb" "P" "sp3" } >> { "pc" "P" "sp3" } >> { "pd" "P" "sp3" } >> { "pe" "P" "sp3" } >> { "pf" "P" "sp3" } >> { "px" "P" "sp3" } >> { "py" "P" "sp3" } >> { "f" "F" "sp3" } >> { "cl" "Cl" "sp3" } >> { "br" "Br" "sp3" } >> { "i" "I" "sp3" } >> } >> # >> # Load the general force field parameter set. >> # >> gaff = loadamberparams gaff.dat Loading parameters: /home/suchetana/amber18/dat/leap/parm/gaff.dat Reading title: AMBER General Force Field for organic molecules (Version 1.81, May 2017) (UNKNOWN ATOM TYPE: oq) (UNKNOWN ATOM TYPE: ni) (UNKNOWN ATOM TYPE: nj) (UNKNOWN ATOM TYPE: nk) (UNKNOWN ATOM TYPE: nl) (UNKNOWN ATOM TYPE: nm) (UNKNOWN ATOM TYPE: sp) (UNKNOWN ATOM TYPE: sq) > > RM = loadpdb try_out.pdb Loading PDB file: ./try_out.pdb /home/suchetana/amber18/bin/teLeap: Warning! Name change in pdb file residue 1 ; this residue is split into MY and ASP. /home/suchetana/amber18/bin/teLeap: Note. 1 residues had naming warnings. Thus, there are split residues; residue sequence numbers will not correspond to those in the pdb. -- residue 1: duplicate [ cD] atoms (total 13) /home/suchetana/amber18/bin/teLeap: Warning! Atom names in each residue should be unique. (Same-name atoms are handled by using the first occurrence and by ignoring the rest. Frequently duplicate atom names stem from alternate conformations in the PDB file.) Matching PDB residue names to LEaP variables. Mapped residue ALA, term: Terminal/last, seq. number: 37 to: CALA. Created a new atom named: cC within residue: .R Created a new atom named: oC within residue: .R Created a new atom named: cD within residue: .R Added missing heavy atom: .R.A Added missing heavy atom: .R.A Added missing heavy atom: .R.A Added missing heavy atom: .R.A Added missing heavy atom: .R.A Added missing heavy atom: .R.A Added missing heavy atom: .R.A Added missing heavy atom: .R.A Added missing heavy atom: .R.A Added missing heavy atom: .R.A Added missing heavy atom: .R.A Added missing heavy atom: .R.A Added missing heavy atom: .R.A total atoms in file: 604 Leap added 44 missing atoms according to residue templates: 13 Heavy 31 H / lone pairs The file contained 3 atoms not in residue templates > savepdb RM try_out_1.pdb Writing pdb file: try_out_1.pdb /home/suchetana/amber18/bin/teLeap: Warning! Converting C-terminal residue name to PDB format: CALA -> ALA > saveAmberParm RM 1.top 1.crd Checking Unit. /home/suchetana/amber18/bin/teLeap: Warning! There is a bond of 14.636845 angstroms between: /home/suchetana/amber18/bin/teLeap: Warning! The unperturbed charge of the unit (-1.000001) is not zero. FATAL: Atom .R.A does not have a type. FATAL: Atom .R.A does not have a type. FATAL: Atom .R.A does not have a type. /home/suchetana/amber18/bin/teLeap: Fatal Error! Failed to generate parameters Exiting LEaP: Errors = 1; Warnings = 5; Notes = 1.