Difference between revisions of "BIN-SEQA-Composition"

From "A B C"
Jump to navigation Jump to search
m
m
 
(13 intermediate revisions by the same user not shown)
Line 1: Line 1:
−
<div id="BIO">
+
<div id="ABC">
−
  <div class="b1">
+
<div style="padding:5px; border:1px solid #000000; background-color:#b3dbce; font-size:300%; font-weight:400; color: #000000; width:100%;">
−
Sequence Analysis: Composition
+
Analysis of Sequence Composition
−
  </div>
+
<div style="padding:5px; margin-top:20px; margin-bottom:10px; background-color:#b3dbce; font-size:30%; font-weight:200; color: #000000; ">
−
 
+
(Analysis of sequence composition.)
−
  {{Vspace}}
+
</div>
−
 
 
−
<div class="keywords">
 
−
<b>Keywords:</b>&nbsp;
 
−
Sequence analysis by composition (with R examples)
 
 
</div>
 
</div>
  
−
{{Vspace}}
+
{{Smallvspace}}
−
 
 
−
 
 
−
__TOC__
 
−
 
 
−
{{Vspace}}
 
−
 
 
−
 
 
−
{{DEV}}
 
−
 
 
−
{{Vspace}}
 
  
  
 +
<div style="padding:5px; border:1px solid #000000; background-color:#b3dbce33; font-size:85%;">
 +
<div style="font-size:118%;">
 +
<b>Abstract:</b><br />
 +
<section begin=abstract />
 +
The most basic step of sequence analysis is to ask: what amino acids does a sequence contain.
 +
<section end=abstract />
 
</div>
 
</div>
−
<div id="ABC-unit-framework">
+
<!-- ============================  -->
−
== Abstract ==
+
<hr>
−
<!-- included from "../components/BIN-SEQA-Composition.components.wtxt", section: "abstract" -->
+
<table>
−
...
+
<tr>
−
 
+
<td style="padding:10px;">
−
{{Vspace}}
+
<b>Objectives:</b><br />
−
 
+
This unit will ...
−
 
+
* ... introduce basic compositional amino acid analysis;
−
== This unit ... ==
+
* ... teach how to prodcue barplots of enrichment ratios;
−
=== Prerequisites ===
+
</td>
−
<!-- included from "../components/BIN-SEQA-Composition.components.wtxt", section: "prerequisites" -->
+
<td style="padding:10px;">
−
<!-- included from "ABC-unit_components.wtxt", section: "notes-external_prerequisites" -->
+
<b>Outcomes:</b><br />
−
You need the following preparation before beginning this unit. If you are not familiar with this material from courses you took previously, you need to prepare yourself from other information sources:
+
After working through this unit you ...
−
<!-- included from "FND-prerequisites.wtxt", section: "biomolecules" -->
+
* ... are familar with comparing observed to reference frequencies;
 +
* ... can use tools from the seqinr package for simple amino acid analysis;
 +
* ... can produce barplots of log(ratio) data.
 +
</td>
 +
</tr>
 +
</table>
 +
<!-- ============================ -->
 +
<hr>
 +
<b>Deliverables:</b><br />
 +
<section begin=deliverables />
 +
<li><b>Time management</b>: Before you begin, estimate how long it will take you to complete this unit. Then, record in your course journal: the number of hours you estimated, the number of hours you worked on the unit, and the amount of time that passed between start and completion of this unit.</li>
 +
<li><b>Journal</b>: Document your progress in your [[FND-Journal|Course Journal]]. Some tasks may ask you to include specific items in your journal. Don't overlook these.</li>
 +
<li><b>Insights</b>: If you find something particularly noteworthy about this unit, make a note in your [[ABC-Insights|'''insights!''' page]].</li>
 +
<section end=deliverables />
 +
<!-- ============================  -->
 +
<hr>
 +
<section begin=prerequisites />
 +
<b>Prerequisites:</b><br />
 +
