name = 'cysteine_with_negative_charge'
symbol = 'Cym'

peptide = Group('peptide')
sidechain = Group('cym_sidechain')

bonds = [Bond(peptide.C_alpha, sidechain.C_beta)]

amber_charge = {peptide.N: -0.463, peptide.H: 0.252, peptide.C_alpha: 0.035, peptide.H_alpha: 0.048, peptide.C: 0.616, peptide.O: -0.504, sidechain.C_beta: -0.736, sidechain.H_beta_3: 0.244, sidechain.H_beta_2: 0.244, sidechain.S_gamma: -0.736, }

chain_links = [peptide.N, peptide.C]
