sidechain = Group('lys_sidechain_neutral')
peptide = Group('peptide')
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.C: 0.59730, peptide.H_alpha: 0.09940, peptide.O: -0.56790,
    peptide.H: 0.27190, peptide.C_alpha: -0.07206,peptide.N: -0.41570,}

name = 'lysine_neutral'
