name ='guanine'

N_9 = Atom('N')
C_8 = Atom('C')
H_8 = Atom('H')
N_7 = Atom('N')
C_5 = Atom('C')
C_6 = Atom('C')
O_6 = Atom('O')
N_1 = Atom('N')
H_1 = Atom('H')
C_2 = Atom('C')
N_2 = Atom('N')
H_21 = Atom('H')
H_22 = 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(O_6, C_6), Bond(N_1, C_6), Bond(H_1, N_1), Bond(C_2, N_1), Bond(N_2, C_2), Bond(H_21, N_2), Bond(H_22, N_2), Bond(N_3, C_2), Bond(C_4, N_3), Bond(C_4, C_5), Bond(C_4, N_9), ]

pdbmap = [('G', {"N2": N_2, "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, "H22": H_22, "H1": H_1, "N9": N_9, "O6": O_6, "H21": H_21, })]

amber_atom_type = {N_2: 'N2', C_8: 'CK', C_6: 'C', C_4: 'CB', C_5: 'CB', C_2: 'CA', N_3: 'NC', H_8: 'H5', N_1: 'NA', N_7: 'NB', H_22: 'H', H_1: 'H', N_9: 'N*', O_6: 'O', H_21: 'H', }

amber12_atom_type = amber_atom_type
