Coder Social home page Coder Social logo

map-conservation's Introduction

map-conservation

Performs mapping of residue conservation into PDB structures

Conservation Mapping Protocol

Objective

To derive residue conservation scoring from a multiple alignment of homologous proteins and map it into the structure of one of the aligned sequences

Requirements

  1. PDB file of the solved structure
  2. Protein sequences of homologs
  3. MUSCLE multiple alignment program
  4. PyMOL molecular graphics program
  5. Conservation mapping program (installation is described below)

Procedure

  1. Installation of the conservation mapping program
    1. Open Terminal (available only in Mac OS and Linux)
    2. Make sure that pip is installed by running this command in the Terminal:
      • easy_install pip
    3. Install the conservation mapping program using pip with the following the Terminal:
    4. Preparation of the structure and derived sequence
      1. Get the amino acid sequence of the solved structure in FASTA format. If you do not have it separately, you can extract the sequence from the PDB file by opening it in PyMOL and saving the molecule in FASTA format.
      2. If you have exported your sequence from PyMOL, check the sequence to make sure the written sequence matches with what is displayed in PyMOL. Correct the exported sequence, if necessary.
      3. Rename the PDB file into structure.pdb.
    5. Multiple sequence alignment
      1. Convert the sequence file of the homologs into FASTA format, if it is not in that format yet.
      2. Append the sequence of the solved structure into the sequence file of the homologs.
      3. Do a multiple alignment of the combined sequence file using MUSCLE. Manually check the alignment and remove odd sequences, if necessary. You may use Jalview for this (www.jalview.org).
      4. Rename the output alignment file into input_alignment.fasta.
    6. Residue conservation scoring and mapping
      1. Put the structure.pdb and input_alignment.fasta file into a single folder / directory.
      2. Open Terminal, change directory to the working folder with the following command:
        • cd [path_to_the_directory]
      3. Execute the conservation mapping script by issuing the following command:
        • map_conservation.py
      4. Drag the resulting conservation_mapping.py file into PyMOL. This will show the color mapping of the residue conservation. By default, the color spectrum is yellow-white-blue in order of increasing conservation. Hence, the least conserved residues are in yellow and the more conserved ones are in blue.

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    ๐Ÿ–– Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. ๐Ÿ“Š๐Ÿ“ˆ๐ŸŽ‰

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google โค๏ธ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.