2 * Jalview - A Sequence Alignment Editor and Viewer (Development Version 2.4.1)
3 * Copyright (C) 2009 AM Waterhouse, J Procter, G Barton, M Clamp, S Searle
5 * This program is free software; you can redistribute it and/or
6 * modify it under the terms of the GNU General Public License
7 * as published by the Free Software Foundation; either version 2
8 * of the License, or (at your option) any later version.
10 * This program is distributed in the hope that it will be useful,
11 * but WITHOUT ANY WARRANTY; without even the implied warranty of
12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13 * GNU General Public License for more details.
15 * You should have received a copy of the GNU General Public License
16 * along with this program; if not, write to the Free Software
17 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA
19 package jalview.schemes;
25 public class ResidueProperties
27 public static Hashtable scoreMatrices = new Hashtable();
29 // Stores residue codes/names and colours and other things
30 public static final int[] aaIndex; // aaHash version 2.1.1 and below
32 public static final int[] nucleotideIndex;
34 public static final Hashtable aa3Hash = new Hashtable();
36 public static final Hashtable aa2Triplet = new Hashtable();
38 public static final Hashtable nucleotideName = new Hashtable();
42 aaIndex = new int[255];
43 for (int i = 0; i < 255; i++)
95 aaIndex['u'] = 22; // TODO: selenocystine triplet and codons needed. also
96 // extend subt. matrices
101 nucleotideIndex = new int[255];
102 for (int i = 0; i < 255; i++)
104 nucleotideIndex[i] = 10; // non-nucleotide symbols are all non-gap gaps.
107 nucleotideIndex['A'] = 0;
108 nucleotideIndex['a'] = 0;
109 nucleotideIndex['C'] = 1;
110 nucleotideIndex['c'] = 1;
111 nucleotideIndex['G'] = 2;
112 nucleotideIndex['g'] = 2;
113 nucleotideIndex['T'] = 3;
114 nucleotideIndex['t'] = 3;
115 nucleotideIndex['U'] = 4;
116 nucleotideIndex['u'] = 4;
117 nucleotideIndex['I'] = 5;
118 nucleotideIndex['i'] = 5;
119 nucleotideIndex['X'] = 6;
120 nucleotideIndex['x'] = 6;
121 nucleotideIndex['R'] = 7;
122 nucleotideIndex['r'] = 7;
123 nucleotideIndex['Y'] = 8;
124 nucleotideIndex['y'] = 8;
125 nucleotideIndex['N'] = 9;
126 nucleotideIndex['n'] = 9;
128 nucleotideName.put("A", "Adenine");
129 nucleotideName.put("a", "Adenine");
130 nucleotideName.put("G", "Guanine");
131 nucleotideName.put("g", "Guanine");
132 nucleotideName.put("C", "Cytosine");
133 nucleotideName.put("c", "Cytosine");
134 nucleotideName.put("T", "Thymine");
135 nucleotideName.put("t", "Thymine");
136 nucleotideName.put("U", "Uracil");
137 nucleotideName.put("u", "Uracil");
138 nucleotideName.put("I", "Inosine");
139 nucleotideName.put("i", "Inosine");
140 nucleotideName.put("X", "Xanthine");
141 nucleotideName.put("x", "Xanthine");
142 nucleotideName.put("R", "Unknown Purine");
143 nucleotideName.put("r", "Unknown Purine");
144 nucleotideName.put("Y", "Unknown Pyrimidine");
145 nucleotideName.put("y", "Unknown Pyrimidine");
146 nucleotideName.put("N", "Unknown");
147 nucleotideName.put("n", "Unknown");
152 aa3Hash.put("ALA", new Integer(0));
153 aa3Hash.put("ARG", new Integer(1));
154 aa3Hash.put("ASN", new Integer(2));
155 aa3Hash.put("ASP", new Integer(3)); // D
156 aa3Hash.put("CYS", new Integer(4));
157 aa3Hash.put("GLN", new Integer(5)); // Q
158 aa3Hash.put("GLU", new Integer(6)); // E
159 aa3Hash.put("GLY", new Integer(7));
160 aa3Hash.put("HIS", new Integer(8));
161 aa3Hash.put("ILE", new Integer(9));
162 aa3Hash.put("LEU", new Integer(10));
163 aa3Hash.put("LYS", new Integer(11));
164 aa3Hash.put("MET", new Integer(12));
165 aa3Hash.put("PHE", new Integer(13));
166 aa3Hash.put("PRO", new Integer(14));
167 aa3Hash.put("SER", new Integer(15));
168 aa3Hash.put("THR", new Integer(16));
169 aa3Hash.put("TRP", new Integer(17));
170 aa3Hash.put("TYR", new Integer(18));
171 aa3Hash.put("VAL", new Integer(19));
172 // IUB Nomenclature for ambiguous peptides
173 aa3Hash.put("ASX", new Integer(20)); // "B";
174 aa3Hash.put("GLX", new Integer(21)); // X
175 aa3Hash.put("XAA", new Integer(22)); // X unknown
176 aa3Hash.put("-", new Integer(23));
177 aa3Hash.put("*", new Integer(23));
178 aa3Hash.put(".", new Integer(23));
179 aa3Hash.put(" ", new Integer(23));
180 aa3Hash.put("Gap", new Integer(23));
185 aa2Triplet.put("A", "ALA");
186 aa2Triplet.put("a", "ALA");
187 aa2Triplet.put("R", "ARG");
188 aa2Triplet.put("r", "ARG");
189 aa2Triplet.put("N", "ASN");
190 aa2Triplet.put("n", "ASN");
191 aa2Triplet.put("D", "ASP");
192 aa2Triplet.put("d", "ASP");
193 aa2Triplet.put("C", "CYS");
194 aa2Triplet.put("c", "CYS");
195 aa2Triplet.put("Q", "GLN");
196 aa2Triplet.put("q", "GLN");
197 aa2Triplet.put("E", "GLU");
198 aa2Triplet.put("e", "GLU");
199 aa2Triplet.put("G", "GLY");
200 aa2Triplet.put("g", "GLY");
201 aa2Triplet.put("H", "HIS");
202 aa2Triplet.put("h", "HIS");
203 aa2Triplet.put("I", "ILE");
204 aa2Triplet.put("i", "ILE");
205 aa2Triplet.put("L", "LEU");
206 aa2Triplet.put("l", "LEU");
207 aa2Triplet.put("K", "LYS");
208 aa2Triplet.put("k", "LYS");
209 aa2Triplet.put("M", "MET");
210 aa2Triplet.put("m", "MET");
211 aa2Triplet.put("F", "PHE");
212 aa2Triplet.put("f", "PHE");
213 aa2Triplet.put("P", "PRO");
214 aa2Triplet.put("p", "PRO");
215 aa2Triplet.put("S", "SER");
216 aa2Triplet.put("s", "SER");
217 aa2Triplet.put("T", "THR");
218 aa2Triplet.put("t", "THR");
219 aa2Triplet.put("W", "TRP");
220 aa2Triplet.put("w", "TRP");
221 aa2Triplet.put("Y", "TYR");
222 aa2Triplet.put("y", "TYR");
223 aa2Triplet.put("V", "VAL");
224 aa2Triplet.put("v", "VAL");
227 public static final String[] aa =
228 { "A", "R", "N", "D", "C", "Q", "E", "G", "H", "I", "L", "K", "M", "F",
229 "P", "S", "T", "W", "Y", "V", "B", "Z", "X", "_", "*", ".", " " };
231 public static final Color midBlue = new Color(100, 100, 255);
233 public static final Vector scaleColours = new Vector();
237 scaleColours.addElement(new Color(114, 0, 147));
238 scaleColours.addElement(new Color(156, 0, 98));
239 scaleColours.addElement(new Color(190, 0, 0));
240 scaleColours.addElement(Color.red);
241 scaleColours.addElement(new Color(255, 125, 0));
242 scaleColours.addElement(Color.orange);
243 scaleColours.addElement(new Color(255, 194, 85));
244 scaleColours.addElement(Color.yellow);
245 scaleColours.addElement(new Color(255, 255, 181));
246 scaleColours.addElement(Color.white);
249 public static final Color[] taylor =
250 { new Color(204, 255, 0), // A Greenish-yellowy-yellow
251 new Color(0, 0, 255), // R Blueish-bluey-blue
252 new Color(204, 0, 255), // N Blueish-reddy-blue
253 new Color(255, 0, 0), // D Reddish-reddy-red
254 new Color(255, 255, 0), // C Yellowish-yellowy-yellow
255 new Color(255, 0, 204), // Q Reddish-bluey-red
256 new Color(255, 0, 102), // E Blueish-reddy-red
257 new Color(255, 153, 0), // G Yellowy-reddy-yellow
258 new Color(0, 102, 255), // H Greenish-bluey-blue
259 new Color(102, 255, 0), // I Greenish-yellowy-green
260 new Color(51, 255, 0), // L Yellowish-greeny-green
261 new Color(102, 0, 255), // K Reddish-bluey-blue
262 new Color(0, 255, 0), // M Greenish-greeny-green
263 new Color(0, 255, 102), // F Blueish-greeny-green
264 new Color(255, 204, 0), // P Reddish-yellowy-yellow
265 new Color(255, 51, 0), // S Yellowish-reddy-red
266 new Color(255, 102, 0), // T Reddish-yellowy-red
267 new Color(0, 204, 255), // W Blueish-greeny-green
268 new Color(0, 255, 204), // Y Greenish-bluey-green
269 new Color(153, 255, 0), // V Yellowish-greeny-yellow
278 public static final Color[] nucleotide =
279 { new Color(100, 247, 63), // A
280 new Color(255, 179, 64), // C
281 new Color(235, 65, 60), // G
282 new Color(60, 136, 238), // T
283 new Color(60, 136, 238), // U
293 public static final Color[] zappo =
302 midBlue,// Color.red, // H
323 // Dunno where I got these numbers from
324 public static final double[] hyd2 =
350 public static final double[] helix =
351 { 1.42, 0.98, 0.67, 1.01, 0.70, 1.11, 1.51, 0.57, 1.00, 1.08, 1.21, 1.16,
352 1.45, 1.13, 0.57, 0.77, 0.83, 1.08, 0.69, 1.06, 0.84, 1.31, 1.00, 0.0 };
354 public static final double helixmin = 0.57;
356 public static final double helixmax = 1.51;
358 public static final double[] strand =
359 { 0.83, 0.93, 0.89, 0.54, 1.19, 1.10, 0.37, 0.75, 0.87, 1.60, 1.30, 0.74,
360 1.05, 1.38, 0.55, 0.75, 1.19, 1.37, 1.47, 1.70, 0.72, 0.74, 1.0, 0.0 };
362 public static final double strandmin = 0.37;
364 public static final double strandmax = 1.7;
366 public static final double[] turn =
367 { 0.66, 0.95, 1.56, 1.46, 1.19, 0.98, 0.74, 1.56, 0.95, 0.47, 0.59, 1.01,
368 0.60, 0.60, 1.52, 1.43, 0.96, 0.96, 1.14, 0.50, 1.51, 0.86, 1.00, 0,
371 public static final double turnmin = 0.47;
373 public static final double turnmax = 1.56;
375 public static final double[] buried =
376 { 1.7, 0.1, 0.4, 0.4, 4.6, 0.3, 0.3, 1.8, 0.8, 3.1, 2.4, 0.05, 1.9, 2.2,
377 0.6, 0.8, 0.7, 1.6, 0.5, 2.9, 0.4, 0.3, 1.358, 0.00 };
379 public static final double buriedmin = 0.05;
381 public static final double buriedmax = 4.6;
383 // This is hydropathy index
384 // Kyte, J., and Doolittle, R.F., J. Mol. Biol.
385 // 1157, 105-132, 1982
386 public static final double[] hyd =
387 { 1.8, -4.5, -3.5, -3.5, 2.5, -3.5, -3.5, -0.4, -3.2, 4.5, 3.8, -3.9,
388 1.9, 2.8, -1.6, -0.8, -0.7, -0.9, -1.3, 4.2, -3.5, -3.5, -0.49, 0.0 };
390 public static final double hydmax = 4.5;
392 public static final double hydmin = -3.9;
394 // public static final double hydmax = 1.38;
395 // public static final double hydmin = -2.53;
396 private static final int[][] BLOSUM62 =
398 { 4, -1, -2, -2, 0, -1, -1, 0, -2, -1, -1, -1, -1, -2, -1, 1, 0, -3,
399 -2, 0, -2, -1, 0, -4 },
400 { -1, 5, 0, -2, -3, 1, 0, -2, 0, -3, -2, 2, -1, -3, -2, -1, -1, -3,
401 -2, -3, -1, 0, -1, -4 },
402 { -2, 0, 6, 1, -3, 0, 0, 0, 1, -3, -3, 0, -2, -3, -2, 1, 0, -4, -2,
404 { -2, -2, 1, 6, -3, 0, 2, -1, -1, -3, -4, -1, -3, -3, -1, 0, -1, -4,
405 -3, -3, 4, 1, -1, -4 },
406 { 0, 3, -3, -3, 9, -3, -4, -3, -3, -1, -1, -3, -1, -2, -3, -1, -1,
407 -2, -2, -1, -3, -3, -2, -4 },
408 { -1, 1, 0, 0, -3, 5, 2, -2, 0, -3, -2, 1, 0, -3, -1, 0, -1, -2, -1,
410 { -1, 0, 0, 2, -4, 2, 5, -2, 0, -3, -3, 1, -2, -3, -1, 0, -1, -3, -2,
412 { 0, -2, 0, -1, -3, -2, -2, 6, -2, -4, -4, -2, -3, -3, -2, 0, -2, -2,
413 -3, -3, -1, -2, -1, -4 },
414 { -2, 0, 1, -1, -3, 0, 0, -2, 8, -3, -3, -1, -2, -1, -2, -1, -2, -2,
415 2, -3, 0, 0, -1, -4 },
416 { -1, -3, -3, -3, -1, -3, -3, -4, -3, 4, 2, -3, 1, 0, -3, -2, -1, -3,
417 -1, 3, -3, -3, -1, -4 },
418 { -1, -2, -3, -4, -1, -2, -3, -4, -3, 2, 4, -2, 2, 0, -3, -2, -1, -2,
419 -1, 1, -4, -3, -1, -4 },
420 { -1, 2, 0, -1, -3, 1, 1, -2, -1, -3, -2, 5, -1, -3, -1, 0, -1, -3,
421 -2, -2, 0, 1, -1, -4 },
422 { -1, -1, -2, -3, -1, 0, -2, -3, -2, 1, 2, -1, 5, 0, -2, -1, -1, -1,
423 -1, 1, -3, -1, -1, -4 },
424 { -2, -3, -3, -3, -2, -3, -3, -3, -1, 0, 0, -3, 0, 6, -4, -2, -2, 1,
425 3, -1, -3, -3, -1, -4 },
426 { -1, -2, -2, -1, -3, -1, -1, -2, -2, -3, -3, -1, -2, -4, 7, -1, -1,
427 -4, -3, -2, -2, -1, -2, -4 },
428 { 1, -1, 1, 0, -1, 0, 0, 0, -1, -2, -2, 0, -1, -2, -1, 4, 1, -3, -2,
430 { 0, -1, 0, -1, -1, -1, -1, -2, -2, -1, -1, -1, -1, -2, -1, 1, 5, -2,
431 -2, 0, -1, -1, 0, -4 },
432 { -3, -3, -4, -4, -2, -2, -3, -2, -2, -3, -2, -3, -1, 1, -4, -3, -2,
433 11, 2, -3, -4, -3, -2, -4 },
434 { -2, -2, -2, -3, -2, -1, -2, -3, 2, -1, -1, -2, -1, 3, -3, -2, -2,
435 2, 7, -1, -3, -2, -1, -4 },
436 { 0, -3, -3, -3, -1, -2, -2, -3, -3, 3, 1, -2, 1, -1, -2, -2, 0, -3,
437 -1, 4, -3, -2, -1, -4 },
438 { -2, -1, 3, 4, -3, 0, 1, -1, 0, -3, -4, 0, -3, -3, -2, 0, -1, -4,
439 -3, -3, 4, 1, -1, -4 },
440 { -1, 0, 0, 1, -3, 3, 4, -2, 0, -3, -3, 1, -1, -3, -1, 0, -1, -3, -2,
442 { 0, -1, -1, -1, -2, -1, -1, -1, -1, -1, -1, -1, -1, -1, -2, 0, 0,
443 -2, -1, -1, -1, -1, -1, -4 },
444 { -4, -4, -4, -4, -4, -4, -4, -4, -4, -4, -4, -4, -4, -4, -4, -4, -4,
445 -4, -4, -4, -4, -4, -4, 1 }, };
447 static final int[][] PAM250 =
449 { 2, -2, 0, 0, -2, 0, 0, 1, -1, -1, -2, -1, -1, -3, 1, 1, 1, -6, -3,
451 { -2, 6, 0, -1, -4, 1, -1, -3, 2, -2, -3, 3, 0, -4, 0, 0, -1, 2, -4,
453 { 0, 0, 2, 2, -4, 1, 1, 0, 2, -2, -3, 1, -2, -3, 0, 1, 0, -4, -2, -2,
455 { 0, -1, 2, 4, -5, 2, 3, 1, 1, -2, -4, 0, -3, -6, -1, 0, 0, -7, -4,
457 { -2, -4, -4, -5, 12, -5, -5, -3, -3, -2, -6, -5, -5, -4, -3, 0, -2,
458 -8, 0, -2, -4, -5, -3, -8 },
459 { 0, 1, 1, 2, -5, 4, 2, -1, 3, -2, -2, 1, -1, -5, 0, -1, -1, -5, -4,
461 { 0, -1, 1, 3, -5, 2, 4, 0, 1, -2, -3, 0, -2, -5, -1, 0, 0, -7, -4,
463 { 1, -3, 0, 1, -3, -1, 0, 5, -2, -3, -4, -2, -3, -5, 0, 1, 0, -7, -5,
465 { -1, 2, 2, 1, -3, 3, 1, -2, 6, -2, -2, 0, -2, -2, 0, -1, -1, -3, 0,
467 { -1, -2, -2, -2, -2, -2, -2, -3, -2, 5, 2, -2, 2, 1, -2, -1, 0, -5,
468 -1, 4, -2, -2, -1, -8 },
469 { -2, -3, -3, -4, -6, -2, -3, -4, -2, 2, 6, -3, 4, 2, -3, -3, -2, -2,
470 -1, 2, -3, -3, -1, -8 },
471 { -1, 3, 1, 0, -5, 1, 0, -2, 0, -2, -3, 5, 0, -5, -1, 0, 0, -3, -4,
473 { -1, 0, -2, -3, -5, -1, -2, -3, -2, 2, 4, 0, 6, 0, -2, -2, -1, -4,
474 -2, 2, -2, -2, -1, -8 },
475 { -3, -4, -3, -6, -4, -5, -5, -5, -2, 1, 2, -5, 0, 9, -5, -3, -3, 0,
476 7, -1, -4, -5, -2, -8 },
477 { 1, 0, 0, -1, -3, 0, -1, 0, 0, -2, -3, -1, -2, -5, 6, 1, 0, -6, -5,
479 { 1, 0, 1, 0, 0, -1, 0, 1, -1, -1, -3, 0, -2, -3, 1, 2, 1, -2, -3,
481 { 1, -1, 0, 0, -2, -1, 0, 0, -1, 0, -2, 0, -1, -3, 0, 1, 3, -5, -3,
483 { -6, 2, -4, -7, -8, -5, -7, -7, -3, -5, -2, -3, -4, 0, -6, -2, -5,
484 17, 0, -6, -5, -6, -4, -8 },
485 { -3, -4, -2, -4, 0, -4, -4, -5, 0, -1, -1, -4, -2, 7, -5, -3, -3, 0,
486 10, -2, -3, -4, -2, -8 },
487 { 0, -2, -2, -2, -2, -2, -2, -1, -2, 4, 2, -2, 2, -1, -1, -1, 0, -6,
488 -2, 4, -2, -2, -1, -8 },
489 { 0, -1, 2, 3, -4, 1, 3, 0, 1, -2, -3, 1, -2, -4, -1, 0, 0, -5, -3,
491 { 0, 0, 1, 3, -5, 3, 3, 0, 2, -2, -3, 0, -2, -5, 0, 0, -1, -6, -4,
493 { 0, -1, 0, -1, -3, -1, -1, -1, -1, -1, -1, -1, -1, -2, -1, 0, 0, -4,
494 -2, -1, -1, -1, -1, -8 },
495 { -8, -8, -8, -8, -8, -8, -8, -8, -8, -8, -8, -8, -8, -8, -8, -8, -8,
496 -8, -8, -8, -8, -8, -8, 1 }, };
498 public static final Hashtable ssHash = new Hashtable(); // stores the number
503 ssHash.put("H", Color.magenta);
504 ssHash.put("E", Color.yellow);
505 ssHash.put("-", Color.white);
506 ssHash.put(".", Color.white);
507 ssHash.put("S", Color.cyan);
508 ssHash.put("T", Color.blue);
509 ssHash.put("G", Color.pink);
510 ssHash.put("I", Color.pink);
511 ssHash.put("B", Color.yellow);
513 /* new Color(60, 136, 238), // U
521 static final int[][] DNA =
523 { 10, -8, -8, -8, 1, 0 ,0 ,0 ,0 ,0 ,1 }, // C
524 { -8, 10, -8, -8, 1, 0 ,0 ,0 ,0 ,0 ,1 }, // T
525 { -8, -8, 10, -8, 1, 0 ,0 ,0 ,0 ,0 ,1 }, // A
526 { -8, -8, -8, 10, 1, 0 ,0 ,0 ,0 ,0 ,1 }, // G
527 { 1, 1, 1, 1, 10,0 ,0 ,0 ,0 ,0 ,1 }, // -
528 { 1, 1, 1, 1, 1 ,10 ,0 ,0 ,0 ,0 ,1 }, // -
529 { 1, 1, 1, 1, 1 ,0 ,10 ,0 ,0 ,0 ,1 }, // -
530 { 1, 1, 1, 1, 1 ,0 ,0 ,10 ,0 ,0 ,1 }, // -
531 { 1, 1, 1, 1, 1 ,0 ,0 ,0 ,10 ,0 ,1 }, // -
532 { 1, 1, 1, 1, 1 ,0 ,0 ,0 ,0 ,10 ,1 }, // -
533 { 1, 1, 1, 1, 1 ,1 ,1 ,1 ,1 ,1 ,1 }, // -
536 * register matrices in list
540 scoreMatrices.put("BLOSUM62", new ScoreMatrix("BLOSUM62", BLOSUM62, 0));
541 scoreMatrices.put("PAM250", new ScoreMatrix("PAM250", PAM250, 0));
542 scoreMatrices.put("DNA", new ScoreMatrix("DNA", DNA, 1));
545 public static final Color[] pidColours =
546 { midBlue, new Color(153, 153, 255),
548 new Color(204, 204, 255), };
550 public static final float[] pidThresholds =
553 public static Hashtable codonHash = new Hashtable();
555 public static Vector Lys = new Vector();
557 public static Vector Asn = new Vector();
559 public static Vector Gln = new Vector();
561 public static Vector His = new Vector();
563 public static Vector Glu = new Vector();
565 public static Vector Asp = new Vector();
567 public static Vector Tyr = new Vector();
569 public static Vector Thr = new Vector();
571 public static Vector Pro = new Vector();
573 public static Vector Ala = new Vector();
575 public static Vector Ser = new Vector();
577 public static Vector Arg = new Vector();
579 public static Vector Gly = new Vector();
581 public static Vector Trp = new Vector();
583 public static Vector Cys = new Vector();
585 public static Vector Ile = new Vector();
587 public static Vector Met = new Vector();
589 public static Vector Leu = new Vector();
591 public static Vector Val = new Vector();
593 public static Vector Phe = new Vector();
595 public static Vector STOP = new Vector();
599 codonHash.put("K", Lys);
600 codonHash.put("N", Asn);
601 codonHash.put("Q", Gln);
602 codonHash.put("H", His);
603 codonHash.put("E", Glu);
604 codonHash.put("D", Asp);
605 codonHash.put("Y", Tyr);
606 codonHash.put("T", Thr);
607 codonHash.put("P", Pro);
608 codonHash.put("A", Ala);
609 codonHash.put("S", Ser);
610 codonHash.put("R", Arg);
611 codonHash.put("G", Gly);
612 codonHash.put("W", Trp);
613 codonHash.put("C", Cys);
614 codonHash.put("I", Ile);
615 codonHash.put("M", Met);
616 codonHash.put("L", Leu);
617 codonHash.put("V", Val);
618 codonHash.put("F", Phe);
619 codonHash.put("STOP", STOP);
622 public static Hashtable codonHash2 = new Hashtable();
626 codonHash2.put("AAA", "K");
627 codonHash2.put("AAG", "K");
628 codonHash2.put("AAC", "N");
629 codonHash2.put("AAT", "N");
631 codonHash2.put("CAA", "E");
632 codonHash2.put("CAG", "E");
633 codonHash2.put("CAC", "H");
634 codonHash2.put("CAT", "H");
636 codonHash2.put("GAA", "Q");
637 codonHash2.put("GAG", "Q");
638 codonHash2.put("GAC", "D");
639 codonHash2.put("GAT", "D");
641 codonHash2.put("TAC", "Y");
642 codonHash2.put("TAT", "Y");
644 codonHash2.put("ACA", "T");
645 codonHash2.put("AAG", "T");
646 codonHash2.put("ACC", "T");
647 codonHash2.put("ACT", "T");
649 codonHash2.put("CCA", "P");
650 codonHash2.put("CCG", "P");
651 codonHash2.put("CCC", "P");
652 codonHash2.put("CCT", "P");
654 codonHash2.put("GCA", "A");
655 codonHash2.put("GCG", "A");
656 codonHash2.put("GCC", "A");
657 codonHash2.put("GCT", "A");
659 codonHash2.put("TCA", "S");
660 codonHash2.put("TCG", "S");
661 codonHash2.put("TCC", "S");
662 codonHash2.put("TCT", "S");
663 codonHash2.put("AGC", "S");
664 codonHash2.put("AGT", "S");
666 codonHash2.put("AGA", "R");
667 codonHash2.put("AGG", "R");
668 codonHash2.put("CGA", "R");
669 codonHash2.put("CGG", "R");
670 codonHash2.put("CGC", "R");
671 codonHash2.put("CGT", "R");
673 codonHash2.put("GGA", "G");
674 codonHash2.put("GGG", "G");
675 codonHash2.put("GGC", "G");
676 codonHash2.put("GGT", "G");
678 codonHash2.put("TGA", "*");
679 codonHash2.put("TAA", "*");
680 codonHash2.put("TAG", "*");
682 codonHash2.put("TGG", "W");
684 codonHash2.put("TGC", "C");
685 codonHash2.put("TGT", "C");
687 codonHash2.put("ATA", "I");
688 codonHash2.put("ATC", "I");
689 codonHash2.put("ATT", "I");
691 codonHash2.put("ATG", "M");
693 codonHash2.put("CTA", "L");
694 codonHash2.put("CTG", "L");
695 codonHash2.put("CTC", "L");
696 codonHash2.put("CTT", "L");
697 codonHash2.put("TTA", "L");
698 codonHash2.put("TTG", "L");
700 codonHash2.put("GTA", "V");
701 codonHash2.put("GTG", "V");
702 codonHash2.put("GTC", "V");
703 codonHash2.put("GTT", "V");
705 codonHash2.put("TTC", "F");
706 codonHash2.put("TTT", "F");
711 Lys.addElement("AAA");
712 Lys.addElement("AAG");
713 Asn.addElement("AAC");
714 Asn.addElement("AAT");
716 Gln.addElement("CAA");
717 Gln.addElement("CAG");
718 His.addElement("CAC");
719 His.addElement("CAT");
721 Glu.addElement("GAA");
722 Glu.addElement("GAG");
723 Asp.addElement("GAC");
724 Asp.addElement("GAT");
726 Tyr.addElement("TAC");
727 Tyr.addElement("TAT");
729 Thr.addElement("ACA");
730 Thr.addElement("ACG");
731 Thr.addElement("ACC");
732 Thr.addElement("ACT");
734 Pro.addElement("CCA");
735 Pro.addElement("CCG");
736 Pro.addElement("CCC");
737 Pro.addElement("CCT");
739 Ala.addElement("GCA");
740 Ala.addElement("GCG");
741 Ala.addElement("GCC");
742 Ala.addElement("GCT");
744 Ser.addElement("TCA");
745 Ser.addElement("TCG");
746 Ser.addElement("TCC");
747 Ser.addElement("TCT");
748 Ser.addElement("AGC");
749 Ser.addElement("AGT");
751 Arg.addElement("AGA");
752 Arg.addElement("AGG");
753 Arg.addElement("CGA");
754 Arg.addElement("CGG");
755 Arg.addElement("CGC");
756 Arg.addElement("CGT");
758 Gly.addElement("GGA");
759 Gly.addElement("GGG");
760 Gly.addElement("GGC");
761 Gly.addElement("GGT");
763 STOP.addElement("TGA");
764 STOP.addElement("TAA");
765 STOP.addElement("TAG");
767 Trp.addElement("TGG");
769 Cys.addElement("TGC");
770 Cys.addElement("TGT");
772 Ile.addElement("ATA");
773 Ile.addElement("ATC");
774 Ile.addElement("ATT");
776 Met.addElement("ATG");
778 Leu.addElement("CTA");
779 Leu.addElement("CTG");
780 Leu.addElement("CTC");
781 Leu.addElement("CTT");
782 Leu.addElement("TTA");
783 Leu.addElement("TTG");
785 Val.addElement("GTA");
786 Val.addElement("GTG");
787 Val.addElement("GTC");
788 Val.addElement("GTT");
790 Phe.addElement("TTC");
791 Phe.addElement("TTT");
794 // Stores residue codes/names and colours and other things
795 public static Hashtable propHash = new Hashtable();
797 public static Hashtable hydrophobic = new Hashtable();
799 public static Hashtable polar = new Hashtable();
801 public static Hashtable small = new Hashtable();
803 public static Hashtable positive = new Hashtable();
805 public static Hashtable negative = new Hashtable();
807 public static Hashtable charged = new Hashtable();
809 public static Hashtable aromatic = new Hashtable();
811 public static Hashtable aliphatic = new Hashtable();
813 public static Hashtable tiny = new Hashtable();
815 public static Hashtable proline = new Hashtable();
819 hydrophobic.put("I", new Integer(1));
820 hydrophobic.put("L", new Integer(1));
821 hydrophobic.put("V", new Integer(1));
822 hydrophobic.put("C", new Integer(1));
823 hydrophobic.put("A", new Integer(1));
824 hydrophobic.put("G", new Integer(1));
825 hydrophobic.put("M", new Integer(1));
826 hydrophobic.put("F", new Integer(1));
827 hydrophobic.put("Y", new Integer(1));
828 hydrophobic.put("W", new Integer(1));
829 hydrophobic.put("H", new Integer(1));
830 hydrophobic.put("K", new Integer(1));
831 hydrophobic.put("X", new Integer(1));
832 hydrophobic.put("-", new Integer(1));
833 hydrophobic.put("*", new Integer(1));
834 hydrophobic.put("R", new Integer(0));
835 hydrophobic.put("E", new Integer(0));
836 hydrophobic.put("Q", new Integer(0));
837 hydrophobic.put("D", new Integer(0));
838 hydrophobic.put("N", new Integer(0));
839 hydrophobic.put("S", new Integer(0));
840 hydrophobic.put("T", new Integer(0));
841 hydrophobic.put("P", new Integer(0));
846 polar.put("Y", new Integer(1));
847 polar.put("W", new Integer(1));
848 polar.put("H", new Integer(1));
849 polar.put("K", new Integer(1));
850 polar.put("R", new Integer(1));
851 polar.put("E", new Integer(1));
852 polar.put("Q", new Integer(1));
853 polar.put("D", new Integer(1));
854 polar.put("N", new Integer(1));
855 polar.put("S", new Integer(1));
856 polar.put("T", new Integer(1));
857 polar.put("X", new Integer(1));
858 polar.put("-", new Integer(1));
859 polar.put("*", new Integer(1));
860 polar.put("I", new Integer(0));
861 polar.put("L", new Integer(0));
862 polar.put("V", new Integer(0));
863 polar.put("C", new Integer(0));
864 polar.put("A", new Integer(0));
865 polar.put("G", new Integer(0));
866 polar.put("M", new Integer(0));
867 polar.put("F", new Integer(0));
868 polar.put("P", new Integer(0));
873 small.put("I", new Integer(0));
874 small.put("L", new Integer(0));
875 small.put("V", new Integer(1));
876 small.put("C", new Integer(1));
877 small.put("A", new Integer(1));
878 small.put("G", new Integer(1));
879 small.put("M", new Integer(0));
880 small.put("F", new Integer(0));
881 small.put("Y", new Integer(0));
882 small.put("W", new Integer(0));
883 small.put("H", new Integer(0));
884 small.put("K", new Integer(0));
885 small.put("R", new Integer(0));
886 small.put("E", new Integer(0));
887 small.put("Q", new Integer(0));
888 small.put("D", new Integer(1));
889 small.put("N", new Integer(1));
890 small.put("S", new Integer(1));
891 small.put("T", new Integer(1));
892 small.put("P", new Integer(1));
893 small.put("-", new Integer(1));
894 small.put("*", new Integer(1));
899 positive.put("I", new Integer(0));
900 positive.put("L", new Integer(0));
901 positive.put("V", new Integer(0));
902 positive.put("C", new Integer(0));
903 positive.put("A", new Integer(0));
904 positive.put("G", new Integer(0));
905 positive.put("M", new Integer(0));
906 positive.put("F", new Integer(0));
907 positive.put("Y", new Integer(0));
908 positive.put("W", new Integer(0));
909 positive.put("H", new Integer(1));
910 positive.put("K", new Integer(1));
911 positive.put("R", new Integer(1));
912 positive.put("E", new Integer(0));
913 positive.put("Q", new Integer(0));
914 positive.put("D", new Integer(0));
915 positive.put("N", new Integer(0));
916 positive.put("S", new Integer(0));
917 positive.put("T", new Integer(0));
918 positive.put("P", new Integer(0));
919 positive.put("-", new Integer(1));
920 positive.put("*", new Integer(1));
925 negative.put("I", new Integer(0));
926 negative.put("L", new Integer(0));
927 negative.put("V", new Integer(0));
928 negative.put("C", new Integer(0));
929 negative.put("A", new Integer(0));
930 negative.put("G", new Integer(0));
931 negative.put("M", new Integer(0));
932 negative.put("F", new Integer(0));
933 negative.put("Y", new Integer(0));
934 negative.put("W", new Integer(0));
935 negative.put("H", new Integer(0));
936 negative.put("K", new Integer(0));
937 negative.put("R", new Integer(0));
938 negative.put("E", new Integer(1));
939 negative.put("Q", new Integer(0));
940 negative.put("D", new Integer(1));
941 negative.put("N", new Integer(0));
942 negative.put("S", new Integer(0));
943 negative.put("T", new Integer(0));
944 negative.put("P", new Integer(0));
945 negative.put("-", new Integer(1));
946 negative.put("*", new Integer(1));
951 charged.put("I", new Integer(0));
952 charged.put("L", new Integer(0));
953 charged.put("V", new Integer(0));
954 charged.put("C", new Integer(0));
955 charged.put("A", new Integer(0));
956 charged.put("G", new Integer(0));
957 charged.put("M", new Integer(0));
958 charged.put("F", new Integer(0));
959 charged.put("Y", new Integer(0));
960 charged.put("W", new Integer(0));
961 charged.put("H", new Integer(1));
962 charged.put("K", new Integer(1));
963 charged.put("R", new Integer(1));
964 charged.put("E", new Integer(1));
965 charged.put("Q", new Integer(0));
966 charged.put("D", new Integer(1));
967 charged.put("N", new Integer(0)); // Asparagine is polar but not charged. Alternative would be charged and negative (in basic form)?
968 charged.put("S", new Integer(0));
969 charged.put("T", new Integer(0));
970 charged.put("P", new Integer(0));
971 charged.put("-", new Integer(1));
972 charged.put("*", new Integer(1));
977 aromatic.put("I", new Integer(0));
978 aromatic.put("L", new Integer(0));
979 aromatic.put("V", new Integer(0));
980 aromatic.put("C", new Integer(0));
981 aromatic.put("A", new Integer(0));
982 aromatic.put("G", new Integer(0));
983 aromatic.put("M", new Integer(0));
984 aromatic.put("F", new Integer(1));
985 aromatic.put("Y", new Integer(1));
986 aromatic.put("W", new Integer(1));
987 aromatic.put("H", new Integer(1));
988 aromatic.put("K", new Integer(0));
989 aromatic.put("R", new Integer(0));
990 aromatic.put("E", new Integer(0));
991 aromatic.put("Q", new Integer(0));
992 aromatic.put("D", new Integer(0));
993 aromatic.put("N", new Integer(0));
994 aromatic.put("S", new Integer(0));
995 aromatic.put("T", new Integer(0));
996 aromatic.put("P", new Integer(0));
997 aromatic.put("-", new Integer(1));
998 aromatic.put("*", new Integer(1));
1003 aliphatic.put("I", new Integer(1));
1004 aliphatic.put("L", new Integer(1));
1005 aliphatic.put("V", new Integer(1));
1006 aliphatic.put("C", new Integer(0));
1007 aliphatic.put("A", new Integer(0));
1008 aliphatic.put("G", new Integer(0));
1009 aliphatic.put("M", new Integer(0));
1010 aliphatic.put("F", new Integer(0));
1011 aliphatic.put("Y", new Integer(0));
1012 aliphatic.put("W", new Integer(0));
1013 aliphatic.put("H", new Integer(0));
1014 aliphatic.put("K", new Integer(0));
1015 aliphatic.put("R", new Integer(0));
1016 aliphatic.put("E", new Integer(0));
1017 aliphatic.put("Q", new Integer(0));
1018 aliphatic.put("D", new Integer(0));
1019 aliphatic.put("N", new Integer(0));
1020 aliphatic.put("S", new Integer(0));
1021 aliphatic.put("T", new Integer(0));
1022 aliphatic.put("P", new Integer(0));
1023 aliphatic.put("-", new Integer(1));
1024 aliphatic.put("*", new Integer(1));
1029 tiny.put("I", new Integer(0));
1030 tiny.put("L", new Integer(0));
1031 tiny.put("V", new Integer(0));
1032 tiny.put("C", new Integer(0));
1033 tiny.put("A", new Integer(1));
1034 tiny.put("G", new Integer(1));
1035 tiny.put("M", new Integer(0));
1036 tiny.put("F", new Integer(0));
1037 tiny.put("Y", new Integer(0));
1038 tiny.put("W", new Integer(0));
1039 tiny.put("H", new Integer(0));
1040 tiny.put("K", new Integer(0));
1041 tiny.put("R", new Integer(0));
1042 tiny.put("E", new Integer(0));
1043 tiny.put("Q", new Integer(0));
1044 tiny.put("D", new Integer(0));
1045 tiny.put("N", new Integer(0));
1046 tiny.put("S", new Integer(1));
1047 tiny.put("T", new Integer(0));
1048 tiny.put("P", new Integer(0));
1049 tiny.put("-", new Integer(1));
1050 tiny.put("*", new Integer(1));
1055 proline.put("I", new Integer(0));
1056 proline.put("L", new Integer(0));
1057 proline.put("V", new Integer(0));
1058 proline.put("C", new Integer(0));
1059 proline.put("A", new Integer(0));
1060 proline.put("G", new Integer(0));
1061 proline.put("M", new Integer(0));
1062 proline.put("F", new Integer(0));
1063 proline.put("Y", new Integer(0));
1064 proline.put("W", new Integer(0));
1065 proline.put("H", new Integer(0));
1066 proline.put("K", new Integer(0));
1067 proline.put("R", new Integer(0));
1068 proline.put("E", new Integer(0));
1069 proline.put("Q", new Integer(0));
1070 proline.put("D", new Integer(0));
1071 proline.put("N", new Integer(0));
1072 proline.put("S", new Integer(0));
1073 proline.put("T", new Integer(0));
1074 proline.put("P", new Integer(1));
1075 proline.put("-", new Integer(1));
1076 proline.put("*", new Integer(1));
1081 propHash.put("hydrophobic", hydrophobic);
1082 propHash.put("small", small);
1083 propHash.put("positive", positive);
1084 propHash.put("negative", negative);
1085 propHash.put("charged", charged);
1086 propHash.put("aromatic", aromatic);
1087 propHash.put("aliphatic", aliphatic);
1088 propHash.put("tiny", tiny);
1089 propHash.put("proline", proline);
1090 propHash.put("polar", polar);
1093 private ResidueProperties()
1097 public static double getHydmax()
1102 public static double getHydmin()
1107 public static double[] getHyd()
1112 public static Hashtable getAA3Hash()
1117 public static int[][] getDNA()
1119 return ResidueProperties.DNA;
1122 public static int[][] getBLOSUM62()
1124 return ResidueProperties.BLOSUM62;
1127 public static int getPAM250(String A1, String A2)
1129 return getPAM250(A1.charAt(0), A2.charAt(0));
1132 public static int getBLOSUM62(char c1, char c2)
1138 int a = aaIndex[c1];
1139 int b = aaIndex[c2];
1141 pog = ResidueProperties.BLOSUM62[a][b];
1142 } catch (Exception e)
1144 // System.out.println("Unknown residue in " + A1 + " " + A2);
1150 public static Vector getCodons(String res)
1152 if (codonHash.containsKey(res))
1154 return (Vector) codonHash.get(res);
1160 public static String codonTranslate(String lccodon)
1162 String codon=lccodon.toUpperCase();
1163 // all base ambiguity codes yield an 'X' amino acid residue
1164 if (codon.indexOf('X')>-1 || codon.indexOf('N')>-1)
1168 Enumeration e = codonHash.keys();
1170 while (e.hasMoreElements())
1172 String key = (String) e.nextElement();
1173 Vector tmp = (Vector) codonHash.get(key);
1175 if (tmp.contains(codon))
1184 public static int[][] getDefaultPeptideMatrix()
1186 return ResidueProperties.getBLOSUM62();
1189 public static int[][] getDefaultDnaMatrix()
1191 return ResidueProperties.getDNA();
1195 * get a ScoreMatrix based on its string name
1198 * @return matrix in scoreMatrices with key pwtype or null
1200 public static ScoreMatrix getScoreMatrix(String pwtype)
1202 Object val = scoreMatrices.get(pwtype);
1205 return (ScoreMatrix) val;
1210 public static int getPAM250(char c, char d)
1215 int pog = ResidueProperties.PAM250[a][b];
1220 public static Hashtable toDssp3State;
1223 toDssp3State = new Hashtable();
1224 toDssp3State.put("H", "H");
1225 toDssp3State.put("E", "E");
1226 toDssp3State.put("C", " ");
1227 toDssp3State.put(" ", " ");
1228 toDssp3State.put("T", " ");
1229 toDssp3State.put("B", "E");
1230 toDssp3State.put("G", "H");
1231 toDssp3State.put("I", "H");
1232 toDssp3State.put("X", " ");
1236 * translate from other dssp secondary structure alphabets to 3-state
1239 * @return ssstring as a three-state secondary structure assignment.
1241 public static String getDssp3state(String ssstring)
1243 if (ssstring == null)
1247 StringBuffer ss = new StringBuffer();
1248 for (int i = 0; i < ssstring.length(); i++)
1250 String ssc = ssstring.substring(i, i + 1);
1251 if (toDssp3State.containsKey(ssc))
1253 ss.append((String) toDssp3State.get(ssc));
1260 return ss.toString();
1263 // main method generates perl representation of residue property hash
1265 public static void main(String[] args)
1267 Hashtable aa = new Hashtable();
1268 System.out.println("my %aa = {");
1269 // invert property hashes
1270 Enumeration prop = propHash.keys();
1271 while (prop.hasMoreElements())
1273 String pname = (String) prop.nextElement();
1274 Hashtable phash = (Hashtable) propHash.get(pname);
1275 Enumeration res = phash.keys();
1276 while (res.hasMoreElements())
1278 String rname = (String) res.nextElement();
1279 Vector aprops = (Vector) aa.get(rname);
1282 aprops = new Vector();
1283 aa.put(rname, aprops);
1285 Integer hasprop = (Integer) phash.get(rname);
1286 if (hasprop.intValue() == 1)
1288 aprops.addElement(pname);
1292 Enumeration res = aa.keys();
1293 while (res.hasMoreElements())
1295 String rname = (String) res.nextElement();
1297 System.out.print("'" + rname + "' => [");
1298 Enumeration props = ((Vector) aa.get(rname)).elements();
1299 while (props.hasMoreElements())
1301 System.out.print("'" + (String) props.nextElement() + "'");
1302 if (props.hasMoreElements())
1304 System.out.println(", ");
1307 System.out.println("]" + (res.hasMoreElements() ? "," : ""));
1309 System.out.println("};");