Word of caution: The pI depends on the sum of the (unknown) local pKa, and various calculators have varying opinion about the 'normal' pKa values. It would be good to state which source of pKa's is used.
BR -----Original Message----- From: CCP4 bulletin board [mailto:[email protected]] On Behalf Of Paul Smith Sent: Monday, May 04, 2009 12:00 PM To: [email protected] Subject: [ccp4bb] Charge vs. pH plot Isoelectric titration curve ascii plotter Hello all, I had the worst time the other day finding an online tool to calculate charge versus pH values for an input protein sequence. There are many programs for calculating the isoelectric point (pI), put not many that output a complete titration curve. Because such information is often useful in protein purification and crystallization, I coded up a quick perl program I would like to share with the CCP4 community. The code is attached. It is pretty self-explanatory. Input a FASTA file as a command line arguement and the program outputs charge vs. pH value pairs (unsorted). The constants defined in the program header can be used to alter the pKa values used, titration starting and ending points, and grid sampling. Please feel free to use, distribute, and modify freely. Also attached is a ASCII plotting tool that takes paired values from a file or STDIN and plots them to STDOUT in ascii using automatically determining ranges and binning values. These two programs work well like so: "./titrate.pl FASTA.seq | ./plot.pl" to produce a nice ascii plot of charge vs. pH. I hope these are helpful to someone. Best, Paul Smith Paladin Scientific paladinscientific.com
