sidechain = Group('lys_sidechain_neutral')
peptide = Group('peptide_ct')
bonds = [Bond(peptide.C_alpha, sidechain.C_beta), ]
symbol = 'Lyp'

amber_charge = {
    sidechain.C_delta: -0.03768,
    sidechain.H_delta_3: 0.01155,
    sidechain.H_delta_2: 0.01155,
    sidechain.H_gamma_2: 0.01041,
    sidechain.H_zeta_3: 0.38604,
    sidechain.N_zeta: -1.03581,
    sidechain.H_epsilon_2: -0.03358,
    sidechain.C_beta: -0.04845,
    sidechain.H_epsilon_3: -0.03358,
    sidechain.H_beta_3: 0.034,
    sidechain.C_epsilon: 0.32604,
    sidechain.C_gamma: 0.06612,
    sidechain.H_beta_2: 0.034,
    sidechain.H_zeta_2: 0.38604,
    sidechain.H_gamma_3: 0.01041,
    peptide.O: -0.8252,
    peptide.H_alpha: 0.1438,
    peptide.N: -0.3481,
    peptide.O_2: -0.8252,
    peptide.H: 0.2764,
    peptide.C: 0.8488,
    peptide.C_alpha: -0.2903,
    }

name = 'lysine_neutral'

chain_links = [peptide.N, peptide.C]
