sidechain = Group('lys_sidechain_neutral')
peptide = Group('peptide_nt')
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_alpha: -0.0015,
    peptide.H_alpha: 0.118,
    peptide.H_1: 0.2165,
    peptide.H_2: 0.2165,
    peptide.C: 0.7214,
    peptide.H_3: 0.2165,
    peptide.O: -0.6013,
    peptide.N: 0.0966,}

name = 'lysine_neutral'

chain_links = [peptide.N, peptide.C]
