log started: Fri Oct 23 10:16:02 2015 Log file: ./leap.log >> # >> # ----- leaprc for loading the ff99SB (Hornak & Simmerling) force field >> # using parmbsc0 for nucleic acid. >> # ----- 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" } >> { "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" "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" } >> { "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" } >> { "EP" "" "sp3" } >> # glycam >> { "OG" "O" "sp3" } >> { "OL" "O" "sp3" } >> { "AC" "C" "sp3" } >> { "EC" "C" "sp3" } >> } >> # >> # Load the main parameter set. >> # >> parm99 = loadamberparams parm99.dat Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/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: IM) (UNKNOWN ATOM TYPE: Li) (UNKNOWN ATOM TYPE: IP) (UNKNOWN ATOM TYPE: K) (UNKNOWN ATOM TYPE: Rb) (UNKNOWN ATOM TYPE: Cs) (UNKNOWN ATOM TYPE: Zn) (UNKNOWN ATOM TYPE: IB) (UNKNOWN ATOM TYPE: LP) >> frcmod99SB = loadamberparams frcmod.ff99SB Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.ff99SB Reading force field modification type file (frcmod) Reading title: Modification/update of parm99.dat (Hornak & Simmerling) >> frcmodbsc0 = loadamberparams frcmod.parmbsc0 Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.parmbsc0 Reading force field modification type file (frcmod) Reading title: modifies parm99 for nucleic acids: see http://mmb.pcb.ub.es/PARMBSC0/ >> # >> # Load DNA/RNA libraries >> # >> loadOff all_nucleic94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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 (i.e. ff94 libs) >> # >> loadOff all_amino94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/all_amino94.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: ILE Loading: LEU Loading: LYN Loading: LYS Loading: MET Loading: PHE Loading: PRO Loading: SER Loading: THR Loading: TRP Loading: TYR Loading: VAL >> loadOff all_aminoct94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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: /oasis/projects/nsf/slc216/droe/GIT/amber/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 DNA_CI.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/DNA_CI.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 >> # >> # Load water and ions >> # >> loadOff ions94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/solvents.lib Loading: CHCL3BOX Loading: DC4 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 >> >> # >> # 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" "RG5" } { 1 "G" "RG3" } { "G" "RG" } { "GN" "RGN" } >> { 0 "A" "RA5" } { 1 "A" "RA3" } { "A" "RA" } { "AN" "RAN" } >> { 0 "C" "RC5" } { 1 "C" "RC3" } { "C" "RC" } { "CN" "RCN" } >> { 0 "U" "RU5" } { 1 "U" "RU3" } { "U" "RU" } { "UN" "RUN" } >> { 0 "DG" "DG5" } { 1 "DG" "DG3" } >> { 0 "DA" "DA5" } { 1 "DA" "DA3" } >> { 0 "DC" "DC5" } { 1 "DC" "DC3" } >> { 0 "DT" "DT5" } { 1 "DT" "DT3" } >> >> } >> >> addPdbAtomMap { >> { "O5*" "O5'" } >> { "C5*" "C5'" } >> { "C4*" "C4'" } >> { "O4*" "O4'" } >> { "C3*" "C3'" } >> { "O3*" "O3'" } >> { "C2*" "C2'" } >> { "O2*" "O2'" } >> { "C1*" "C1'" } >> { "C5M" "C7" } >> { "H1*" "H1'" } >> { "H2*1" "H2'1" } >> { "H2*2" "H2'2" } >> { "H2'" "H2'1" } >> { "H2''" "H2'2" } >> { "H3*" "H3'" } >> { "H4*" "H4'" } >> { "H5*1" "H5'1" } >> { "H5*2" "H5'2" } >> { "H5'" "H5'1" } >> { "H5''" "H5'2" } >> { "HO2'" "HO'2" } >> { "HO5'" "H5T" } >> { "HO3'" "H3T" } >> # old ff atom names -> new >> { "O1'" "O4'" } >> { "OA" "O1P" } >> { "OB" "O2P" } >> { "OP1" "O1P" } >> { "OP2" "O2P" } >> } >> >> >> # >> # assumed that most often proteins use HIE >> # >> NHIS = NHIE >> HIS = HIE >> CHIS = CHIE > > loadOff RNA_CI.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/RNA_CI.lib 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 > loadamberparams frcmod.tip3p Could not open file frcmod.tip3p: not found > loadamberparams frcmod.ionsjc_tip3p Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.ionsjc_tip3p Reading force field modification type file (frcmod) Reading title: Monovalent ion parameters for Ewald and TIP3P water from Joung & Cheatham JPCB (2008) (UNKNOWN ATOM TYPE: Li+) (UNKNOWN ATOM TYPE: Na+) (UNKNOWN ATOM TYPE: K+) (UNKNOWN ATOM TYPE: Rb+) (UNKNOWN ATOM TYPE: Cs+) (UNKNOWN ATOM TYPE: F-) (UNKNOWN ATOM TYPE: Cl-) (UNKNOWN ATOM TYPE: Br-) (UNKNOWN ATOM TYPE: I-) > loadamberparams frcmod.ff99Chen-Garcia Loading parameters: ./frcmod.ff99Chen-Garcia Reading force field modification type file (frcmod) Reading title: glycosdic bond torsion modifications to parm99 from Chen and Garcia (PNAS, 2013, 110, 16820) (UNKNOWN ATOM TYPE: CS) > > m = sequence {RA5 RA3} Sequence: RA5 Sequence: RA3 Joining RA5 - RA3 > > check m Checking 'm'.... WARNING: The unperturbed charge of the unit: -1.000000 is not zero. Warning: Close contact of 1.378936 angstroms between .R.A and .R.A Warning: Close contact of 0.666938 angstroms between .R.A and .R.A Checking parameters for unit 'm'. Checking for bond parameters. Checking for angle parameters. check: Warnings: 3 Unit is OK. > > solvatebox m TIP3PBOX 7.35 iso 2 Solute vdw bounding box: 11.977 20.213 9.591 Total bounding box for atom centers: 34.913 34.913 34.913 (box expansion for 'iso' is 88.1%) Solvent unit box: 18.774 18.774 18.774 The number of boxes: x= 2 y= 2 z= 2 Adding box at: x=0 y=0 z=0 Center of solvent box is: 9.387175, 9.387175, 9.387175 Adding box at: x=0 y=0 z=1 Center of solvent box is: 9.387175, 9.387175, -9.387175 Adding box at: x=0 y=1 z=0 Center of solvent box is: 9.387175, -9.387175, 9.387175 Adding box at: x=0 y=1 z=1 Center of solvent box is: 9.387175, -9.387175, -9.387175 Adding box at: x=1 y=0 z=0 Center of solvent box is: -9.387175, 9.387175, 9.387175 Adding box at: x=1 y=0 z=1 Center of solvent box is: -9.387175, 9.387175, -9.387175 Adding box at: x=1 y=1 z=0 Center of solvent box is: -9.387175, -9.387175, 9.387175 Adding box at: x=1 y=1 z=1 Center of solvent box is: -9.387175, -9.387175, -9.387175 Volume: 54302.015 A^3 Total mass 20034.726 amu, Density 0.613 g/cc Added 1079 residues. > > addions m Na+ 0 1 Na+ ion required to neutralize. Adding 1 counter ions to "m" using 1A grid Total solute charge: -1.00 Max atom radius: 2.10 Grid extends from solute vdw + 1.87 to 7.97 Box: enclosing: -12.94 -17.23 -11.75 12.74 16.75 11.12 sized: 51.06 46.77 52.25 edge: 64.00 Resolution: 1.00 Angstrom. Tree depth: 6 grid build: 0 sec Volume = 2.80% of box, grid points 7335 Solvent present: replacing closest with ion when steric overlaps occur Calculating grid charges charges: 0 sec (No solvent overlap) Placed Na+ in m at (-0.94, -1.23, 5.25). Done adding ions. > addions m Na+ 5 Cl- 5 Adding 10 counter ions to "m" using 1A grid Total solute charge: -0.00 Max atom radius: 2.10 Grid extends from solute vdw + 5.49 to 11.59 Box: enclosing: -16.56 -20.85 -15.37 16.36 20.37 16.85 sized: 47.44 43.15 48.63 edge: 64.00 Resolution: 1.00 Angstrom. Tree depth: 6 grid build: 0 sec Volume = 5.09% of box, grid points 13347 Solvent present: replacing closest with ion when steric overlaps occur Calculating grid charges charges: 0 sec (Replacing solvent molecule) Placed Na+ in m at (4.21, 10.71, -8.22). (Replacing solvent molecule) Placed Cl- in m at (1.86, 13.42, -11.23). (Replacing solvent molecule) Placed Na+ in m at (-2.00, 16.13, -9.38). (Replacing solvent molecule) Placed Cl- in m at (-2.47, 14.93, -6.89). (No solvent overlap) Placed Na+ in m at (-4.56, 18.15, -4.37). (Replacing solvent molecule) Placed Cl- in m at (-5.42, 15.36, 2.44). (Replacing solvent molecule) Placed Na+ in m at (-4.56, 14.81, 8.46). (Replacing solvent molecule) Placed Cl- in m at (-1.84, 7.30, 12.31). (Replacing solvent molecule) Placed Na+ in m at (-7.70, 5.66, 7.06). (Replacing solvent molecule) Placed Cl- in m at (-10.49, 1.06, 8.30). Done adding ions. > > saveamberparm m full-ff99.topo non-randomized.crds Checking Unit. Building topology. Building atom parameters. Building bond parameters. Building angle parameters. Building proper torsion parameters. Building improper torsion parameters. ** Warning: No sp2 improper torsion term for N*-CB-CB-NC atoms are: N9 C5 C4 N3 ** Warning: No sp2 improper torsion term for NB-CA-CB-CB atoms are: N7 C6 C5 C4 ** Warning: No sp2 improper torsion term for OH-H1-CI-H1 atoms are: O5' H5'1 C5' H5'2 ** Warning: No sp2 improper torsion term for OH-H1-CI-CT atoms are: O5' H5'1 C5' C4' ** Warning: No sp2 improper torsion term for OH-H1-CI-CT atoms are: O5' H5'2 C5' C4' ** Warning: No sp2 improper torsion term for H1-H1-CI-CT atoms are: H5'1 H5'2 C5' C4' ** Warning: No sp2 improper torsion term for N*-CB-CB-NC atoms are: N9 C5 C4 N3 ** Warning: No sp2 improper torsion term for NB-CA-CB-CB atoms are: N7 C6 C5 C4 ** Warning: No sp2 improper torsion term for OS-H1-CI-H1 atoms are: O5' H5'1 C5' H5'2 ** Warning: No sp2 improper torsion term for OS-H1-CI-CT atoms are: O5' H5'1 C5' C4' ** Warning: No sp2 improper torsion term for OS-H1-CI-CT atoms are: O5' H5'2 C5' C4' ** Warning: No sp2 improper torsion term for H1-H1-CI-CT atoms are: H5'1 H5'2 C5' C4' total 10 improper torsions applied Building H-Bond parameters. Incorporating Non-Bonded adjustments. Not Marking per-residue atom chain types. Marking per-residue atom chain types. (Residues lacking connect0/connect1 - these don't have chain types marked: res total affected WAT 1070 ) (no restraints) > check m Checking 'm'.... Warning: Close contact of 1.378936 angstroms between .R.A and .R.A Warning: Close contact of 0.666938 angstroms between .R.A and .R.A Checking parameters for unit 'm'. Checking for bond parameters. Checking for angle parameters. check: Warnings: 2 Unit is OK. > quit Quit log started: Tue Oct 27 12:49:30 2015 Log file: ./leap.log >> # >> # ----- leaprc for loading the ff99SB (Hornak & Simmerling) force field >> # using parmbsc0 for nucleic acid. >> # ----- 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" } >> { "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" "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" } >> { "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" } >> { "EP" "" "sp3" } >> # glycam >> { "OG" "O" "sp3" } >> { "OL" "O" "sp3" } >> { "AC" "C" "sp3" } >> { "EC" "C" "sp3" } >> } >> # >> # Load the main parameter set. >> # >> parm99 = loadamberparams parm99.dat Loading parameters: /uufs/chpc.utah.edu/common/home/u0827715/Amber/GIT/amber/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: IM) (UNKNOWN ATOM TYPE: Li) (UNKNOWN ATOM TYPE: IP) (UNKNOWN ATOM TYPE: K) (UNKNOWN ATOM TYPE: Rb) (UNKNOWN ATOM TYPE: Cs) (UNKNOWN ATOM TYPE: Zn) (UNKNOWN ATOM TYPE: IB) (UNKNOWN ATOM TYPE: LP) >> frcmod99SB = loadamberparams frcmod.ff99SB Loading parameters: /uufs/chpc.utah.edu/common/home/u0827715/Amber/GIT/amber/dat/leap/parm/frcmod.ff99SB Reading force field modification type file (frcmod) Reading title: Modification/update of parm99.dat (Hornak & Simmerling) >> frcmodbsc0 = loadamberparams frcmod.parmbsc0 Loading parameters: /uufs/chpc.utah.edu/common/home/u0827715/Amber/GIT/amber/dat/leap/parm/frcmod.parmbsc0 Reading force field modification type file (frcmod) Reading title: modifies parm99 for nucleic acids: see http://mmb.pcb.ub.es/PARMBSC0/ >> # >> # Load DNA/RNA libraries >> # >> loadOff all_nucleic94.lib Loading library: /uufs/chpc.utah.edu/common/home/u0827715/Amber/GIT/amber/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 (i.e. ff94 libs) >> # >> loadOff all_amino94.lib Loading library: /uufs/chpc.utah.edu/common/home/u0827715/Amber/GIT/amber/dat/leap/lib/all_amino94.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: ILE Loading: LEU Loading: LYN Loading: LYS Loading: MET Loading: PHE Loading: PRO Loading: SER Loading: THR Loading: TRP Loading: TYR Loading: VAL >> loadOff all_aminoct94.lib Loading library: /uufs/chpc.utah.edu/common/home/u0827715/Amber/GIT/amber/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: /uufs/chpc.utah.edu/common/home/u0827715/Amber/GIT/amber/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 DNA_CI.lib Loading library: /uufs/chpc.utah.edu/common/home/u0827715/Amber/GIT/amber/dat/leap/lib/DNA_CI.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 >> # >> # Load water and ions >> # >> loadOff ions94.lib Loading library: /uufs/chpc.utah.edu/common/home/u0827715/Amber/GIT/amber/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: /uufs/chpc.utah.edu/common/home/u0827715/Amber/GIT/amber/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 >> >> # >> # 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" "RG5" } { 1 "G" "RG3" } { "G" "RG" } { "GN" "RGN" } >> { 0 "A" "RA5" } { 1 "A" "RA3" } { "A" "RA" } { "AN" "RAN" } >> { 0 "C" "RC5" } { 1 "C" "RC3" } { "C" "RC" } { "CN" "RCN" } >> { 0 "U" "RU5" } { 1 "U" "RU3" } { "U" "RU" } { "UN" "RUN" } >> { 0 "DG" "DG5" } { 1 "DG" "DG3" } >> { 0 "DA" "DA5" } { 1 "DA" "DA3" } >> { 0 "DC" "DC5" } { 1 "DC" "DC3" } >> { 0 "DT" "DT5" } { 1 "DT" "DT3" } >> >> } >> >> addPdbAtomMap { >> { "O5*" "O5'" } >> { "C5*" "C5'" } >> { "C4*" "C4'" } >> { "O4*" "O4'" } >> { "C3*" "C3'" } >> { "O3*" "O3'" } >> { "C2*" "C2'" } >> { "O2*" "O2'" } >> { "C1*" "C1'" } >> { "C5M" "C7" } >> { "H1*" "H1'" } >> { "H2*1" "H2'1" } >> { "H2*2" "H2'2" } >> { "H2'" "H2'1" } >> { "H2''" "H2'2" } >> { "H3*" "H3'" } >> { "H4*" "H4'" } >> { "H5*1" "H5'1" } >> { "H5*2" "H5'2" } >> { "H5'" "H5'1" } >> { "H5''" "H5'2" } >> { "HO2'" "HO'2" } >> { "HO5'" "H5T" } >> { "HO3'" "H3T" } >> # old ff atom names -> new >> { "O1'" "O4'" } >> { "OA" "O1P" } >> { "OB" "O2P" } >> { "OP1" "O1P" } >> { "OP2" "O2P" } >> } >> >> >> # >> # assumed that most often proteins use HIE >> # >> NHIS = NHIE >> HIS = HIE >> CHIS = CHIE > > loadOff RNA_CI.lib Loading library: /uufs/chpc.utah.edu/common/home/u0827715/Amber/GIT/amber/dat/leap/lib/RNA_CI.lib 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 > loadamberparams frcmod.tip3p Could not open file frcmod.tip3p: not found > loadamberparams frcmod.ionsjc_tip3p Loading parameters: /uufs/chpc.utah.edu/common/home/u0827715/Amber/GIT/amber/dat/leap/parm/frcmod.ionsjc_tip3p Reading force field modification type file (frcmod) Reading title: Monovalent ion parameters for Ewald and TIP3P water from Joung & Cheatham JPCB (2008) (UNKNOWN ATOM TYPE: Li+) (UNKNOWN ATOM TYPE: Na+) (UNKNOWN ATOM TYPE: K+) (UNKNOWN ATOM TYPE: Rb+) (UNKNOWN ATOM TYPE: Cs+) (UNKNOWN ATOM TYPE: F-) (UNKNOWN ATOM TYPE: Cl-) (UNKNOWN ATOM TYPE: Br-) (UNKNOWN ATOM TYPE: I-) > loadamberparams frcmod.ff99Chen-Garcia Loading parameters: ./frcmod.ff99Chen-Garcia Reading force field modification type file (frcmod) Reading title: glycosdic bond torsion modifications to parm99 from Chen and Garcia (PNAS, 2013, 110, 16820) (UNKNOWN ATOM TYPE: CS) > > m = sequence {RA5 RC3} Sequence: RA5 Sequence: RC3 Joining RA5 - RC3 > > check m Checking 'm'.... WARNING: The unperturbed charge of the unit: -1.000000 is not zero. Warning: Close contact of 1.378936 angstroms between .R.A and .R.A Warning: Close contact of 0.666938 angstroms between .R.A and .R.A Checking parameters for unit 'm'. Checking for bond parameters. Checking for angle parameters. check: Warnings: 3 Unit is OK. > > solvatebox m TIP3PBOX 7.35 iso 2 Solute vdw bounding box: 11.847 18.599 9.384 Total bounding box for atom centers: 33.299 33.299 33.299 (box expansion for 'iso' is 73.4%) Solvent unit box: 18.774 18.774 18.774 The number of boxes: x= 2 y= 2 z= 2 Adding box at: x=0 y=0 z=0 Center of solvent box is: 9.387175, 9.387175, 9.387175 Adding box at: x=0 y=0 z=1 Center of solvent box is: 9.387175, 9.387175, -9.387175 Adding box at: x=0 y=1 z=0 Center of solvent box is: 9.387175, -9.387175, 9.387175 Adding box at: x=0 y=1 z=1 Center of solvent box is: 9.387175, -9.387175, -9.387175 Adding box at: x=1 y=0 z=0 Center of solvent box is: -9.387175, 9.387175, 9.387175 Adding box at: x=1 y=0 z=1 Center of solvent box is: -9.387175, 9.387175, -9.387175 Adding box at: x=1 y=1 z=0 Center of solvent box is: -9.387175, -9.387175, 9.387175 Adding box at: x=1 y=1 z=1 Center of solvent box is: -9.387175, -9.387175, -9.387175 Volume: 47483.242 A^3 Total mass 16569.640 amu, Density 0.579 g/cc Added 888 residues. > > addions m Na+ 0 1 Na+ ion required to neutralize. Adding 1 counter ions to "m" using 1A grid Total solute charge: -1.00 Max atom radius: 2.10 Grid extends from solute vdw + 1.87 to 7.97 Box: enclosing: -12.68 -16.60 -11.44 12.61 15.49 10.91 sized: 51.32 47.40 52.56 edge: 64.00 Resolution: 1.00 Angstrom. Tree depth: 6 grid build: 0 sec Volume = 2.66% of box, grid points 6977 Solvent present: replacing closest with ion when steric overlaps occur Calculating grid charges charges: 0 sec (No solvent overlap) Placed Na+ in m at (-2.68, -0.60, 3.56). Done adding ions. > addions m Na+ 5 Cl- 5 Na+ has a charge of 0. Adding the ions anyway. Adding 10 counter ions to "m" using 1A grid Total solute charge: 0.00 Max atom radius: 2.10 Grid extends from solute vdw + 5.49 to 11.59 Box: enclosing: -16.31 -20.23 -15.06 16.24 19.12 15.16 sized: 47.69 43.77 48.94 edge: 64.00 Resolution: 1.00 Angstrom. Tree depth: 6 grid build: 0 sec Volume = 4.86% of box, grid points 12749 Solvent present: replacing closest with ion when steric overlaps occur Calculating grid charges charges: 0 sec (Replacing solvent molecule) Placed Na+ in m at (6.29, 11.32, -6.73). (Replacing solvent molecule) Placed Cl- in m at (10.91, 13.93, -3.23). (Replacing solvent molecule) Placed Na+ in m at (12.21, 13.33, -0.68). (Replacing solvent molecule) Placed Cl- in m at (10.79, 8.94, 4.93). (Replacing solvent molecule) Placed Na+ in m at (11.19, 5.47, 6.95). (Replacing solvent molecule) Placed Cl- in m at (13.40, 9.14, -5.13). (Replacing solvent molecule) Placed Na+ in m at (13.77, 6.81, 2.35). (Replacing solvent molecule) Placed Cl- in m at (15.43, 3.18, -0.99). (Replacing solvent molecule) Placed Na+ in m at (12.43, 4.66, -5.68). (Replacing solvent molecule) Placed Cl- in m at (4.33, 10.52, 10.45). Done adding ions. > > saveamberparm m full-ff99.topo non-randomized.crds Checking Unit. Building topology. Building atom parameters. Building bond parameters. Building angle parameters. Building proper torsion parameters. Building improper torsion parameters. ** Warning: No sp2 improper torsion term for N*-CB-CB-NC atoms are: N9 C5 C4 N3 ** Warning: No sp2 improper torsion term for NB-CA-CB-CB atoms are: N7 C6 C5 C4 ** Warning: No sp2 improper torsion term for OH-H1-CI-H1 atoms are: O5' H5'1 C5' H5'2 ** Warning: No sp2 improper torsion term for OH-H1-CI-CT atoms are: O5' H5'1 C5' C4' ** Warning: No sp2 improper torsion term for OH-H1-CI-CT atoms are: O5' H5'2 C5' C4' ** Warning: No sp2 improper torsion term for H1-H1-CI-CT atoms are: H5'1 H5'2 C5' C4' ** Warning: No sp2 improper torsion term for OS-H1-CI-H1 atoms are: O5' H5'1 C5' H5'2 ** Warning: No sp2 improper torsion term for OS-H1-CI-CT atoms are: O5' H5'1 C5' C4' ** Warning: No sp2 improper torsion term for OS-H1-CI-CT atoms are: O5' H5'2 C5' C4' ** Warning: No sp2 improper torsion term for H1-H1-CI-CT atoms are: H5'1 H5'2 C5' C4' total 11 improper torsions applied Building H-Bond parameters. Incorporating Non-Bonded adjustments. Not Marking per-residue atom chain types. Marking per-residue atom chain types. (Residues lacking connect0/connect1 - these don't have chain types marked: res total affected WAT 878 ) (no restraints) > check m Checking 'm'.... Warning: Close contact of 1.378936 angstroms between .R.A and .R.A Warning: Close contact of 0.666938 angstroms between .R.A and .R.A Checking parameters for unit 'm'. Checking for bond parameters. Checking for angle parameters. check: Warnings: 2 Unit is OK. > quit Quit log started: Mon Nov 2 10:16:35 2015 Log file: ./leap.log >> # >> # ----- leaprc for loading the ff99SB (Hornak & Simmerling) force field >> # using parmbsc0 for nucleic acid. >> # ----- 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" } >> { "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" "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" } >> { "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" } >> { "EP" "" "sp3" } >> # glycam >> { "OG" "O" "sp3" } >> { "OL" "O" "sp3" } >> { "AC" "C" "sp3" } >> { "EC" "C" "sp3" } >> } >> # >> # Load the main parameter set. >> # >> parm99 = loadamberparams parm99.dat Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/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: IM) (UNKNOWN ATOM TYPE: Li) (UNKNOWN ATOM TYPE: IP) (UNKNOWN ATOM TYPE: K) (UNKNOWN ATOM TYPE: Rb) (UNKNOWN ATOM TYPE: Cs) (UNKNOWN ATOM TYPE: Zn) (UNKNOWN ATOM TYPE: IB) (UNKNOWN ATOM TYPE: LP) >> frcmod99SB = loadamberparams frcmod.ff99SB Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.ff99SB Reading force field modification type file (frcmod) Reading title: Modification/update of parm99.dat (Hornak & Simmerling) >> frcmodbsc0 = loadamberparams frcmod.parmbsc0 Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.parmbsc0 Reading force field modification type file (frcmod) Reading title: modifies parm99 for nucleic acids: see http://mmb.pcb.ub.es/PARMBSC0/ >> # >> # Load DNA/RNA libraries >> # >> loadOff all_nucleic94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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 (i.e. ff94 libs) >> # >> loadOff all_amino94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/all_amino94.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: ILE Loading: LEU Loading: LYN Loading: LYS Loading: MET Loading: PHE Loading: PRO Loading: SER Loading: THR Loading: TRP Loading: TYR Loading: VAL >> loadOff all_aminoct94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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: /oasis/projects/nsf/slc216/droe/GIT/amber/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 DNA_CI.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/DNA_CI.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 >> # >> # Load water and ions >> # >> loadOff ions94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/solvents.lib Loading: CHCL3BOX Loading: DC4 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 >> >> # >> # 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" "RG5" } { 1 "G" "RG3" } { "G" "RG" } { "GN" "RGN" } >> { 0 "A" "RA5" } { 1 "A" "RA3" } { "A" "RA" } { "AN" "RAN" } >> { 0 "C" "RC5" } { 1 "C" "RC3" } { "C" "RC" } { "CN" "RCN" } >> { 0 "U" "RU5" } { 1 "U" "RU3" } { "U" "RU" } { "UN" "RUN" } >> { 0 "DG" "DG5" } { 1 "DG" "DG3" } >> { 0 "DA" "DA5" } { 1 "DA" "DA3" } >> { 0 "DC" "DC5" } { 1 "DC" "DC3" } >> { 0 "DT" "DT5" } { 1 "DT" "DT3" } >> >> } >> >> addPdbAtomMap { >> { "O5*" "O5'" } >> { "C5*" "C5'" } >> { "C4*" "C4'" } >> { "O4*" "O4'" } >> { "C3*" "C3'" } >> { "O3*" "O3'" } >> { "C2*" "C2'" } >> { "O2*" "O2'" } >> { "C1*" "C1'" } >> { "C5M" "C7" } >> { "H1*" "H1'" } >> { "H2*1" "H2'1" } >> { "H2*2" "H2'2" } >> { "H2'" "H2'1" } >> { "H2''" "H2'2" } >> { "H3*" "H3'" } >> { "H4*" "H4'" } >> { "H5*1" "H5'1" } >> { "H5*2" "H5'2" } >> { "H5'" "H5'1" } >> { "H5''" "H5'2" } >> { "HO2'" "HO'2" } >> { "HO5'" "H5T" } >> { "HO3'" "H3T" } >> # old ff atom names -> new >> { "O1'" "O4'" } >> { "OA" "O1P" } >> { "OB" "O2P" } >> { "OP1" "O1P" } >> { "OP2" "O2P" } >> } >> >> >> # >> # assumed that most often proteins use HIE >> # >> NHIS = NHIE >> HIS = HIE >> CHIS = CHIE > > loadOff RNA_CI.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/RNA_CI.lib 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 > loadamberparams frcmod.tip3p Could not open file frcmod.tip3p: not found > loadamberparams frcmod.ionsjc_tip3p Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.ionsjc_tip3p Reading force field modification type file (frcmod) Reading title: Monovalent ion parameters for Ewald and TIP3P water from Joung & Cheatham JPCB (2008) (UNKNOWN ATOM TYPE: Li+) (UNKNOWN ATOM TYPE: Na+) (UNKNOWN ATOM TYPE: K+) (UNKNOWN ATOM TYPE: Rb+) (UNKNOWN ATOM TYPE: Cs+) (UNKNOWN ATOM TYPE: F-) (UNKNOWN ATOM TYPE: Cl-) (UNKNOWN ATOM TYPE: Br-) (UNKNOWN ATOM TYPE: I-) > loadamberparams frcmod.ff99Chen-Garcia Loading parameters: ./frcmod.ff99Chen-Garcia Reading force field modification type file (frcmod) Reading title: glycosdic bond torsion modifications to parm99 from Chen and Garcia (PNAS, 2013, 110, 16820) (UNKNOWN ATOM TYPE: CS) > > m = sequence {RA5 RC3} Sequence: RA5 Sequence: RC3 Joining RA5 - RC3 > > check m Checking 'm'.... WARNING: The unperturbed charge of the unit: -1.000000 is not zero. Warning: Close contact of 1.378936 angstroms between .R.A and .R.A Warning: Close contact of 0.666938 angstroms between .R.A and .R.A Checking parameters for unit 'm'. Checking for bond parameters. Checking for angle parameters. check: Warnings: 3 Unit is OK. > > savepdb m init.pdb Writing pdb file: init.pdb > quit Quit log started: Mon Nov 2 10:17:04 2015 Log file: ./leap.log >> # >> # ----- leaprc for loading the ff99SB (Hornak & Simmerling) force field >> # using parmbsc0 for nucleic acid. >> # ----- 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" } >> { "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" "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" } >> { "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" } >> { "EP" "" "sp3" } >> # glycam >> { "OG" "O" "sp3" } >> { "OL" "O" "sp3" } >> { "AC" "C" "sp3" } >> { "EC" "C" "sp3" } >> } >> # >> # Load the main parameter set. >> # >> parm99 = loadamberparams parm99.dat Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/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: IM) (UNKNOWN ATOM TYPE: Li) (UNKNOWN ATOM TYPE: IP) (UNKNOWN ATOM TYPE: K) (UNKNOWN ATOM TYPE: Rb) (UNKNOWN ATOM TYPE: Cs) (UNKNOWN ATOM TYPE: Zn) (UNKNOWN ATOM TYPE: IB) (UNKNOWN ATOM TYPE: LP) >> frcmod99SB = loadamberparams frcmod.ff99SB Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.ff99SB Reading force field modification type file (frcmod) Reading title: Modification/update of parm99.dat (Hornak & Simmerling) >> frcmodbsc0 = loadamberparams frcmod.parmbsc0 Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.parmbsc0 Reading force field modification type file (frcmod) Reading title: modifies parm99 for nucleic acids: see http://mmb.pcb.ub.es/PARMBSC0/ >> # >> # Load DNA/RNA libraries >> # >> loadOff all_nucleic94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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 (i.e. ff94 libs) >> # >> loadOff all_amino94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/all_amino94.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: ILE Loading: LEU Loading: LYN Loading: LYS Loading: MET Loading: PHE Loading: PRO Loading: SER Loading: THR Loading: TRP Loading: TYR Loading: VAL >> loadOff all_aminoct94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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: /oasis/projects/nsf/slc216/droe/GIT/amber/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 DNA_CI.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/DNA_CI.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 >> # >> # Load water and ions >> # >> loadOff ions94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/solvents.lib Loading: CHCL3BOX Loading: DC4 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 >> >> # >> # 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" "RG5" } { 1 "G" "RG3" } { "G" "RG" } { "GN" "RGN" } >> { 0 "A" "RA5" } { 1 "A" "RA3" } { "A" "RA" } { "AN" "RAN" } >> { 0 "C" "RC5" } { 1 "C" "RC3" } { "C" "RC" } { "CN" "RCN" } >> { 0 "U" "RU5" } { 1 "U" "RU3" } { "U" "RU" } { "UN" "RUN" } >> { 0 "DG" "DG5" } { 1 "DG" "DG3" } >> { 0 "DA" "DA5" } { 1 "DA" "DA3" } >> { 0 "DC" "DC5" } { 1 "DC" "DC3" } >> { 0 "DT" "DT5" } { 1 "DT" "DT3" } >> >> } >> >> addPdbAtomMap { >> { "O5*" "O5'" } >> { "C5*" "C5'" } >> { "C4*" "C4'" } >> { "O4*" "O4'" } >> { "C3*" "C3'" } >> { "O3*" "O3'" } >> { "C2*" "C2'" } >> { "O2*" "O2'" } >> { "C1*" "C1'" } >> { "C5M" "C7" } >> { "H1*" "H1'" } >> { "H2*1" "H2'1" } >> { "H2*2" "H2'2" } >> { "H2'" "H2'1" } >> { "H2''" "H2'2" } >> { "H3*" "H3'" } >> { "H4*" "H4'" } >> { "H5*1" "H5'1" } >> { "H5*2" "H5'2" } >> { "H5'" "H5'1" } >> { "H5''" "H5'2" } >> { "HO2'" "HO'2" } >> { "HO5'" "H5T" } >> { "HO3'" "H3T" } >> # old ff atom names -> new >> { "O1'" "O4'" } >> { "OA" "O1P" } >> { "OB" "O2P" } >> { "OP1" "O1P" } >> { "OP2" "O2P" } >> } >> >> >> # >> # assumed that most often proteins use HIE >> # >> NHIS = NHIE >> HIS = HIE >> CHIS = CHIE > > loadOff RNA_CI.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/RNA_CI.lib 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 > loadamberparams frcmod.ionsjc_tip3p Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.ionsjc_tip3p Reading force field modification type file (frcmod) Reading title: Monovalent ion parameters for Ewald and TIP3P water from Joung & Cheatham JPCB (2008) (UNKNOWN ATOM TYPE: Li+) (UNKNOWN ATOM TYPE: Na+) (UNKNOWN ATOM TYPE: K+) (UNKNOWN ATOM TYPE: Rb+) (UNKNOWN ATOM TYPE: Cs+) (UNKNOWN ATOM TYPE: F-) (UNKNOWN ATOM TYPE: Cl-) (UNKNOWN ATOM TYPE: Br-) (UNKNOWN ATOM TYPE: I-) > loadamberparams ../../../../../scripts/frcmod.ff99Chen-Garcia Loading parameters: ../../../../../scripts/frcmod.ff99Chen-Garcia Reading force field modification type file (frcmod) Reading title: glycosdic bond torsion modifications to parm99 from Chen and Garcia (PNAS, 2013, 110, 16820) (UNKNOWN ATOM TYPE: CS) > > m = loadpdb init.pdb Loading PDB file: ./init.pdb Enter zPdbReadScan from call depth 0. Exit zPdbReadScan from call depth 0. Matching PDB residue names to LEaP variables. (Residue 0: RA5, Terminal/beginning, was not found in name map.) (Residue 1: RC3, Terminal/last, was not found in name map.) Joining RA5 - RC3 total atoms in file: 63 > > > check m Checking 'm'.... WARNING: The unperturbed charge of the unit: -1.000000 is not zero. Warning: Close contact of 1.379097 angstroms between .R.A and .R.A Warning: Close contact of 0.666498 angstroms between .R.A and .R.A Checking parameters for unit 'm'. Checking for bond parameters. Checking for angle parameters. check: Warnings: 3 Unit is OK. > > solvatebox m TIP3PBOX 8 iso Solute vdw bounding box: 11.847 18.599 9.384 Total bounding box for atom centers: 34.599 34.599 34.599 (box expansion for 'iso' is 69.3%) Solvent unit box: 18.774 18.774 18.774 The number of boxes: x= 2 y= 2 z= 2 Adding box at: x=0 y=0 z=0 Center of solvent box is: 9.387175, 9.387175, 9.387175 Adding box at: x=0 y=0 z=1 Center of solvent box is: 9.387175, 9.387175, -9.387175 Adding box at: x=0 y=1 z=0 Center of solvent box is: 9.387175, -9.387175, 9.387175 Adding box at: x=0 y=1 z=1 Center of solvent box is: 9.387175, -9.387175, -9.387175 Adding box at: x=1 y=0 z=0 Center of solvent box is: -9.387175, 9.387175, 9.387175 Adding box at: x=1 y=0 z=1 Center of solvent box is: -9.387175, 9.387175, -9.387175 Adding box at: x=1 y=1 z=0 Center of solvent box is: -9.387175, -9.387175, 9.387175 Adding box at: x=1 y=1 z=1 Center of solvent box is: -9.387175, -9.387175, -9.387175 Volume: 53253.502 A^3 Total mass 21866.344 amu, Density 0.682 g/cc Added 1182 residues. > quit Quit log started: Mon Nov 2 10:17:04 2015 Log file: ./leap.log >> # >> # ----- leaprc for loading the ff99SB (Hornak & Simmerling) force field >> # using parmbsc0 for nucleic acid. >> # ----- 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" } >> { "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" "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" } >> { "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" } >> { "EP" "" "sp3" } >> # glycam >> { "OG" "O" "sp3" } >> { "OL" "O" "sp3" } >> { "AC" "C" "sp3" } >> { "EC" "C" "sp3" } >> } >> # >> # Load the main parameter set. >> # >> parm99 = loadamberparams parm99.dat Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/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: IM) (UNKNOWN ATOM TYPE: Li) (UNKNOWN ATOM TYPE: IP) (UNKNOWN ATOM TYPE: K) (UNKNOWN ATOM TYPE: Rb) (UNKNOWN ATOM TYPE: Cs) (UNKNOWN ATOM TYPE: Zn) (UNKNOWN ATOM TYPE: IB) (UNKNOWN ATOM TYPE: LP) >> frcmod99SB = loadamberparams frcmod.ff99SB Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.ff99SB Reading force field modification type file (frcmod) Reading title: Modification/update of parm99.dat (Hornak & Simmerling) >> frcmodbsc0 = loadamberparams frcmod.parmbsc0 Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.parmbsc0 Reading force field modification type file (frcmod) Reading title: modifies parm99 for nucleic acids: see http://mmb.pcb.ub.es/PARMBSC0/ >> # >> # Load DNA/RNA libraries >> # >> loadOff all_nucleic94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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 (i.e. ff94 libs) >> # >> loadOff all_amino94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/all_amino94.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: ILE Loading: LEU Loading: LYN Loading: LYS Loading: MET Loading: PHE Loading: PRO Loading: SER Loading: THR Loading: TRP Loading: TYR Loading: VAL >> loadOff all_aminoct94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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: /oasis/projects/nsf/slc216/droe/GIT/amber/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 DNA_CI.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/DNA_CI.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 >> # >> # Load water and ions >> # >> loadOff ions94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/solvents.lib Loading: CHCL3BOX Loading: DC4 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 >> >> # >> # 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" "RG5" } { 1 "G" "RG3" } { "G" "RG" } { "GN" "RGN" } >> { 0 "A" "RA5" } { 1 "A" "RA3" } { "A" "RA" } { "AN" "RAN" } >> { 0 "C" "RC5" } { 1 "C" "RC3" } { "C" "RC" } { "CN" "RCN" } >> { 0 "U" "RU5" } { 1 "U" "RU3" } { "U" "RU" } { "UN" "RUN" } >> { 0 "DG" "DG5" } { 1 "DG" "DG3" } >> { 0 "DA" "DA5" } { 1 "DA" "DA3" } >> { 0 "DC" "DC5" } { 1 "DC" "DC3" } >> { 0 "DT" "DT5" } { 1 "DT" "DT3" } >> >> } >> >> addPdbAtomMap { >> { "O5*" "O5'" } >> { "C5*" "C5'" } >> { "C4*" "C4'" } >> { "O4*" "O4'" } >> { "C3*" "C3'" } >> { "O3*" "O3'" } >> { "C2*" "C2'" } >> { "O2*" "O2'" } >> { "C1*" "C1'" } >> { "C5M" "C7" } >> { "H1*" "H1'" } >> { "H2*1" "H2'1" } >> { "H2*2" "H2'2" } >> { "H2'" "H2'1" } >> { "H2''" "H2'2" } >> { "H3*" "H3'" } >> { "H4*" "H4'" } >> { "H5*1" "H5'1" } >> { "H5*2" "H5'2" } >> { "H5'" "H5'1" } >> { "H5''" "H5'2" } >> { "HO2'" "HO'2" } >> { "HO5'" "H5T" } >> { "HO3'" "H3T" } >> # old ff atom names -> new >> { "O1'" "O4'" } >> { "OA" "O1P" } >> { "OB" "O2P" } >> { "OP1" "O1P" } >> { "OP2" "O2P" } >> } >> >> >> # >> # assumed that most often proteins use HIE >> # >> NHIS = NHIE >> HIS = HIE >> CHIS = CHIE > > loadOff RNA_CI.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/RNA_CI.lib 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 > loadamberparams frcmod.ionsjc_tip3p Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.ionsjc_tip3p Reading force field modification type file (frcmod) Reading title: Monovalent ion parameters for Ewald and TIP3P water from Joung & Cheatham JPCB (2008) (UNKNOWN ATOM TYPE: Li+) (UNKNOWN ATOM TYPE: Na+) (UNKNOWN ATOM TYPE: K+) (UNKNOWN ATOM TYPE: Rb+) (UNKNOWN ATOM TYPE: Cs+) (UNKNOWN ATOM TYPE: F-) (UNKNOWN ATOM TYPE: Cl-) (UNKNOWN ATOM TYPE: Br-) (UNKNOWN ATOM TYPE: I-) > loadamberparams ../../../../../scripts/frcmod.ff99Chen-Garcia Loading parameters: ../../../../../scripts/frcmod.ff99Chen-Garcia Reading force field modification type file (frcmod) Reading title: glycosdic bond torsion modifications to parm99 from Chen and Garcia (PNAS, 2013, 110, 16820) (UNKNOWN ATOM TYPE: CS) > > m = loadpdb init.pdb Loading PDB file: ./init.pdb Enter zPdbReadScan from call depth 0. Exit zPdbReadScan from call depth 0. Matching PDB residue names to LEaP variables. (Residue 0: RA5, Terminal/beginning, was not found in name map.) (Residue 1: RC3, Terminal/last, was not found in name map.) Joining RA5 - RC3 total atoms in file: 63 > > > check m Checking 'm'.... WARNING: The unperturbed charge of the unit: -1.000000 is not zero. Warning: Close contact of 1.379097 angstroms between .R.A and .R.A Warning: Close contact of 0.666498 angstroms between .R.A and .R.A Checking parameters for unit 'm'. Checking for bond parameters. Checking for angle parameters. check: Warnings: 3 Unit is OK. > > solvatebox m TIP3PBOX 7 iso Solute vdw bounding box: 11.847 18.599 9.384 Total bounding box for atom centers: 32.599 32.599 32.599 (box expansion for 'iso' is 75.8%) Solvent unit box: 18.774 18.774 18.774 The number of boxes: x= 2 y= 2 z= 2 Adding box at: x=0 y=0 z=0 Center of solvent box is: 9.387175, 9.387175, 9.387175 Adding box at: x=0 y=0 z=1 Center of solvent box is: 9.387175, 9.387175, -9.387175 Adding box at: x=0 y=1 z=0 Center of solvent box is: 9.387175, -9.387175, 9.387175 Adding box at: x=0 y=1 z=1 Center of solvent box is: 9.387175, -9.387175, -9.387175 Adding box at: x=1 y=0 z=0 Center of solvent box is: -9.387175, 9.387175, 9.387175 Adding box at: x=1 y=0 z=1 Center of solvent box is: -9.387175, 9.387175, -9.387175 Adding box at: x=1 y=1 z=0 Center of solvent box is: -9.387175, -9.387175, 9.387175 Adding box at: x=1 y=1 z=1 Center of solvent box is: -9.387175, -9.387175, -9.387175 Volume: 45545.493 A^3 Total mass 18028.936 amu, Density 0.657 g/cc Added 969 residues. > quit Quit log started: Mon Nov 2 10:17:04 2015 Log file: ./leap.log >> # >> # ----- leaprc for loading the ff99SB (Hornak & Simmerling) force field >> # using parmbsc0 for nucleic acid. >> # ----- 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" } >> { "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" "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" } >> { "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" } >> { "EP" "" "sp3" } >> # glycam >> { "OG" "O" "sp3" } >> { "OL" "O" "sp3" } >> { "AC" "C" "sp3" } >> { "EC" "C" "sp3" } >> } >> # >> # Load the main parameter set. >> # >> parm99 = loadamberparams parm99.dat Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/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: IM) (UNKNOWN ATOM TYPE: Li) (UNKNOWN ATOM TYPE: IP) (UNKNOWN ATOM TYPE: K) (UNKNOWN ATOM TYPE: Rb) (UNKNOWN ATOM TYPE: Cs) (UNKNOWN ATOM TYPE: Zn) (UNKNOWN ATOM TYPE: IB) (UNKNOWN ATOM TYPE: LP) >> frcmod99SB = loadamberparams frcmod.ff99SB Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.ff99SB Reading force field modification type file (frcmod) Reading title: Modification/update of parm99.dat (Hornak & Simmerling) >> frcmodbsc0 = loadamberparams frcmod.parmbsc0 Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.parmbsc0 Reading force field modification type file (frcmod) Reading title: modifies parm99 for nucleic acids: see http://mmb.pcb.ub.es/PARMBSC0/ >> # >> # Load DNA/RNA libraries >> # >> loadOff all_nucleic94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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 (i.e. ff94 libs) >> # >> loadOff all_amino94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/all_amino94.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: ILE Loading: LEU Loading: LYN Loading: LYS Loading: MET Loading: PHE Loading: PRO Loading: SER Loading: THR Loading: TRP Loading: TYR Loading: VAL >> loadOff all_aminoct94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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: /oasis/projects/nsf/slc216/droe/GIT/amber/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 DNA_CI.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/DNA_CI.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 >> # >> # Load water and ions >> # >> loadOff ions94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/solvents.lib Loading: CHCL3BOX Loading: DC4 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 >> >> # >> # 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" "RG5" } { 1 "G" "RG3" } { "G" "RG" } { "GN" "RGN" } >> { 0 "A" "RA5" } { 1 "A" "RA3" } { "A" "RA" } { "AN" "RAN" } >> { 0 "C" "RC5" } { 1 "C" "RC3" } { "C" "RC" } { "CN" "RCN" } >> { 0 "U" "RU5" } { 1 "U" "RU3" } { "U" "RU" } { "UN" "RUN" } >> { 0 "DG" "DG5" } { 1 "DG" "DG3" } >> { 0 "DA" "DA5" } { 1 "DA" "DA3" } >> { 0 "DC" "DC5" } { 1 "DC" "DC3" } >> { 0 "DT" "DT5" } { 1 "DT" "DT3" } >> >> } >> >> addPdbAtomMap { >> { "O5*" "O5'" } >> { "C5*" "C5'" } >> { "C4*" "C4'" } >> { "O4*" "O4'" } >> { "C3*" "C3'" } >> { "O3*" "O3'" } >> { "C2*" "C2'" } >> { "O2*" "O2'" } >> { "C1*" "C1'" } >> { "C5M" "C7" } >> { "H1*" "H1'" } >> { "H2*1" "H2'1" } >> { "H2*2" "H2'2" } >> { "H2'" "H2'1" } >> { "H2''" "H2'2" } >> { "H3*" "H3'" } >> { "H4*" "H4'" } >> { "H5*1" "H5'1" } >> { "H5*2" "H5'2" } >> { "H5'" "H5'1" } >> { "H5''" "H5'2" } >> { "HO2'" "HO'2" } >> { "HO5'" "H5T" } >> { "HO3'" "H3T" } >> # old ff atom names -> new >> { "O1'" "O4'" } >> { "OA" "O1P" } >> { "OB" "O2P" } >> { "OP1" "O1P" } >> { "OP2" "O2P" } >> } >> >> >> # >> # assumed that most often proteins use HIE >> # >> NHIS = NHIE >> HIS = HIE >> CHIS = CHIE > > loadOff RNA_CI.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/RNA_CI.lib 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 > loadamberparams frcmod.ionsjc_tip3p Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.ionsjc_tip3p Reading force field modification type file (frcmod) Reading title: Monovalent ion parameters for Ewald and TIP3P water from Joung & Cheatham JPCB (2008) (UNKNOWN ATOM TYPE: Li+) (UNKNOWN ATOM TYPE: Na+) (UNKNOWN ATOM TYPE: K+) (UNKNOWN ATOM TYPE: Rb+) (UNKNOWN ATOM TYPE: Cs+) (UNKNOWN ATOM TYPE: F-) (UNKNOWN ATOM TYPE: Cl-) (UNKNOWN ATOM TYPE: Br-) (UNKNOWN ATOM TYPE: I-) > loadamberparams ../../../../../scripts/frcmod.ff99Chen-Garcia Loading parameters: ../../../../../scripts/frcmod.ff99Chen-Garcia Reading force field modification type file (frcmod) Reading title: glycosdic bond torsion modifications to parm99 from Chen and Garcia (PNAS, 2013, 110, 16820) (UNKNOWN ATOM TYPE: CS) > > m = loadpdb init.pdb Loading PDB file: ./init.pdb Enter zPdbReadScan from call depth 0. Exit zPdbReadScan from call depth 0. Matching PDB residue names to LEaP variables. (Residue 0: RA5, Terminal/beginning, was not found in name map.) (Residue 1: RC3, Terminal/last, was not found in name map.) Joining RA5 - RC3 total atoms in file: 63 > > > check m Checking 'm'.... WARNING: The unperturbed charge of the unit: -1.000000 is not zero. Warning: Close contact of 1.379097 angstroms between .R.A and .R.A Warning: Close contact of 0.666498 angstroms between .R.A and .R.A Checking parameters for unit 'm'. Checking for bond parameters. Checking for angle parameters. check: Warnings: 3 Unit is OK. > > solvatebox m TIP3PBOX 7.5 iso Solute vdw bounding box: 11.847 18.599 9.384 Total bounding box for atom centers: 33.599 33.599 33.599 (box expansion for 'iso' is 72.4%) Solvent unit box: 18.774 18.774 18.774 The number of boxes: x= 2 y= 2 z= 2 Adding box at: x=0 y=0 z=0 Center of solvent box is: 9.387175, 9.387175, 9.387175 Adding box at: x=0 y=0 z=1 Center of solvent box is: 9.387175, 9.387175, -9.387175 Adding box at: x=0 y=1 z=0 Center of solvent box is: 9.387175, -9.387175, 9.387175 Adding box at: x=0 y=1 z=1 Center of solvent box is: 9.387175, -9.387175, -9.387175 Adding box at: x=1 y=0 z=0 Center of solvent box is: -9.387175, 9.387175, 9.387175 Adding box at: x=1 y=0 z=1 Center of solvent box is: -9.387175, 9.387175, -9.387175 Adding box at: x=1 y=1 z=0 Center of solvent box is: -9.387175, -9.387175, 9.387175 Adding box at: x=1 y=1 z=1 Center of solvent box is: -9.387175, -9.387175, -9.387175 Volume: 48794.245 A^3 Total mass 19596.328 amu, Density 0.667 g/cc Added 1056 residues. > quit Quit log started: Mon Nov 2 10:17:05 2015 Log file: ./leap.log >> # >> # ----- leaprc for loading the ff99SB (Hornak & Simmerling) force field >> # using parmbsc0 for nucleic acid. >> # ----- 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" } >> { "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" "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" } >> { "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" } >> { "EP" "" "sp3" } >> # glycam >> { "OG" "O" "sp3" } >> { "OL" "O" "sp3" } >> { "AC" "C" "sp3" } >> { "EC" "C" "sp3" } >> } >> # >> # Load the main parameter set. >> # >> parm99 = loadamberparams parm99.dat Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/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: IM) (UNKNOWN ATOM TYPE: Li) (UNKNOWN ATOM TYPE: IP) (UNKNOWN ATOM TYPE: K) (UNKNOWN ATOM TYPE: Rb) (UNKNOWN ATOM TYPE: Cs) (UNKNOWN ATOM TYPE: Zn) (UNKNOWN ATOM TYPE: IB) (UNKNOWN ATOM TYPE: LP) >> frcmod99SB = loadamberparams frcmod.ff99SB Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.ff99SB Reading force field modification type file (frcmod) Reading title: Modification/update of parm99.dat (Hornak & Simmerling) >> frcmodbsc0 = loadamberparams frcmod.parmbsc0 Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.parmbsc0 Reading force field modification type file (frcmod) Reading title: modifies parm99 for nucleic acids: see http://mmb.pcb.ub.es/PARMBSC0/ >> # >> # Load DNA/RNA libraries >> # >> loadOff all_nucleic94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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 (i.e. ff94 libs) >> # >> loadOff all_amino94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/all_amino94.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: ILE Loading: LEU Loading: LYN Loading: LYS Loading: MET Loading: PHE Loading: PRO Loading: SER Loading: THR Loading: TRP Loading: TYR Loading: VAL >> loadOff all_aminoct94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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: /oasis/projects/nsf/slc216/droe/GIT/amber/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 DNA_CI.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/DNA_CI.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 >> # >> # Load water and ions >> # >> loadOff ions94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/solvents.lib Loading: CHCL3BOX Loading: DC4 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 >> >> # >> # 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" "RG5" } { 1 "G" "RG3" } { "G" "RG" } { "GN" "RGN" } >> { 0 "A" "RA5" } { 1 "A" "RA3" } { "A" "RA" } { "AN" "RAN" } >> { 0 "C" "RC5" } { 1 "C" "RC3" } { "C" "RC" } { "CN" "RCN" } >> { 0 "U" "RU5" } { 1 "U" "RU3" } { "U" "RU" } { "UN" "RUN" } >> { 0 "DG" "DG5" } { 1 "DG" "DG3" } >> { 0 "DA" "DA5" } { 1 "DA" "DA3" } >> { 0 "DC" "DC5" } { 1 "DC" "DC3" } >> { 0 "DT" "DT5" } { 1 "DT" "DT3" } >> >> } >> >> addPdbAtomMap { >> { "O5*" "O5'" } >> { "C5*" "C5'" } >> { "C4*" "C4'" } >> { "O4*" "O4'" } >> { "C3*" "C3'" } >> { "O3*" "O3'" } >> { "C2*" "C2'" } >> { "O2*" "O2'" } >> { "C1*" "C1'" } >> { "C5M" "C7" } >> { "H1*" "H1'" } >> { "H2*1" "H2'1" } >> { "H2*2" "H2'2" } >> { "H2'" "H2'1" } >> { "H2''" "H2'2" } >> { "H3*" "H3'" } >> { "H4*" "H4'" } >> { "H5*1" "H5'1" } >> { "H5*2" "H5'2" } >> { "H5'" "H5'1" } >> { "H5''" "H5'2" } >> { "HO2'" "HO'2" } >> { "HO5'" "H5T" } >> { "HO3'" "H3T" } >> # old ff atom names -> new >> { "O1'" "O4'" } >> { "OA" "O1P" } >> { "OB" "O2P" } >> { "OP1" "O1P" } >> { "OP2" "O2P" } >> } >> >> >> # >> # assumed that most often proteins use HIE >> # >> NHIS = NHIE >> HIS = HIE >> CHIS = CHIE > > loadOff RNA_CI.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/RNA_CI.lib 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 > loadamberparams frcmod.ionsjc_tip3p Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.ionsjc_tip3p Reading force field modification type file (frcmod) Reading title: Monovalent ion parameters for Ewald and TIP3P water from Joung & Cheatham JPCB (2008) (UNKNOWN ATOM TYPE: Li+) (UNKNOWN ATOM TYPE: Na+) (UNKNOWN ATOM TYPE: K+) (UNKNOWN ATOM TYPE: Rb+) (UNKNOWN ATOM TYPE: Cs+) (UNKNOWN ATOM TYPE: F-) (UNKNOWN ATOM TYPE: Cl-) (UNKNOWN ATOM TYPE: Br-) (UNKNOWN ATOM TYPE: I-) > loadamberparams ../../../../../scripts/frcmod.ff99Chen-Garcia Loading parameters: ../../../../../scripts/frcmod.ff99Chen-Garcia Reading force field modification type file (frcmod) Reading title: glycosdic bond torsion modifications to parm99 from Chen and Garcia (PNAS, 2013, 110, 16820) (UNKNOWN ATOM TYPE: CS) > > m = loadpdb init.pdb Loading PDB file: ./init.pdb Enter zPdbReadScan from call depth 0. Exit zPdbReadScan from call depth 0. Matching PDB residue names to LEaP variables. (Residue 0: RA5, Terminal/beginning, was not found in name map.) (Residue 1: RC3, Terminal/last, was not found in name map.) Joining RA5 - RC3 total atoms in file: 63 > > > check m Checking 'm'.... WARNING: The unperturbed charge of the unit: -1.000000 is not zero. Warning: Close contact of 1.379097 angstroms between .R.A and .R.A Warning: Close contact of 0.666498 angstroms between .R.A and .R.A Checking parameters for unit 'm'. Checking for bond parameters. Checking for angle parameters. check: Warnings: 3 Unit is OK. > > solvatebox m TIP3PBOX 7.6 iso Solute vdw bounding box: 11.847 18.599 9.384 Total bounding box for atom centers: 33.799 33.799 33.799 (box expansion for 'iso' is 71.8%) Solvent unit box: 18.774 18.774 18.774 The number of boxes: x= 2 y= 2 z= 2 Adding box at: x=0 y=0 z=0 Center of solvent box is: 9.387175, 9.387175, 9.387175 Adding box at: x=0 y=0 z=1 Center of solvent box is: 9.387175, 9.387175, -9.387175 Adding box at: x=0 y=1 z=0 Center of solvent box is: 9.387175, -9.387175, 9.387175 Adding box at: x=0 y=1 z=1 Center of solvent box is: 9.387175, -9.387175, -9.387175 Adding box at: x=1 y=0 z=0 Center of solvent box is: -9.387175, 9.387175, 9.387175 Adding box at: x=1 y=0 z=1 Center of solvent box is: -9.387175, 9.387175, -9.387175 Adding box at: x=1 y=1 z=0 Center of solvent box is: -9.387175, -9.387175, 9.387175 Adding box at: x=1 y=1 z=1 Center of solvent box is: -9.387175, -9.387175, -9.387175 Volume: 49828.855 A^3 Total mass 20208.872 amu, Density 0.673 g/cc Added 1090 residues. > quit Quit log started: Mon Nov 2 10:17:05 2015 Log file: ./leap.log >> # >> # ----- leaprc for loading the ff99SB (Hornak & Simmerling) force field >> # using parmbsc0 for nucleic acid. >> # ----- 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" } >> { "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" "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" } >> { "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" } >> { "EP" "" "sp3" } >> # glycam >> { "OG" "O" "sp3" } >> { "OL" "O" "sp3" } >> { "AC" "C" "sp3" } >> { "EC" "C" "sp3" } >> } >> # >> # Load the main parameter set. >> # >> parm99 = loadamberparams parm99.dat Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/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: IM) (UNKNOWN ATOM TYPE: Li) (UNKNOWN ATOM TYPE: IP) (UNKNOWN ATOM TYPE: K) (UNKNOWN ATOM TYPE: Rb) (UNKNOWN ATOM TYPE: Cs) (UNKNOWN ATOM TYPE: Zn) (UNKNOWN ATOM TYPE: IB) (UNKNOWN ATOM TYPE: LP) >> frcmod99SB = loadamberparams frcmod.ff99SB Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.ff99SB Reading force field modification type file (frcmod) Reading title: Modification/update of parm99.dat (Hornak & Simmerling) >> frcmodbsc0 = loadamberparams frcmod.parmbsc0 Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.parmbsc0 Reading force field modification type file (frcmod) Reading title: modifies parm99 for nucleic acids: see http://mmb.pcb.ub.es/PARMBSC0/ >> # >> # Load DNA/RNA libraries >> # >> loadOff all_nucleic94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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 (i.e. ff94 libs) >> # >> loadOff all_amino94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/all_amino94.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: ILE Loading: LEU Loading: LYN Loading: LYS Loading: MET Loading: PHE Loading: PRO Loading: SER Loading: THR Loading: TRP Loading: TYR Loading: VAL >> loadOff all_aminoct94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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: /oasis/projects/nsf/slc216/droe/GIT/amber/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 DNA_CI.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/DNA_CI.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 >> # >> # Load water and ions >> # >> loadOff ions94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/solvents.lib Loading: CHCL3BOX Loading: DC4 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 >> >> # >> # 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" "RG5" } { 1 "G" "RG3" } { "G" "RG" } { "GN" "RGN" } >> { 0 "A" "RA5" } { 1 "A" "RA3" } { "A" "RA" } { "AN" "RAN" } >> { 0 "C" "RC5" } { 1 "C" "RC3" } { "C" "RC" } { "CN" "RCN" } >> { 0 "U" "RU5" } { 1 "U" "RU3" } { "U" "RU" } { "UN" "RUN" } >> { 0 "DG" "DG5" } { 1 "DG" "DG3" } >> { 0 "DA" "DA5" } { 1 "DA" "DA3" } >> { 0 "DC" "DC5" } { 1 "DC" "DC3" } >> { 0 "DT" "DT5" } { 1 "DT" "DT3" } >> >> } >> >> addPdbAtomMap { >> { "O5*" "O5'" } >> { "C5*" "C5'" } >> { "C4*" "C4'" } >> { "O4*" "O4'" } >> { "C3*" "C3'" } >> { "O3*" "O3'" } >> { "C2*" "C2'" } >> { "O2*" "O2'" } >> { "C1*" "C1'" } >> { "C5M" "C7" } >> { "H1*" "H1'" } >> { "H2*1" "H2'1" } >> { "H2*2" "H2'2" } >> { "H2'" "H2'1" } >> { "H2''" "H2'2" } >> { "H3*" "H3'" } >> { "H4*" "H4'" } >> { "H5*1" "H5'1" } >> { "H5*2" "H5'2" } >> { "H5'" "H5'1" } >> { "H5''" "H5'2" } >> { "HO2'" "HO'2" } >> { "HO5'" "H5T" } >> { "HO3'" "H3T" } >> # old ff atom names -> new >> { "O1'" "O4'" } >> { "OA" "O1P" } >> { "OB" "O2P" } >> { "OP1" "O1P" } >> { "OP2" "O2P" } >> } >> >> >> # >> # assumed that most often proteins use HIE >> # >> NHIS = NHIE >> HIS = HIE >> CHIS = CHIE > > loadOff RNA_CI.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/RNA_CI.lib 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 > loadamberparams frcmod.ionsjc_tip3p Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.ionsjc_tip3p Reading force field modification type file (frcmod) Reading title: Monovalent ion parameters for Ewald and TIP3P water from Joung & Cheatham JPCB (2008) (UNKNOWN ATOM TYPE: Li+) (UNKNOWN ATOM TYPE: Na+) (UNKNOWN ATOM TYPE: K+) (UNKNOWN ATOM TYPE: Rb+) (UNKNOWN ATOM TYPE: Cs+) (UNKNOWN ATOM TYPE: F-) (UNKNOWN ATOM TYPE: Cl-) (UNKNOWN ATOM TYPE: Br-) (UNKNOWN ATOM TYPE: I-) > loadamberparams ../../../../../scripts/frcmod.ff99Chen-Garcia Loading parameters: ../../../../../scripts/frcmod.ff99Chen-Garcia Reading force field modification type file (frcmod) Reading title: glycosdic bond torsion modifications to parm99 from Chen and Garcia (PNAS, 2013, 110, 16820) (UNKNOWN ATOM TYPE: CS) > > m = loadpdb init.pdb Loading PDB file: ./init.pdb Enter zPdbReadScan from call depth 0. Exit zPdbReadScan from call depth 0. Matching PDB residue names to LEaP variables. (Residue 0: RA5, Terminal/beginning, was not found in name map.) (Residue 1: RC3, Terminal/last, was not found in name map.) Joining RA5 - RC3 total atoms in file: 63 > > > check m Checking 'm'.... WARNING: The unperturbed charge of the unit: -1.000000 is not zero. Warning: Close contact of 1.379097 angstroms between .R.A and .R.A Warning: Close contact of 0.666498 angstroms between .R.A and .R.A Checking parameters for unit 'm'. Checking for bond parameters. Checking for angle parameters. check: Warnings: 3 Unit is OK. > > solvateoct m TIP3PBOX 7.6 iso Scaling up box by a factor of 1.260118 to meet diagonal cut criterion Solute vdw bounding box: 11.847 18.599 9.384 Total bounding box for atom centers: 37.753 37.753 37.753 (box expansion for 'iso' is 61.1%) 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 Volume: 29099.711 A^3 (oct) Total mass 14335.656 amu, Density 0.818 g/cc Added 764 residues. > remove m m.1092 remove: Argument #2 is type String must be of type: [molecule residue atom] usage: remove > remove m m.1091 remove: Argument #2 is type String must be of type: [molecule residue atom] usage: remove > remove m m.1090 remove: Argument #2 is type String must be of type: [molecule residue atom] usage: remove > remove m m.1089 remove: Argument #2 is type String must be of type: [molecule residue atom] usage: remove > > addions m Na+ 0 1 Na+ ion required to neutralize. Adding 1 counter ions to "m" using 1A grid Total solute charge: -1.00 Max atom radius: 2.10 Grid extends from solute vdw + 1.87 to 7.97 Box: enclosing: -14.64 -12.14 -13.13 13.53 13.29 15.99 sized: 17.36 19.86 18.87 edge: 32.00 Resolution: 1.00 Angstrom. Tree depth: 5 grid build: 0 sec Volume = 21.39% of box, grid points 7010 Solvent present: replacing closest with ion when steric overlaps occur Calculating grid charges charges: 0 sec (Replacing solvent molecule) Placed Na+ in m at (-3.57, -3.04, -2.74). Done adding ions. > addions m Na+ 5 Cl- 5 Na+ has a charge of 0. Adding the ions anyway. Adding 10 counter ions to "m" using 1A grid Total solute charge: 0.00 Max atom radius: 2.10 Grid extends from solute vdw + 5.49 to 11.59 Box: enclosing: -18.26 -15.76 -16.76 17.16 16.91 19.62 sized: 45.74 48.24 47.24 edge: 64.00 Resolution: 1.00 Angstrom. Tree depth: 6 grid build: 0 sec Volume = 4.89% of box, grid points 12810 Solvent present: replacing closest with ion when steric overlaps occur Calculating grid charges charges: 0 sec (Replacing solvent molecule) Placed Na+ in m at (-0.33, 8.03, 10.87). (Replacing solvent molecule) Placed Cl- in m at (1.84, 12.90, 10.31). (Replacing solvent molecule) Placed Na+ in m at (2.09, 14.07, 6.93). (Replacing solvent molecule) Placed Cl- in m at (1.41, 8.36, 5.35). (Replacing solvent molecule) Placed Na+ in m at (6.76, 9.86, 2.88). (Replacing solvent molecule) Placed Cl- in m at (9.31, 5.71, 3.29). (Replacing solvent molecule) Placed Na+ in m at (10.97, 5.96, -0.45). (Replacing solvent molecule) Placed Cl- in m at (9.47, 10.98, -0.93). (Replacing solvent molecule) Placed Na+ in m at (12.51, 1.68, 5.07). (Replacing solvent molecule) Placed Cl- in m at (12.77, 3.18, 0.42). Done adding ions. > > saveamberparm m full.topo non-randomized.crds Checking Unit. Building topology. Building atom parameters. Building bond parameters. Building angle parameters. Building proper torsion parameters. Building improper torsion parameters. ** Warning: No sp2 improper torsion term for N*-CB-CB-NC atoms are: N9 C5 C4 N3 ** Warning: No sp2 improper torsion term for NB-CA-CB-CB atoms are: N7 C6 C5 C4 ** Warning: No sp2 improper torsion term for OH-H1-CI-H1 atoms are: O5' H5'1 C5' H5'2 ** Warning: No sp2 improper torsion term for OH-H1-CI-CT atoms are: O5' H5'1 C5' C4' ** Warning: No sp2 improper torsion term for OH-H1-CI-CT atoms are: O5' H5'2 C5' C4' ** Warning: No sp2 improper torsion term for H1-H1-CI-CT atoms are: H5'1 H5'2 C5' C4' ** Warning: No sp2 improper torsion term for OS-H1-CI-H1 atoms are: O5' H5'1 C5' H5'2 ** Warning: No sp2 improper torsion term for OS-H1-CI-CT atoms are: O5' H5'1 C5' C4' ** Warning: No sp2 improper torsion term for OS-H1-CI-CT atoms are: O5' H5'2 C5' C4' ** Warning: No sp2 improper torsion term for H1-H1-CI-CT atoms are: H5'1 H5'2 C5' C4' total 11 improper torsions applied Building H-Bond parameters. Incorporating Non-Bonded adjustments. Not Marking per-residue atom chain types. Marking per-residue atom chain types. (Residues lacking connect0/connect1 - these don't have chain types marked: res total affected WAT 753 ) (no restraints) > quit Quit log started: Mon Nov 2 10:17:05 2015 Log file: ./leap.log >> # >> # ----- leaprc for loading the ff99SB (Hornak & Simmerling) force field >> # using parmbsc0 for nucleic acid. >> # ----- 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" } >> { "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" "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" } >> { "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" } >> { "EP" "" "sp3" } >> # glycam >> { "OG" "O" "sp3" } >> { "OL" "O" "sp3" } >> { "AC" "C" "sp3" } >> { "EC" "C" "sp3" } >> } >> # >> # Load the main parameter set. >> # >> parm99 = loadamberparams parm99.dat Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/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: IM) (UNKNOWN ATOM TYPE: Li) (UNKNOWN ATOM TYPE: IP) (UNKNOWN ATOM TYPE: K) (UNKNOWN ATOM TYPE: Rb) (UNKNOWN ATOM TYPE: Cs) (UNKNOWN ATOM TYPE: Zn) (UNKNOWN ATOM TYPE: IB) (UNKNOWN ATOM TYPE: LP) >> frcmod99SB = loadamberparams frcmod.ff99SB Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.ff99SB Reading force field modification type file (frcmod) Reading title: Modification/update of parm99.dat (Hornak & Simmerling) >> frcmodbsc0 = loadamberparams frcmod.parmbsc0 Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.parmbsc0 Reading force field modification type file (frcmod) Reading title: modifies parm99 for nucleic acids: see http://mmb.pcb.ub.es/PARMBSC0/ >> # >> # Load DNA/RNA libraries >> # >> loadOff all_nucleic94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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 (i.e. ff94 libs) >> # >> loadOff all_amino94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/all_amino94.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: ILE Loading: LEU Loading: LYN Loading: LYS Loading: MET Loading: PHE Loading: PRO Loading: SER Loading: THR Loading: TRP Loading: TYR Loading: VAL >> loadOff all_aminoct94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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: /oasis/projects/nsf/slc216/droe/GIT/amber/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 DNA_CI.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/DNA_CI.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 >> # >> # Load water and ions >> # >> loadOff ions94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/solvents.lib Loading: CHCL3BOX Loading: DC4 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 >> >> # >> # 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" "RG5" } { 1 "G" "RG3" } { "G" "RG" } { "GN" "RGN" } >> { 0 "A" "RA5" } { 1 "A" "RA3" } { "A" "RA" } { "AN" "RAN" } >> { 0 "C" "RC5" } { 1 "C" "RC3" } { "C" "RC" } { "CN" "RCN" } >> { 0 "U" "RU5" } { 1 "U" "RU3" } { "U" "RU" } { "UN" "RUN" } >> { 0 "DG" "DG5" } { 1 "DG" "DG3" } >> { 0 "DA" "DA5" } { 1 "DA" "DA3" } >> { 0 "DC" "DC5" } { 1 "DC" "DC3" } >> { 0 "DT" "DT5" } { 1 "DT" "DT3" } >> >> } >> >> addPdbAtomMap { >> { "O5*" "O5'" } >> { "C5*" "C5'" } >> { "C4*" "C4'" } >> { "O4*" "O4'" } >> { "C3*" "C3'" } >> { "O3*" "O3'" } >> { "C2*" "C2'" } >> { "O2*" "O2'" } >> { "C1*" "C1'" } >> { "C5M" "C7" } >> { "H1*" "H1'" } >> { "H2*1" "H2'1" } >> { "H2*2" "H2'2" } >> { "H2'" "H2'1" } >> { "H2''" "H2'2" } >> { "H3*" "H3'" } >> { "H4*" "H4'" } >> { "H5*1" "H5'1" } >> { "H5*2" "H5'2" } >> { "H5'" "H5'1" } >> { "H5''" "H5'2" } >> { "HO2'" "HO'2" } >> { "HO5'" "H5T" } >> { "HO3'" "H3T" } >> # old ff atom names -> new >> { "O1'" "O4'" } >> { "OA" "O1P" } >> { "OB" "O2P" } >> { "OP1" "O1P" } >> { "OP2" "O2P" } >> } >> >> >> # >> # assumed that most often proteins use HIE >> # >> NHIS = NHIE >> HIS = HIE >> CHIS = CHIE > > loadOff RNA_CI.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/RNA_CI.lib 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 > loadamberparams frcmod.ionsjc_tip3p Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.ionsjc_tip3p Reading force field modification type file (frcmod) Reading title: Monovalent ion parameters for Ewald and TIP3P water from Joung & Cheatham JPCB (2008) (UNKNOWN ATOM TYPE: Li+) (UNKNOWN ATOM TYPE: Na+) (UNKNOWN ATOM TYPE: K+) (UNKNOWN ATOM TYPE: Rb+) (UNKNOWN ATOM TYPE: Cs+) (UNKNOWN ATOM TYPE: F-) (UNKNOWN ATOM TYPE: Cl-) (UNKNOWN ATOM TYPE: Br-) (UNKNOWN ATOM TYPE: I-) > loadamberparams ../../../../../scripts/frcmod.ff99Chen-Garcia Loading parameters: ../../../../../scripts/frcmod.ff99Chen-Garcia Reading force field modification type file (frcmod) Reading title: glycosdic bond torsion modifications to parm99 from Chen and Garcia (PNAS, 2013, 110, 16820) (UNKNOWN ATOM TYPE: CS) > > m = loadpdb init.pdb Loading PDB file: ./init.pdb Enter zPdbReadScan from call depth 0. Exit zPdbReadScan from call depth 0. Matching PDB residue names to LEaP variables. (Residue 0: RA5, Terminal/beginning, was not found in name map.) (Residue 1: RC3, Terminal/last, was not found in name map.) Joining RA5 - RC3 total atoms in file: 63 > > > check m Checking 'm'.... WARNING: The unperturbed charge of the unit: -1.000000 is not zero. Warning: Close contact of 1.379097 angstroms between .R.A and .R.A Warning: Close contact of 0.666498 angstroms between .R.A and .R.A Checking parameters for unit 'm'. Checking for bond parameters. Checking for angle parameters. check: Warnings: 3 Unit is OK. > > solvatebox m TIP3PBOX 7.6 iso Solute vdw bounding box: 11.847 18.599 9.384 Total bounding box for atom centers: 33.799 33.799 33.799 (box expansion for 'iso' is 71.8%) Solvent unit box: 18.774 18.774 18.774 The number of boxes: x= 2 y= 2 z= 2 Adding box at: x=0 y=0 z=0 Center of solvent box is: 9.387175, 9.387175, 9.387175 Adding box at: x=0 y=0 z=1 Center of solvent box is: 9.387175, 9.387175, -9.387175 Adding box at: x=0 y=1 z=0 Center of solvent box is: 9.387175, -9.387175, 9.387175 Adding box at: x=0 y=1 z=1 Center of solvent box is: 9.387175, -9.387175, -9.387175 Adding box at: x=1 y=0 z=0 Center of solvent box is: -9.387175, 9.387175, 9.387175 Adding box at: x=1 y=0 z=1 Center of solvent box is: -9.387175, 9.387175, -9.387175 Adding box at: x=1 y=1 z=0 Center of solvent box is: -9.387175, -9.387175, 9.387175 Adding box at: x=1 y=1 z=1 Center of solvent box is: -9.387175, -9.387175, -9.387175 Volume: 49828.855 A^3 Total mass 20208.872 amu, Density 0.673 g/cc Added 1090 residues. > remove m m.1092 > remove m m.1091 > remove m m.1090 > remove m m.1089 > > addions m Na+ 0 1 Na+ ion required to neutralize. Adding 1 counter ions to "m" using 1A grid Total solute charge: -1.00 Max atom radius: 2.10 Grid extends from solute vdw + 1.87 to 7.97 Box: enclosing: -12.69 -16.22 -11.34 12.61 15.88 11.00 sized: 51.31 47.78 52.66 edge: 64.00 Resolution: 1.00 Angstrom. Tree depth: 6 grid build: 0 sec Volume = 2.66% of box, grid points 6976 Solvent present: replacing closest with ion when steric overlaps occur Calculating grid charges charges: 0 sec (Replacing solvent molecule) Placed Na+ in m at (-1.05, -0.51, 5.30). Done adding ions. > addions m Na+ 5 Cl- 5 Na+ has a charge of 0. Adding the ions anyway. Adding 10 counter ions to "m" using 1A grid Total solute charge: 0.00 Max atom radius: 2.10 Grid extends from solute vdw + 5.49 to 11.59 Box: enclosing: -16.31 -19.85 -14.97 16.24 19.50 16.89 sized: 47.69 44.15 49.03 edge: 64.00 Resolution: 1.00 Angstrom. Tree depth: 6 grid build: 0 sec Volume = 4.92% of box, grid points 12896 Solvent present: replacing closest with ion when steric overlaps occur Calculating grid charges charges: 0 sec (Replacing solvent molecule) Placed Na+ in m at (6.29, 11.70, -6.64). (Replacing solvent molecule) Placed Cl- in m at (10.91, 14.32, -3.13). (Replacing solvent molecule) Placed Na+ in m at (12.21, 13.71, -0.59). (Replacing solvent molecule) Placed Cl- in m at (10.79, 9.32, 5.03). (Replacing solvent molecule) Placed Na+ in m at (11.19, 5.85, 7.05). (Replacing solvent molecule) Placed Cl- in m at (4.72, 3.77, 11.97). (Replacing solvent molecule) Placed Na+ in m at (6.76, 7.42, 9.69). (Replacing solvent molecule) Placed Cl- in m at (13.40, 9.52, -5.04). (Replacing solvent molecule) Placed Na+ in m at (13.77, 7.19, 2.44). (Replacing solvent molecule) Placed Cl- in m at (15.43, 3.56, -0.89). Done adding ions. > > saveamberparm m full.topo non-randomized.crds Checking Unit. Building topology. Building atom parameters. Building bond parameters. Building angle parameters. Building proper torsion parameters. Building improper torsion parameters. ** Warning: No sp2 improper torsion term for N*-CB-CB-NC atoms are: N9 C5 C4 N3 ** Warning: No sp2 improper torsion term for NB-CA-CB-CB atoms are: N7 C6 C5 C4 ** Warning: No sp2 improper torsion term for OH-H1-CI-H1 atoms are: O5' H5'1 C5' H5'2 ** Warning: No sp2 improper torsion term for OH-H1-CI-CT atoms are: O5' H5'1 C5' C4' ** Warning: No sp2 improper torsion term for OH-H1-CI-CT atoms are: O5' H5'2 C5' C4' ** Warning: No sp2 improper torsion term for H1-H1-CI-CT atoms are: H5'1 H5'2 C5' C4' ** Warning: No sp2 improper torsion term for OS-H1-CI-H1 atoms are: O5' H5'1 C5' H5'2 ** Warning: No sp2 improper torsion term for OS-H1-CI-CT atoms are: O5' H5'1 C5' C4' ** Warning: No sp2 improper torsion term for OS-H1-CI-CT atoms are: O5' H5'2 C5' C4' ** Warning: No sp2 improper torsion term for H1-H1-CI-CT atoms are: H5'1 H5'2 C5' C4' total 11 improper torsions applied Building H-Bond parameters. Incorporating Non-Bonded adjustments. Not Marking per-residue atom chain types. Marking per-residue atom chain types. (Residues lacking connect0/connect1 - these don't have chain types marked: res total affected WAT 1075 ) (no restraints) > quit Quit log started: Mon Nov 2 10:17:06 2015 Log file: ./leap.log >> # >> # ----- leaprc for loading the ff99SB (Hornak & Simmerling) force field >> # using parmbsc0 for nucleic acid. >> # ----- 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" } >> { "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" "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" } >> { "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" } >> { "EP" "" "sp3" } >> # glycam >> { "OG" "O" "sp3" } >> { "OL" "O" "sp3" } >> { "AC" "C" "sp3" } >> { "EC" "C" "sp3" } >> } >> # >> # Load the main parameter set. >> # >> parm99 = loadamberparams parm99.dat Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/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: IM) (UNKNOWN ATOM TYPE: Li) (UNKNOWN ATOM TYPE: IP) (UNKNOWN ATOM TYPE: K) (UNKNOWN ATOM TYPE: Rb) (UNKNOWN ATOM TYPE: Cs) (UNKNOWN ATOM TYPE: Zn) (UNKNOWN ATOM TYPE: IB) (UNKNOWN ATOM TYPE: LP) >> frcmod99SB = loadamberparams frcmod.ff99SB Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.ff99SB Reading force field modification type file (frcmod) Reading title: Modification/update of parm99.dat (Hornak & Simmerling) >> frcmodbsc0 = loadamberparams frcmod.parmbsc0 Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.parmbsc0 Reading force field modification type file (frcmod) Reading title: modifies parm99 for nucleic acids: see http://mmb.pcb.ub.es/PARMBSC0/ >> # >> # Load DNA/RNA libraries >> # >> loadOff all_nucleic94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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 (i.e. ff94 libs) >> # >> loadOff all_amino94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/all_amino94.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: ILE Loading: LEU Loading: LYN Loading: LYS Loading: MET Loading: PHE Loading: PRO Loading: SER Loading: THR Loading: TRP Loading: TYR Loading: VAL >> loadOff all_aminoct94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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: /oasis/projects/nsf/slc216/droe/GIT/amber/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 DNA_CI.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/DNA_CI.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 >> # >> # Load water and ions >> # >> loadOff ions94.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/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: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/solvents.lib Loading: CHCL3BOX Loading: DC4 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 >> >> # >> # 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" "RG5" } { 1 "G" "RG3" } { "G" "RG" } { "GN" "RGN" } >> { 0 "A" "RA5" } { 1 "A" "RA3" } { "A" "RA" } { "AN" "RAN" } >> { 0 "C" "RC5" } { 1 "C" "RC3" } { "C" "RC" } { "CN" "RCN" } >> { 0 "U" "RU5" } { 1 "U" "RU3" } { "U" "RU" } { "UN" "RUN" } >> { 0 "DG" "DG5" } { 1 "DG" "DG3" } >> { 0 "DA" "DA5" } { 1 "DA" "DA3" } >> { 0 "DC" "DC5" } { 1 "DC" "DC3" } >> { 0 "DT" "DT5" } { 1 "DT" "DT3" } >> >> } >> >> addPdbAtomMap { >> { "O5*" "O5'" } >> { "C5*" "C5'" } >> { "C4*" "C4'" } >> { "O4*" "O4'" } >> { "C3*" "C3'" } >> { "O3*" "O3'" } >> { "C2*" "C2'" } >> { "O2*" "O2'" } >> { "C1*" "C1'" } >> { "C5M" "C7" } >> { "H1*" "H1'" } >> { "H2*1" "H2'1" } >> { "H2*2" "H2'2" } >> { "H2'" "H2'1" } >> { "H2''" "H2'2" } >> { "H3*" "H3'" } >> { "H4*" "H4'" } >> { "H5*1" "H5'1" } >> { "H5*2" "H5'2" } >> { "H5'" "H5'1" } >> { "H5''" "H5'2" } >> { "HO2'" "HO'2" } >> { "HO5'" "H5T" } >> { "HO3'" "H3T" } >> # old ff atom names -> new >> { "O1'" "O4'" } >> { "OA" "O1P" } >> { "OB" "O2P" } >> { "OP1" "O1P" } >> { "OP2" "O2P" } >> } >> >> >> # >> # assumed that most often proteins use HIE >> # >> NHIS = NHIE >> HIS = HIE >> CHIS = CHIE > > loadOff RNA_CI.lib Loading library: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/lib/RNA_CI.lib 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 > loadamberparams frcmod.ionsjc_tip3p Loading parameters: /oasis/projects/nsf/slc216/droe/GIT/amber/dat/leap/parm/frcmod.ionsjc_tip3p Reading force field modification type file (frcmod) Reading title: Monovalent ion parameters for Ewald and TIP3P water from Joung & Cheatham JPCB (2008) (UNKNOWN ATOM TYPE: Li+) (UNKNOWN ATOM TYPE: Na+) (UNKNOWN ATOM TYPE: K+) (UNKNOWN ATOM TYPE: Rb+) (UNKNOWN ATOM TYPE: Cs+) (UNKNOWN ATOM TYPE: F-) (UNKNOWN ATOM TYPE: Cl-) (UNKNOWN ATOM TYPE: Br-) (UNKNOWN ATOM TYPE: I-) > loadamberparams ../../../../../scripts/frcmod.ff99Chen-Garcia Loading parameters: ../../../../../scripts/frcmod.ff99Chen-Garcia Reading force field modification type file (frcmod) Reading title: glycosdic bond torsion modifications to parm99 from Chen and Garcia (PNAS, 2013, 110, 16820) (UNKNOWN ATOM TYPE: CS) > > m = loadpdb init.pdb Loading PDB file: ./init.pdb Enter zPdbReadScan from call depth 0. Exit zPdbReadScan from call depth 0. Matching PDB residue names to LEaP variables. (Residue 0: RA5, Terminal/beginning, was not found in name map.) (Residue 1: RC3, Terminal/last, was not found in name map.) Joining RA5 - RC3 total atoms in file: 63 > > > check m Checking 'm'.... WARNING: The unperturbed charge of the unit: -1.000000 is not zero. Warning: Close contact of 1.379097 angstroms between .R.A and .R.A Warning: Close contact of 0.666498 angstroms between .R.A and .R.A Checking parameters for unit 'm'. Checking for bond parameters. Checking for angle parameters. check: Warnings: 3 Unit is OK. > > solvatebox m TIP3PBOX 7.6 iso Solute vdw bounding box: 11.847 18.599 9.384 Total bounding box for atom centers: 33.799 33.799 33.799 (box expansion for 'iso' is 71.8%) Solvent unit box: 18.774 18.774 18.774 The number of boxes: x= 2 y= 2 z= 2 Adding box at: x=0 y=0 z=0 Center of solvent box is: 9.387175, 9.387175, 9.387175 Adding box at: x=0 y=0 z=1 Center of solvent box is: 9.387175, 9.387175, -9.387175 Adding box at: x=0 y=1 z=0 Center of solvent box is: 9.387175, -9.387175, 9.387175 Adding box at: x=0 y=1 z=1 Center of solvent box is: 9.387175, -9.387175, -9.387175 Adding box at: x=1 y=0 z=0 Center of solvent box is: -9.387175, 9.387175, 9.387175 Adding box at: x=1 y=0 z=1 Center of solvent box is: -9.387175, 9.387175, -9.387175 Adding box at: x=1 y=1 z=0 Center of solvent box is: -9.387175, -9.387175, 9.387175 Adding box at: x=1 y=1 z=1 Center of solvent box is: -9.387175, -9.387175, -9.387175 Volume: 49828.855 A^3 Total mass 20208.872 amu, Density 0.673 g/cc Added 1090 residues. > remove m m.1092 > remove m m.1091 > remove m m.1090 > remove m m.1089 > > addions m Na+ 0 1 Na+ ion required to neutralize. Adding 1 counter ions to "m" using 1A grid Total solute charge: -1.00 Max atom radius: 2.10 Grid extends from solute vdw + 1.87 to 7.97 Box: enclosing: -12.69 -16.22 -11.34 12.61 15.88 11.00 sized: 51.31 47.78 52.66 edge: 64.00 Resolution: 1.00 Angstrom. Tree depth: 6 grid build: 0 sec Volume = 2.66% of box, grid points 6976 Solvent present: replacing closest with ion when steric overlaps occur Calculating grid charges charges: 0 sec (Replacing solvent molecule) Placed Na+ in m at (-1.05, -0.51, 5.30). Done adding ions. > addions m Na+ 5 Cl- 5 Na+ has a charge of 0. Adding the ions anyway. Adding 10 counter ions to "m" using 1A grid Total solute charge: 0.00 Max atom radius: 2.10 Grid extends from solute vdw + 5.49 to 11.59 Box: enclosing: -16.31 -19.85 -14.97 16.24 19.50 16.89 sized: 47.69 44.15 49.03 edge: 64.00 Resolution: 1.00 Angstrom. Tree depth: 6 grid build: 0 sec Volume = 4.92% of box, grid points 12896 Solvent present: replacing closest with ion when steric overlaps occur Calculating grid charges charges: 0 sec (Replacing solvent molecule) Placed Na+ in m at (6.29, 11.70, -6.64). (Replacing solvent molecule) Placed Cl- in m at (10.91, 14.32, -3.13). (Replacing solvent molecule) Placed Na+ in m at (12.21, 13.71, -0.59). (Replacing solvent molecule) Placed Cl- in m at (10.79, 9.32, 5.03). (Replacing solvent molecule) Placed Na+ in m at (11.19, 5.85, 7.05). (Replacing solvent molecule) Placed Cl- in m at (4.72, 3.77, 11.97). (Replacing solvent molecule) Placed Na+ in m at (6.76, 7.42, 9.69). (Replacing solvent molecule) Placed Cl- in m at (13.40, 9.52, -5.04). (Replacing solvent molecule) Placed Na+ in m at (13.77, 7.19, 2.44). (Replacing solvent molecule) Placed Cl- in m at (15.43, 3.56, -0.89). Done adding ions. > > saveamberparm m full.topo non-randomized.crds Checking Unit. Building topology. Building atom parameters. Building bond parameters. Building angle parameters. Building proper torsion parameters. Building improper torsion parameters. ** Warning: No sp2 improper torsion term for N*-CB-CB-NC atoms are: N9 C5 C4 N3 ** Warning: No sp2 improper torsion term for NB-CA-CB-CB atoms are: N7 C6 C5 C4 ** Warning: No sp2 improper torsion term for OH-H1-CI-H1 atoms are: O5' H5'1 C5' H5'2 ** Warning: No sp2 improper torsion term for OH-H1-CI-CT atoms are: O5' H5'1 C5' C4' ** Warning: No sp2 improper torsion term for OH-H1-CI-CT atoms are: O5' H5'2 C5' C4' ** Warning: No sp2 improper torsion term for H1-H1-CI-CT atoms are: H5'1 H5'2 C5' C4' ** Warning: No sp2 improper torsion term for OS-H1-CI-H1 atoms are: O5' H5'1 C5' H5'2 ** Warning: No sp2 improper torsion term for OS-H1-CI-CT atoms are: O5' H5'1 C5' C4' ** Warning: No sp2 improper torsion term for OS-H1-CI-CT atoms are: O5' H5'2 C5' C4' ** Warning: No sp2 improper torsion term for H1-H1-CI-CT atoms are: H5'1 H5'2 C5' C4' total 11 improper torsions applied Building H-Bond parameters. Incorporating Non-Bonded adjustments. Not Marking per-residue atom chain types. Marking per-residue atom chain types. (Residues lacking connect0/connect1 - these don't have chain types marked: res total affected WAT 1075 ) (no restraints) > quit Quit