Log file: ./leap.log >> # >> # ----- leaprc for loading the ff03 (Duan et al.) force field >> # >> # 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" } >> { "OH" "O" "sp3" } >> { "OS" "O" "sp3" } >> { "O" "O" "sp2" } >> { "O2" "O" "sp2" } >> { "OW" "O" "sp3" } >> { "CT" "C" "sp3" } >> { "CH" "C" "sp3" } >> { "C2" "C" "sp3" } >> { "C3" "C" "sp3" } >> { "C" "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" } >> { "CE" "C" "sp2" } >> { "CF" "C" "sp2" } >> { "CP" "C" "sp2" } >> { "CI" "C" "sp2" } >> { "CJ" "C" "sp2" } >> { "CW" "C" "sp2" } >> { "CV" "C" "sp2" } >> { "CR" "C" "sp2" } >> { "CA" "C" "sp2" } >> { "CY" "C" "sp2" } >> { "C0" "C" "sp2" } >> { "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" } >> { "N3" "N" "sp3" } >> { "S" "S" "sp3" } >> { "SH" "S" "sp3" } >> { "P" "P" "sp3" } >> { "LP" "" "sp3" } >> { "F" "F" "sp3" } >> { "CL" "Cl" "sp3" } >> { "BR" "Br" "sp3" } >> { "I" "I" "sp3" } >> { "FE" "Fe" "sp3" } >> # things should be there >> { "IM" "Cl" "sp3" } >> { "IP" "Na" "sp3" } >> { "Li" "Li" "sp3" } >> { "K" "K" "sp3" } >> { "Rb" "Rb" "sp3" } >> { "Cs" "Cs" "sp3" } >> { "Zn" "Zn" "sp3" } >> { "IB" "Na" "sp3" } >> # "new" types >> { "H0" "H" "sp3" } >> >> } >> # >> # Load the main parameter set. >> # >> parm99 = loadamberparams parm99.dat Loading parameters: /home/sbl/Documents/amber16/dat/leap/parm/parm99.dat Reading title: PARM99 for DNA,RNA,AA, organic molecules, TIP3P wat. Polariz.& LP incl.02/04/99 (UNKNOWN ATOM TYPE: HZ) (UNKNOWN ATOM TYPE: LP) >> frcmod03 = loadamberparams frcmod.ff03 Loading parameters: /home/sbl/Documents/amber16/dat/leap/parm/frcmod.ff03 Reading force field modification type file (frcmod) Reading title: Duan et al ff03 phi psi torsions >> >> # >> # Load water and ions >> # >> loadOff ions94.lib Loading library: /home/sbl/Documents/amber16/dat/leap/lib/ions94.lib Loading: CIO Loading: Cl- Loading: Cs+ Loading: IB Loading: K+ Loading: Li+ Loading: MG2 Loading: Na+ Loading: Rb+ >> loadOff solvents.lib Loading library: /home/sbl/Documents/amber16/dat/leap/lib/solvents.lib Loading: CHCL3BOX Loading: DC4 Loading: FB3 Loading: FB3BOX Loading: FB4 Loading: FB4BOX Loading: MEOHBOX Loading: NMABOX Loading: OPC Loading: OPCBOX Loading: PL3 Loading: POL3BOX Loading: QSPCFWBOX Loading: SPC Loading: SPCBOX Loading: SPCFWBOX Loading: SPF Loading: SPG Loading: T4E Loading: TIP3PBOX Loading: TIP3PFBOX Loading: TIP4PBOX Loading: TIP4PEWBOX Loading: TIP5PBOX Loading: TP3 Loading: TP4 Loading: TP5 Loading: TPF >> HOH = TP3 >> WAT = TP3 >> >> # >> # Load the existing 94 lib for nucleic acids and C- and N- terminal groups >> # >> # >> # Load DNA/RNA libraries >> # >> loadOff all_nucleic94.lib Loading library: /home/sbl/Documents/amber16/dat/leap/lib/all_nucleic94.lib Loading: DA Loading: DA3 Loading: DA5 Loading: DAN Loading: DC Loading: DC3 Loading: DC5 Loading: DCN Loading: DG Loading: DG3 Loading: DG5 Loading: DGN Loading: DT Loading: DT3 Loading: DT5 Loading: DTN Loading: RA Loading: RA3 Loading: RA5 Loading: RAN Loading: RC Loading: RC3 Loading: RC5 Loading: RCN Loading: RG Loading: RG3 Loading: RG5 Loading: RGN Loading: RU Loading: RU3 Loading: RU5 Loading: RUN >> # >> # Load main chain and terminating >> # amino acid libraries. >> # >> loadOff all_aminoct94.lib Loading library: /home/sbl/Documents/amber16/dat/leap/lib/all_aminoct94.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: 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 all_aminont94.lib Loading library: /home/sbl/Documents/amber16/dat/leap/lib/all_aminont94.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 >> loadOff all_amino03.lib Loading library: /home/sbl/Documents/amber16/dat/leap/lib/all_amino03.lib Loading: ACE 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: ILE Loading: LEU Loading: LYN Loading: LYS Loading: MET Loading: NHE Loading: NME Loading: PHE Loading: PRO Loading: SER Loading: THR Loading: TRP Loading: TYR Loading: VAL >> >> # >> # Define the PDB name map for the amino acids and DNA. >> # >> addPdbResMap { >> { 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" } >> { 0 "GUA" "DG5" } { 1 "GUA" "DG3" } { "GUA" "DG" } >> { 0 "ADE" "DA5" } { 1 "ADE" "DA3" } { "ADE" "DA" } >> { 0 "CYT" "DC5" } { 1 "CYT" "DC3" } { "CYT" "DC" } >> { 0 "THY" "DT5" } { 1 "THY" "DT3" } { "THY" "DT" } >> { 0 "G" "DG5" } { 1 "G" "DG3" } { "G" "DG" } { "GN" "DGN" } >> { 0 "A" "DA5" } { 1 "A" "DA3" } { "A" "DA" } { "AN" "DAN" } >> { 0 "C" "DC5" } { 1 "C" "DC3" } { "C" "DC" } { "CN" "DCN" } >> { 0 "T" "DT5" } { 1 "T" "DT3" } { "T" "DT" } { "TN" "DTN" } >> { 0 "C5" "DC5" } >> { 0 "G5" "DG5" } >> { 0 "A5" "DA5" } >> { 0 "T5" "DT5" } >> { 1 "C3" "DC3" } >> { 1 "G3" "DG3" } >> { 1 "A3" "DA3" } >> { 1 "T3" "DT3" } >> >> } >> >> addPdbAtomMap { >> { "O5*" "O5'" } >> { "C5*" "C5'" } >> { "C4*" "C4'" } >> { "O4*" "O4'" } >> { "C3*" "C3'" } >> { "O3*" "O3'" } >> { "C2*" "C2'" } >> { "C1*" "C1'" } >> { "C5M" "C7" } >> { "H1*" "H1'" } >> { "H2*1" "H2'1" } >> { "H2*2" "H2'2" } >> { "H3*" "H3'" } >> { "H4*" "H4'" } >> { "H5*1" "H5'1" } >> { "H5*2" "H5'2" } >> # old ff atom names -> new >> { "O1'" "O4'" } >> { "OA" "O1P" } >> { "OB" "O2P" } >> } >> >> >> # >> # assumed that most often proteins use HIE >> # >> NHIS = NHIE >> HIS = HIE >> CHIS = CHIE > > loadoff all_prot_nucleic10.lib Loading library: /home/sbl/Documents/amber16/dat/leap/lib/all_prot_nucleic10.lib Loading: A3E Loading: A3P Loading: A5E Loading: A5P Loading: AE Loading: AP Loading: C3E Loading: C3P Loading: C5E Loading: C5P Loading: CE Loading: CP Loading: D3A Loading: D3AE Loading: D3C Loading: D3CE Loading: D3G Loading: D3GE Loading: D3T Loading: D3TE Loading: D5A Loading: D5AE Loading: D5C Loading: D5CE Loading: D5G Loading: D5GE Loading: D5T Loading: D5TE Loading: DAE Loading: DAP Loading: DCE Loading: DCP Loading: DGD Loading: DGE Loading: DTD Loading: DTE Loading: G3D Loading: G3E Loading: G5D Loading: G5E Loading: GD Loading: GE Loading: U3D Loading: U3E Loading: U5D Loading: U5E Loading: UD Loading: UE > loadAmberParams frcmod.protonated_nucleic Loading parameters: /home/sbl/Documents/amber16/dat/leap/parm/frcmod.protonated_nucleic Reading force field modification type file (frcmod) Reading title: Force field modifications for protonated nucleic acids (UNKNOWN ATOM TYPE: NX) (UNKNOWN ATOM TYPE: OX) (UNKNOWN ATOM TYPE: HD) > dna1=loadpdb motif.pdb Loading PDB file: ./motif.pdb Enter zPdbReadScan from call depth 0. Exit zPdbReadScan from call depth 0. Matching PDB residue names to LEaP variables. (Residue 0: DC, Terminal/beginning, was not found in name map.) (Residue 1: DC, Nonterminal, was not found in name map.) (Residue 2: CP, Terminal/last, was not found in name map.) (Residue 3: DC, Terminal/beginning, was not found in name map.) (Residue 4: DC, Nonterminal, was not found in name map.) (Residue 5: DC, Terminal/last, was not found in name map.) (Residue 6: DC, Terminal/beginning, was not found in name map.) (Residue 7: DC, Nonterminal, was not found in name map.) (Residue 8: DC, Terminal/last, was not found in name map.) (Residue 9: DA, Terminal/beginning, was not found in name map.) (Residue 10: DA, Nonterminal, was not found in name map.) (Residue 11: AP, Terminal/last, was not found in name map.) Joining DC - DC Joining DC - CP Added missing heavy atom: .R.A Joining DC - DC Joining DC - DC Joining DC - DC Joining DC - DC Joining DA - DA Joining DA - AP Added missing heavy atom: .R.A total atoms in file: 256 Leap added 114 missing atoms according to residue templates: 2 Heavy 112 H / lone pairs > addions dna1 Na+ 0 10 Na+ ions required to neutralize. Adding 10 counter ions to "dna1" using 1A grid Used default radius 1.50 for 5 atoms Total solute charge: -10.00 Max atom radius: 2.10 Grid extends from solute vdw + 4.15 to 10.25 Box: enclosing: 11.90 37.25 -19.72 54.81 82.66 23.21 sized: 75.90 101.25 44.28 edge: 64.00 Resolution: 1.00 Angstrom. Tree depth: 6 grid build: 0 sec Volume = 7.35% of box, grid points 19279 (no solvent present) Calculating grid charges charges: 0 sec Placed Na+ in dna1 at (31.90, 68.25, -6.72). Placed Na+ in dna1 at (23.90, 66.25, 4.28). Placed Na+ in dna1 at (44.90, 56.25, -9.72). Placed Na+ in dna1 at (31.90, 67.25, 12.28). Placed Na+ in dna1 at (41.90, 66.25, -11.72). Placed Na+ in dna1 at (29.90, 47.25, -6.72). Placed Na+ in dna1 at (42.90, 73.25, 3.28). Placed Na+ in dna1 at (23.90, 56.25, 6.28). Placed Na+ in dna1 at (20.90, 66.25, -6.72). Placed Na+ in dna1 at (40.90, 49.25, 0.28). Done adding ions. > solvatebox dna1 TIP3PBOX 10 (using default radius 1.500000 for C5') (using default radius 1.500000 for C6) (using default radius 1.500000 for C5) (using default radius 1.500000 for C5') (using default radius 1.500000 for C8) Solute vdw bounding box: 27.736 29.736 28.195 Total bounding box for atom centers: 47.736 49.736 48.195 Solvent unit box: 18.774 18.774 18.774 The number of boxes: x= 3 y= 3 z= 3 Adding box at: x=0 y=0 z=0 Center of solvent box is: 18.774349, 18.774349, 18.774349 Adding box at: x=0 y=0 z=1 Center of solvent box is: 18.774349, 18.774349, 0.000000 Adding box at: x=0 y=0 z=2 Center of solvent box is: 18.774349, 18.774349, -18.774349 Adding box at: x=0 y=1 z=0 Center of solvent box is: 18.774349, 0.000000, 18.774349 Adding box at: x=0 y=1 z=1 Center of solvent box is: 18.774349, 0.000000, 0.000000 Adding box at: x=0 y=1 z=2 Center of solvent box is: 18.774349, 0.000000, -18.774349 Adding box at: x=0 y=2 z=0 Center of solvent box is: 18.774349, -18.774349, 18.774349 Adding box at: x=0 y=2 z=1 Center of solvent box is: 18.774349, -18.774349, 0.000000 Adding box at: x=0 y=2 z=2 Center of solvent box is: 18.774349, -18.774349, -18.774349 Adding box at: x=1 y=0 z=0 Center of solvent box is: 0.000000, 18.774349, 18.774349 Adding box at: x=1 y=0 z=1 Center of solvent box is: 0.000000, 18.774349, 0.000000 Adding box at: x=1 y=0 z=2 Center of solvent box is: 0.000000, 18.774349, -18.774349 Adding box at: x=1 y=1 z=0 Center of solvent box is: 0.000000, 0.000000, 18.774349 Adding box at: x=1 y=1 z=1 Center of solvent box is: 0.000000, 0.000000, 0.000000 Adding box at: x=1 y=1 z=2 Center of solvent box is: 0.000000, 0.000000, -18.774349 Adding box at: x=1 y=2 z=0 Center of solvent box is: 0.000000, -18.774349, 18.774349 Adding box at: x=1 y=2 z=1 Center of solvent box is: 0.000000, -18.774349, 0.000000 Adding box at: x=1 y=2 z=2 Center of solvent box is: 0.000000, -18.774349, -18.774349 Adding box at: x=2 y=0 z=0 Center of solvent box is: -18.774349, 18.774349, 18.774349 Adding box at: x=2 y=0 z=1 Center of solvent box is: -18.774349, 18.774349, 0.000000 Adding box at: x=2 y=0 z=2 Center of solvent box is: -18.774349, 18.774349, -18.774349 Adding box at: x=2 y=1 z=0 Center of solvent box is: -18.774349, 0.000000, 18.774349 Adding box at: x=2 y=1 z=1 Center of solvent box is: -18.774349, 0.000000, 0.000000 Adding box at: x=2 y=1 z=2 Center of solvent box is: -18.774349, 0.000000, -18.774349 Adding box at: x=2 y=2 z=0 Center of solvent box is: -18.774349, -18.774349, 18.774349 Adding box at: x=2 y=2 z=1 Center of solvent box is: -18.774349, -18.774349, 0.000000 Adding box at: x=2 y=2 z=2 Center of solvent box is: -18.774349, -18.774349, -18.774349 (using default radius 1.500000 for C5') (using default radius 1.500000 for C6) (using default radius 1.500000 for C5) (using default radius 1.500000 for C5') (using default radius 1.500000 for C8) Total vdw box size: 50.985 52.874 50.629 angstroms. Volume: 136486.948 A^3 Mass > 61331.244 amu, Density > 0.746 g/cc (type - hence mass - of one or more atoms could not be found) Added 3197 residues. > saveamberparm dna1 motif.prmtop motif.inpcrd Checking Unit. WARNING: There is a bond of 6.592593 angstroms between: ------- .R.A and .R.A WARNING: There is a bond of 6.140893 angstroms between: ------- .R.A and .R.A WARNING: There is a bond of 5.742315 angstroms between: ------- .R.A and .R.A WARNING: There is a bond of 8.164783 angstroms between: ------- .R.A and .R.A WARNING: There is a bond of 7.927619 angstroms between: ------- .R.A and .R.A WARNING: There is a bond of 7.738285 angstroms between: ------- .R.A and .R.A -- ignoring the warnings. Building topology. Building atom parameters. For atom: .R.A Could not find vdW (or other) parameters for type: CI For atom: .R.A Could not find vdW (or other) parameters for type: C4 For atom: .R.A Could not find vdW (or other) parameters for type: C4 For atom: .R.A Could not find vdW (or other) parameters for type: CI For atom: .R.A Could not find vdW (or other) parameters for type: C5 Parameter file was not saved.