You need the following preparation before beginning this unit. If you are not familiar with this material from courses you took previously, you need to prepare yourself from other information sources:<br />
 
*<b>Biomolecules</b>: The molecules of life; nucleic acids and amino acids; the genetic code; protein folding; post-translational modifications and protein biochemistry; membrane proteins; biological function.
 
*<b>Biomolecules</b>: The molecules of life; nucleic acids and amino acids; the genetic code; protein folding; post-translational modifications and protein biochemistry; membrane proteins; biological function.
−
<!-- included from "ABC-unit_components.wtxt", section: "notes-prerequisites" -->
+
This unit builds on material covered in the following prerequisite units:<br />
−
You need to complete the following units before beginning this one:
+
*[[BIN-SEQA-Concepts|BIN-SEQA-Concepts (Concepts of Sequence Analysis)]]
−
*[[BIN-SEQA-Concepts]]
+
<section end=prerequisites />
 +
<!-- ============================  -->
 +
</div>
  
−
{{Vspace}}
+
{{Smallvspace}}
  
  
−
=== Objectives ===
 
−
<!-- included from "../components/BIN-SEQA-Composition.components.wtxt", section: "objectives" -->
 
−
...
 
  
−
{{Vspace}}
+
{{Smallvspace}}
  
  
−
=== Outcomes ===
+
__TOC__
−
<!-- included from "../components/BIN-SEQA-Composition.components.wtxt", section: "outcomes" -->
 
−
...
 
−
 
 
−
{{Vspace}}
 
−
 
 
−
 
 
−
=== Deliverables ===
 
−
<!-- included from "../components/BIN-SEQA-Composition.components.wtxt", section: "deliverables" -->
 
−
<!-- included from "ABC-unit_components.wtxt", section: "deliverables-time_management" -->
 
−
*<b>Time management</b>: Before you begin, estimate how long it will take you to complete this unit. Then, record in your course journal: the number of hours you estimated, the number of hours you worked on the unit, and the amount of time that passed between start and completion of this unit.
 
−
<!-- included from "ABC-unit_components.wtxt", section: "deliverables-journal" -->
 
−
*<b>Journal</b>: Document your progress in your [[FND-Journal|course journal]].
 
−
<!-- included from "ABC-unit_components.wtxt", section: "deliverables-insights" -->
 
−
*<b>Insights</b>: If you find something particularly noteworthy about this unit, make a note in your [[ABC-Insights|insights! page]].
 
  
 
{{Vspace}}
 
{{Vspace}}
Line 74: Line 69:
  
 
=== Evaluation ===
 
=== Evaluation ===
−
<!-- included from "../components/BIN-SEQA-Composition.components.wtxt", section: "evaluation" -->
 
−
<!-- included from "ABC-unit_components.wtxt", section: "eval-none" -->
 
 
<b>Evaluation: NA</b><br />
 
<b>Evaluation: NA</b><br />
−
:This unit is not evaluated for course marks.
+
<div style="margin-left: 2rem;">This unit is not evaluated for course marks.</div>
−
 
 
−
{{Vspace}}
 
−
 
 
−
 
 
−
</div>
 
−
<div id="BIO">
 
 
== Contents ==
 
