CS68 Lab 0: Databases and Central Dogma

Due by 11:59 p.m., Saturday January 28, 2017
Overview

The goal of this week's lab is to reinforce basic concepts of biology and bioinformatics. Specifically, Part 1 will have you explore existing bioinformatic databases, discovering what information is available for a particular gene of interest and answering questions along the way. Part 2 will have you implement a short program that will simulate many aspects of the central dogma and explore the mechanisms more in detail.

This lab will be done in pairs, which I have pre-assigned for this lab. In future labs, you will have freedom to choose a partner. Your starting point files (and submissions) will be handled using the GitHub Enterprise interface (see below). This is similar to what most of you have seen in the CS31, CS35, and other upper-level courses. Your partner for this week will be listed in the name of your Lab 0 repo (see below). You may discuss concepts with a fellow classmate, especially if you are having difficulty with the details of transcription or translation. You may not share code, however, with students that are not your lab partner.

The submission is broken into two parts:

The programming for this week's lab may seem basic upon first read, but that is partially because we haven't covered any algorithms in class yet! It is designed to you get you back in the practice of using Python and to get you to see how transcription and translation work.

Getting Started

To get started, first create a cs68 directory in your home directory, and add a labs subdirectory to it:

mkdir cs68
cd cs68
mkdir labs
cd labs
pwd
We will be using git repos hosted on the college's GitHub server for labs in this class. If you have not used git or the college's GitHub server before, here are some instructions: Using Git page (follow the instructions for repos on Swarthmore's GitHub Enterprise server).

Next find your git repo for this lab assignment off the GitHub server for our class: CS68-S17

Clone your git repo with the lab 0 starting point files into your labs directory:

cd ~/cs68/labs
git clone [the ssh url to your your repo)
Then cd into your Lab0-id1_id2 subdirectory (id1 being replaced by the user id of you and your partner).
Part 1: Working with Databases

One of the most well studied proteins in molecular biology is the green fluorescent protein (GFP). It's discovery was recently awarded a Nobel Prize in Chemistry in 2008 for redefining how fluorescent microscopy is utilized in biology. It's also being used to create a breed of glow in the dark pets that may give you nightmares.

In this portion of the lab, you will learn about GFP using three well known databases for genomics: GenBank (for nucleotide sequences), UniProt (for protein sequences), and the Protein Data Bank (PDB; for protein structures). Along the way, you will answer questions that you will submit in your README file.

Genbank:

  1. Genbank is a database of nucleotide sequences. It can be accessed at the NCBI website (National Center for Biotechnology Information) at http://www.ncbi.nlm.nih.gov/genbank.   In the search pull down menu at the top, make sure Nucleotide is selected.   In the text box at the top of the screen where it solicits input for searching, type "GFP" and hit the Go button.
  2. This search will bring up over 1000 results.  To narrow the search, click on Advanced just below the search box. Select Gene Name from the Field Menu and type gfp in the text box. For field two, select Organism and type Aequorea victoria (this is the species name). Click Search and then change Sort by to be Date Release.
  3. These last two entries (8 and 9), M62653 and M62654, are from a seminal 1992 paper.  Click on M62653.1 to look over the Genbank record; answer the questions below in README
Before continuing, answer these questions in your README file:
  1. How long is the nucleotide sequence? Amino acid sequence? (HINT: the field CDS details the encoding, or protein sequence. Either do the math or click on the protein_id for the number of amino acids.)
  2. How many variations of the gene are found in the species population (HINT: read the abstract of the original paper. You should find a link about 12 lines down labeled PUBMED)
  3. GFP was isolated from the genome of Aequorea victoria. What is a more common name for this creature?
  4. How many bases at the beginning of the sequence are not involved in encoding the protein? At the end of the sequence? (HINT: the FEATURES table describes the segments of the gene, including the CDS tag which states which portions are encoded. The Graphics view at the top is also helpful.)

UniProt:

  1. UniProt is a database of amino acid sequences that can be accessed at UniProt.  At the UniProt homepage, type gene:GFP in the search box and click the Search button.  The first link should be GFP_AEQVI (P42212). Click on the link for the Entry ID.
  2. Examine the web page for this protein, noting the wide range of information made available. You will need this information to answer the questions below.
Before continuing, answer these questions in your README file:
  1. Ontologies are categorizations, or labelings, for functions of a protein. In particular, Gene Ontologies (GO) are a standardized set of functions that can be assigned a protein to help organize the large amounts of data about genes. What is the primary GO term on the listing for GFP? Who is the first author of the paper cited for detailing this ontology?
  2. GFP is commonly used by biomedical engineers for various purposes. If the fluorescence is not strong enough, what site on the amino acid sequence can be mutated to increase fluorescence?

Protein Data Bank:

  1. The PDB (Protein Data Bank) is a database of protein structures at http://www.rcsb.org/pdb. Type GFP into the search text box and click the Search button.
  2. Note that the GFP was once the molecule of the month! Click the Tab News & PDB-101 Articles atop the search results and read the story for GBP molecule of the month; it contains a nice history of the protein.
  3. Back to the search results. Sort by release data (oldest to newest) and click at the result 1EMA (should be first or second).
  4. Notice that a lot of the sequence and annotations in the other databases are also accessible here. This is a recent modification to the PDB, making it a great resource for known structures (genes with no known structures will not be here).
  5. If you have Java applets enabled, you can view the molecule. For now, you can at least see a static image of the molecule. This is known as a ribbon representation where instead of atoms, the shape of the protein indicates the type of secondary structure.
  6. On the far right at the top, click on Display File, then click on the link to display the structure file in PDB file format.
  7. In this file the majority of lines are ATOM lines.  Scroll down until you see those lines and note how the atoms are numbered (in this case, 1 to 1771).  In many cases, the full structure is not known so some atoms will be missing. Answer the questions for this section.
Before continuing, answer these questions in your README file:
  1. On the original page (before click on the PDF file format), there is a section titled Structure Validation. Why is this information important? Pick one of the metrics provided and describe what it means.
  2. For atom #16 (in the PDB file), what type of atom is it (HINT: look at column three. This is an abbreviation, e.g., O is oxygen, N is nitrogen, anything beginning with a C is a carbon)? What type of amino acid is it (column 4 has a three letter abbreviation; look up the full name)? What are the (X,Y,Z) coordinates for this atom (columns 7,8,9)?
  3. Note where this amino acid is in the sequence (column 6). On the original page for GFP, click the Sequence tab at the top. What type of secondary structure is the amino acid from question #2 in?

Part 2: Central Dogma

In this portion of the lab, you will create a Python library and main program to simulate operations described in the central dogma in order to better understand the link between a DNA sequence and resulting protein sequence(s).

First, you will construct 3 class definitions, one each for DNA, RNA, and Protein. I will describe the main functionalities that are expected, you can feel free to add additional information/methods. All three should be defined in a file sequences.py.

Sequence classes

First, define a DNA class. Your class should have, at a minimum, the following functionality:

Next, define an RNA class. Your class should have the following methods:

Lastly, you should create a Protein class. This class will look exactly the same as the RNA (e.g., have an amino acid strand, start, and stop) class minus the translate method. The constructor will take in an amino acid sequence and a start and stop index for finding the original encoding region in the DNA sequence. You do not need to print out directions for a protein sequence (i.e., there is no 5' to 3' designation).

Main program

You will define your main program in dogma.py. At a high level, your program should: The main loop can be as creative as you like. At a minimum, you should define behavior for the following options:
Hints and Tips

Reading FASTA file

FASTA is a standardized format used across the field to represent DNA and/or protein sequences. You can read in detail about the format at the NCBI manual page. For this lab, you only need to know that there are two types of lines in the file: description lines and sequence lines. For example:
>gi|129295|sp|P01013|OVAX_CHICK GENE X PROTEIN (OVALBUMIN-RELATED)
QIKDLLVSSSTDLDTTLVLVNAIYFKGMWKTAFNAEDTREMPFHVTKQESKPVQMMCMNNSFNVATLPAE
KMKILELPFASGDLSMLVLLPDEVSDLERIEKTINFEKLTEWTNPNTMEKRRVKVYLPQMKIEEKYNLTS
VLMALGMTDLFIPSANLTGISSAESLKISQAVHGAFMELSEDGIEMAGSTGVIEDIKHSPESEQFRADHP
FLFLIKHNPTNTIVYFGRYWSP
The first line describes the gene and can be ignored for this lab. The next four lines are the gene's protein sequence. When loading your file, you can ignore description lines. The first character on a description line will be the greater than symbol ">". Each line below the description line is part of the sequence, with 80 characters per line. Simply finish reading the file line-by-line concatenating the lines together to create one large string for the sequence.

Reading Codon Table

A codon table maps three-letter RNA codons to a single-letter amino acid that it produces. Look at the codon.txt file and note that each line contains the amino-acid abbreviation first, and then a list of all codons that map to that amino acid. You should load this file into a dictionary data structure (go here to read up on using the built-in dictionary class in Python). You should map codons to their amino acid equivalent. E.g., codonTable["AUG"] = 'M'

Program Requirements

In addition to the requirements listed above, you should ensure your code satisfies these general guidelines

Sample Runs

In your labs directory, I have placed two sample sequence files, test.fasta and gfp.fasta. The latter is the sequence for the green flourescent protein, while the former is a toy example for which I have results below. Try your code on the test file first, and then see what happens with your GFP gene (can you recover the protein sequence you find in Part 1?). If you want to try a large example, try running your code on the E. coli UTI89 genome in ecoli_uti89.fasta. It is located at /home/soni/public/cs68/ecoli_uti89.fasta. DO NOT COPY this file, it is quite large. Note that your program will take awhile to run for certain operations since it is a large sequence.
Welcome to the gene translator

Enter FASTA file name: test.fasta
Enter Codon Table file name: codon.txt

DNA sequence of length 126 successfully loaded:
5' TTAATAGCGTGGAAT...CATTTTATTTTAAAA 3'

Options:
0) Exit
1) Print raw DNA sequence
2) Invert DNA sequence
3) Transcribe DNA sequence and print summary
4) Print raw RNA sequences
5) Translate RNA sequences and print summary
6) Print raw sequences to file

