Question: Protein Sequence Analyses
gravatar for l.souza
3.0 years ago by
Brasilia, Brazil
l.souza70 wrote:

Hello everybody,

I have ~2000 sequences of a capsid protein from different serotypes of a virus. I aligned them to check the identity, and now I wanna study the protein properties, such aa composition, isoelectric point, molecular weight, GRAVY, instability index, etc. Should I analyze each of the sequences I have, or I can use a consensus one? (If using the consensus, what is the best software to get it?)

I really appreciate any help!

protein sequence consensus • 1.1k views
ADD COMMENTlink modified 3.0 years ago by Elisabeth Gasteiger1.7k • written 3.0 years ago by l.souza70
gravatar for Elisabeth Gasteiger
3.0 years ago by
Elisabeth Gasteiger1.7k wrote:

Have you tried the ProtParam tool, ?

ADD COMMENTlink written 3.0 years ago by Elisabeth Gasteiger1.7k

Yes! I usually use this when analyzing physical properties of proteins...

ADD REPLYlink written 3.0 years ago by l.souza70
gravatar for Joe
3.0 years ago by
United Kingdom
Joe17k wrote:

I don't think there's a simple answer to this, and I'd encourage you to think about what you actually want to know about the protein. It's all well and good to say that the pI of your capsid protein is x, but is that useful? Do you care about minute fluctuations in pI which you would almost certainly lose by taking a consensus?

A consensus is useful in certain situations but not others.


Propka can give you pI for your sequences, and if you have that many, you might want to run just a consensus through - I'd say this only makes sense if your 2000 sequences are already pretty similar however.

GRAVY etc:

I believe CodonW will output a lot of parameters, including GRAVY and AA composition, but I'm not super familiar with the program.

Consensus sequences:

Consensus sequences can be gotten out of BioPython (I've actually been writing this myself recently. It's what they call a 'dumb' consensus though, it doesn't have many configurable parameters. See

from Bio.Align import AlignInfo
from Bio import AlignIO

alignment ='~/path/to/alignment.aln', 'format') # where format is a supported type (see BioPython docs)
summary_align = AlignInfo.SummaryInfo(alignment)
consensus = summary_align.dumb_consensus()
ADD COMMENTlink modified 3.0 years ago • written 3.0 years ago by Joe17k
Please log in to add an answer.


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