Performing Needleman-Wunsch Global Alignment for multiple mutant proteins against wild-type
2
0
Entering edit mode
29 days ago
VolEr ▴ 40

Hi all,

I have a list of 60 different insertion and deletion mutant proteins, and I’d like to perform a multiple sequence alignment comparing each mutant with the wild-type sequence. The Needleman-Wunsch global alignment tool on NCBI works well, but it only allows aligning a single mutant with the wild type at a time.

Do you know of a way to perform this global alignment on multiple proteins from the mutant list in one go?

Note that standard local alignment tools don’t work properly, as they tend to mismatch the insertion or deletion mutations (stretches of 5 to 10 amino acids) with the wild-type sequence.

Thanks

Global-alignment Needleman-Wunsch • 565 views
ADD COMMENT
0
Entering edit mode

If you're comfortable using Python, you can try this module: https://github.com/aziele/pairwise-sequence-alignment.

ADD REPLY
0
Entering edit mode

Thanks, you guys helped a lot.

ADD REPLY
0
Entering edit mode
29 days ago

If you don't want to program, you can use exonerate or needle from the EMBOSS suite.

ADD COMMENT
0
Entering edit mode
15 days ago

You can use the bio with the align tool. Install with pip install bio then

bio align query.fa target.fa

You can put all the targets in the target.fa, and it will generate alignments of the query to each target.

The tools support local, global, and semi-global alignments; you can pass scoring matrices. Works seamlessly for nucleotide and protein alignments.

It does have limitations since it relies on an optimal alignment algorithm from BioPython—realistically, it tops out at around 30K bases. I would not use it for sequences longer than that.

Docs: https://www.bioinfo.help/bio-align.html

ADD COMMENT

Login before adding your answer.

Traffic: 1208 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6