:>"^:@sdZddddddddiZd d d d d dddddiZd dd diZddddddddddddd d iZd$ZGd!d"d"eZd#S)%aQCalculate isoelectric points of polypeptides using methods of Bjellqvist. pK values and the methos are taken from:: * Bjellqvist, B.,Hughes, G.J., Pasquali, Ch., Paquet, N., Ravier, F., Sanchez, J.-Ch., Frutiger, S. & Hochstrasser, D.F. The focusing positions of polypeptides in immobilized pH gradients can be predicted from their amino acid sequences. Electrophoresis 1993, 14, 1023-1031. * Bjellqvist, B., Basse, B., Olsen, E. and Celis, J.E. Reference points for comparisons of two-dimensional maps of proteins from different human cell types defined in a pH scale where isoelectric points correlate with polypeptide compositions. Electrophoresis 1994, 15, 529-539. I designed the algorithm according to a note by David L. Tabb, available at: http://fields.scripps.edu/DTASelect/20010710-pI-Algorithm.pdf Ntermg@Kg$@Rg(@HgQ@Ctermgffffff @Dg333333@Eg@Cg"@Yg333333@g@Ag\(\@Mg@SgQ@PgQ @TgHzG@Vg(\@g@c@saeZdZdZdddZddZddZd d Zd d Zd dZ dS)IsoelectricPointaA class for calculating the IEP or charge at given pH of a protein. Parameters ---------- :protein_sequence: A ``Bio.Seq`` or string object containing a protein sequence. :aa_content: A dictionary with amino acid letters as keys and it's occurences as integers, e.g. ``{"A": 3, "C": 0, ...}``. Default: ``None``. If ``None``, the dic will be calculated from the given sequence. Methods ------- :charge_at_pH(pH): Calculates the charge of the protein for a given pH :pi(): Calculates the isoelectric point Examples -------- The methods of this class can either be accessed from the class itself or from a ``ProtParam.ProteinAnalysis`` object (with partially different names): >>> from Bio.SeqUtils.IsoelectricPoint import IsoelectricPoint as IP >>> protein = IP("INGAR") >>> print("IEP of peptide {} is {:.2f}" ... .format(protein.sequence, protein.pi())) IEP of peptide INGAR is 9.75 >>> print("It's charge at pH 7 is {:.2f}" ... .format(protein.charge_at_pH(7.0))) It's charge at pH 7 is 0.76 >>> from Bio.SeqUtils.ProtParam import ProteinAnalysis as PA >>> protein = PA("PETER") >>> print("IEP of {}: {:.2f}".format(protein.sequence, ... protein.isoelectric_point())) IEP of PETER: 4.53 >>> print("Charge at pH 4.53: {:.2f}" ... .format(protein.charge_at_pH(4.53))) Charge at pH 4.53: 0.00 NcCsnt|j|_|s@ddlm}||jj}|j||_|j\|_ |_ dS)zInitialize the class.)ProteinAnalysisN) struppersequenceZBio.SeqUtils.ProtParamrZcount_amino_acids_select_chargedcharged_aas_content_update_pKs_tablespos_pKsneg_pKs)selfZprotein_sequence aa_contentZ_PArB/tmp/pip-build-ww9dw3qa/biopython/Bio/SeqUtils/IsoelectricPoint.py__init__Us zIsoelectricPoint.__init__cCsCi}x"tD]}t||||s$