C_beta = Atom('CH2')
C_delta_2 = Atom('CH')
C_epsilon_1 = Atom('CH')
C_gamma = Atom('C')
N_delta_1 = Atom('N')
H_delta_1 = Atom('H')
H_delta_2 = Atom('H')
N_epsilon_2 = Atom('N')
H_epsilon_1 = Atom('H')

bonds = [Bond(C_gamma, C_beta), Bond(N_delta_1, C_gamma), Bond(H_delta_1, N_delta_1), Bond(C_epsilon_1, N_delta_1), Bond(N_epsilon_2, C_epsilon_1), Bond(C_delta_2, N_epsilon_2), Bond(C_gamma, C_delta_2), Bond(C_delta_2,H_delta_2), Bond(H_epsilon_1,C_epsilon_1)]

pdbmap = [('HID', {'CG': C_gamma, 'CE1': C_epsilon_1, 'CB': C_beta, 'ND1': N_delta_1, '1HD': H_delta_1, 'CD2': C_delta_2, 'NE2': N_epsilon_2, '2HD': H_delta_2, '1HE': H_epsilon_1 }, ), ]

pdb_alternative = {'HD1': '1HD', 'HB1': '3HB', 'HND': '1HD'}

name = 'hid_sidechain'

opls_atom_type = { C_beta: 'C2', C_gamma: 'CK', N_delta_1: 'NA', H_delta_1: 'H', C_epsilon_1: 'CK', N_epsilon_2: 'NB', C_delta_2: 'CK', H_epsilon_1: 'HK', H_delta_2: 'HK'  }

opls_charge = { C_beta: .115, C_gamma: 0.015, N_delta_1: -0.57, H_delta_1: 0.42, C_epsilon_1: 0.295, N_epsilon_2: -0.49, C_delta_2: -.015, H_epsilon_1: .115, H_delta_2: .115 }
