o RŀgI4@sxdZddlZddlmZddlmZddlmZddlmZddlm Z Gdd d Z e d kr:dd l m Z e dSdS) aSimple protein analysis. Examples -------- >>> from Bio.SeqUtils.ProtParam import ProteinAnalysis >>> X = ProteinAnalysis("MAEGEITTFTALTEKFNLPPGNYKKPKLLYCSNGGHFLRILPDGTVDGT" ... "RDRSDQHIQLQLSAESVGEVYIKSTETGQYLAMDTSGLLYGSQTPSEEC" ... "LFLERLEENHYNTYTSKKHAEKNWFVGLKKNGSCKRGPRTHYGQKAILF" ... "LPLPV") >>> print(X.count_amino_acids()['A']) 6 >>> print(X.count_amino_acids()['E']) 12 >>> print("%0.2f" % X.get_amino_acids_percent()['A']) 0.04 >>> print("%0.2f" % X.get_amino_acids_percent()['L']) 0.12 >>> print("%0.2f" % X.molecular_weight()) 17103.16 >>> print("%0.2f" % X.aromaticity()) 0.10 >>> print("%0.2f" % X.instability_index()) 41.98 >>> print("%0.2f" % X.isoelectric_point()) 7.72 >>> sec_struc = X.secondary_structure_fraction() # [helix, turn, sheet] >>> print("%0.2f" % sec_struc[0]) # helix 0.33 >>> print("%0.2f" % sec_struc[1]) # turn 0.29 >>> print("%0.2f" % sec_struc[2]) # sheet 0.37 >>> epsilon_prot = X.molar_extinction_coefficient() # [reduced, oxidized] >>> print(epsilon_prot[0]) # with reduced cysteines 17420 >>> print(epsilon_prot[1]) # with disulfid bridges 17545 Other public methods are: - gravy - protein_scale - flexibility - charge_at_pH N) IUPACData)Seq)IsoelectricPoint)molecular_weight) ProtParamDatac@seZdZdZd"ddZddZddZd d Zd d Zd dZ ddZ d#ddZ ddZ d$ddZ ddZddZddZdd Zd!S)%ProteinAnalysisaPClass containing methods for protein analysis. The constructor takes two arguments. The first is the protein sequence as a string or a Seq object. The second argument is optional. If set to True, the weight of the amino acids will be calculated using their monoisotopic mass (the weight of the most abundant isotopes for each element), instead of the average molecular mass (the averaged weight of all stable isotopes for each element). If set to false (the default value) or left out, the IUPAC average molecular mass will be used for the calculation. FcCs,||_d|_d|_t|j|_||_dS)zInitialize the class.N)uppersequenceamino_acids_contentamino_acids_percentlenlength monoisotopic)self prot_sequencerrJ/var/www/html/myenv/lib/python3.10/site-packages/Bio/SeqUtils/ProtParam.py__init__Ls   zProteinAnalysis.__init__cCs@|jdurddtjD}|D] }|j|||<q||_|jS)a(Count standard amino acids, return a dict. Counts the number times each amino acid is in the protein sequence. Returns a dictionary {AminoAcid:Number}. The return value is cached in self.amino_acids_content. It is not recalculated upon subsequent calls. NcSsi|]}|dqS)rr).0krrr ^sz5ProteinAnalysis.count_amino_acids..)r rprotein_lettersr count)rprot_dicaarrrcount_amino_acidsTs z!ProteinAnalysis.count_amino_acidscs4jdur}fdd|D}|_jS)aCalculate the amino acid content in percentages. The same as count_amino_acids only returns the Number in percentage of entire sequence. Returns a dictionary of {AminoAcid:percentage}. The return value is cached in self.amino_acids_percent. input is the dictionary self.amino_acids_content. output is a dictionary with amino acids as keys. Ncsi|] \}}||jqSr)r )rrrrrrrtsz;ProteinAnalysis.get_amino_acids_percent..)r ritems)r aa_counts percentagesrrrget_amino_acids_percentfs z'ProteinAnalysis.get_amino_acids_percentcCst|jd|jdS)z#Calculate MW from Protein sequence.protein)seq_typer)rr rrrrrrzs z ProteinAnalysis.molecular_weightcs&d}|tfdd|D}|S)zCalculate the aromaticity according to Lobry, 1994. Calculates the aromaticity value of a protein according to Lobry, 1994. It is simply the relative frequency of Phe+Trp+Tyr. YWFc3|]}|VqdSNrrraa_percentagesrr z.ProteinAnalysis.aromaticity..r sum)r aromatic_aas aromaticityrr'rr.szProteinAnalysis.aromaticitycCsVtj}d}t|jdD]}|j||d\}}|||}||7}q d|j|S)akCalculate the instability index according to Guruprasad et al 1990. Implementation of the method of Guruprasad et al. 1990 to test a protein for stability. Any value above 40 means the protein is unstable (has a short half life). See: Guruprasad K., Reddy B.V.B., Pandit M.W. Protein Engineering 4:155-161(1990). g$@)rDIWVranger r )rindexscoreithisnextdipeptide_valuerrrinstability_indexs   z!ProteinAnalysis.instability_indexc Cstj}d}gd}g}t|j|D]E}|j|||}d}t|dD]}||} |||d} ||| || ||7}q%||dd} ||| 7}||dq|S)zCalculate the flexibility according to Vihinen, 1994. No argument to change window size because parameters are specific for a window=9. The parameters used are optimized for determining the flexibility. )g?g?g?g?r0r/r1r0g@)rFlexr3r r append) r flexibilities window_sizeweightsscoresr6 subsequencer5jfrontbackmiddlerrr flexibilitys zProteinAnalysis.flexibility KyteDoolitlecsHtj|ddkrtd|dtfdd|jD}||jS)aBCalculate the GRAVY (Grand Average of Hydropathy) according to Kyte and Doolitle, 1982. Utilizes the given Hydrophobicity scale, by default uses the original proposed by Kyte and Doolittle (KyteDoolitle). Other options are: Aboderin, AbrahamLeo, Argos, BlackMould, BullBreese, Casari, Cid, Cowan3.4, Cowan7.5, Eisenberg, Engelman, Fasman, Fauchere, GoldSack, Guy, Jones, Juretic, Kidera, Miyazawa, Parker,Ponnuswamy, Rose, Roseman, Sweet, Tanford, Wilson and Zimmerman. New scales can be added in ProtParamData. zscale: z not knownc3r$r%rr&selected_scalerrr)r*z(ProteinAnalysis.gravy..)r gravy_scalesget ValueErrorr,r r )rscale total_gravyrrJrgravys   zProteinAnalysis.gravycCsHdd||d}dg|d}t|dD] }|||||<q|S)a1Make list of relative weight of window edges (PRIVATE). The relative weight of window edges are compared to the window center. The weights are linear. It actually generates half a list. For a window of size 9 and edge 0.4 you get a list of [0.4, 0.55, 0.7, 0.85]. r1?r0r/)r3)rwindowedgeunitr@r6rrr _weight_lists zProteinAnalysis._weight_listrRc Cs|||}g}t|dd}t|j|dD]s}|j|||}d} t|dD]=} z ||| } |||| d} | || | || | 7} Wq,tyitjd|| ||| dfYq,w||d} | |vr{| || 7} n tjd| d| | |q|S)aCompute a profile by any amino acid scale. An amino acid scale is defined by a numerical value assigned to each type of amino acid. The most frequently used scales are the hydrophobicity or hydrophilicity scales and the secondary structure conformational parameters scales, but many other scales exist which are based on different chemical and physical properties of the amino acids. You can set several parameters that control the computation of a scale profile, such as the window size and the window edge relative weight value. WindowSize: The window size is the length of the interval to use for the profile computation. For a window size n, we use the i-(n-1)/2 neighboring residues on each side to compute the score for residue i. The score for residue i is the sum of the scaled values for these amino acids, optionally weighted according to their position in the window. Edge: The central amino acid of the window always has a weight of 1. By default, the amino acids at the remaining window positions have the same weight, but you can make the residue at the center of the window have a larger weight than the others by setting the edge value for the residues at the beginning and end of the interval to a value between 0 and 1. For instance, for Edge=0.4 and a window size of 5 the weights will be: 0.4, 0.7, 1.0, 0.7, 0.4. The method returns a list of values which can be plotted to view the change along a protein sequence. Many scales exist. Just add your favorites to the ProtParamData modules. Similar to expasy's ProtScale: http://www.expasy.org/cgi-bin/protscale.pl r1r0r/z0warning: %s or %s is not a standard amino acid. z warning: z is not a standard amino acid. ) rVr,r3r r KeyErrorsysstderrwriter=)r param_dictrSrTr@rAsum_of_weightsr6rBr5rCrDrErFrrr protein_scales0 &    zProteinAnalysis.protein_scalecCs|}t|j|}|S)zuCalculate the isoelectric point. Uses the module IsoelectricPoint to calculate the pI of a protein. )rrr pi)r aa_contentie_pointrrrisoelectric_point,sz!ProteinAnalysis.isoelectric_pointcCs |}t|j|}||S)z.Calculate the charge of a protein at given pH.)rrr charge_at_pH)rpHr_chargerrrrb6s zProteinAnalysis.charge_at_pHcsT|tfdddD}tfdddD}tfdddD}|||fS)akCalculate fraction of helix, turn and sheet. Returns a list of the fraction of amino acids which tend to be in Helix, Turn or Sheet, according to Haimov and Srebnik, 2016; Hutchinson and Thornton, 1994; and Kim and Berg, 1993, respectively. Amino acids in helix: E, M, A, L, K. Amino acids in turn: N, P, G, S, D. Amino acids in sheet: V, I, Y, F, W, L, T. Note that, prior to v1.82, this method wrongly returned (Sheet, Turn, Helix) while claiming to return (Helix, Turn, Sheet). Returns a tuple of three floats (Helix, Turn, Sheet). c3r$r%rrrr'rrr)Nr*z?ProteinAnalysis.secondary_structure_fraction..EMALKc3r$r%rrer'rrr)Or*NPGSDc3r$r%rrer'rrr)Pr*VIYFWLTr+)rhelixturnsheetrr'rsecondary_structure_fraction<s  z,ProteinAnalysis.secondary_structure_fractioncCs<|}|dd|dd}||ddd}||fS)zCalculate the molar extinction coefficient. Calculates the molar extinction coefficient assuming cysteines (reduced) and cystines residues (Cys-Cys-bond) Wi|YiCr1})r)rnum_aa mec_reduced mec_cystinesrrrmolar_extinction_coefficientTsz,ProteinAnalysis.molar_extinction_coefficientN)F)rH)rR)__name__ __module__ __qualname____doc__rrr rr.r:rGrQrVr]rarbrmrurrrrr=s     J  r__main__) run_doctest)ryrXBio.DatarBio.Seqr Bio.SeqUtilsrrrrrv Bio._utilsr{rrrrs.     %