How many did hydropathy predict but not tmhmm

Assignment Help Computer Engineering

Reference no: EM131905239

Practical Computing for Biologists Assignment: Integrating Tools: Python + UNIX

The last assignment asked you to build a script that calculated the hydropathy scores across a protein sequence using a "sliding window" approach. Highly positive regions for hydropathy scores are proposed to indicate regions of a protein predicted to be associated with transmembrane domains, providing both information on the function of the protein and its cellular location. However, a "sliding window" approach is not the only way to estimate whether a protein is predicted to have a transmembrane domain. One of the other more well--known methods is implemented in a program called TMHMM (Transmembrane Hidden Markov Model). Without going into the details of this method (see iLearn site if you're really interested - I'll post a paper there), this approach does a global analysis to predict, based on amino acid sequence, if there are regions of a protein that are highly likely to be associated with a secondary α--helix structure with certain characteristics; those regions are then predicted to be transmembrane domains. In this assignment, you compare the predictions of using hydropathy scores to predict transmembrane proteins with those of a program like TMHMM. To do this you will examine a subset of proteins from the ~20,000 reference proteins for humans in the SwissProt database and will do the following:

1) Modify the provided template script from the previous assignment to allow arguments to define the input file name (changes each time script is run), a threshold value for determining "transmembrane" prediction and window size to be used in the python script from UNIX command line interface. Also, change the code such that instead of outputting the scores for each bin, only output the name of the gene (it's the base name of the fasta file) if it has any score for a bin greater than an argument define threshold. (35 points)

2) Provide a UNIX script that can run your python script for the following 9 different parameter combinations on the ~1000 sequences you are assigned and generate individual output files for each parameter set with names of genes that had hydropathy scores above the assigned thresholds: window size=9, threshold=30; window size=15, threshold=30; window size=21, threshold=30; window size=9, threshold=40; window size=15, threshold=40; window size=21, threshold=40; window size=9, threshold=50; window size=15, threshold=50; window size=21, threshold=50. (35 points)

3) Compare your results with the TMHMM predictions for all human reference proteins (Hint: Think "regular expressions"). I only want responses to questions here. Not code.

a. How many of your hydropathy--score predicted transmembrane proteins are also predicted by TMHMM.
b. How many did hydropathy predict, but not TMHMM?
c. How many did TMHMM predict, but not hydropathy--scrore?
d. Assume that the TMHMM analysis is the standard of accuracy, let's calculate the type--1 (false positive) and type--2 (false negative) errors associated with using the hydropathy method. Type--1 errors would be where hydropathy values predicted transmembrane domains but TMHMM does not, while type--2 errors would be where TMHMM predicted transmembrane domains but hydropathy scores do not. Compare the proportions of false positives and negatives of all runs and describe which conditions (window size and threshold value) are best for using hydropathy values to predict transmembrane proteins. Make sure use evidence for this.

Turn in the revised python script for part 1, the UNIX script that does the tasks in part 2, and a BRIEF (a few succinct, well--written sentences here people) summary your results associated with part. Your name must be in all three of those files and all of your code is carefully commented.

Attachment:- Attachments.rar

Reference no: EM131905239

Questions Cloud

What is the current value of your call option : What is your net gain or loss on this investment if the price of WXO is $33.60 on the option expiration date?

Why is cross-border cybercrime expanding rapidly : Why are authentication and authorization important to e-commerce? Why is cross-border cybercrime expanding rapidly? Dis-cuss possible solutions.

What is the process whereby a company goes public : Are there sufficient safeguards in the law to ensure the veracity of the information companies provide in the prospectus?

Why is the office development no longer an attractive : What is the payback period? If the interest rate is 6 percent, what is the project NPV? Should the project be accepted?

How many did hydropathy predict but not tmhmm : How many of your hydropathy--score predicted transmembrane proteins are also predicted by TMHMM. How many did hydropathy predict, but not TMHMM?

How many pins must it process each year to break even : The Rubenstein Blood Bank has a variable cost of $10 per pint of blood obtained and annual fixed cost are $ 300,000.

What priniciple of value best describes this situation : Seller smith purchased a lot for $150,000 and spent an additional $400,000 on the consrtuction of a new home.

Identify the method for validating that policy is followed : Identify the method for validating that policy is followed (ex. Review process of DDR-101 system, or procedure for handling CDs place in the general trash).

Show the respective changes for balance sheet accounts : Recently the Fed has been engaging in a contraction policy by selling Treasurysecurities to the non-bank public .

User Account

All Pages