Enter choice: 1

Entire DNA sequence:
TTAATAGCGTGGAATGATCCTTATTAAAGAGTGTCACGAAGAGTCGGAATAGAATATGGAGGCGACAGTCGAGGGTGGGATAGAGTCCTAAAGATAACATTAAGTGTTAATCATTTTATTTTAAAA

Options:
0) Exit
1) Print raw DNA sequence
2) Invert DNA sequence
3) Transcribe DNA sequence and print summary
4) Print raw RNA sequences
5) Translate RNA sequences and print summary
6) Print raw sequences to file

Enter choice: 3

2 Resulting mRNA sequences:
16-48: 5' AUCCUUAUUAAAGAG...CACGAAGAGUCGGAA 3'
58-87: 5' GAGGCGACAGUCGAGGGUGGGAUAGAGUCC 3'

Options:
0) Exit
1) Print raw DNA sequence
2) Invert DNA sequence
3) Transcribe DNA sequence and print summary
4) Print raw RNA sequences
5) Translate RNA sequences and print summary
6) Print raw sequences to file

Enter choice: 4

mRNA Sequence 0
AUCCUUAUUAAAGAGUGUCACGAAGAGUCGGAA
mRNA Sequence 1
GAGGCGACAGUCGAGGGUGGGAUAGAGUCC

