Biol425 2011

Computational Molecular Biology

BIOL 425/790.49, Spring 2011

Hunter College of the City University of New York

Course information

Instructors: Che Martin and Yozen Hernandez

Class Hours: Room 1000G HN; Saturday 11am-2pm

Office Hours: Room 830 HN; Thursdays 12-2pm or by appointment

Contact information:

• Che: cmartin@gc.cuny.edu, 917-684-0864
• Yozen: yzhernand@gmail.com, 347-829-6936

Course Description

Background

Biomedical research is becoming a high-throughput science. As a result, information technology plays an increasingly important role in biomedical discovery. Bioinformatics is a new interdisciplinary field formed between molecular biology and computer science.

Contents

This course will introduce both bioinformatics theories and practices. Topics include: database searching, sequence alignment, molecular phylogenetics, structure prediction, and microarray analysis. The course is held in a UNIX-based instructional lab specifically configured for bioinformatics applications. Each session consists of a first-half instruction on bioinformatics theories and a second-half session of hands-on exercises.

Learning Goals

Students are expected to be able to:

• Approach biological questions evolutionarily ("Tree-thinking")
• Evaluate and interpret computational results statistically ("Statistical-thinking")
• Formulate informatics questions quantitatively and precisely ("Abstraction")
• Design efficient procedures to solve problems ("Algorithm-thinking")
• Manipulate high-volume textual data using UNIX tools, Perl/BioPerl, R, and Relational Database ("Data Visualization")

Pre-requisites

This 3-credit course is designed for upper-level undergraduates and graduate students. Prior experiences in the UNIX Operating System and at least one programming language are required. Hunter pre-requisites are CSCI132 (Practical Unix and Perl Programming) and BIOL300 (Biochemistry) or BIOL302 (Molecular Genetics), or permission by the instructor.

Textbook

Krane & Raymer (2003). Fundamental Concepts of Bioinformatics. Pearson Education, Inc. (ISBN 0-8053-4633-3)

This book should be available in the Hunter Bookstore, as well as through several popular retailers and resellers online.

Grading & Academic Honesty

Hunter College regards acts of academic dishonesty (e.g., plagiarism, cheating on examinations, obtaining unfair advantage, and falsification of records and official documents) as serious offenses against the values of intellectual honesty. The College is committed to enforcing the CUNY Policy on Academic Integrity and will pursue cases of academic dishonesty according to the Hunter College Academic Integrity Procedures.

Student performance will be evaluated by weekly assignments and projects. While these are take-home projects and students are allowed to work in groups and answers to some of the questions are provided in the back of the textbook, students are expected to compose the final short answers, computer commands, and code independently. There are virtually an unlimited number of ways to solve a computational problem, as are ways and personal styles to implement an algorithm. Writings and blocks of codes that are virtually exact copies between individual students will be investigated as possible cases of plagiarism (e.g., copies from the Internet, text book, or each other). In such a case, the instructor will hold closed-door exams for involved individuals. Zero credits will be given to ALL involved individuals if the instructor considers there is enough evidence for plagiarism. To avoid being investigated for plagiarism, Do Not Copy from Others & Do Not Let Others Copy Your Work.

Submit assignments in Printed Hard Copies. Email attachments will NOT be accepted. Each assignment will be graded based on timeliness (10%), completeness (30%), whether executable or having major errors (20%), correctness of the final output (20%), algorithm efficiency (10%), and cleanness and readability in programming styles (10%).

The grading scheme for the course, is as follows (Subject to some change. You will be notified with sufficient time):

• Assignments (50%): 10 exercises.
• Mid-term (20%): In class Assignment + Take home to be collected on the same day.
• Final exam (20%)
• Classroom Q & A (5%): Read the chapters before lecture.
• Attendance (5%): 1-2 absences = -2.5%. More than 2 = -5%.

Course Schedule (All Saturdays)

Dates and assignments below are subject to some change

"Lecture slides" links will be available either during or before each lecture, in PDF.

January 29

• Course Overview
• Tutorial: UNIX Account, Tools, & Emacs [Lecture Slides]
• UNIX Tutorial: Please check the new #Useful Links section below
• How to connect remotely: (Windows) (Mac)
• Homework:

Assignment #1
Linux Proficiency Display the absolute path of your home directory List files in your home directory in long format & ordered by their time stamps List files in the "/data/yoda/b/student.accounts/bio425_2011/" directory from your home directory Copy of the file "/data/yoda/b/student.accounts/bio425_2011/data/GBB.seq" into your home directory Count the number of lines in the file "GBB.seq" Show the first five lines of the file "GBB.seq" & save it to a file with arbitrary name Show your last ten commands using "history"
Read Chapter 1

February 5

• Chapter 1. Central Dogma & Wet Lab Tools [Lecture Slides]
• Homework:

Assignment #2
Linux and Molecular Sequences Stay in your home directory (without changing into the "bio425_2011" directory), count the number of ORFs in the sequence file "GBB.seq" using a single "grep" command and one option Achieve the same purpose as above, by using "grep" and "wc" Again, stay in your home directory, filter the DNA sequences in the file "/data/yoda/b/student.accounts/bio425_2011/data/TenSeq.nuc" into upper cases Filter out the "a"s in the above sequence Using "sed", filter out the string "plasmid" in the sequence names from the file "/data/yoda/b/student.accounts/bio425_2011/data/GBB.1con" Using "cut", extract the third field delimited by "_" Compose regular expressions matching the Thr, Ser, and Glu codons, respectively Write a BASH loop to count the number of sequences in every *.pep file in the directory "/data/yoda/b/student.accounts/bio425_2011/data/"
Problems (pg.31-32): 1.2, 1.3, 1.5, 1.7, 1.9, 1.10, 1.11

February 12

NO CLASS

(Read Chapter 6 for next class)

February 19

Yozen will not be lecturing

• Chapter 6. Gene and Genome Structures [Lecture Slides]
• Tutorial: ORF Prediction using GLIMMER

February 26

• Appendix 1. Basic PERL [Lecture Slides]

March 5

• Object-Oriented PERL & BioPerl [Lecture Slides]

March 12

• Information Theory
• Tutorial: Sequence Logo

March 19

• Chapter 2. Data Search and Alignments
• Tutorial: Pairwise Alignment using BLAST & NUCMER
• Tutorial: Multiple Alignment using CLUSTALW

March 26

• Chapter 3. Molecular Evolution

April 2

NO CLASSES

April 9

• Chapter 4. Phylogenetics I. Distance Methods
• Tutorial: PROTDIST and NEIGHBOR using Mobyle Pasteur

April 16

• Chapter 5. Phylogenetics II. Character-Based Methods
• Tutorial: DNAML and bootstrap analysis using Mobyle Pasteur

April 23

• Relational Database and SQL
• Tutorial: the Borrelia Genome Database

April 30

• Statistics
• Tutorial: Statistical Visualization using R

May 7

• Chapter 6 (Gene Expression) & Chapter 8 (Proteomics)
• Tutorial: Array Data Visualization and Analysis

May 14

• Chapter 7. Protein Structure Prediction

May 21

• Final Project Due (TBA)

Useful Links

• A very nice UNIX tutorial (you will only need up to, and including, tutorial 4).
• FOSSWire's Unix/Linux command reference (PDF). Of use to you: "File commands", "SSH", "Searching" and "Shortcuts".

© Weigang Qiu, Hunter College, Last Update Jan 2011