@datagrok/peptides 1.3.3 → 1.3.5

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@datagrok/peptides",
-  "version": "1.3.3",
+  "version": "1.3.5",
   "author": {
     "name": "Volodymyr Dyma",
     "email": "vdyma@datagrok.ai"
@@ -12,10 +12,10 @@
     "directory": "packages/Peptides"
   },
   "dependencies": {
-    "@datagrok-libraries/bio": "^5.1.0",
+    "@datagrok-libraries/bio": "^5.5.0",
     "@datagrok-libraries/ml": "^2.0.1",
     "@datagrok-libraries/statistics": "^0.1.6",
-    "@datagrok-libraries/utils": "^1.10.1",
+    "@datagrok-libraries/utils": "^1.11.1",
     "cash-dom": "latest",
     "datagrok-api": "^1.6.0",
     "file-loader": "^6.2.0",

package/setup-unlink-clean.cmd

@@ -5,6 +5,8 @@ set dirs=^
 \..\..\libraries\utils\ ^
 \..\..\libraries\ml\ ^
 \..\..\libraries\bio\ ^
+\..\..\packages\Bio\ ^
+\..\..\packages\Helm\ ^
 \
 
 call npm uninstall -g ^

package/setup.cmd

@@ -7,6 +7,8 @@ set dirs=^
 \..\..\libraries\utils\ ^
 \..\..\libraries\ml\ ^
 \..\..\libraries\bio\ ^
+\..\..\packages\Bio\ ^
+\..\..\packages\Helm\ ^
 \
 
 for %%p in (%dirs%) do cd %package_dir%\%%p & call npm install

package/setup.sh

@@ -1,15 +1,38 @@
-npm unlink datagrok-api
-npm unlink @datagrok-libraries/utils
-npm unlink @datagrok-libraries/ml
-cd ../../js-api
-npm install
-npm link
-cd ../libraries/utils
-npm install
-npm link
-cd ../../libraries/ml
-npm install
-npm link datagrok-api @datagrok-libraries/utils
-cd ../../packages/Peptides
-npm install
-npm link datagrok-api @datagrok-libraries/utils @datagrok-libraries/ml 
+#!/bin/bash
+
+./setup-unlink-clean.sh
+
+GREEN='\e[0;32m'
+NO_COLOR='\e[0m'
+
+package_dir=$(pwd)
+
+dirs=(
+  "../../js-api/"
+  "../../libraries/utils/"
+  "../../libraries/ml/"
+  "../../libraries/bio/"
+)
+
+for dir in ${dirs[@]}; do
+  cd $package_dir
+  cd $dir 
+  echo -e $GREEN npm install in $(pwd) $NO_COLOR
+  npm install 
+  echo -e $GREEN npm link in $(pwd) $NO_COLOR
+  npm link
+done
+
+for dir in ${dirs[@]}; do
+  cd $package_dir
+  cd $dir 
+  if [ $dir != "../../js-api/" ]; then
+    echo -e $GREEN npm link-all in $(pwd) $NO_COLOR
+    npm run link-all
+  fi
+  echo -e $GREEN npm run build in$(pwd) $NO_COLOR
+  npm run build || exit
+done
+
+cd $package_dir
+npm run link-all

package/src/model.ts

@@ -1,6 +1,7 @@
 import * as ui from 'datagrok-api/ui';
 import * as grok from 'datagrok-api/grok';
 import * as DG from 'datagrok-api/dg';
+import * as bio from '@datagrok-libraries/bio';
 
 import {splitAlignedSequences} from '@datagrok-libraries/bio/src/utils/splitter';
 
@@ -36,9 +37,9 @@ export class PeptidesModel {
   edf: DG.DataFrame | null = null;
   monomerPositionStatsDf!: DG.DataFrame;
   clusterStatsDf!: DG.DataFrame;
-  _mutationCliffsSelection: type.PositionToAARList = {};
-  _invariantMapSelection: type.PositionToAARList = {};
-  _logoSummarySelection: number[] = [];
+  _mutationCliffsSelection!: type.PositionToAARList;
+  _invariantMapSelection!: type.PositionToAARList;
+  _logoSummarySelection!: number[];
   substitutionsInfo: type.SubstitutionsInfo = new Map();
   isInitialized = false;
   currentView!: DG.TableView;
@@ -49,12 +50,12 @@ export class PeptidesModel {
   mutationCliffsViewer!: MutationCliffsViewer;
   mostPotentResiduesViewer!: MostPotentResiduesViewer;
 
-  _usedProperties: {[propName: string]: string | number | boolean} = {};
-  monomerMap: {[key: string]: {molfile: string, fullName: string}} = {};
+  _usedProperties: { [propName: string]: string | number | boolean } = {};
+  monomerMap: { [key: string]: { molfile: string, fullName: string } } = {};
   barData: type.MonomerDfStats = {};
-  barsBounds: {[position: string]: type.BarCoordinates} = {};
-  cachedBarchartTooltip: {bar: string, tooltip: null | HTMLDivElement} = {bar: '', tooltip: null};
-  monomerLib: any;
+  barsBounds: { [position: string]: type.BarCoordinates } = {};
+  cachedBarchartTooltip: { bar: string, tooltip: null | HTMLDivElement } = {bar: '', tooltip: null};
+  monomerLib: bio.IMonomerLib | null = null;
 
   private constructor(dataFrame: DG.DataFrame) {
     this.df = dataFrame;
@@ -82,6 +83,7 @@ export class PeptidesModel {
     this._mutationCliffsSelection ??= JSON.parse(this.df.tags[C.TAGS.SELECTION] || '{}');
     return this._mutationCliffsSelection;
   }
+
   set mutationCliffsSelection(selection: type.PositionToAARList) {
     this._mutationCliffsSelection = selection;
     this.df.tags[C.TAGS.SELECTION] = JSON.stringify(selection);
@@ -93,6 +95,7 @@ export class PeptidesModel {
     this._invariantMapSelection ??= JSON.parse(this.df.tags[C.TAGS.FILTER] || '{}');
     return this._invariantMapSelection;
   }
+
   set invariantMapSelection(selection: type.PositionToAARList) {
     this._invariantMapSelection = selection;
     this.df.tags[C.TAGS.FILTER] = JSON.stringify(selection);
@@ -104,6 +107,7 @@ export class PeptidesModel {
     this._logoSummarySelection ??= JSON.parse(this.df.tags[C.TAGS.CLUSTER_SELECTION] || '[]');
     return this._logoSummarySelection;
   }
+
   set logoSummarySelection(selection: number[]) {
     this._logoSummarySelection = selection;
     this.df.tags[C.TAGS.CLUSTER_SELECTION] = JSON.stringify(selection);
@@ -111,11 +115,12 @@ export class PeptidesModel {
     this.invalidateGrids();
   }
 
-  get usedProperties(): {[propName: string]: string | number | boolean} {
+  get usedProperties(): { [propName: string]: string | number | boolean } {
     this._usedProperties = JSON.parse(this.df.tags['sarProperties'] ?? '{}');
     return this._usedProperties;
   }
-  set usedProperties(properties: {[propName: string]: string | number | boolean}) {
+
+  set usedProperties(properties: { [propName: string]: string | number | boolean }) {
     this.df.tags['sarProperties'] = JSON.stringify(properties);
     this._usedProperties = properties;
   }
@@ -124,6 +129,7 @@ export class PeptidesModel {
     const splitByPosFlag = (this.df.tags['distributionSplit'] ?? '00')[0];
     return splitByPosFlag == '1' ? true : false;
   }
+
   set splitByPos(flag: boolean) {
     const splitByAARFlag = (this.df.tags['distributionSplit'] ?? '00')[1];
     this.df.tags['distributionSplit'] = `${flag ? 1 : 0}${splitByAARFlag}`;
@@ -133,6 +139,7 @@ export class PeptidesModel {
     const splitByPosFlag = (this.df.tags['distributionSplit'] ?? '00')[1];
     return splitByPosFlag == '1' ? true : false;
   }
+
   set splitByAAR(flag: boolean) {
     const splitByAARFlag = (this.df.tags['distributionSplit'] ?? '00')[0];
     this.df.tags['distributionSplit'] = `${splitByAARFlag}${flag ? 1 : 0}`;
@@ -141,10 +148,22 @@ export class PeptidesModel {
   get isInvariantMap(): boolean {
     return this.df.getTag('isInvariantMap') === '1';
   }
+
   set isInvariantMap(x: boolean) {
     this.df.setTag('isInvariantMap', x ? '1' : '0');
   }
 
+  get isMutationCliffSelectionEmpty(): boolean {
+    for (const aarList of Object.values(this.mutationCliffsSelection))
+      if (aarList.length !== 0)
+        return false;
+    return true;
+  }
+
+  get isLogoSummarySelectionEmpty(): boolean {
+    return this.logoSummarySelection.length === 0;
+  }
+
   createAccordion(): DG.Accordion {
     const acc = ui.accordion();
     acc.root.style.width = '100%';
@@ -295,7 +314,7 @@ export class PeptidesModel {
           continue;
 
         let substCounterFlag = false;
-        const tempData: {pos: string, seq1monomer: string, seq2monomer: string, seq1Idx: number, seq2Idx: number}[] =
+        const tempData: { pos: string, seq1monomer: string, seq2monomer: string, seq1Idx: number, seq2Idx: number }[] =
           [];
         for (const currentPosCol of columnList) {
           const seq1monomer = currentPosCol.get(seq1Idx)!;
@@ -375,14 +394,15 @@ export class PeptidesModel {
   joinDataFrames(positionColumns: string[], splitSeqDf: DG.DataFrame, alphabet: string): void {
     // append splitSeqDf columns to source table and make sure columns are not added more than once
     const name = this.df.name;
+    const cols = this.df.columns;
     for (const colName of positionColumns) {
       const col = this.df.col(colName);
       const newCol = splitSeqDf.getCol(colName);
       if (col === null)
-        this.df.columns.add(newCol);
+        cols.add(newCol);
       else {
-        this.df.columns.remove(colName);
-        this.df.columns.add(newCol);
+        cols.remove(colName);
+        cols.add(newCol);
       }
       CR.setAARRenderer(newCol, alphabet, this.sourceGrid);
     }
@@ -395,7 +415,7 @@ export class PeptidesModel {
     for (let i = 1; i < this.sourceGrid.columns.length; i++)
       colNames.push(this.sourceGrid.columns.byIndex(i)!);
 
-    colNames.sort((a, b)=>{
+    colNames.sort((a, b) => {
       if (a.column!.semType == C.SEM_TYPES.MONOMER) {
         if (b.column!.semType == C.SEM_TYPES.MONOMER)
           return 0;
@@ -409,21 +429,15 @@ export class PeptidesModel {
   }
 
   createScaledCol(activityScaling: string, splitSeqDf: DG.DataFrame): void {
-    const [scaledDf, _, newColName] =
-      scaleActivity(activityScaling, this.df.getCol(C.COLUMNS_NAMES.ACTIVITY));
+    const scaledCol = scaleActivity(activityScaling, this.df.getCol(C.COLUMNS_NAMES.ACTIVITY));
     //TODO: make another func
-    const scaledCol = scaledDf.getCol(C.COLUMNS_NAMES.ACTIVITY_SCALED);
-    scaledCol.semType = C.SEM_TYPES.ACTIVITY_SCALED;
     splitSeqDf.columns.add(scaledCol);
-    const oldScaledCol = this.df.getCol(C.COLUMNS_NAMES.ACTIVITY_SCALED);
-    this.df.columns.replace(oldScaledCol, scaledCol);
+    this.df.columns.replace(C.COLUMNS_NAMES.ACTIVITY_SCALED, scaledCol);
     const gridCol = this.sourceGrid.col(C.COLUMNS_NAMES.ACTIVITY_SCALED);
-    if (gridCol !== null) {
-      gridCol.name = newColName;
-      this.df.tags[C.COLUMNS_NAMES.ACTIVITY_SCALED] = newColName;
-    }
+    if (gridCol)
+      gridCol.name = scaledCol.getTag('gridName');
 
-    this.sourceGrid.columns.setOrder([newColName]);
+    this.sourceGrid.columns.setOrder([scaledCol.getTag('gridName')]);
   }
 
   calculateMonomerPositionStatistics(matrixDf: DG.DataFrame): DG.DataFrame {
@@ -431,7 +445,7 @@ export class PeptidesModel {
 
     //calculate p-values based on t-test
     const matrixCols = matrixDf.columns;
-    const mdCol= matrixCols.addNewFloat(C.COLUMNS_NAMES.MEAN_DIFFERENCE);
+    const mdCol = matrixCols.addNewFloat(C.COLUMNS_NAMES.MEAN_DIFFERENCE);
     const pValCol = matrixCols.addNewFloat(C.COLUMNS_NAMES.P_VALUE);
     const countCol = matrixCols.addNewInt(C.COLUMNS_NAMES.COUNT);
     const ratioCol = matrixCols.addNewFloat(C.COLUMNS_NAMES.RATIO);
@@ -462,7 +476,7 @@ export class PeptidesModel {
     const statsDf = this.df.groupBy([C.COLUMNS_NAMES.CLUSTERS]).aggregate();
     const clustersCol = statsDf.getCol(C.COLUMNS_NAMES.CLUSTERS);
     const statsDfCols = statsDf.columns;
-    const mdCol= statsDfCols.addNewFloat(C.COLUMNS_NAMES.MEAN_DIFFERENCE);
+    const mdCol = statsDfCols.addNewFloat(C.COLUMNS_NAMES.MEAN_DIFFERENCE);
     const pValCol = statsDfCols.addNewFloat(C.COLUMNS_NAMES.P_VALUE);
     const countCol = statsDfCols.addNewInt(C.COLUMNS_NAMES.COUNT);
     const ratioCol = statsDfCols.addNewFloat(C.COLUMNS_NAMES.RATIO);
@@ -512,7 +526,7 @@ export class PeptidesModel {
       .aggregate();
 
     let tempStats: DG.Stats;
-    const maxAtPos: {[index: string]: number} = {};
+    const maxAtPos: { [index: string]: number } = {};
     const posColCategories = sequenceDf.getCol(C.COLUMNS_NAMES.POSITION).categories;
     const mdCol = sequenceDf.getCol(C.COLUMNS_NAMES.MEAN_DIFFERENCE);
     const posCol = sequenceDf.getCol(C.COLUMNS_NAMES.POSITION);
@@ -660,7 +674,7 @@ export class PeptidesModel {
       .subscribe((mouseMove: MouseEvent) => eventAction(mouseMove));
   }
 
-  findAARandPosition(cell: DG.GridCell, ev: MouseEvent): {monomer: string, position: string} | null {
+  findAARandPosition(cell: DG.GridCell, ev: MouseEvent): { monomer: string, position: string } | null {
     const barCoords = this.barsBounds[cell.tableColumn!.name];
     for (const [monomer, coords] of Object.entries(barCoords)) {
       const isIntersectingX = ev.offsetX >= coords.x && ev.offsetX <= coords.x + coords.width;
@@ -672,7 +686,7 @@ export class PeptidesModel {
     return null;
   }
 
-  requestBarchartAction(ev: MouseEvent, barPart: {position: string, monomer: string} | null): void {
+  requestBarchartAction(ev: MouseEvent, barPart: { position: string, monomer: string } | null): void {
     if (!barPart)
       return;
     const monomer = barPart.monomer;
@@ -798,7 +812,9 @@ export class PeptidesModel {
 
   showMonomerTooltip(aar: string, x: number, y: number): void {
     const tooltipElements: HTMLDivElement[] = [];
-    const monomer: type.HELMMonomer = this.monomerLib[aar.toLowerCase()];
+    const monomerName = aar.toLowerCase();
+    const monomer: bio.Monomer = this.monomerLib!.get('HELM_AA', monomerName) ??
+      this.monomerLib!.get('HELM_CHEM', monomerName)!;
 
     if (monomer) {
       tooltipElements.push(ui.div(monomer.n));
@@ -964,10 +980,10 @@ export class PeptidesModel {
       const getBitAt = (i: number): boolean => {
         for (const position of positionList) {
           const positionCol: DG.Column<string> = this.df.getCol(position);
-          if (this._mutationCliffsSelection[position].includes(positionCol.get(i)!))
+          if (this.mutationCliffsSelection[position].includes(positionCol.get(i)!))
             return true;
         }
-        if (this._logoSummarySelection.includes(clusterCol?.get(i)!))
+        if (this.logoSummarySelection.includes(clusterCol?.get(i)!))
           return true;
         return false;
       };
@@ -979,14 +995,16 @@ export class PeptidesModel {
     selection.onChanged.subscribe(() => changeSelectionBitset(selection));
     filter.onChanged.subscribe(() => {
       const positionList = Object.keys(this.invariantMapSelection);
-      filter.init((i) => {
+      for (let index = 0; index < this.df.rowCount; ++index) {
         let result = true;
         for (const position of positionList) {
-          const aarList = this._invariantMapSelection[position];
-          result &&= aarList.length == 0 || aarList.includes(this.df.get(position, i));
+          const aarList = this.invariantMapSelection[position];
+          result &&= aarList.length === 0 || aarList.includes(this.df.get(position, index));
+          if (!result)
+            break;
         }
-        return result;
-      }, false);
+        filter.set(index, filter.get(index) && result, false);
+      }
     });
     this.isBitsetChangedInitialized = true;
   }
@@ -1049,7 +1067,7 @@ export class PeptidesModel {
       const [sourceViewer, targetViewer] = isSourceSAR ? [this.mutationCliffsViewer, this.mostPotentResiduesViewer] :
         [this.mostPotentResiduesViewer, this.mutationCliffsViewer];
       const properties = sourceViewer.props.getProperties();
-      const newProps: {[propName: string]: string | number | boolean} = {};
+      const newProps: { [propName: string]: string | number | boolean } = {};
       for (const property of properties) {
         const propName = property.name;
         const propVal = property.get(sourceViewer);
@@ -1066,8 +1084,11 @@ export class PeptidesModel {
     if (this.isInitialized)
       return;
 
-    this.monomerLib =
-      JSON.parse(await grok.functions.call('Helm:getMonomerLib', {type: this.df.getTag('monomerType')}) as string);
+    // Get monomer library through bio library
+    this.monomerLib = (await bio.getMonomerLib()) as bio.IMonomerLib;
+    this.monomerLib.onChanged.subscribe(() => {
+      this.sourceGrid.invalidate();
+    });
 
     this.currentView = this.df.tags[C.PEPTIDES_ANALYSIS] == 'true' ? grok.shell.v as DG.TableView :
       grok.shell.addTableView(this.df);
@@ -1077,9 +1098,9 @@ export class PeptidesModel {
 
     this.df.tags[C.PEPTIDES_ANALYSIS] = 'true';
     const scaledGridCol = this.sourceGrid.col(C.COLUMNS_NAMES.ACTIVITY_SCALED)!;
-    scaledGridCol.name = this.df.tags[C.COLUMNS_NAMES.ACTIVITY_SCALED];
+    scaledGridCol.name = scaledGridCol.column!.getTag('gridName');
     scaledGridCol.format = '#.000';
-    this.sourceGrid.columns.setOrder([this.df.tags[C.COLUMNS_NAMES.ACTIVITY_SCALED]]);
+    this.sourceGrid.columns.setOrder([scaledGridCol.name]);
     this.sourceGrid.props.allowColSelection = false;
 
     this.df.temp[C.EMBEDDING_STATUS] = false;

package/src/package-test.ts

@@ -1,5 +1,5 @@
 import * as DG from 'datagrok-api/dg';
-import {runTests, tests} from '@datagrok-libraries/utils/src/test';
+import {runTests, tests, TestContext} from '@datagrok-libraries/utils/src/test';
 
 import './tests/core';
 import './tests/peptide-space-test';
@@ -10,8 +10,9 @@ export {tests};
 //name: test
 //input: string category {optional: true}
 //input: string test {optional: true}
+//input: object testContext {optional: true}
 //output: dataframe result
-export async function test(category: string, test: string): Promise<DG.DataFrame> {
-  const data = await runTests({category, test: test});
+export async function test(category: string, test: string, testContext: TestContext): Promise<DG.DataFrame> {
+  const data = await runTests({category, test, testContext});
   return DG.DataFrame.fromObjects(data)!;
 }

package/src/tests/core.ts

@@ -1,5 +1,6 @@
 import * as grok from 'datagrok-api/grok';
 import * as DG from 'datagrok-api/dg';
+import * as bio from '@datagrok-libraries/bio';
 
 import {category, test, expect, delay} from '@datagrok-libraries/utils/src/test';
 
@@ -13,14 +14,14 @@ category('Core', () => {
   let simpleTable: DG.DataFrame;
   let simpleActivityCol: DG.Column<number>;
   let simpleAlignedSeqCol: DG.Column<string>;
-  let simpleScaledDf: DG.DataFrame;
+  let simpleScaledCol: DG.Column<number>;
   let scalingFormula: (x: number) => number;
   let simpleScaledColName: string;
 
   let complexTable: DG.DataFrame;
   let complexActivityCol: DG.Column<number>;
   let complexAlignedSeqCol: DG.Column<string>;
-  let complexScaledDf: DG.DataFrame;
+  let complexScaledCol: DG.Column<number>;
   let complexScaledColName: string;
   const alignedSequenceCol = 'AlignedSequence';
 
@@ -32,13 +33,12 @@ category('Core', () => {
     simpleActivityCol = simpleTable.getCol(simpleActivityColName);
     simpleAlignedSeqCol = simpleTable.getCol(alignedSequenceCol);
     simpleAlignedSeqCol.semType = C.SEM_TYPES.MACROMOLECULE;
-    simpleAlignedSeqCol.tags[C.TAGS.ALPHABET] = 'PT';
-    simpleAlignedSeqCol.tags[DG.TAGS.UNITS] = 'fasta';
-    simpleAlignedSeqCol.tags['aligned'] = 'SEQ.MSA';
-    [simpleScaledDf, scalingFormula, simpleScaledColName] = scaleActivity('-lg', simpleActivityCol);
+    simpleAlignedSeqCol.setTag(C.TAGS.ALPHABET, bio.ALPHABET.PT);
+    simpleAlignedSeqCol.setTag(DG.TAGS.UNITS, bio.NOTATION.FASTA);
+    simpleAlignedSeqCol.setTag(bio.TAGS.aligned, bio.ALIGNMENT.SEQ_MSA);
+    simpleScaledCol = scaleActivity('-lg', simpleActivityCol);
 
-    model = await startAnalysis(
-      simpleActivityCol, simpleAlignedSeqCol, null, simpleTable, scalingFormula, simpleScaledColName, '-lg', []);
+    model = await startAnalysis(simpleActivityCol, simpleAlignedSeqCol, null, simpleTable, simpleScaledCol, '-lg');
     expect(model instanceof PeptidesModel, true);
 
     if (model != null) {
@@ -53,14 +53,14 @@ category('Core', () => {
     complexActivityCol = complexTable.getCol(complexActivityColName);
     complexAlignedSeqCol = complexTable.getCol('MSA');
     complexAlignedSeqCol.semType = C.SEM_TYPES.MACROMOLECULE;
-    complexAlignedSeqCol.tags[C.TAGS.ALPHABET] = 'UN';
-    complexAlignedSeqCol.tags[DG.TAGS.UNITS] = 'separator';
-    complexAlignedSeqCol.tags['aligned'] = 'SEQ.MSA';
+    complexAlignedSeqCol.setTag(C.TAGS.ALPHABET, bio.ALPHABET.UN);
+    complexAlignedSeqCol.setTag(DG.TAGS.UNITS, bio.NOTATION.SEPARATOR);
+    complexAlignedSeqCol.setTag(bio.TAGS.aligned, bio.ALIGNMENT.SEQ_MSA);
     complexAlignedSeqCol.tags[C.TAGS.SEPARATOR] = '/';
-    [complexScaledDf, scalingFormula, complexScaledColName] = scaleActivity('-lg', complexActivityCol);
+    complexScaledCol = scaleActivity('-lg', complexActivityCol);
 
     model = await startAnalysis(
-      complexActivityCol, complexAlignedSeqCol, null, complexTable, scalingFormula, complexScaledColName, '-lg', []);
+      complexActivityCol, complexAlignedSeqCol, null, complexTable, complexScaledCol, '-lg');
     expect(model instanceof PeptidesModel, true);
 
     if (model != null) {
@@ -75,13 +75,12 @@ category('Core', () => {
     simpleActivityCol = simpleTable.getCol(simpleActivityColName);
     simpleAlignedSeqCol = simpleTable.getCol(alignedSequenceCol);
     simpleAlignedSeqCol.semType = C.SEM_TYPES.MACROMOLECULE;
-    simpleAlignedSeqCol.tags[C.TAGS.ALPHABET] = 'PT';
-    simpleAlignedSeqCol.tags[DG.TAGS.UNITS] = 'fasta';
-    simpleAlignedSeqCol.tags['aligned'] = 'SEQ.MSA';
-    [simpleScaledDf, scalingFormula, simpleScaledColName] = scaleActivity('-lg', simpleActivityCol);
+    simpleAlignedSeqCol.setTag(C.TAGS.ALPHABET, bio.ALPHABET.PT);
+    simpleAlignedSeqCol.setTag(DG.TAGS.UNITS, bio.NOTATION.FASTA);
+    simpleAlignedSeqCol.setTag(bio.TAGS.aligned, bio.ALIGNMENT.SEQ_MSA);
+    simpleScaledCol = scaleActivity('-lg', simpleActivityCol);
 
-    model = await startAnalysis(
-      simpleActivityCol, simpleAlignedSeqCol, null, simpleTable, scalingFormula, simpleScaledColName, '-lg', []);
+    model = await startAnalysis(simpleActivityCol, simpleAlignedSeqCol, null, simpleTable, simpleScaledCol, '-lg');
     let v = grok.shell.getTableView('Peptides analysis');
     const d = v.dataFrame;
     const layout = v.saveLayout();

package/src/utils/misc.ts

@@ -18,13 +18,10 @@ export function getSeparator(col: DG.Column<string>): string {
   return col.getTag(C.TAGS.SEPARATOR) ?? '';
 }
 
-export function scaleActivity(activityScaling: string, activityCol: DG.Column<number>, indexes?: number[],
-): [DG.DataFrame, (x: number) => number, string] {
-  const tempDf = DG.DataFrame.create(activityCol.length);
-
+export function scaleActivity(scaling: string, activityCol: DG.Column<number>): DG.Column<number> {
   let formula = (x: number): number => x;
   let newColName = 'activity';
-  switch (activityScaling) {
+  switch (scaling) {
   case 'none':
     break;
   case 'lg':
@@ -36,15 +33,16 @@ export function scaleActivity(activityScaling: string, activityCol: DG.Column<nu
     newColName = `-Log10(${newColName})`;
     break;
   default:
-    throw new Error(`ScalingError: method \`${activityScaling}\` is not available.`);
+    throw new Error(`ScalingError: method \`${scaling}\` is not available.`);
   }
-  tempDf.columns.addNewVirtual(
-    C.COLUMNS_NAMES.ACTIVITY_SCALED, (i) => {
-      const val = activityCol.get(indexes ? indexes[i] : i);
-      return val ? formula(val) : val;
-    }, DG.TYPE.FLOAT);
+  const scaledCol = DG.Column.float(C.COLUMNS_NAMES.ACTIVITY_SCALED, activityCol.length).init((i) => {
+    const val = activityCol.get(i);
+    return val ? formula(val) : val;
+  });
+  scaledCol.semType = C.SEM_TYPES.ACTIVITY_SCALED;
+  scaledCol.setTag('gridName', newColName);
 
-  return [tempDf, formula, newColName];
+  return scaledCol;
 }
 
 export function calculateBarsData(columns: DG.Column<string>[], selection: DG.BitSet): type.MonomerDfStats {

package/src/utils/types.ts

@@ -1,4 +1,5 @@
 import * as DG from 'datagrok-api/dg';
+import * as bio from '@datagrok-libraries/bio';
 
 export type DataFrameDict = {[key: string]: DG.DataFrame};
 
@@ -7,7 +8,7 @@ export type UTypedArray = Uint8Array | Uint16Array | Uint32Array;
 export type SubstitutionsInfo = Map<string, Map<string, Map<number, number[] | UTypedArray>>>;
 export type PositionToAARList = {[postiton: string]: string[]};
 
-export type HELMMonomer = {at: {[key: string]: string}, id: string, m: string, na: string, n: string, rs: number};
+export type HELMMonomer = bio.Monomer;
 
 export type MonomerColStats = {[monomer: string]: {count: number, selected: number}};
 export type MonomerDfStats = {[position: string]: MonomerColStats};

package/src/viewers/sar-viewer.ts

@@ -53,10 +53,10 @@ export class SARViewerBase extends DG.JsViewer {
 
   detach(): void {this.subs.forEach((sub) => sub.unsubscribe());}
 
-  get state(): string {
-    return this.dataFrame.getTag(C.TAGS.SAR_MODE) ?? '10';
+  get isMutationCliffsMode(): string {
+    return this.dataFrame.getTag(C.TAGS.SAR_MODE) ?? '1';
   }
-  set state(s: string) {
+  set isMutationCliffsMode(s: string) {
     this.dataFrame.setTag(C.TAGS.SAR_MODE, s);
   }
 
@@ -67,22 +67,24 @@ export class SARViewerBase extends DG.JsViewer {
       $(this.root).empty();
       let switchHost = ui.div();
       if (this.name == 'MC') {
-        const mutationCliffsMode = ui.boolInput('', this.state[0] === '1', () => {
+        const mutationCliffsMode = ui.boolInput('', this.isMutationCliffsMode === '1', () => {
           if (this.isModeChanging)
             return;
           this.isModeChanging = true;
           invariantMapMode.value = !invariantMapMode.value;
+          this.isMutationCliffsMode = '1';
           this.isModeChanging = false;
           this._titleHost.innerText = 'Mutation Cliffs';
           this.model.isInvariantMap = false;
           this.viewerGrid.invalidate();
         });
         mutationCliffsMode.addPostfix('Mutation Cliffs');
-        const invariantMapMode = ui.boolInput('', this.state[1] === '1', () => {
+        const invariantMapMode = ui.boolInput('', this.isMutationCliffsMode === '0', () => {
           if (this.isModeChanging)
             return;
           this.isModeChanging = true;
           mutationCliffsMode.value = !mutationCliffsMode.value;
+          this.isMutationCliffsMode = '0';
           this.isModeChanging = false;
           this._titleHost.innerText = 'Invariant Map';
           this.model.isInvariantMap = true;

package/src/widgets/distribution.ts

@@ -15,9 +15,6 @@ export function getDistributionWidget(table: DG.DataFrame, model: PeptidesModel)
   const activityScaledCol = table.columns.bySemType(C.SEM_TYPES.ACTIVITY_SCALED)!;
   const rowCount = activityScaledCol.length;
   const selectionObject = model.mutationCliffsSelection;
-  let isMutationCliffsSelectionEmpty = true;
-  for (const aarList of Object.values(selectionObject))
-    isMutationCliffsSelectionEmpty &&= aarList.length === 0;
   const clustersObject = model.logoSummarySelection;
   const positions = Object.keys(selectionObject);
   const positionsLen = positions.length;
@@ -145,16 +142,23 @@ export function getDistributionWidget(table: DG.DataFrame, model: PeptidesModel)
   };
 
   const setDefaultProperties = (input: DG.InputBase): void => {
-    input.enabled = !isMutationCliffsSelectionEmpty;
+    input.enabled = !model.isMutationCliffSelectionEmpty;
     $(input.root).find('.ui-input-editor').css('margin', '0px');
     $(input.root).find('.ui-input-description').css('padding', '0px').css('padding-left', '5px');
   };
 
-  const splitByPosition = ui.boolInput('', model.splitByPos, updateDistributionHost);
+  let defaultValuePos = model.splitByPos;
+  let defaultValueAAR = model.splitByAAR;
+  if (!model.isLogoSummarySelectionEmpty && model.isMutationCliffSelectionEmpty) {
+    defaultValuePos = false;
+    defaultValueAAR = false;
+  }    
+
+  const splitByPosition = ui.boolInput('', defaultValuePos, updateDistributionHost);
   splitByPosition.addPostfix('Split by position');
   setDefaultProperties(splitByPosition);
   $(splitByPosition.root).css('margin-right', '10px');
-  const splitByAAR = ui.boolInput('', model.splitByAAR, updateDistributionHost);
+  const splitByAAR = ui.boolInput('', defaultValueAAR, updateDistributionHost);
   splitByAAR.addPostfix('Split by monomer');
   setDefaultProperties(splitByAAR);