Bioinformatics Sequence Alignment
Introduction to Sequence Alignment
Finding the Similarities
Imagine you have two slightly different versions of the same sentence. To see how they differ, you'd line them up character by character. Bioinformatics does the same thing with biological sequences, which are long strings of letters representing DNA, RNA, or proteins. This process is called sequence alignment.
The goal is to arrange two or more sequences to identify regions of similarity. This similarity isn't just a curiosity; it can be a clue to a shared evolutionary past, a similar function, or a related three-dimensional structure. By comparing a newly discovered gene to a database of known genes, scientists can make educated guesses about what it does without spending months or years in a lab.
Think of it as recognizing that "colour" and "color" are variations of the same word. Sequence alignment helps us spot these kinds of relationships in the language of life.
The Language of Alignment
When we align two sequences, we're looking for how well the letters, or residues, line up. There are three basic possibilities at each position.
A match occurs when the letters at a position are identical. This suggests the position has been conserved over time. A mismatch means the letters are different. This indicates a substitution mutation has occurred at some point in the past. Here’s a simple alignment showing both:
Sequence 1: A C T G G A
Sequence 2: A C T C G A
| | | | |
Matches: A C T G A
Mismatch: C vs G
But what if one sequence is longer than the other? Evolution doesn't just substitute letters; it can also add or remove them through insertions or deletions. To account for this, we introduce gaps into one of the sequences.
gap
noun
A space introduced into a sequence alignment to compensate for an insertion or deletion in one of the sequences.
Introducing a gap in one sequence allows us to better align the remaining parts of both sequences. For example, to align ACGT and ACT, we can add a gap:
Sequence 1: A C G T
Sequence 2: A C - T
This shows a better relationship than forcing a mismatch between G and T. Together, matches, mismatches, and gaps help us construct the most plausible story of how two sequences might be related.
Keeping Score
There are many ways to align two sequences, especially long ones. How do we know which alignment is the best? We use a scoring system to grade them. A scoring system is a simple set of rules that assigns points to different events in an alignment.
Generally, the rules look something like this:
- Match: Positive score (e.g., +1)
- Mismatch: Negative score (e.g., -1)
- Gap: Negative score, often called a gap penalty (e.g., -2)
The goal is to find the alignment with the highest possible total score. This score represents the most likely evolutionary relationship between the sequences, based on our rules.
| Alignment Event | Typical Score | Rationale |
|---|---|---|
| Match | Positive | Rewards similarity, suggesting conservation. |
| Mismatch | Negative | Penalizes differences, suggesting mutation. |
| Gap | Negative | Penalizes insertions/deletions, which are often rare. |
More complex scoring systems exist. For instance, in protein sequences, some amino acid substitutions are more likely than others from a biochemical perspective. A mismatch between two chemically similar amino acids might receive a smaller penalty than a mismatch between two very different ones. These nuanced rules help produce more biologically meaningful alignments.
Now, let's test your understanding of these core concepts.
What is the primary goal of sequence alignment in bioinformatics?
In sequence alignment, what does a gap represent?
Understanding these fundamentals of alignment is the first step in using powerful bioinformatics tools to decode the stories hidden within DNA and protein sequences.