@datagrok-libraries/bio 0.0.3 → 0.0.4

package/package.json

@@ -5,7 +5,7 @@
   },
   "beta": true,
   "friendlyName": "Datagrok bio library",
-  "version": "0.0.3",
+  "version": "0.0.4",
   "description": "",
   "dependencies": {
     "datagrok-api": ">=0.108.0",

package/src/sequence-encoder.js

@@ -0,0 +1,155 @@
+import { assert } from '@datagrok-libraries/utils/src/operations';
+class SideChainScales {
+    static getAvailableScales() {
+        return Object.entries(this.scales).map(([k, _]) => k);
+    }
+    static getScale(name) {
+        assert(!(this.scales[name] === undefined), `Scale '${name}' was not found.`);
+        return this.scales[name];
+    }
+}
+SideChainScales.scales = {
+    // Wimley-White interfacial hydrophobicity scale
+    'WimleyWhite': {
+        '-': 0,
+        'A': 0.17,
+        'C': -0.24,
+        'D': -0.07,
+        'E': -0.01,
+        'F': -1.13,
+        'G': 0.01,
+        'H': 0.17,
+        'I': -0.31,
+        'K': 0.99,
+        'L': -0.56,
+        'M': -0.23,
+        'N': 0.42,
+        'P': 0.45,
+        'Q': 0.58,
+        'R': 0.81,
+        'S': 0.13,
+        'T': 0.14,
+        'V': 0.07,
+        'W': -1.85,
+        'Y': -0.94,
+    },
+    'categorial': {
+        '-': 0,
+        'A': 1,
+        'C': 2,
+        'D': 3,
+        'E': 4,
+        'F': 5,
+        'G': 6,
+        'H': 7,
+        'I': 8,
+        'K': 9,
+        'L': 10,
+        'M': 11,
+        'N': 12,
+        'P': 13,
+        'Q': 14,
+        'R': 15,
+        'S': 16,
+        'T': 17,
+        'V': 18,
+        'W': 19,
+        'Y': 20,
+    },
+};
+/**
+ * Class to categorial encode/decode aligned amino acid residues sequence.
+ *
+ * @export
+ * @class AlignedSequenceEncoder
+ */
+export class AlignedSequenceEncoder {
+    constructor(scale = 'categorial') {
+        this.aa2num = SideChainScales.getScale(scale);
+        this.num2aa = {};
+        Object.entries(this.aa2num).forEach(([k, v]) => (this.num2aa[v] = k));
+    }
+    /**
+       * Truncate NH2 and -COOH terminals of the given sequence.
+       *
+       * @static
+       * @param {string} seq The sequence provided.
+       * @return {string} Truncated sequence.
+       * @memberof AlignedSequenceEncoder
+       */
+    static _truncateSequence(seq) {
+        let start = 0;
+        let end = seq.length;
+        const termina = ['NH2', 'COOH'];
+        if (seq.startsWith(termina[0])) {
+            const l = termina[0].length; // Cut only 'NH2' without following '-'.
+            assert(seq[l] == '-', `Wrong sequence format: ${termina[0]} without following '-' in '${seq}'.`);
+            start = l;
+        }
+        if (seq.endsWith(termina[1])) {
+            const l = termina[1].length + 1; // Cut both 'COOH' and precending '-'.
+            assert(seq[end - l] == '-', `Wrong sequence format: ${termina[1]} without '-' precending in '${seq}'.`);
+            end -= l;
+        }
+        return seq.substring(start, end);
+    }
+    /**
+       * Cuts auxiliary defises before a residue.
+       *
+       * @static
+       * @param {string} seq The sequence to process.
+       * @return {string} Processed sequence.
+       * @memberof AlignedSequenceEncoder
+       */
+    static _dropDefises(seq) {
+        return seq.replace(/(-)([^-]+)/g, '$2');
+    }
+    /**
+       * Performs truncation and cutting auxiliary defises.
+       *
+       * @static
+       * @param {string} sequence The sequence work under process.
+       * @return {string} Result of cleaning.
+       * @memberof AlignedSequenceEncoder
+       */
+    static clean(sequence) {
+        return AlignedSequenceEncoder._dropDefises(AlignedSequenceEncoder._truncateSequence(sequence));
+    }
+    /**
+       * Categorial encode of the sequence provided.
+       *
+       * @param {string} sequence The sequence.
+       * @return {number[]} Encoded vector.
+       * @memberof AlignedSequenceEncoder
+       */
+    encode(sequence) {
+        const nItems = sequence.length;
+        const values = new Array(nItems).fill(0);
+        for (let i = 0; i < nItems; ++i) {
+            const char = sequence[i];
+            assert(char in this.aa2num, `Unknown char '${char}' found in sequence '${sequence}'`);
+            values[i] = this.encodeLettter(char);
+        }
+        return values;
+    }
+    encodeLettter(letter) {
+        return this.aa2num[letter];
+    }
+    /**
+       * Decode the encoded vector into the sequence back.
+       *
+       * @param {number[]} value The vector encoded.
+       * @return {string} Decoded sequence.
+       * @memberof AlignedSequenceEncoder
+       */
+    decode(value) {
+        let s = '';
+        for (let i = 0; i < value.length; ++i) {
+            const code = value[i];
+            assert(code in this.num2aa, `Unknown code '${code}' found in vector '${value}'`);
+            s += this.num2aa[code];
+        }
+        return s;
+    }
+}
+//# sourceMappingURL=data:application/json;base64,eyJ2ZXJzaW9uIjozLCJmaWxlIjoic2VxdWVuY2UtZW5jb2Rlci5qcyIsInNvdXJjZVJvb3QiOiIiLCJzb3VyY2VzIjpbInNlcXVlbmNlLWVuY29kZXIudHMiXSwibmFtZXMiOltdLCJtYXBwaW5ncyI6IkFBQUEsT0FBTyxFQUFDLE1BQU0sRUFBQyxNQUFNLDBDQUEwQyxDQUFDO0FBS2hFLE1BQU0sZUFBZTtJQW1EbkIsTUFBTSxDQUFDLGtCQUFrQjtRQUN2QixPQUFPLE1BQU0sQ0FBQyxPQUFPLENBQUMsSUFBSSxDQUFDLE1BQU0sQ0FBQyxDQUFDLEdBQUcsQ0FBQyxDQUFDLENBQUMsQ0FBQyxFQUFFLENBQUMsQ0FBQyxFQUFFLEVBQUUsQ0FBQyxDQUFDLENBQUMsQ0FBQztJQUN4RCxDQUFDO0lBRUQsTUFBTSxDQUFDLFFBQVEsQ0FBQyxJQUFZO1FBQzFCLE1BQU0sQ0FBQyxDQUFDLENBQUMsSUFBSSxDQUFDLE1BQU0sQ0FBQyxJQUFJLENBQUMsS0FBSyxTQUFTLENBQUMsRUFBRSxVQUFVLElBQUksa0JBQWtCLENBQUMsQ0FBQztRQUM3RSxPQUFPLElBQUksQ0FBQyxNQUFNLENBQUMsSUFBSSxDQUFDLENBQUM7SUFDM0IsQ0FBQzs7QUF6RE0sc0JBQU0sR0FBNkI7SUFDeEMsZ0RBQWdEO0lBQ2hELGFBQWEsRUFBRTtRQUNiLEdBQUcsRUFBRSxDQUFDO1FBQ04sR0FBRyxFQUFFLElBQUk7UUFDVCxHQUFHLEVBQUUsQ0FBQyxJQUFJO1FBQ1YsR0FBRyxFQUFFLENBQUMsSUFBSTtRQUNWLEdBQUcsRUFBRSxDQUFDLElBQUk7UUFDVixHQUFHLEVBQUUsQ0FBQyxJQUFJO1FBQ1YsR0FBRyxFQUFFLElBQUk7UUFDVCxHQUFHLEVBQUUsSUFBSTtRQUNULEdBQUcsRUFBRSxDQUFDLElBQUk7UUFDVixHQUFHLEVBQUUsSUFBSTtRQUNULEdBQUcsRUFBRSxDQUFDLElBQUk7UUFDVixHQUFHLEVBQUUsQ0FBQyxJQUFJO1FBQ1YsR0FBRyxFQUFFLElBQUk7UUFDVCxHQUFHLEVBQUUsSUFBSTtRQUNULEdBQUcsRUFBRSxJQUFJO1FBQ1QsR0FBRyxFQUFFLElBQUk7UUFDVCxHQUFHLEVBQUUsSUFBSTtRQUNULEdBQUcsRUFBRSxJQUFJO1FBQ1QsR0FBRyxFQUFFLElBQUk7UUFDVCxHQUFHLEVBQUUsQ0FBQyxJQUFJO1FBQ1YsR0FBRyxFQUFFLENBQUMsSUFBSTtLQUNYO0lBQ0QsWUFBWSxFQUFFO1FBQ1osR0FBRyxFQUFFLENBQUM7UUFDTixHQUFHLEVBQUUsQ0FBQztRQUNOLEdBQUcsRUFBRSxDQUFDO1FBQ04sR0FBRyxFQUFFLENBQUM7UUFDTixHQUFHLEVBQUUsQ0FBQztRQUNOLEdBQUcsRUFBRSxDQUFDO1FBQ04sR0FBRyxFQUFFLENBQUM7UUFDTixHQUFHLEVBQUUsQ0FBQztRQUNOLEdBQUcsRUFBRSxDQUFDO1FBQ04sR0FBRyxFQUFFLENBQUM7UUFDTixHQUFHLEVBQUUsRUFBRTtRQUNQLEdBQUcsRUFBRSxFQUFFO1FBQ1AsR0FBRyxFQUFFLEVBQUU7UUFDUCxHQUFHLEVBQUUsRUFBRTtRQUNQLEdBQUcsRUFBRSxFQUFFO1FBQ1AsR0FBRyxFQUFFLEVBQUU7UUFDUCxHQUFHLEVBQUUsRUFBRTtRQUNQLEdBQUcsRUFBRSxFQUFFO1FBQ1AsR0FBRyxFQUFFLEVBQUU7UUFDUCxHQUFHLEVBQUUsRUFBRTtRQUNQLEdBQUcsRUFBRSxFQUFFO0tBQ1I7Q0FDRixDQUFDO0FBWUo7Ozs7O0dBS0c7QUFDSCxNQUFNLE9BQU8sc0JBQXNCO0lBSWpDLFlBQVksUUFBZ0IsWUFBWTtRQUN0QyxJQUFJLENBQUMsTUFBTSxHQUFHLGVBQWUsQ0FBQyxRQUFRLENBQUMsS0FBSyxDQUFDLENBQUM7UUFDOUMsSUFBSSxDQUFDLE1BQU0sR0FBRyxFQUFFLENBQUM7UUFDakIsTUFBTSxDQUFDLE9BQU8sQ0FBQyxJQUFJLENBQUMsTUFBTSxDQUFDLENBQUMsT0FBTyxDQUFDLENBQUMsQ0FBQyxDQUFDLEVBQUUsQ0FBQyxDQUFDLEVBQUUsRUFBRSxDQUFDLENBQUMsSUFBSSxDQUFDLE1BQU0sQ0FBQyxDQUFDLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQyxDQUFDO0lBQ3hFLENBQUM7SUFFRDs7Ozs7OztTQU9LO0lBQ0wsTUFBTSxDQUFDLGlCQUFpQixDQUFDLEdBQVc7UUFDbEMsSUFBSSxLQUFLLEdBQUcsQ0FBQyxDQUFDO1FBQ2QsSUFBSSxHQUFHLEdBQUcsR0FBRyxDQUFDLE1BQU0sQ0FBQztRQUNyQixNQUFNLE9BQU8sR0FBRyxDQUFDLEtBQUssRUFBRSxNQUFNLENBQUMsQ0FBQztRQUVoQyxJQUFJLEdBQUcsQ0FBQyxVQUFVLENBQUMsT0FBTyxDQUFDLENBQUMsQ0FBQyxDQUFDLEVBQUU7WUFDOUIsTUFBTSxDQUFDLEdBQUcsT0FBTyxDQUFDLENBQUMsQ0FBQyxDQUFDLE1BQU0sQ0FBQyxDQUFDLHdDQUF3QztZQUNyRSxNQUFNLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQyxJQUFJLEdBQUcsRUFBRSwwQkFBMEIsT0FBTyxDQUFDLENBQUMsQ0FBQyw4QkFBOEIsR0FBRyxJQUFJLENBQUMsQ0FBQztZQUNqRyxLQUFLLEdBQUcsQ0FBQyxDQUFDO1NBQ1g7UUFDRCxJQUFJLEdBQUcsQ0FBQyxRQUFRLENBQUMsT0FBTyxDQUFDLENBQUMsQ0FBQyxDQUFDLEVBQUU7WUFDNUIsTUFBTSxDQUFDLEdBQUcsT0FBTyxDQUFDLENBQUMsQ0FBQyxDQUFDLE1BQU0sR0FBQyxDQUFDLENBQUMsQ0FBQyxzQ0FBc0M7WUFDckUsTUFBTSxDQUFDLEdBQUcsQ0FBQyxHQUFHLEdBQUMsQ0FBQyxDQUFDLElBQUksR0FBRyxFQUFFLDBCQUEwQixPQUFPLENBQUMsQ0FBQyxDQUFDLCtCQUErQixHQUFHLElBQUksQ0FBQyxDQUFDO1lBQ3RHLEdBQUcsSUFBSSxDQUFDLENBQUM7U0FDVjtRQUNELE9BQU8sR0FBRyxDQUFDLFNBQVMsQ0FBQyxLQUFLLEVBQUUsR0FBRyxDQUFDLENBQUM7SUFDbkMsQ0FBQztJQUVEOzs7Ozs7O1NBT0s7SUFDTCxNQUFNLENBQUMsWUFBWSxDQUFDLEdBQVc7UUFDN0IsT0FBTyxHQUFHLENBQUMsT0FBTyxDQUFDLGFBQWEsRUFBRSxJQUFJLENBQUMsQ0FBQztJQUMxQyxDQUFDO0lBRUQ7Ozs7Ozs7U0FPSztJQUNMLE1BQU0sQ0FBQyxLQUFLLENBQUMsUUFBZ0I7UUFDM0IsT0FBTyxzQkFBc0IsQ0FBQyxZQUFZLENBQUMsc0JBQXNCLENBQUMsaUJBQWlCLENBQUMsUUFBUSxDQUFDLENBQUMsQ0FBQztJQUNqRyxDQUFDO0lBRUQ7Ozs7OztTQU1LO0lBQ0UsTUFBTSxDQUFDLFFBQWdCO1FBQzVCLE1BQU0sTUFBTSxHQUFHLFFBQVEsQ0FBQyxNQUFNLENBQUM7UUFDL0IsTUFBTSxNQUFNLEdBQUcsSUFBSSxLQUFLLENBQUMsTUFBTSxDQUFDLENBQUMsSUFBSSxDQUFDLENBQUMsQ0FBQyxDQUFDO1FBRXpDLEtBQUssSUFBSSxDQUFDLEdBQUcsQ0FBQyxFQUFFLENBQUMsR0FBRyxNQUFNLEVBQUUsRUFBRSxDQUFDLEVBQUU7WUFDL0IsTUFBTSxJQUFJLEdBQUcsUUFBUSxDQUFDLENBQUMsQ0FBQyxDQUFDO1lBRXpCLE1BQU0sQ0FBQyxJQUFJLElBQUksSUFBSSxDQUFDLE1BQU0sRUFBRSxpQkFBaUIsSUFBSSx3QkFBd0IsUUFBUSxHQUFHLENBQUMsQ0FBQztZQUV0RixNQUFNLENBQUMsQ0FBQyxDQUFDLEdBQUcsSUFBSSxDQUFDLGFBQWEsQ0FBQyxJQUFJLENBQUMsQ0FBQztTQUN0QztRQUNELE9BQU8sTUFBTSxDQUFDO0lBQ2hCLENBQUM7SUFFTSxhQUFhLENBQUMsTUFBYztRQUNqQyxPQUFPLElBQUksQ0FBQyxNQUFNLENBQUMsTUFBTSxDQUFDLENBQUM7SUFDN0IsQ0FBQztJQUVEOzs7Ozs7U0FNSztJQUNFLE1BQU0sQ0FBQyxLQUFlO1FBQzNCLElBQUksQ0FBQyxHQUFXLEVBQUUsQ0FBQztRQUVuQixLQUFLLElBQUksQ0FBQyxHQUFHLENBQUMsRUFBRSxDQUFDLEdBQUcsS0FBSyxDQUFDLE1BQU0sRUFBRSxFQUFFLENBQUMsRUFBRTtZQUNyQyxNQUFNLElBQUksR0FBRyxLQUFLLENBQUMsQ0FBQyxDQUFDLENBQUM7WUFFdEIsTUFBTSxDQUFDLElBQUksSUFBSSxJQUFJLENBQUMsTUFBTSxFQUFFLGlCQUFpQixJQUFJLHNCQUFzQixLQUFLLEdBQUcsQ0FBQyxDQUFDO1lBRWpGLENBQUMsSUFBSSxJQUFJLENBQUMsTUFBTSxDQUFDLElBQUksQ0FBQyxDQUFDO1NBQ3hCO1FBQ0QsT0FBTyxDQUFDLENBQUM7SUFDWCxDQUFDO0NBQ0YiLCJzb3VyY2VzQ29udGVudCI6WyJpbXBvcnQge2Fzc2VydH0gZnJvbSAnQGRhdGFncm9rLWxpYnJhcmllcy91dGlscy9zcmMvb3BlcmF0aW9ucyc7XHJcblxyXG50eXBlIFNpZGVDaGFpblNjYWxlID0ge1tuYW1lOiBzdHJpbmddOiBudW1iZXJ9O1xyXG50eXBlIFNpZGVDaGFpblNjYWxlQ29sbGVjdGlvbiA9IHtbbmFtZTogc3RyaW5nXTogU2lkZUNoYWluU2NhbGV9O1xyXG5cclxuY2xhc3MgU2lkZUNoYWluU2NhbGVzIHtcclxuICBzdGF0aWMgc2NhbGVzOiBTaWRlQ2hhaW5TY2FsZUNvbGxlY3Rpb24gPSB7XHJcbiAgICAvLyBXaW1sZXktV2hpdGUgaW50ZXJmYWNpYWwgaHlkcm9waG9iaWNpdHkgc2NhbGVcclxuICAgICdXaW1sZXlXaGl0ZSc6IHtcclxuICAgICAgJy0nOiAwLFxyXG4gICAgICAnQSc6IDAuMTcsXHJcbiAgICAgICdDJzogLTAuMjQsXHJcbiAgICAgICdEJzogLTAuMDcsIC8vIEFzcC06IDEuMjNcclxuICAgICAgJ0UnOiAtMC4wMSwgLy8gR2x1LTogMi4wMlxyXG4gICAgICAnRic6IC0xLjEzLCAvL1xyXG4gICAgICAnRyc6IDAuMDEsXHJcbiAgICAgICdIJzogMC4xNywgLy8gSGlzKzogMC45NlxyXG4gICAgICAnSSc6IC0wLjMxLFxyXG4gICAgICAnSyc6IDAuOTksIC8vIEx5cytcclxuICAgICAgJ0wnOiAtMC41NixcclxuICAgICAgJ00nOiAtMC4yMyxcclxuICAgICAgJ04nOiAwLjQyLFxyXG4gICAgICAnUCc6IDAuNDUsXHJcbiAgICAgICdRJzogMC41OCxcclxuICAgICAgJ1InOiAwLjgxLCAvLyBBcmcrXHJcbiAgICAgICdTJzogMC4xMyxcclxuICAgICAgJ1QnOiAwLjE0LFxyXG4gICAgICAnVic6IDAuMDcsXHJcbiAgICAgICdXJzogLTEuODUsXHJcbiAgICAgICdZJzogLTAuOTQsXHJcbiAgICB9LFxyXG4gICAgJ2NhdGVnb3JpYWwnOiB7XHJcbiAgICAgICctJzogMCxcclxuICAgICAgJ0EnOiAxLFxyXG4gICAgICAnQyc6IDIsXHJcbiAgICAgICdEJzogMyxcclxuICAgICAgJ0UnOiA0LFxyXG4gICAgICAnRic6IDUsXHJcbiAgICAgICdHJzogNixcclxuICAgICAgJ0gnOiA3LFxyXG4gICAgICAnSSc6IDgsXHJcbiAgICAgICdLJzogOSxcclxuICAgICAgJ0wnOiAxMCxcclxuICAgICAgJ00nOiAxMSxcclxuICAgICAgJ04nOiAxMixcclxuICAgICAgJ1AnOiAxMyxcclxuICAgICAgJ1EnOiAxNCxcclxuICAgICAgJ1InOiAxNSxcclxuICAgICAgJ1MnOiAxNixcclxuICAgICAgJ1QnOiAxNyxcclxuICAgICAgJ1YnOiAxOCxcclxuICAgICAgJ1cnOiAxOSxcclxuICAgICAgJ1knOiAyMCxcclxuICAgIH0sXHJcbiAgfTtcclxuXHJcbiAgc3RhdGljIGdldEF2YWlsYWJsZVNjYWxlcygpOiBzdHJpbmdbXSB7XHJcbiAgICByZXR1cm4gT2JqZWN0LmVudHJpZXModGhpcy5zY2FsZXMpLm1hcCgoW2ssIF9dKSA9PiBrKTtcclxuICB9XHJcblxyXG4gIHN0YXRpYyBnZXRTY2FsZShuYW1lOiBzdHJpbmcpOiBTaWRlQ2hhaW5TY2FsZSB7XHJcbiAgICBhc3NlcnQoISh0aGlzLnNjYWxlc1tuYW1lXSA9PT0gdW5kZWZpbmVkKSwgYFNjYWxlICcke25hbWV9JyB3YXMgbm90IGZvdW5kLmApO1xyXG4gICAgcmV0dXJuIHRoaXMuc2NhbGVzW25hbWVdO1xyXG4gIH1cclxufVxyXG5cclxuLyoqXHJcbiAqIENsYXNzIHRvIGNhdGVnb3JpYWwgZW5jb2RlL2RlY29kZSBhbGlnbmVkIGFtaW5vIGFjaWQgcmVzaWR1ZXMgc2VxdWVuY2UuXHJcbiAqXHJcbiAqIEBleHBvcnRcclxuICogQGNsYXNzIEFsaWduZWRTZXF1ZW5jZUVuY29kZXJcclxuICovXHJcbmV4cG9ydCBjbGFzcyBBbGlnbmVkU2VxdWVuY2VFbmNvZGVyIHtcclxuICBwcm90ZWN0ZWQgYWEybnVtOiBTaWRlQ2hhaW5TY2FsZTtcclxuICBwcm90ZWN0ZWQgbnVtMmFhOiB7W2NvZGU6IG51bWJlcl06IHN0cmluZ307XHJcblxyXG4gIGNvbnN0cnVjdG9yKHNjYWxlOiBzdHJpbmcgPSAnY2F0ZWdvcmlhbCcpIHtcclxuICAgIHRoaXMuYWEybnVtID0gU2lkZUNoYWluU2NhbGVzLmdldFNjYWxlKHNjYWxlKTtcclxuICAgIHRoaXMubnVtMmFhID0ge307XHJcbiAgICBPYmplY3QuZW50cmllcyh0aGlzLmFhMm51bSkuZm9yRWFjaCgoW2ssIHZdKSA9PiAodGhpcy5udW0yYWFbdl0gPSBrKSk7XHJcbiAgfVxyXG5cclxuICAvKipcclxuICAgICAqIFRydW5jYXRlIE5IMiBhbmQgLUNPT0ggdGVybWluYWxzIG9mIHRoZSBnaXZlbiBzZXF1ZW5jZS5cclxuICAgICAqXHJcbiAgICAgKiBAc3RhdGljXHJcbiAgICAgKiBAcGFyYW0ge3N0cmluZ30gc2VxIFRoZSBzZXF1ZW5jZSBwcm92aWRlZC5cclxuICAgICAqIEByZXR1cm4ge3N0cmluZ30gVHJ1bmNhdGVkIHNlcXVlbmNlLlxyXG4gICAgICogQG1lbWJlcm9mIEFsaWduZWRTZXF1ZW5jZUVuY29kZXJcclxuICAgICAqL1xyXG4gIHN0YXRpYyBfdHJ1bmNhdGVTZXF1ZW5jZShzZXE6IHN0cmluZyk6IHN0cmluZyB7XHJcbiAgICBsZXQgc3RhcnQgPSAwO1xyXG4gICAgbGV0IGVuZCA9IHNlcS5sZW5ndGg7XHJcbiAgICBjb25zdCB0ZXJtaW5hID0gWydOSDInLCAnQ09PSCddO1xyXG5cclxuICAgIGlmIChzZXEuc3RhcnRzV2l0aCh0ZXJtaW5hWzBdKSkge1xyXG4gICAgICBjb25zdCBsID0gdGVybWluYVswXS5sZW5ndGg7IC8vIEN1dCBvbmx5ICdOSDInIHdpdGhvdXQgZm9sbG93aW5nICctJy5cclxuICAgICAgYXNzZXJ0KHNlcVtsXSA9PSAnLScsIGBXcm9uZyBzZXF1ZW5jZSBmb3JtYXQ6ICR7dGVybWluYVswXX0gd2l0aG91dCBmb2xsb3dpbmcgJy0nIGluICcke3NlcX0nLmApO1xyXG4gICAgICBzdGFydCA9IGw7XHJcbiAgICB9XHJcbiAgICBpZiAoc2VxLmVuZHNXaXRoKHRlcm1pbmFbMV0pKSB7XHJcbiAgICAgIGNvbnN0IGwgPSB0ZXJtaW5hWzFdLmxlbmd0aCsxOyAvLyBDdXQgYm90aCAnQ09PSCcgYW5kIHByZWNlbmRpbmcgJy0nLlxyXG4gICAgICBhc3NlcnQoc2VxW2VuZC1sXSA9PSAnLScsIGBXcm9uZyBzZXF1ZW5jZSBmb3JtYXQ6ICR7dGVybWluYVsxXX0gd2l0aG91dCAnLScgcHJlY2VuZGluZyBpbiAnJHtzZXF9Jy5gKTtcclxuICAgICAgZW5kIC09IGw7XHJcbiAgICB9XHJcbiAgICByZXR1cm4gc2VxLnN1YnN0cmluZyhzdGFydCwgZW5kKTtcclxuICB9XHJcblxyXG4gIC8qKlxyXG4gICAgICogQ3V0cyBhdXhpbGlhcnkgZGVmaXNlcyBiZWZvcmUgYSByZXNpZHVlLlxyXG4gICAgICpcclxuICAgICAqIEBzdGF0aWNcclxuICAgICAqIEBwYXJhbSB7c3RyaW5nfSBzZXEgVGhlIHNlcXVlbmNlIHRvIHByb2Nlc3MuXHJcbiAgICAgKiBAcmV0dXJuIHtzdHJpbmd9IFByb2Nlc3NlZCBzZXF1ZW5jZS5cclxuICAgICAqIEBtZW1iZXJvZiBBbGlnbmVkU2VxdWVuY2VFbmNvZGVyXHJcbiAgICAgKi9cclxuICBzdGF0aWMgX2Ryb3BEZWZpc2VzKHNlcTogc3RyaW5nKTogc3RyaW5nIHtcclxuICAgIHJldHVybiBzZXEucmVwbGFjZSgvKC0pKFteLV0rKS9nLCAnJDInKTtcclxuICB9XHJcblxyXG4gIC8qKlxyXG4gICAgICogUGVyZm9ybXMgdHJ1bmNhdGlvbiBhbmQgY3V0dGluZyBhdXhpbGlhcnkgZGVmaXNlcy5cclxuICAgICAqXHJcbiAgICAgKiBAc3RhdGljXHJcbiAgICAgKiBAcGFyYW0ge3N0cmluZ30gc2VxdWVuY2UgVGhlIHNlcXVlbmNlIHdvcmsgdW5kZXIgcHJvY2Vzcy5cclxuICAgICAqIEByZXR1cm4ge3N0cmluZ30gUmVzdWx0IG9mIGNsZWFuaW5nLlxyXG4gICAgICogQG1lbWJlcm9mIEFsaWduZWRTZXF1ZW5jZUVuY29kZXJcclxuICAgICAqL1xyXG4gIHN0YXRpYyBjbGVhbihzZXF1ZW5jZTogc3RyaW5nKTogc3RyaW5nIHtcclxuICAgIHJldHVybiBBbGlnbmVkU2VxdWVuY2VFbmNvZGVyLl9kcm9wRGVmaXNlcyhBbGlnbmVkU2VxdWVuY2VFbmNvZGVyLl90cnVuY2F0ZVNlcXVlbmNlKHNlcXVlbmNlKSk7XHJcbiAgfVxyXG5cclxuICAvKipcclxuICAgICAqIENhdGVnb3JpYWwgZW5jb2RlIG9mIHRoZSBzZXF1ZW5jZSBwcm92aWRlZC5cclxuICAgICAqXHJcbiAgICAgKiBAcGFyYW0ge3N0cmluZ30gc2VxdWVuY2UgVGhlIHNlcXVlbmNlLlxyXG4gICAgICogQHJldHVybiB7bnVtYmVyW119IEVuY29kZWQgdmVjdG9yLlxyXG4gICAgICogQG1lbWJlcm9mIEFsaWduZWRTZXF1ZW5jZUVuY29kZXJcclxuICAgICAqL1xyXG4gIHB1YmxpYyBlbmNvZGUoc2VxdWVuY2U6IHN0cmluZyk6IG51bWJlcltdIHtcclxuICAgIGNvbnN0IG5JdGVtcyA9IHNlcXVlbmNlLmxlbmd0aDtcclxuICAgIGNvbnN0IHZhbHVlcyA9IG5ldyBBcnJheShuSXRlbXMpLmZpbGwoMCk7XHJcblxyXG4gICAgZm9yIChsZXQgaSA9IDA7IGkgPCBuSXRlbXM7ICsraSkge1xyXG4gICAgICBjb25zdCBjaGFyID0gc2VxdWVuY2VbaV07XHJcblxyXG4gICAgICBhc3NlcnQoY2hhciBpbiB0aGlzLmFhMm51bSwgYFVua25vd24gY2hhciAnJHtjaGFyfScgZm91bmQgaW4gc2VxdWVuY2UgJyR7c2VxdWVuY2V9J2ApO1xyXG5cclxuICAgICAgdmFsdWVzW2ldID0gdGhpcy5lbmNvZGVMZXR0dGVyKGNoYXIpO1xyXG4gICAgfVxyXG4gICAgcmV0dXJuIHZhbHVlcztcclxuICB9XHJcblxyXG4gIHB1YmxpYyBlbmNvZGVMZXR0dGVyKGxldHRlcjogc3RyaW5nKTogbnVtYmVyIHtcclxuICAgIHJldHVybiB0aGlzLmFhMm51bVtsZXR0ZXJdO1xyXG4gIH1cclxuXHJcbiAgLyoqXHJcbiAgICAgKiBEZWNvZGUgdGhlIGVuY29kZWQgdmVjdG9yIGludG8gdGhlIHNlcXVlbmNlIGJhY2suXHJcbiAgICAgKlxyXG4gICAgICogQHBhcmFtIHtudW1iZXJbXX0gdmFsdWUgVGhlIHZlY3RvciBlbmNvZGVkLlxyXG4gICAgICogQHJldHVybiB7c3RyaW5nfSBEZWNvZGVkIHNlcXVlbmNlLlxyXG4gICAgICogQG1lbWJlcm9mIEFsaWduZWRTZXF1ZW5jZUVuY29kZXJcclxuICAgICAqL1xyXG4gIHB1YmxpYyBkZWNvZGUodmFsdWU6IG51bWJlcltdKTogc3RyaW5nIHtcclxuICAgIGxldCBzOiBzdHJpbmcgPSAnJztcclxuXHJcbiAgICBmb3IgKGxldCBpID0gMDsgaSA8IHZhbHVlLmxlbmd0aDsgKytpKSB7XHJcbiAgICAgIGNvbnN0IGNvZGUgPSB2YWx1ZVtpXTtcclxuXHJcbiAgICAgIGFzc2VydChjb2RlIGluIHRoaXMubnVtMmFhLCBgVW5rbm93biBjb2RlICcke2NvZGV9JyBmb3VuZCBpbiB2ZWN0b3IgJyR7dmFsdWV9J2ApO1xyXG5cclxuICAgICAgcyArPSB0aGlzLm51bTJhYVtjb2RlXTtcclxuICAgIH1cclxuICAgIHJldHVybiBzO1xyXG4gIH1cclxufVxyXG4iXX0=

package/src/sequence-encoder.ts

@@ -1,177 +0,0 @@
-import {assert} from '@datagrok-libraries/utils/src/operations';
-
-type SideChainScale = {[name: string]: number};
-type SideChainScaleCollection = {[name: string]: SideChainScale};
-
-class SideChainScales {
-  static scales: SideChainScaleCollection = {
-    // Wimley-White interfacial hydrophobicity scale
-    'WimleyWhite': {
-      '-': 0,
-      'A': 0.17,
-      'C': -0.24,
-      'D': -0.07, // Asp-: 1.23
-      'E': -0.01, // Glu-: 2.02
-      'F': -1.13, //
-      'G': 0.01,
-      'H': 0.17, // His+: 0.96
-      'I': -0.31,
-      'K': 0.99, // Lys+
-      'L': -0.56,
-      'M': -0.23,
-      'N': 0.42,
-      'P': 0.45,
-      'Q': 0.58,
-      'R': 0.81, // Arg+
-      'S': 0.13,
-      'T': 0.14,
-      'V': 0.07,
-      'W': -1.85,
-      'Y': -0.94,
-    },
-    'categorial': {
-      '-': 0,
-      'A': 1,
-      'C': 2,
-      'D': 3,
-      'E': 4,
-      'F': 5,
-      'G': 6,
-      'H': 7,
-      'I': 8,
-      'K': 9,
-      'L': 10,
-      'M': 11,
-      'N': 12,
-      'P': 13,
-      'Q': 14,
-      'R': 15,
-      'S': 16,
-      'T': 17,
-      'V': 18,
-      'W': 19,
-      'Y': 20,
-    },
-  };
-
-  static getAvailableScales(): string[] {
-    return Object.entries(this.scales).map(([k, _]) => k);
-  }
-
-  static getScale(name: string): SideChainScale {
-    assert(!(this.scales[name] === undefined), `Scale '${name}' was not found.`);
-    return this.scales[name];
-  }
-}
-
-/**
- * Class to categorial encode/decode aligned amino acid residues sequence.
- *
- * @export
- * @class AlignedSequenceEncoder
- */
-export class AlignedSequenceEncoder {
-  protected aa2num: SideChainScale;
-  protected num2aa: {[code: number]: string};
-
-  constructor(scale: string = 'categorial') {
-    this.aa2num = SideChainScales.getScale(scale);
-    this.num2aa = {};
-    Object.entries(this.aa2num).forEach(([k, v]) => (this.num2aa[v] = k));
-  }
-
-  /**
-     * Truncate NH2 and -COOH terminals of the given sequence.
-     *
-     * @static
-     * @param {string} seq The sequence provided.
-     * @return {string} Truncated sequence.
-     * @memberof AlignedSequenceEncoder
-     */
-  static _truncateSequence(seq: string): string {
-    let start = 0;
-    let end = seq.length;
-    const termina = ['NH2', 'COOH'];
-
-    if (seq.startsWith(termina[0])) {
-      const l = termina[0].length; // Cut only 'NH2' without following '-'.
-      assert(seq[l] == '-', `Wrong sequence format: ${termina[0]} without following '-' in '${seq}'.`);
-      start = l;
-    }
-    if (seq.endsWith(termina[1])) {
-      const l = termina[1].length+1; // Cut both 'COOH' and precending '-'.
-      assert(seq[end-l] == '-', `Wrong sequence format: ${termina[1]} without '-' precending in '${seq}'.`);
-      end -= l;
-    }
-    return seq.substring(start, end);
-  }
-
-  /**
-     * Cuts auxiliary defises before a residue.
-     *
-     * @static
-     * @param {string} seq The sequence to process.
-     * @return {string} Processed sequence.
-     * @memberof AlignedSequenceEncoder
-     */
-  static _dropDefises(seq: string): string {
-    return seq.replace(/(-)([^-]+)/g, '$2');
-  }
-
-  /**
-     * Performs truncation and cutting auxiliary defises.
-     *
-     * @static
-     * @param {string} sequence The sequence work under process.
-     * @return {string} Result of cleaning.
-     * @memberof AlignedSequenceEncoder
-     */
-  static clean(sequence: string): string {
-    return AlignedSequenceEncoder._dropDefises(AlignedSequenceEncoder._truncateSequence(sequence));
-  }
-
-  /**
-     * Categorial encode of the sequence provided.
-     *
-     * @param {string} sequence The sequence.
-     * @return {number[]} Encoded vector.
-     * @memberof AlignedSequenceEncoder
-     */
-  public encode(sequence: string): number[] {
-    const nItems = sequence.length;
-    const values = new Array(nItems).fill(0);
-
-    for (let i = 0; i < nItems; ++i) {
-      const char = sequence[i];
-
-      assert(char in this.aa2num, `Unknown char '${char}' found in sequence '${sequence}'`);
-
-      values[i] = this.encodeLettter(char);
-    }
-    return values;
-  }
-
-  public encodeLettter(letter: string): number {
-    return this.aa2num[letter];
-  }
-
-  /**
-     * Decode the encoded vector into the sequence back.
-     *
-     * @param {number[]} value The vector encoded.
-     * @return {string} Decoded sequence.
-     * @memberof AlignedSequenceEncoder
-     */
-  public decode(value: number[]): string {
-    let s: string = '';
-
-    for (let i = 0; i < value.length; ++i) {
-      const code = value[i];
-
-      assert(code in this.num2aa, `Unknown code '${code}' found in vector '${value}'`);
-
-      s += this.num2aa[code];
-    }
-    return s;
-  }
-}