C_beta = Atom('CH2')
C_delta_1 = Atom('CH')
C_delta_2 = Atom('CH')
C_epsilon_1 = Atom('CH')
C_epsilon_2 = Atom('CH')
C_gamma = Atom('C')
C_zeta = Atom('C')
H_eta = Atom('H')
O_eta = Atom('O')
H_delta_1 = Atom('H')
H_delta_2 = Atom('H')
H_epsilon_1 = Atom('H')
H_epsilon_2 = Atom('H')

bonds = [Bond(C_gamma, C_beta), Bond(C_delta_1, C_gamma), Bond(C_epsilon_1, C_delta_1), Bond(C_zeta, C_epsilon_1), Bond(O_eta, C_zeta), Bond(H_eta, O_eta), Bond(C_epsilon_2, C_zeta), Bond(C_delta_2, C_epsilon_2), Bond(C_delta_2, C_gamma), Bond(C_delta_1,H_delta_1), Bond(C_delta_2,H_delta_2), Bond(C_epsilon_1,H_epsilon_1), Bond(C_epsilon_2,H_epsilon_2)]
opls_atom_type = {C_gamma: 'CK', C_zeta: 'CK', C_epsilon_2: 'CK', C_delta_2: 'CK', C_delta_1: 'CK', H_eta: 'HO', C_beta: 'C2', O_eta: 'OH', C_epsilon_1: 'CK', H_epsilon_1: 'HK', H_epsilon_2: 'HK', H_delta_1: 'HK', H_delta_2: 'HK' }
opls_charge = {C_gamma: -.115, C_zeta: .150, C_epsilon_2: -.115, C_delta_2: -.115, C_delta_1: -.115, H_eta: .435, C_beta: .115, O_eta: -.585, C_epsilon_1: -.115, H_epsilon_1: .115, H_epsilon_2: .115, H_delta_1: .115, H_delta_2: .115 }
name = 'tyr_sidechain'
pdbmap = [('TYR', {'2HD': H_delta_2, '2HE': H_epsilon_2, '1HD': H_delta_1, '1HE': H_epsilon_1, 'CE2': C_epsilon_2, 'CE1': C_epsilon_1, 'CD1': C_delta_1, 'HH': H_eta, 'CZ': C_zeta, 'CG': C_gamma, 'CB': C_beta, 'OH': O_eta, 'CD2': C_delta_2, }, ), ]
pdb_alternative = {'HD1': '1HD', 'HD2': '2HD', 'HE1': '1HE', 'HE2': '2HE',
                   'HOH': 'HH'}
