C_beta = Atom('CH2')
C_delta_1 = Atom('C')
C_delta_2 = Atom('C')
C_epsilon_1 = Atom('C')
C_epsilon_2 = Atom('C')
C_gamma = Atom('C')
C_zeta = Atom('C')
H_delta_1 = Atom('H')
H_delta_2 = Atom('H')
H_epsilon_1 = Atom('H')
H_epsilon_2 = Atom('H')
H_zeta = 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(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), Bond(C_zeta,H_zeta)]

pdbmap = [('PHE', {'2HD': H_delta_2, '2HE': H_epsilon_2, '1HD': H_delta_1, '1HE': H_epsilon_1, 'CE1': C_epsilon_1, 'CE2': C_epsilon_2, 'CZ': C_zeta, 'CG': C_gamma, 'CD1': C_delta_1, 'CB': C_beta, 'CD2': C_delta_2, 'HZ': H_zeta, }, ), ]
pdb_alternative = {'HD1': '1HD', 'HD2': '2HD', 'HE1': '1HE', 'HE2': '2HE', }

name = 'phe_sidechain'

opls_atom_type = { C_beta: 'C2', C_gamma: 'CK', C_delta_1: 'CK', C_epsilon_1: 'CK', C_zeta: 'CK', C_epsilon_2: 'CK', C_delta_2: 'CK', H_delta_1: 'HK', H_delta_2: 'HK', H_epsilon_1: 'HK', H_epsilon_2: 'HK', H_zeta: 'HK'  }

opls_charge = { C_beta: .115, C_gamma: -.115, C_delta_1: -.115, C_epsilon_1: -.115, C_zeta: -.115, C_epsilon_2: -.115, C_delta_2: -.115, H_delta_1: .115, H_delta_2: .115, H_epsilon_1: .115, H_epsilon_2: .115, H_zeta: .115  }
