def chargeandprop(AAseq):
    """ Returns the net charge of a protein sequence
    and proportion of charged amino acids 
    """ 
    protseq = AAseq.upper()
    charge = -0.002
    cp = 0
    AACharge={"C":-.045,"D":-.999,"E":-.998,"H":.091,
              "K":1,"R":1,"Y":-.001}
    for aa in protseq:
        charge += AACharge.get(aa,0)
        if aa in AACharge:
            cp += 1
    prop = 100.*cp/len(AAseq)
    return (charge,prop)

This code is part of the book "Python for Bioinformatics", by Sebastian Bassi (sbassi@genesdigitales.com). Return to home page.