Multiple sequence alignment (MSA) in bioinformatics involves aligning three or more biological
sequences (DNA, RNA, or protein) to identify conserved regions and infer evolutionary and
functional relationships.
Applications:
Phylogenetic Analysis: MSAs are crucial for constructing phylogenetic trees, which
represent the evolutionary history of a group of organisms.
Identification of Conserved Regions: By aligning multiple sequences, conserved regions
(motifs, domains, patterns) that are functionally important can be identified.
Structure and Function Prediction: MSAs can be used to predict protein structure and
identify residues crucial for protein function.
Single nucleotide polymorphism (SNP) analysis
Extrapolation (to know the protein family of sequences)
DNA regulatory elements prediction (finding potentially similar binding sites)
Gene annotation
Steps of MSA
1. Gather sequences
2. Compute MSA
3. Evaluate the quality of the alignment: Remove insertions/deletions, keep trimming to
interpret, remove gap-rich regions
4. Interpret the MSA
5. Save the sequences for further analysis
Input Required
FASTA format files with multiple sequences.
All sequences should be homologous (same gene/protein family).
Software for computing MSAs
-ClustalW
-MUSCLE
-MEGA
-3D-Coffee
-BioEdit
-TCOFFEE
FASTA format: In bioinformatics and biochemistry, the FASTA format is a text-based format
for representing either nucleotide sequences or amino acid (protein) sequences, in which
nucleotides or amino acids are represented using single-letter codes. The format allows for
sequence names and comments to precede the sequences. It originated from the FASTA software
package and has since become a near-universal standard in bioinformatics. T he simplicity of
FASTA format makes it easy to manipulate and parse sequences using text-processing tools and
scripting languages. A sequence begins with a greater-than character (">") followed by a
description of the sequence (all in a single line). The lines immediately following the description
line are the sequence representation, with one letter per amino acid or nucleic acid, and are
typically no more than 80 characters in length
Example:
>MCHU - Calmodulin - Human, rabbit, bovine, rat, and chicken
MADQLTEEQIAEFKEAFSLFDKDGDGTITTKELGTVMRSLGQNPTEAELQDMINEVDAD
GNGTIDFPEFLTMMARKMKDTDSEEEIREAFRVFDKDGNGYISAAELRHVMTNLGEKLT
DEEVDEMIREADIDGDGQVNYEEFVQMMTAK