name ='uracil'

N_1 = Atom('N')
C_6 = Atom('C')
H_6 = Atom('H')
C_5 = Atom('C')
H_5 = Atom('H')
C_4 = Atom('C')
O_4 = Atom('O')
N_3 = Atom('N')
H_3 = Atom('H')
C_2 = Atom('C')
O_2 = Atom('O')

bonds = [Bond(C_6, N_1), Bond(H_6, C_6), Bond(C_5, C_6), Bond(H_5, C_5), Bond(C_4, C_5), Bond(O_4, C_4), Bond(N_3, C_4), Bond(H_3, N_3), Bond(C_2, N_3), Bond(O_2, C_2), Bond(C_2, N_1), ]

pdbmap = [('U', {"N3": N_3, "C2": C_2, "N1": N_1, "H3": H_3, "C6": C_6, "C4": C_4, "C5": C_5, "H5": H_5, "O2": O_2, "O4": O_4, "H6": H_6, })]

amber_atom_type = {N_3: 'NA', C_2: 'C', N_1: 'N*', H_3: 'H', C_6: 'CM', C_4: 'C', C_5: 'CM', H_5: 'HA', O_2: 'O', O_4: 'O', H_6: 'H4', }
amber12_atom_type = {N_3: 'NA', C_2: 'C', N_1: 'N*', H_3: 'H', C_6: 'CS', C_4: 'C', C_5: 'CS', H_5: 'HA', O_2: 'O', O_4: 'O', H_6: 'H4', }

