Question: How To Calculate The Net Charge Of Amino Acids *Sequences/Peptides* At Ph 7 Using R?
0
gravatar for bioinfo12
21 months ago by
bioinfo120
bioinfo120 wrote:

Hi,

I would like to calculate the net charge of amino acids at pH 7. Is there any R package available for this calculation? Edit: op wants the charge of a sequence, not single AA

ADD COMMENTlink modified 21 months ago by Neilfws41k • written 21 months ago by bioinfo120
3
gravatar for Michael Dondrup
21 months ago by
Bergen
Michael Dondrup27k wrote:

There is nothing to calculate, this is textbook knowledge about amino acids, there are amino acids with charged and uncharged side chains. Just look up the charges in Wikipedia http://en.wikipedia.org/wiki/Amino_acid

Edit: you can simply count the positively and negatively charged AAs in the sequence, the difference is the charge of the sequence.

ADD COMMENTlink modified 21 months ago • written 21 months ago by Michael Dondrup27k

@Michael I am really sorry for my incorrect question.

ADD REPLYlink written 21 months ago by bioinfo120
2
gravatar for seidel
21 months ago by
seidel4.5k
United States
seidel4.5k wrote:

I think you're looking for either of two things. I disagree with Michael that there is nothing to calculate. On one hand, you could make a sort of lookup table using a named array of the amino acids and their pK values, then at any given pH you could access the value and assess the direction of charge. However you might actually be asking something else. At any given pH, only a fraction of any amino acids exist in a given charged state, so perhaps by "net charge" you're looking for the relatives proportion of charge at that pH? In that case I'm not aware of a package (I haven't looked), but it wouldn't be that hard to write a function for the Henderson-Hasselbalch equation (complete with amino acid pK lookup table), that returns this value for any pH. It would only be a few lines of simple code.

ADD COMMENTlink written 21 months ago by seidel4.5k

Thank you for your answer. I would like to calculate the net charge of amino acid sequence. I am not a programmer. Is there any software available to calculate the net charge from the amino acid sequences?

ADD REPLYlink written 21 months ago by bioinfo120

You should have written that immediately in your question. That is ofc something very different.

ADD REPLYlink written 21 months ago by Michael Dondrup27k

I thought as well too that bionfo12 meant net charge of an peptide. Pathetic with poorly framed questions!! What is really sad is that there are people out here trying to help others out and in the end it turns out that it has been a waste of time for them just because of a poorly framed question!!

ADD REPLYlink modified 21 months ago • written 21 months ago by abi70

I don't bother too much, people are themselves responsible to pose their question correctly. If they don't the answers are not going to help them, if they really need a good answer, they will come back, if not, no damage done except for the poster.

I only answer if I have time to spare and I know misunderstanding can happen, so for me it's not a waste of time.

ADD REPLYlink modified 21 months ago • written 21 months ago by Michael Dondrup27k

You could speculate a lot about what a question really means. I prefer to not do that and answer the literal question if possible. In this case the original question was quite clearly phrased, but unfortunately, the terms were used incorrectly. So, bad question lead to incorrect answer, that's what happened.

ADD REPLYlink written 21 months ago by Michael Dondrup27k
2
gravatar for Neilfws
21 months ago by
Neilfws41k
Sydney, Australia
Neilfws41k wrote:

You can calculate net charge of a protein sequence using the seqinr package. See the example function computeCharge() on page 133 of the manual (PDF).

ADD COMMENTlink modified 21 months ago • written 21 months ago by Neilfws41k

Good, that's more reliable than trying to code this yourself! There should be functions in bioperl and bio python too.

ADD REPLYlink written 21 months ago by Michael Dondrup27k

I expect so, but OP asked for an R package.

ADD REPLYlink written 21 months ago by Neilfws41k
Please log in to add an answer.

Help
Access
  • RSS
  • Stats
  • API

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.0.0
Traffic: 404 users visited in the last hour