~ubuntu-branches/ubuntu/raring/libbpp-seq/raring

« back to all changes in this revision

Viewing changes to src/Bpp/Seq/StateProperties/AAChenGuHuangHydrophobicityIndex.h

  • Committer: Bazaar Package Importer
  • Author(s): Julien Dutheil
  • Date: 2011-06-09 11:00:00 UTC
  • Revision ID: james.westby@ubuntu.com-20110609110000-hnfrd9it0np58l54
Tags: upstream-2.0.2
Import upstream version 2.0.2

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
//
 
2
// File: AAChenGuHuangHydrophobicityIndex.h
 
3
// Created by: Bastien Boussau
 
4
// Created on: Fri Jan 14 10:31 2011
 
5
//
 
6
 
 
7
/*
 
8
 Copyright or © or Copr. CNRS, (November 17, 2004)
 
9
 
 
10
 This software is a computer program whose purpose is to provide classes
 
11
 for sequences analysis.
 
12
 
 
13
 This software is governed by the CeCILL  license under French law and
 
14
 abiding by the rules of distribution of free software.  You can  use, 
 
15
 modify and/ or redistribute the software under the terms of the CeCILL
 
16
 license as circulated by CEA, CNRS and INRIA at the following URL
 
17
 "http://www.cecill.info". 
 
18
 
 
19
 As a counterpart to the access to the source code and  rights to copy,
 
20
 modify and redistribute granted by the license, users are provided only
 
21
 with a limited warranty  and the software's author,  the holder of the
 
22
 economic rights,  and the successive licensors  have only  limited
 
23
 liability. 
 
24
 
 
25
 In this respect, the user's attention is drawn to the risks associated
 
26
 with loading,  using,  modifying and/or developing or reproducing the
 
27
 software by the user in light of its specific status of free software,
 
28
 that may mean  that it is complicated to manipulate,  and  that  also
 
29
 therefore means  that it is reserved for developers  and  experienced
 
30
 professionals having in-depth computer knowledge. Users are therefore
 
31
 encouraged to load and test the software's suitability as regards their
 
32
 requirements in conditions enabling the security of their systems and/or 
 
33
 data to be ensured and,  more generally, to use and operate it in the 
 
34
 same conditions as regards security. 
 
35
 
 
36
 The fact that you are presently reading this means that you have had
 
37
 knowledge of the CeCILL license and that you accept its terms.
 
38
 */
 
39
 
 
40
#ifndef _AACHENGUHUANGHYDROPHOBICITYINDEX_H_
 
41
#define _AACHENGUHUANGHYDROPHOBICITYINDEX_H_
 
42
 
 
43
#include "AlphabetIndex1.h"
 
44
#include "../Alphabet/ProteicAlphabet.h"
 
45
 
 
46
namespace bpp
 
47
{
 
48
  
 
49
  /**
 
50
   * @brief Hydrophobicity of each amino acid, according to Table 1 in Chen, Gu and Huang, BMC Bioinformatics 2006.
 
51
   *
 
52
   * @code
 
53
   * Database: AAChenGuHuangHydrophobicity
 
54
   * Entry: CHENH06
 
55
   *
 
56
   * H FAUJ880111
 
57
   * D Hydrophobicity (Chen et al., 2006)
 
58
   * R PMCID:PMC1780123 PMID:17217506
 
59
   * A Hang Chen, Fei Gu, and Zhengge Huang.
 
60
   * T Improved Chou-Fasman method for protein secondary structure prediction
 
61
   * J BMC Bioinformatics. 2006; 7(Suppl 4): S14.  (2006)
 
62
   * //
 
63
   * @endcode
 
64
   * 
 
65
   */
 
66
  class AAChenGuHuangHydrophobicityIndex:
 
67
  public AlphabetIndex1<double>
 
68
  {
 
69
  private:
 
70
  std::vector<double> hydrophobicity_;
 
71
  
 
72
  public:
 
73
  AAChenGuHuangHydrophobicityIndex() :
 
74
  hydrophobicity_()
 
75
  {
 
76
  hydrophobicity_.resize(20);
 
77
  hydrophobicity_[ 0] =  0.87; //A
 
78
  hydrophobicity_[ 1] =  0.85; //R
 
79
  hydrophobicity_[ 2] =  0.09; //N
 
80
  hydrophobicity_[ 3] =  0.66; //D
 
81
  hydrophobicity_[ 4] =  1.52; //C
 
82
  hydrophobicity_[ 5] =  0.00; //Q
 
83
  hydrophobicity_[ 6] =  0.67; //E
 
84
  hydrophobicity_[ 7] =  0.00; //G
 
85
  hydrophobicity_[ 8] =  0.87; //H
 
86
  hydrophobicity_[ 9] =  3.15; //I
 
87
  hydrophobicity_[10] =  2.17; //L
 
88
  hydrophobicity_[11] =  1.64; //K
 
89
  hydrophobicity_[12] =  1.67; //M
 
90
  hydrophobicity_[13] =  2.87; //F
 
91
  hydrophobicity_[14] =  2.77; //P
 
92
  hydrophobicity_[15] =  0.07; //S
 
93
  hydrophobicity_[16] =  0.07; //T
 
94
  hydrophobicity_[17] =  3.77; //W
 
95
  hydrophobicity_[18] =  2.76; //Y
 
96
  hydrophobicity_[19] =  1.87; //V
 
97
  }
 
98
  
 
99
  virtual ~AAChenGuHuangHydrophobicityIndex() {}
 
100
  
 
101
  AAChenGuHuangHydrophobicityIndex* clone() const { return new AAChenGuHuangHydrophobicityIndex(); }
 
102
  
 
103
  public:
 
104
  double getIndex(int state) const throw (BadIntException)
 
105
  {
 
106
  if (state < 0 || state > 19) throw BadIntException(state, "AAChenGuHuangHydrophobicityIndex::getIndex(). Invalid state.", &AlphabetTools::PROTEIN_ALPHABET);
 
107
  return hydrophobicity_[state];
 
108
  }
 
109
  
 
110
  double getIndex(const std::string& state) const throw (BadCharException)
 
111
  {
 
112
  return hydrophobicity_[AlphabetTools::PROTEIN_ALPHABET.charToInt(state)];
 
113
  }
 
114
  
 
115
  std::vector<double>* getIndexVector() const { return new std::vector<double>(hydrophobicity_); }
 
116
  
 
117
  const Alphabet* getAlphabet() const { return &AlphabetTools::PROTEIN_ALPHABET; }
 
118
  
 
119
  };
 
120
  
 
121
} //end of namespace bpp.
 
122
 
 
123
#endif //_AACHENGUHUANGHYDROPHOBICITYINDEX_H_
 
124