Re: [AMBER] How to solve FATAL: Atom .R<*>.A<*> does not have a type.

From: Rodrigo Galindo-Murillo via AMBER <amber.ambermd.org>
Date: Tue, 18 Mar 2025 19:57:27 -0600

Hello, you need to use (in your tleap file)

loadOff terminal_monophosphate.lib

This will load residues in the 5’ end with one or more phosphate groups.

Hope that helps!
Rodrigo.

On Tue, Mar 18, 2025 at 7:49 PM 1464903098 via AMBER <amber.ambermd.org> wrote:
>
> Hello everyone,
> I have a problem with the topology file and coordinate file generated by tleap at the beginning of using AMBER. I imported a single DNA strand, but I don't know why the four phosphate atoms at the 5' end of the first nucleic acid are not recognized. But these four atoms are very important for my subsequent calculations. I would like to ask you how to deal with this?
>
>
> Thank you and wish you all good luck.
>
>
>
>
>
>
>
>
>
>
> Below is the content of my leap.log file
> log started: Tue Mar 18 10:22:28 2025
>
>
> Log file: ./leap.log
> &gt;&gt; #
> &gt;&gt; # ----- NOTE: this is designed for PDB format 3!
> &gt;&gt; # &nbsp; &nbsp; &nbsp; &nbsp;Uses ff99bsc0+OL15 for DNA
> &gt;&gt; # &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;(ff99): Wang et al. J. Comp. Chem. 21: 1049, 2000
> &gt;&gt; # &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;(bsc0): Perez et al. Biophys J. 92: 3817, 2007
> &gt;&gt; # &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;(OL15): Zgarbova et al. JCTC 11: 5723, 2015
> &gt;&gt; #  load atom type hybridizations
> &gt;&gt; #
> &gt;&gt; addAtomTypes {
> &gt;&gt; { "H" &nbsp; "H" "sp3" }
> &gt;&gt; { "HO" &nbsp;"H" "sp3" }
> &gt;&gt; { "HS" &nbsp;"H" "sp3" }
> &gt;&gt; { "H1" &nbsp;"H" "sp3" }
> &gt;&gt; { "H2" &nbsp;"H" "sp3" }
> &gt;&gt; { "H3" &nbsp;"H" "sp3" }
> &gt;&gt; { "H4" &nbsp;"H" "sp3" }
> &gt;&gt; { "H5" &nbsp;"H" "sp3" }
> &gt;&gt; { "HW" &nbsp;"H" "sp3" }
> &gt;&gt; { "HC" &nbsp;"H" "sp3" }
> &gt;&gt; { "HA" &nbsp;"H" "sp3" }
> &gt;&gt; { "HP" &nbsp;"H" "sp3" }
> &gt;&gt; { "HZ" &nbsp;"H" "sp3" }
> &gt;&gt; { "OH" &nbsp;"O" "sp3" }
> &gt;&gt; { "OS" &nbsp;"O" "sp3" }
> &gt;&gt; { "O" &nbsp; "O" "sp2" }
> &gt;&gt; { "O2" &nbsp;"O" "sp2" }
> &gt;&gt; { "OP" &nbsp;"O" "sp2" }
> &gt;&gt; { "OW" &nbsp;"O" "sp3" }
> &gt;&gt; { "CT" &nbsp;"C" "sp3" }
> &gt;&gt; { "CX" &nbsp;"C" "sp3" }
> &gt;&gt; { "C8" &nbsp;"C" "sp3" }
> &gt;&gt; { "2C" &nbsp;"C" "sp3" }
> &gt;&gt; { "3C" &nbsp;"C" "sp3" }
> &gt;&gt; { "CH" &nbsp;"C" "sp3" }
> &gt;&gt; { "CS" &nbsp;"C" "sp2" }
> &gt;&gt; { "C" &nbsp; "C" "sp2" }
> &gt;&gt; { "CO" &nbsp;"C" "sp2" }
> &gt;&gt; { "C*" &nbsp;"C" "sp2" }
> &gt;&gt; { "CA" &nbsp;"C" "sp2" }
> &gt;&gt; { "CB" &nbsp;"C" "sp2" }
> &gt;&gt; { "CC" &nbsp;"C" "sp2" }
> &gt;&gt; { "CN" &nbsp;"C" "sp2" }
> &gt;&gt; { "CM" &nbsp;"C" "sp2" }
> &gt;&gt; { "CK" &nbsp;"C" "sp2" }
> &gt;&gt; { "CQ" &nbsp;"C" "sp2" }
> &gt;&gt; { "CD" &nbsp;"C" "sp2" }
> &gt;&gt; { "C5" &nbsp;"C" "sp2" }
> &gt;&gt; { "C4" &nbsp;"C" "sp2" }
> &gt;&gt; { "CP" &nbsp;"C" "sp2" }
> &gt;&gt; { "CI" &nbsp;"C" "sp3" }
> &gt;&gt; { "CJ" &nbsp;"C" "sp2" }
> &gt;&gt; { "CW" &nbsp;"C" "sp2" }
> &gt;&gt; { "CV" &nbsp;"C" "sp2" }
> &gt;&gt; { "CR" &nbsp;"C" "sp2" }
> &gt;&gt; { "CA" &nbsp;"C" "sp2" }
> &gt;&gt; { "CY" &nbsp;"C" "sp2" }
> &gt;&gt; { "C0" &nbsp;"Ca" "sp3" }
> &gt;&gt; { "MG" &nbsp;"Mg" "sp3" }
> &gt;&gt; { "N" &nbsp; "N" "sp2" }
> &gt;&gt; { "NA" &nbsp;"N" "sp2" }
> &gt;&gt; { "N2" &nbsp;"N" "sp2" }
> &gt;&gt; { "N*" &nbsp;"N" "sp2" }
> &gt;&gt; { "NP" &nbsp;"N" "sp2" }
> &gt;&gt; { "NQ" &nbsp;"N" "sp2" }
> &gt;&gt; { "NB" &nbsp;"N" "sp2" }
> &gt;&gt; { "NC" &nbsp;"N" "sp2" }
> &gt;&gt; { "NT" &nbsp;"N" "sp3" }
> &gt;&gt; { "NY" &nbsp;"N" "sp2" }
> &gt;&gt; { "N3" &nbsp;"N" "sp3" }
> &gt;&gt; { "S" &nbsp; "S" "sp3" }
> &gt;&gt; { "SH" &nbsp;"S" "sp3" }
> &gt;&gt; { "P" &nbsp; "P" "sp3" }
> &gt;&gt; { "LP" &nbsp;"" &nbsp;"sp3" }
> &gt;&gt; { "EP" &nbsp;"" &nbsp;"sp3" }
> &gt;&gt; { "F" &nbsp; "F" "sp3" }
> &gt;&gt; { "Cl" &nbsp;"Cl" "sp3" }
> &gt;&gt; { "Br" &nbsp;"Br" "sp3" }
> &gt;&gt; { "I" &nbsp; "I" &nbsp;"sp3" }
> &gt;&gt; { "F-" &nbsp; "F" "sp3" }
> &gt;&gt; { "Cl-" &nbsp;"Cl" "sp3" }
> &gt;&gt; { "Br-" &nbsp;"Br" "sp3" }
> &gt;&gt; { "I-" &nbsp; "I" &nbsp;"sp3" }
> &gt;&gt; { "Li+" &nbsp;"Li" &nbsp;"sp3" }
> &gt;&gt; { "Na+" &nbsp;"Na" &nbsp;"sp3" }
> &gt;&gt; { "K+" &nbsp;"K" &nbsp;"sp3" }
> &gt;&gt; { "Rb+" &nbsp;"Rb" &nbsp;"sp3" }
> &gt;&gt; { "Cs+" &nbsp;"Cs" &nbsp;"sp3" }
> &gt;&gt; { "Mg+" &nbsp;"Mg" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Tl+" &nbsp;"Tl" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Cu+" &nbsp;"Cu" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Ag+" &nbsp;"Ag" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Be2+" "Be" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Cu2+" "Cu" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Ni2+" "Ni" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Pt2+" "Pt" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Zn2+" "Zn" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Co2+" "Co" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Pd2+" "Pd" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Ag2+" "Ag" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Cr2+" "Cr" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Fe2+" "Fe" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Mg2+" "Mg" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "V2+" &nbsp;"V" &nbsp; "sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Mn2+" "Mn" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Hg2+" "Hg" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Cd2+" "Cd" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Yb2+" "Yb" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Ca2+" "Ca" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Sn2+" "Sn" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Pb2+" "Pb" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Eu2+" "Eu" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Sr2+" "Sr" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Sm2+" "Sm" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Ba2+" "Ba" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Ra2+" "Ra" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Al3+" "Al" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Fe3+" "Fe" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Cr3+" "Cr" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "In3+" "In" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Tl3+" "Tl" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Y3+" &nbsp;"Y" &nbsp; "sp3" }
> &gt;&gt; &nbsp; &nbsp; { "La3+" "La" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Ce3+" "Ce" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Pr3+" "Pr" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Nd3+" "Nd" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Sm3+" "Sm" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Eu3+" "Eu" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Gd3+" "Gd" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Tb3+" "Tb" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Dy3+" "Dy" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Er3+" "Er" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Tm3+" "Tm" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Lu3+" "Lu" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Hf4+" "Hf" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Zr4+" "Zr" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Ce4+" "Ce" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "U4+" &nbsp;"U" &nbsp; "sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Pu4+" "Pu" &nbsp;"sp3" }
> &gt;&gt; &nbsp; &nbsp; { "Th4+" "Th" &nbsp;"sp3" }
> &gt;&gt; }
> &gt;&gt; #
> &gt;&gt; #  Define the PDB name map for the nucleic acids
> &gt;&gt; #
> &gt;&gt; addPdbResMap {
> &gt;&gt; &nbsp; { 0 "DG" "DG5" &nbsp;} { 1 "DG" "DG3" &nbsp;} &nbsp;
> &gt;&gt; &nbsp; { 0 "DA" "DA5" &nbsp;} { 1 "DA" "DA3" &nbsp;} &nbsp;
> &gt;&gt; &nbsp; { 0 "DC" "DC5" &nbsp;} { 1 "DC" "DC3" &nbsp;} &nbsp;
> &gt;&gt; &nbsp; { 0 "DT" "DT5" &nbsp;} { 1 "DT" "DT3" &nbsp;} &nbsp;
> &gt;&gt; # &nbsp;some really old Amber residue names, assuming DNA:
> &gt;&gt; &nbsp; { 0 "GUA" "DG5" &nbsp;} { 1 "GUA" "DG3" &nbsp;} { "GUA" "DG" }
> &gt;&gt; &nbsp; { 0 "ADE" "DA5" &nbsp;} { 1 "ADE" "DA3" &nbsp;} { "ADE" "DA" }
> &gt;&gt; &nbsp; { 0 "CYT" "DC5" &nbsp;} { 1 "CYT" "DC3" &nbsp;} { "CYT" "DC" }
> &gt;&gt; &nbsp; { 0 "THY" "DT5" &nbsp;} { 1 "THY" "DT3" &nbsp;} { "THY" "DT" }
> &gt;&gt;&nbsp;
> &gt;&gt; }
> &gt;&gt;&nbsp;
> &gt;&gt; # &nbsp;try to be good about reading in really old atom names as well:
> &gt;&gt; addPdbAtomMap {
> &gt;&gt; &nbsp; { "O5*" "O5'" }
> &gt;&gt; &nbsp; { "C5*" "C5'" }
> &gt;&gt; &nbsp; { "C4*" "C4'" }
> &gt;&gt; &nbsp; { "O4*" "O4'" }
> &gt;&gt; &nbsp; { "C3*" "C3'" }
> &gt;&gt; &nbsp; { "O3*" "O3'" }
> &gt;&gt; &nbsp; { "C2*" "C2'" }
> &gt;&gt; &nbsp; { "O2*" "O2'" }
> &gt;&gt; &nbsp; { "C1*" "C1'" }
> &gt;&gt; &nbsp; { "C5M" "C7" &nbsp;}
> &gt;&gt; &nbsp; { "H1*" "H1'" }
> &gt;&gt; &nbsp; { "H2*1" "H2'" }
> &gt;&gt; &nbsp; { "H2*2" "H2''" }
> &gt;&gt; &nbsp; { "H2'1" "H2'" }
> &gt;&gt; &nbsp; { "H2'2" "H2''" }
> &gt;&gt; &nbsp; { "'H2'" "H2''" }
> &gt;&gt; &nbsp; { "H3*" "H3'" }
> &gt;&gt; &nbsp; { "H4*" "H4'" }
> &gt;&gt; &nbsp; { "H5*1" "H5'" }
> &gt;&gt; &nbsp; { "H5*2" "H5''" }
> &gt;&gt; &nbsp; { "H5'1" "H5'" }
> &gt;&gt; &nbsp; { "H5'2" "H5''" }
> &gt;&gt; &nbsp; { "'H5'" "H5''" }
> &gt;&gt; &nbsp; { "HO'2" "HO2'" }
> &gt;&gt; &nbsp; { "H5T" &nbsp;"HO5'" }
> &gt;&gt; &nbsp; { "H3T" &nbsp;"HO3'" }
> &gt;&gt; &nbsp; { "O1'" "O4'" }
> &gt;&gt; &nbsp; { "OA" &nbsp;"OP1" }
> &gt;&gt; &nbsp; { "OB" &nbsp;"OP2" }
> &gt;&gt; &nbsp; { "O1P" "OP1" }
> &gt;&gt; &nbsp; { "O2P" "OP2" }
> &gt;&gt; }
> &gt;&gt;&nbsp;
> &gt;&gt; #
> &gt;&gt; # Load the OL15 force field for DNA
> &gt;&gt; #&nbsp;
> &gt;&gt;&nbsp;
> &gt;&gt; addAtomTypes {
> &gt;&gt; &nbsp; &nbsp; { "CJ" &nbsp;"C" "sp3" }
> &gt;&gt; &nbsp; &nbsp; { "C1" &nbsp;"C" "sp2" }
> &gt;&gt; &nbsp; &nbsp; { "C2" &nbsp;"C" "sp2" }
> &gt;&gt; &nbsp; &nbsp; { "C7" &nbsp;"C" "sp3" }
> &gt;&gt; }
> &gt;&gt;&nbsp;
> &gt;&gt; loadoff DNA.OL15.lib&nbsp;
> Loading library: /home/xu/anaconda3/envs/ambertools23/dat/leap/lib/DNA.OL15.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
> &gt;&gt; parm10 = loadamberparams parm10.dat
> Loading parameters: /home/xu/anaconda3/envs/ambertools23/dat/leap/parm/parm10.dat
> Reading title:
> PARM99 + frcmod.ff99SB + frcmod.parmbsc0 + OL3 for RNA
> &gt;&gt; loadamberparams frcmod.DNA.OL15
> Loading parameters: /home/xu/anaconda3/envs/ambertools23/dat/leap/parm/frcmod.DNA.OL15
> Reading force field modification type file (frcmod)
> Reading title:
> OL15 force field for DNA (99bsc0-betaOL1-eps-zetaOL1-chiOL4) see http://ffol.upol.cz
> &gt;&nbsp;
> &gt; lig = loadpdb &nbsp; lig1.pdb
> Loading PDB file: ./lig1.pdb
> Matching PDB residue names to LEaP variables.
> Mapped residue DT, term: Terminal/beginning, seq. number: 0 to: DT5.
> Mapped residue DC, term: Terminal/last, seq. number: 19 to: DC3.
> Created a new atom named: P within residue: .R<DT5 1&gt;
> Created a new atom named: OP1 within residue: .R<DT5 1&gt;
> Created a new atom named: OP2 within residue: .R<DT5 1&gt;
> Created a new atom named: OP3 within residue: .R<DT5 1&gt;
> &nbsp; total atoms in file: 412
> &nbsp; Leap added 229 missing atoms according to residue templates:
> &nbsp; &nbsp; &nbsp; &nbsp;229 H / lone pairs
> &nbsp; The file contained 4 atoms not in residue templates
> &gt; charge lig&nbsp;
> Total unperturbed charge: -19.000000
> Total perturbed charge: &nbsp; -19.000000
> &gt; saveamberparm lig lig.parm7 &nbsp;lig.rst7
> Checking Unit.
>
>
> /home/xu/anaconda3/envs/ambertools23/bin/teLeap: Warning!
> The unperturbed charge of the unit (-19.000000) is not zero.
> FATAL: &nbsp;Atom .R<DT5 1&gt;.A<P 31&gt; does not have a type.
> FATAL: &nbsp;Atom .R<DT5 1&gt;.A<OP1 32&gt; does not have a type.
> FATAL: &nbsp;Atom .R<DT5 1&gt;.A<OP2 33&gt; does not have a type.
> FATAL: &nbsp;Atom .R<DT5 1&gt;.A<OP3 34&gt; does not have a type.
>
>
> /home/xu/anaconda3/envs/ambertools23/bin/teLeap: Fatal Error!
> Failed to generate parameters
>
>
> Exiting LEaP: Errors = 1; Warnings = 1; Notes = 0.
>
>
>
>
>
>
> 1464903098
> xx1464903098.qq.com_______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber

_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue Mar 18 2025 - 19:00:03 PDT
Custom Search