name ='adenine'

N_9 = Atom('N')
C_8 = Atom('C')
H_8 = Atom('H')
N_7 = Atom('N')
C_5 = Atom('C')
C_6 = Atom('C')
N_6 = Atom('N')
H_61 = Atom('H')
H_62 = Atom('H')
N_1 = Atom('N')
C_2 = Atom('C')
H_2 = Atom('H')
N_3 = Atom('N')
C_4 = Atom('C')

bonds = [Bond(C_8, N_9), Bond(H_8, C_8), Bond(N_7, C_8), Bond(C_5, N_7), Bond(C_6, C_5), Bond(N_6, C_6), Bond(H_61, N_6), Bond(H_62, N_6), Bond(N_1, C_6), Bond(C_2, N_1), Bond(H_2, C_2), Bond(N_3, C_2), Bond(C_4, N_3), Bond(C_4, C_5), Bond(C_4, N_9), ]

pdbmap = [('A', {"H62": H_62, "H61": H_61, "C8": C_8, "C6": C_6, "C4": C_4, "C5": C_5, "C2": C_2, "N3": N_3, "H8": H_8, "N1": N_1, "N7": N_7, "N6": N_6, "N9": N_9, "H2": H_2, })]

amber_atom_type = {H_62: 'H', H_61: 'H', C_8: 'CK', C_6: 'CA', C_4: 'CB', C_5: 'CB', C_2: 'CQ', N_3: 'NC', H_8: 'H5', N_1: 'NC', N_7: 'NB', N_6: 'N2', N_9: 'N*', H_2: 'H5', }