Options:
0) Exit
1) Print raw DNA sequence
2) Invert DNA sequence
3) Transcribe DNA sequence and print summary
4) Print raw RNA sequences
5) Translate RNA sequences and print summary
6) Print raw sequences to file

Enter choice: 5

2 Resulting protein sequences:
16-48: ILIKECHEESE
58-87: EATVEGGIES

Options:
0) Exit
1) Print raw DNA sequence
2) Invert DNA sequence
3) Transcribe DNA sequence and print summary
4) Print raw RNA sequences
5) Translate RNA sequences and print summary
6) Print raw sequences to file

Enter choice: 6
Enter output filename: test.pro

File output complete

Options:
0) Exit
1) Print raw DNA sequence
2) Invert DNA sequence
3) Transcribe DNA sequence and print summary
4) Print raw RNA sequences
5) Translate RNA sequences and print summary
6) Print raw sequences to file

Enter choice: 2

DNA sequence successfully inverted:
5' TTTTAAAATAAAATG...ATTCCACGCTATTAA 3'

Options:
0) Exit
1) Print raw DNA sequence
2) Invert DNA sequence
3) Transcribe DNA sequence and print summary
4) Print raw RNA sequences
5) Translate RNA sequences and print summary
6) Print raw sequences to file

Enter choice: 3

1 Resulting mRNA sequences:
15-23: 5' AUUAACACU 3'

Options:
0) Exit
1) Print raw DNA sequence
2) Invert DNA sequence
3) Transcribe DNA sequence and print summary
4) Print raw RNA sequences
5) Translate RNA sequences and print summary
6) Print raw sequences to file

Enter choice: 5
1 Resulting protein sequences:
15-23: INT

Options:
0) Exit
1) Print raw DNA sequence
2) Invert DNA sequence
3) Transcribe DNA sequence and print summary
4) Print raw RNA sequences
5) Translate RNA sequences and print summary
6) Print raw sequences to file

Enter choice: 0

The output protein files for the other test cases are available as well:
Submitting your work

Be sure to commit your work often to prevent lost data. Only your final pushed solution will be graded.

About the Data
Thanks to Mark Goadrich for sharing his test example sequence for part 2.