== Contents ==
−
<!-- included from "../components/BIN-SEQA-Composition.components.wtxt", section: "contents" -->
 
 
{{Task|1=
 
{{Task|1=
 
*Read the introductory notes on {{ABC-PDF|BIN-SEQA-Composition|"Composition" as a paradigm for sequence analysis}}.
 
*Read the introductory notes on {{ABC-PDF|BIN-SEQA-Composition|"Composition" as a paradigm for sequence analysis}}.
Line 91: Line 77:
  
  
−
{{Vspace}}
+
===EMBOSS Tools ===
−
 
+
{{Smallvspace}}
 +
{{task|1=
  
 +
*Navigate to an instance of the "EMBOSS explorer" e.g. at [http://www.bioinformatics.nl/emboss-explorer/ '''bioinformatics.nl''']
  
 +
;Local composition
 +
# Find <code>pepinfo</code> under the '''PROTEIN COMPOSITION''' heading.
 +
# Retrieve the MYSPE Mbp1 related sequence from your '''R''' database, e.g. with something like <br /><code>  cat(myDB$protein[myDB$protein$name == "PHD1_SACCE" , "sequence"])</code>
 +
# Copy and paste the sequence into the input field.
 +
# Run with default parameters.
 +
# Scroll to the figures <ul>all the way at the bottom</ul>.
 +
# Do the same in a separate window for yeast Mbp1.
 +
# Try to compare ... <small>(kind of hard without reference, overlay and expectation, isn't it?)</small>
 
{{Vspace}}
 
{{Vspace}}
  
−
 
+
;Global composition
−
== Further reading, links and resources ==
+
# Find <code>pepstats</code> under the '''PROTEIN COMPOSITION''' heading.
−
<!-- {{#pmid: 19957275}} -->
+
# Paste the MYSPE Mbp1 sequence into the input field.
−
<!-- {{WWW|WWW_GMOD}} -->
+
# Run with default parameters.
−
<!-- <div class="reference-box">[http://www.ncbi.nlm.nih.gov]</div> -->
+
# Do the same in a separate window for yeast Mbp1.
−
 
+
# Try to compare ... are there significant and unexpected differences?
−
{{Vspace}}
+
}}
−
 
 
−
 
 
−
== Notes ==
 
−
<!-- included from "../components/BIN-SEQA-Composition.components.wtxt", section: "notes" -->
 
−
<!-- included from "ABC-unit_components.wtxt", section: "notes" -->
 
−
<references />
 
−
 
 
−
{{Vspace}}
 
−
 
 
−
 
 
−
</div>
 
−
<div id="ABC-unit-framework">
 
−
== Self-evaluation ==
 
−
<!-- included from "../components/BIN-SEQA-Composition.components.wtxt", section: "self-evaluation" -->
 
−
<!--
 
−
=== Question 1===
 
−
 
 
−
Question ...
 
−
 
 
−
<div class="toccolours mw-collapsible mw-collapsed" style="width:800px">
 
−
Answer ...
 
−
<div class="mw-collapsible-content">
 
−
Answer ...
 
−
 
 
−
</div>
 
−
  </div>
 
−
 
 
−
  {{Vspace}}
 
−
 
 
−
-->
 
  
 
{{Vspace}}
 
{{Vspace}}
  
−
 
+
{{ABC-unit|BIN-SEQA-Composition.R}}
  
 
{{Vspace}}
 
{{Vspace}}
  
−
 
−
<!-- included from "ABC-unit_components.wtxt", section: "ABC-unit_ask" -->
 
−
 
−
----
 
  
 
{{Vspace}}
 
{{Vspace}}
  
−
<b>If in doubt, ask!</b> If anything about this learning unit is not clear to you, do not proceed blindly but ask for clarification. Post your question on the course mailing list: others are likely to have similar problems. Or send an email to your instructor.
 
−
 
−
----
 
−
 
−
{{Vspace}}
 
  
 
<div class="about">
 
<div class="about">
Line 162: Line 119:
 
:2017-08-05
 
:2017-08-05
 
<b>Modified:</b><br />
 
<b>Modified:</b><br />
−
:2017-08-05
+
:2020-09-24
 
<b>Version:</b><br />
 
<b>Version:</b><br />
−
:0.1
+
:1.1
 
<b>Version history:</b><br />
 
<b>Version history:</b><br />
 +
*1.1 2020 Maintenance
 +
*1.0 First live
 
*0.1 First stub
 
*0.1 First stub
 
</div>
 
</div>
−
[[Category:ABC-units]]
 
−
<!-- included from "ABC-unit_components.wtxt", section: "ABC-unit_footer" -->
 
  
 
{{CC-BY}}
 
{{CC-BY}}
  
 +
[[Category:ABC-units]]
 +
{{UNIT}}
 +
{{LIVE}}
 
</div>
 
</div>
 
<!-- [END] -->
 
<!-- [END] -->

Latest revision as of 15:51, 24 September 2020

Analysis of Sequence Composition

(Analysis of sequence composition.)


 


Abstract:

The most basic step of sequence analysis is to ask: what amino acids does a sequence contain.

Objectives:
This unit will ...

  • ... introduce basic compositional amino acid analysis;
  • ... teach how to prodcue barplots of enrichment ratios;

Outcomes:
After working through this unit you ...

  • ... are familar with comparing observed to reference frequencies;
  • ... can use tools from the seqinr package for simple amino acid analysis;
  • ... can produce barplots of log(ratio) data.

Deliverables:

  • Time management: Before you begin, estimate how long it will take you to complete this unit. Then, record in your course journal: the number of hours you estimated, the number of hours you worked on the unit, and the amount of time that passed between start and completion of this unit.
  • Journal: Document your progress in your Course Journal. Some tasks may ask you to include specific items in your journal. Don't overlook these.
  • Insights: If you find something particularly noteworthy about this unit, make a note in your insights! page.

    • Prerequisites:
    You need the following preparation before beginning this unit. If you are not familiar with this material from courses you took previously, you need to prepare yourself from other information sources:

    • Biomolecules: The molecules of life; nucleic acids and amino acids; the genetic code; protein folding; post-translational modifications and protein biochemistry; membrane proteins; biological function.

    This unit builds on material covered in the following prerequisite units:


     



     



     


    Evaluation

    Evaluation: NA

    This unit is not evaluated for course marks.

    Contents

    Task:


    EMBOSS Tools

     

    Task:

    Local composition
    1. Find pepinfo under the PROTEIN COMPOSITION heading.
    2. Retrieve the MYSPE Mbp1 related sequence from your R database, e.g. with something like
      cat(myDB$protein[myDB$protein$name == "PHD1_SACCE" , "sequence"])
    3. Copy and paste the sequence into the input field.
    4. Run with default parameters.
    5. Scroll to the figures
        all the way at the bottom
      .
    6. Do the same in a separate window for yeast Mbp1.
    7. Try to compare ... (kind of hard without reference, overlay and expectation, isn't it?)
     
    Global composition
    1. Find pepstats under the PROTEIN COMPOSITION heading.
    2. Paste the MYSPE Mbp1 sequence into the input field.
    3. Run with default parameters.
    4. Do the same in a separate window for yeast Mbp1.
    5. Try to compare ... are there significant and unexpected differences?


     

    Task:

     
    • Open RStudio and load the ABC-units R project. If you have loaded it before, choose File → Recent projects → ABC-Units. If you have not loaded it before, follow the instructions in the RPR-Introduction unit.
    • Choose Tools → Version Control → Pull Branches to fetch the most recent version of the project from its GitHub repository with all changes and bug fixes included.
    • Type init() if requested.
    • Open the file BIN-SEQA-Composition.R and follow the instructions.


     

    Note: take care that you understand all of the code in the script. Evaluation in this course is cumulative and you may be asked to explain any part of code.


     


     


     


    About ...
     
    Author:

    Boris Steipe <boris.steipe@utoronto.ca>

    Created:

    2017-08-05

    Modified:

    2020-09-24

    Version:

    1.1

    Version history:

    • 1.1 2020 Maintenance
    • 1.0 First live
    • 0.1 First stub

    CreativeCommonsBy.png This copyrighted material is licensed under a Creative Commons Attribution 4.0 International License. Follow the link to learn more.

    Retrieved from "http://steipe.biochemistry.utoronto.ca/abc/index.php?title=BIN-SEQA-Composition&oldid=11640"