sidechain = Group('met_sidechain_uni')
peptide = Group('peptide_uni')
bonds = [Bond(peptide.C_alpha, sidechain.C_beta), ]
symbol = 'Met'
name = 'methionine'
chain_links = [peptide.N, peptide.C]
# amber91_charge = {sidechain.LP_1 ???: -0.381, peptide.O: -0.5, peptide.N: -0.52, peptide.C_alpha: 0.137, sidechain.C_beta: 0.037, sidechain.S_delta: 0.737, sidechain.LP_2 ???: -0.381, sidechain.C_epsilon: 0.007, sidechain.C_gamma: 0.09, peptide.H: 0.248, peptide.C: 0.526, }

