@datagrok/sequence-translator 1.4.7 → 1.4.10

package/package.json

@@ -1,10 +1,10 @@
 {
   "name": "@datagrok/sequence-translator",
   "friendlyName": "Sequence Translator",
-  "version": "1.4.7",
+  "version": "1.4.10",
   "author": {
-    "name": "Alexey Chopovsky",
-    "email": "achopovsky@datagrok.ai"
+    "name": "Leonid Stolbov",
+    "email": "lstolbov@datagrok.ai"
   },
   "description": "SequenceTranslator translates [oligonucleotide](https://en.wikipedia.org/wiki/Oligonucleotide) sequences between [different representations](https://github.com/datagrok-ai/public/tree/master/packages/SequenceTranslator#sequence-representations).",
   "repository": {
@@ -22,10 +22,10 @@
     }
   ],
   "dependencies": {
-    "@datagrok-libraries/bio": "^5.45.5",
+    "@datagrok-libraries/bio": "^5.45.6",
     "@datagrok-libraries/chem-meta": "^1.2.7",
     "@datagrok-libraries/tutorials": "^1.4.3",
-    "@datagrok-libraries/utils": "^4.3.6",
+    "@datagrok-libraries/utils": "^4.3.7",
     "@types/react": "^18.0.15",
     "cash-dom": "^8.1.0",
     "datagrok-api": "^1.21.1",
@@ -39,10 +39,10 @@
     "wu": "^2.1.0"
   },
   "devDependencies": {
-    "@datagrok-libraries/helm-web-editor": "^1.1.12",
-    "@datagrok-libraries/js-draw-lite": "^0.0.9",
-    "@datagrok/bio": "^2.16.6",
-    "@datagrok/helm": "^2.5.7",
+    "@datagrok-libraries/helm-web-editor": "^1.1.13",
+    "@datagrok-libraries/js-draw-lite": "^0.0.10",
+    "@datagrok/bio": "^2.16.8",
+    "@datagrok/helm": "^2.5.8",
     "@datagrok/chem": "^1.12.4",
     "@types/jquery": "^3.5.14",
     "@types/js-yaml": "^4.0.5",

package/src/apps/common/model/oligo-toolkit-package.ts

@@ -5,6 +5,8 @@ import * as DG from 'datagrok-api/dg';
 import {getMonomerLibHelper, IMonomerLibHelper} from '@datagrok-libraries/bio/src/monomer-works/monomer-utils';
 import {IMonomerLib, Monomer} from '@datagrok-libraries/bio/src/types';
 import {LoggerWrapper} from '@datagrok-libraries/bio/src/utils/logger';
+import {ISeqHelper} from '@datagrok-libraries/bio/src/utils/seq-helper';
+import {IHelmHelper} from '@datagrok-libraries/bio/src/helm/helm-helper';
 
 import {APP_NAME} from '../view/const';
 import {DEFAULT_LIB_FILENAME, FALLBACK_LIB_PATH} from './data-loader/const';
@@ -16,14 +18,18 @@ import {MonomerLibWrapper} from './monomer-lib/lib-wrapper';
 import {FormatConverter} from '../../translator/model/format-converter';
 import {FormatDetector} from './parsing-validation/format-detector';
 import {highlightInvalidSubsequence} from '../view/components/colored-input/input-painters';
-import {ISeqHelper} from '@datagrok-libraries/bio/src/utils/seq-helper';
 
 export class OligoToolkitPackage extends DG.Package implements ITranslationHelper {
-  private _seqHelper: ISeqHelper;
+
+  private _helmHelper: IHelmHelper;
+  public get helmHelper(): IHelmHelper {
+    if (!this._helmHelper)
+      throw new Error('Package SequenceTranslator .helmHelper is not initialized');
+    return this._helmHelper;
+  }
+
   public get seqHelper(): ISeqHelper {
-    if (!this._seqHelper)
-      throw new Error('Package SequenceTranslator .seqHelper is not initialized');
-    return this._seqHelper;
+    return this._helmHelper.seqHelper;
   }
 
   private _monomerLib?: IMonomerLib;
@@ -60,8 +66,8 @@ export class OligoToolkitPackage extends DG.Package implements ITranslationHelpe
     this._initPromise = initPromise;
   }
 
-  completeInit(seqHelper: ISeqHelper): void {
-    this._seqHelper = seqHelper;
+  completeInit(helmHelper: IHelmHelper): void {
+    this._helmHelper = helmHelper;
   }
 
   private initLibDataPromise?: Promise<void>;

package/src/package.ts

@@ -18,7 +18,7 @@ import {defaultErrorHandler} from './utils/err-info';
 //polytool specific
 import {polyToolConvert, polyToolConvertUI} from './polytool/pt-dialog';
 import {polyToolEnumerateChemUI} from './polytool/pt-dialog';
-import {polyToolEnumerateHelmUI} from './polytool/pt-enumeration-helm-dialog';
+import {polyToolEnumerateHelmUI} from './polytool/pt-enumerate-seq-dialog';
 import {_setPeptideColumn} from './polytool/utils';
 import {PolyToolCsvLibHandler} from './polytool/csv-to-json-monomer-lib-converter';
 import {ITranslationHelper} from './types';
@@ -27,6 +27,7 @@ import {PolyToolConvertFuncEditor} from './polytool/pt-convert-editor';
 import {CyclizedNotationProvider} from './utils/cyclized';
 import {getSeqHelper} from '@datagrok-libraries/bio/src/utils/seq-helper';
 import {PolyToolTags} from './consts';
+import {getHelmHelper} from '@datagrok-libraries/bio/src/helm/helm-helper';
 
 export const _package: OligoToolkitPackage = new OligoToolkitPackage({debug: true}/**/);
 
@@ -34,17 +35,17 @@ let initSequenceTranslatorPromise: Promise<void> | null = null;
 
 //tags: init
 export async function init(): Promise<void> {
-  if (initSequenceTranslatorPromise === null) {
+  if (initSequenceTranslatorPromise === null)
     _package.startInit(initSequenceTranslatorPromise = initSequenceTranslatorInt());
-  }
+
   return initSequenceTranslatorPromise;
 }
 
 async function initSequenceTranslatorInt(): Promise<void> {
-  const [seqHelper] = await Promise.all([
-    getSeqHelper(),
+  const [helmHelper] = await Promise.all([
+    getHelmHelper(),
   ]);
-  _package.completeInit(seqHelper);
+  _package.completeInit(helmHelper);
 }
 
 //name: Oligo Toolkit
@@ -302,5 +303,5 @@ export async function ptEnumeratorChemApp(): Promise<void> {
 export function applyNotationProviderForCyclized(col: DG.Column<string>, separator: string) {
   col.meta.units = NOTATION.CUSTOM;
   col.tags[PolyToolTags.dataRole] = 'template';
-  col.temp[SeqTemps.notationProvider] = new CyclizedNotationProvider(separator, _package.seqHelper);
+  col.temp[SeqTemps.notationProvider] = new CyclizedNotationProvider(separator, _package.helmHelper);
 }

package/src/polytool/const.ts

@@ -33,6 +33,7 @@ export const PT_WARNING_COLUMN = 'No marcomolecule column chosen!';
 export const PT_UI_GET_HELM = 'Get HELM';
 export const PT_UI_ADD_HELM = 'Add HELM column';
 export const PT_UI_USE_CHIRALITY = 'Chirality engine';
+export const PT_UI_HIGHLIGHT_MONOMERS = 'Highlight monomers';
 export const PT_UI_DIALOG_CONVERSION = 'Poly Tool Conversion';
 export const PT_UI_DIALOG_UNRULE = 'Poly Tool Unrule';
 export const PT_UI_DIALOG_ENUMERATION = 'Poly Tool Enumeration';

package/src/polytool/{pt-conversion.ts → conversion/pt-chain.ts}

@@ -1,24 +1,10 @@
-import * as grok from 'datagrok-api/grok';
-import * as ui from 'datagrok-api/ui';
-import * as DG from 'datagrok-api/dg';
-
-import wu from 'wu';
-
+import {Atom, IHelmBio, HelmMol, HelmType, JSDraw2ModuleType, OrgType,} from '@datagrok-libraries/bio/src/helm/types';
+import {getHelmHelper, IHelmHelper} from '@datagrok-libraries/bio/src/helm/helm-helper';
+import {HelmTypes} from '@datagrok-libraries/bio/src/helm/consts';
 import {cleanupHelmSymbol} from '@datagrok-libraries/bio/src/helm/utils';
-import {HelmTypes, PolymerTypes} from '@datagrok-libraries/bio/src/helm/consts';
-import {getMonomerLibHelper} from '@datagrok-libraries/bio/src/monomer-works/monomer-utils';
-import {IMonomerLib, IMonomerLibBase, Monomer, MonomerLibData, RGroup} from '@datagrok-libraries/bio/src/types';
-import {RDModule, RDMol, RDReaction, MolList, RDReactionResult} from '@datagrok-libraries/chem-meta/src/rdkit-api';
-import {HELM_REQUIRED_FIELD as REQ, HELM_OPTIONAL_FIELDS as OPT, HELM_RGROUP_FIELDS, HELM_OPTIONAL_FIELDS} from '@datagrok-libraries/bio/src/utils/const';
-import {getRdKitModule} from '@datagrok-libraries/bio/src/chem/rdkit-module';
-import {errInfo} from '@datagrok-libraries/bio/src/utils/err-info';
-import {HelmAtom, HelmBio, HelmMol, HelmType, JSDraw2ModuleType, OrgType} from '@datagrok-libraries/bio/src/helm/types';
-import {IHelmHelper} from '@datagrok-libraries/bio/src/helm/helm-helper';
-
-import {Rules, RuleLink, RuleReaction} from './pt-rules';
-import {InvalidReactionError, MonomerNotFoundError} from './types';
 
-import {_package} from '../package';
+import {getOuterIdx, getInnerIdx} from './pt-misc';
+import {Rules, RuleLink, RuleReaction} from '../pt-rules';
 
 declare const JSDraw2: JSDraw2ModuleType;
 declare const org: OrgType;
@@ -31,8 +17,9 @@ type Linkage = {
   fR: number,
   sR: number
 }
+export type PtBio = IHelmBio & { i: number, j: number };
 
-type PolyToolBio = HelmBio & { i: number, j: number };
+type PtAtom = Atom<HelmType, PtBio>
 
 export class Chain {
   linkages: Linkage[];
@@ -40,9 +27,9 @@ export class Chain {
   mol: HelmMol;
 
   constructor(
-    monomers: string[][],
-    linkages: Linkage[],
-    mol: HelmMol) {
+    monomers: string[][], linkages: Linkage[], mol: HelmMol,
+    protected helmHelper: IHelmHelper,
+  ) {
     this.linkages = linkages;
     this.monomers = monomers;
     this.mol = mol;
@@ -73,14 +60,14 @@ export class Chain {
       monomers[i] = rawFragments[i].slice(idxStart + 1, idxEnd).split('.').map((s) => cleanupHelmSymbol(s));
       for (let j = 0; j < monomers[i].length; j++) {
         const elem = monomers[i][j];
-        const bio: PolyToolBio = {type: HelmTypes.AA, i: i, j: j, continuousId: counter};
-        const atom = new JSDraw2.Atom<HelmType>(p, elem, bio);
+        const bio: PtBio = {type: HelmTypes.AA, i: i, j: j, continuousId: counter};
+        const atom = new JSDraw2.Atom<HelmType, IHelmBio>(p, elem, bio);
         resMol.addAtom(atom);
 
         if (j !== 0) {
           const atom1 = resMol.atoms[counter - 1];
           const atom2 = resMol.atoms[counter];
-          const bond = new JSDraw2.Bond<HelmType>(atom1, atom2);
+          const bond = new JSDraw2.Bond<HelmType, IHelmBio>(atom1, atom2);
           bond.r1 = 2;
           bond.r2 = 1;
           resMol.addBond(bond);
@@ -116,29 +103,21 @@ export class Chain {
     for (let i = 0; i < linkages.length; i++) {
       const atom1 = resMol.atoms[linkages[i].fMonomer - 1];
       const atom2 = resMol.atoms[linkages[i].sMonomer - 1];
-      const bond = new JSDraw2.Bond<HelmType>(atom1, atom2);
+      const bond = new JSDraw2.Bond<HelmType, IHelmBio>(atom1, atom2);
       bond.r1 = linkages[i].fR;
       bond.r2 = linkages[i].sR;
       resMol.addBond(bond);
     }
 
-    return new Chain(monomers, linkages, resMol);
+    return new Chain(monomers, linkages, resMol, helmHelper);
   }
 
   /** Get macromolecule from harmonized sequence (template) */
   applyRules(rules: Rules): Chain {
-    // Clone this
-    const resMonomers: string[][] = this.monomers.map((mL) => [...mL]);
-    const resLinkages: Linkage[] = [...this.linkages];
-    const resMol: HelmMol = this.mol.clone();
-
-    throw new Error('not implemented');
-
-    const chain = new Chain(resMonomers, resLinkages, resMol);
-    return chain;
+    const newNotation = this.getNotation();
+    return Chain.fromNotation(newNotation, rules, this.helmHelper);
   }
 
-  /** @deprecated Use {@link parseNotation} and {@link applyRules} instead. */
   static fromNotation(sequence: string, rules: Rules, helmHelper: IHelmHelper): Chain {
     const heterodimerCode = rules.heterodimerCode;
     const homodimerCode = rules.homodimerCode;
@@ -230,14 +209,14 @@ export class Chain {
         const ch2 = new Array<string>(monomersCycled.length - allPos2[0]);
         for (let j = 0; j < allPos2[0] - 1; j++) {
           const elem = ch1[j] = monomersCycled[j];
-          const bio: PolyToolBio = {type: HelmTypes.AA, i: i, j: j, continuousId: counter};
-          const atom: HelmAtom = new JSDraw2.Atom<HelmType>(p, elem, bio);
+          const bio: PtBio = {type: HelmTypes.AA, i: i, j: j, continuousId: counter};
+          const atom: PtAtom = new JSDraw2.Atom<HelmType, PtBio>(p, elem, bio);
           resMol.addAtom(atom);
 
           if (j > 0) {
             const atom1 = resMol.atoms[counter - 1];
             const atom2 = resMol.atoms[counter];
-            const bond = new JSDraw2.Bond<HelmType>(atom1, atom2);
+            const bond = new JSDraw2.Bond<HelmType, IHelmBio>(atom1, atom2);
             bond.r1 = 2;
             bond.r2 = 1;
             resMol.addBond(bond);
@@ -247,14 +226,14 @@ export class Chain {
 
         for (let j = allPos2[0]; j < monomersCycled.length; j++) {
           const elem = ch2[j - allPos2[0]] = monomersCycled[j];
-          const bio: PolyToolBio = {type: HelmTypes.AA, i: i, j: j, continuousId: counter};
-          const atom: HelmAtom = new JSDraw2.Atom<HelmType>(p, elem, bio);
+          const bio: PtBio = {type: HelmTypes.AA, i: i, j: j, continuousId: counter};
+          const atom: PtAtom = new JSDraw2.Atom<HelmType, PtBio>(p, elem, bio);
           resMol.addAtom(atom);
 
           if (j > allPos2[0]) {
             const atom1 = resMol.atoms[counter - 1];
             const atom2 = resMol.atoms[counter];
-            const bond = new JSDraw2.Bond<HelmType>(atom1, atom2);
+            const bond = new JSDraw2.Bond<HelmType, IHelmBio>(atom1, atom2);
             bond.r1 = 2;
             bond.r2 = 1;
             resMol.addBond(bond);
@@ -304,14 +283,14 @@ export class Chain {
       } else {
         for (let j = 0; j < monomers[i].length; j++) {
           const elem = monomers[i][j];
-          const bio: PolyToolBio = {type: HelmTypes.AA, i: i, j: j, continuousId: counter};
-          const atom: HelmAtom = new JSDraw2.Atom<HelmType>(p, elem, bio);
+          const bio: PtBio = {type: HelmTypes.AA, i: i, j: j, continuousId: counter};
+          const atom: PtAtom = new JSDraw2.Atom<HelmType, PtBio>(p, elem, bio);
           resMol.addAtom(atom);
 
           if (j > 0) {
             const atom1 = resMol.atoms[counter - 1];
             const atom2 = resMol.atoms[counter];
-            const bond = new JSDraw2.Bond<HelmType>(atom1, atom2);
+            const bond = new JSDraw2.Bond<HelmType, IHelmBio>(atom1, atom2);
             bond.r1 = 2;
             bond.r2 = 1;
             resMol.addBond(bond);
@@ -325,13 +304,13 @@ export class Chain {
     for (const l of linkages) {
       const atom1 = resMol.atoms[l.fMonomer - 1];
       const atom2 = resMol.atoms[l.sMonomer - 1];
-      const bond = new JSDraw2.Bond<HelmType>(atom1, atom2);
+      const bond = new JSDraw2.Bond<HelmType, IHelmBio>(atom1, atom2);
       bond.r1 = l.fR;
       bond.r2 = l.sR;
       resMol.addBond(bond);
     }
 
-    const chain = new Chain(monomersAll, linkages, resMol);
+    const chain = new Chain(monomersAll, linkages, resMol, helmHelper);
     return chain;
   }
 
@@ -382,14 +361,14 @@ export class Chain {
       for (let j = 0; j < mainFragments[i].length; j++) {
         if (!!mainFragments[i][j]) {
           const elem = mainFragments[i][j];
-          const bio: PolyToolBio = {type: HelmTypes.AA, i: i, j: j, continuousId: counter};
-          const atom = new JSDraw2.Atom<HelmType>(p, elem, bio);
+          const bio: PtBio = {type: HelmTypes.AA, i: i, j: j, continuousId: counter};
+          const atom = new JSDraw2.Atom<HelmType, IHelmBio>(p, elem, bio);
           resMol.addAtom(atom);
 
           if (j !== 0) {
             const atom1 = resMol.atoms[counter - 1];
             const atom2 = resMol.atoms[counter];
-            const bond = new JSDraw2.Bond<HelmType>(atom1, atom2);
+            const bond = new JSDraw2.Bond<HelmType, IHelmBio>(atom1, atom2);
             bond.r1 = 2;
             bond.r2 = 1;
             resMol.addBond(bond);
@@ -405,13 +384,13 @@ export class Chain {
     for (let i = 0; i < linkages.length; i++) {
       const atom1 = resMol.atoms[linkages[i].fMonomer - 1];
       const atom2 = resMol.atoms[linkages[i].sMonomer - 1];
-      const bond = new JSDraw2.Bond<HelmType>(atom1, atom2);
+      const bond = new JSDraw2.Bond<HelmType, IHelmBio>(atom1, atom2);
       bond.r1 = linkages[i].fR;
       bond.r2 = linkages[i].sR;
       resMol.addBond(bond);
     }
 
-    const chain = new Chain(mainFragments, linkages, resMol);
+    const chain = new Chain(mainFragments, linkages, resMol, helmHelper);
     return chain;
   }
 
@@ -648,15 +627,17 @@ export class Chain {
     const errors: string[] = [];
 
     const chainsMonomerCount = this.monomers.map((ch) => ch.length).reduce((acc, curr) => acc + curr, 0);
-    if (this.mol.atoms.length !== chainsMonomerCount)
+    if (this.mol.atoms.length !== chainsMonomerCount) {
       errors.push(`The mol atoms count ${this.mol.atoms.length} does not match ` +
         `the total number ${chainsMonomerCount} of chains' monomers.`);
+    }
 
     const internalBondsCount = this.monomers.map((ch) => ch.length - 1).reduce((acc, curr) => acc + curr, 0);
     const chainsBondCount = internalBondsCount + this.linkages.length;
-    if (this.mol.bonds.length !== chainsBondCount)
+    if (this.mol.bonds.length !== chainsBondCount) {
       errors.push(`The mol bonds count ${this.mol.bonds.length} does not match ` +
         `the total number ${chainsBondCount} in- and inter-chain linkages.`);
+    }
 
     let counter: number = 0;
     for (let spIdx = 0; spIdx < this.monomers.length; ++spIdx) {
@@ -677,188 +658,3 @@ export class Chain {
     return errors;
   }
 }
-
-/** The main PolyTool convert engine. Returns list of Helms. Covered with tests. */
-export function doPolyToolConvert(sequences: string[], rules: Rules, helmHelper: IHelmHelper): string[] {
-  const helms = new Array<string>(sequences.length);
-  for (let i = 0; i < sequences.length; i++) {
-    try {
-      if (sequences[i] == null) { helms[i] = ''; } else {
-        const chain = Chain.fromNotation(sequences[i], rules, helmHelper);
-        helms[i] = chain.getHelm();
-      }
-    } catch (err: any) {
-      const [errMsg, errStack] = errInfo(err);
-      _package.logger.error(errMsg, undefined, errStack);
-      helms[i] = '';
-    }
-  }
-  return helms;
-}
-
-function getMonomersMolBlocks(monomer1: Monomer, monomer2: Monomer): [string, string] {
-  const mb1 = monomer1.molfile;
-  let mb2 = monomer2.molfile;
-  const addGroups = monomer1.rgroups.length;
-
-  //mol v2000 monomer
-  const rgpIdx = mb2.indexOf('M  RGP');
-  if (rgpIdx !== -1) {
-    const groupsCountStr = mb2.substring(rgpIdx + 6, rgpIdx + 9);
-    const groupsCount = Number(groupsCountStr);
-
-    for (let i = 0; i < groupsCount; i++) {
-      const start = rgpIdx + 9 + 4 + i * 8;
-      const end = rgpIdx + 9 + 8 + i * 8;
-      const rGroupSpecifier = mb2.substring(start, end);
-      const groupPosition = Number(rGroupSpecifier) + addGroups;
-      const digits = Math.floor(Math.log10(groupPosition) + 1);
-      const newSpecifier = ' '.repeat(4 - digits) + String(groupPosition);
-      mb2 = mb2.substring(0, start) + newSpecifier + mb2.substring(end, mb2.length);
-    }
-  }
-
-  //TODO: same for v3000 monomer
-
-  return [mb1, mb2];
-}
-
-function getSyntheticMolBlock(rdkit: RDModule, reaction: string,
-  mb1: string, mb2: string, monomerName: string): string {
-  let rxn: RDReaction | null = null;
-  let mols: MolList | null = null;
-  let mol1: RDMol | null = null;
-  let mol2: RDMol | null = null;
-  let rctns: RDReactionResult | null = null;
-  let molP: RDMol | null = null;
-  let molBlock = '';
-
-  try {
-    rxn = rdkit.get_rxn(reaction);
-    if (!rxn) throw new InvalidReactionError(reaction);
-    mols = new rdkit.MolList();
-    mol1 = rdkit.get_mol(mb1!);
-    mol2 = rdkit.get_mol(mb2!);
-    mols.append(mol1!);
-    mols.append(mol2!);
-
-    rctns = rxn.run_reactants(mols, 1);
-    //const size = rctns.size();
-    const element = rctns.get(0);
-
-    molP = element.next();
-    molBlock = molP?.get_molblock();//molP?.get_v3Kmolblock();//
-  } catch (err: any) {
-    const [errMsg, _errStack] = errInfo(err);
-    grok.shell.error(`Can not assemble monomer '${monomerName}': ${errMsg}.`);
-    throw err;
-  } finally {
-    rxn?.delete();
-    mols?.delete();
-    mol1?.delete();
-    mol2?.delete();
-    rctns?.delete();
-    molP?.delete();
-  }
-
-  return molBlock;
-}
-
-function getNewGroups(monomer1: Monomer, monomer2: Monomer): RGroup[] {
-  const groups = new Array<RGroup>(monomer1?.rgroups.length! + monomer2?.rgroups.length!);
-  const length1 = monomer1?.rgroups.length!;
-  const length2 = monomer2?.rgroups.length!;
-
-  for (let i = 0; i < length1; i++)
-    groups[i] = monomer1?.rgroups[i]!;
-
-  for (let i = 0; i < length2; i++) {
-    const rGroupSpecifier = monomer2?.rgroups[i]!.label.replace('R', '');
-    const groupPosition = Number(rGroupSpecifier) + length1;
-    const group: RGroup = {
-      //@ts-ignore
-      [HELM_RGROUP_FIELDS.CAP_GROUP_SMILES_UPPERCASE]: monomer2?.rgroups[i].capGroupSMILES.replace(rGroupSpecifier, String(groupPosition)),
-      [HELM_RGROUP_FIELDS.ALTERNATE_ID]: monomer2?.rgroups[i].alternateId.replace(rGroupSpecifier, String(groupPosition)),
-      [HELM_RGROUP_FIELDS.CAP_GROUP_NAME]: monomer2?.rgroups[i].capGroupName,
-      [HELM_RGROUP_FIELDS.LABEL]: monomer2?.rgroups[i].label.replace(rGroupSpecifier, String(groupPosition)),
-    };
-
-    groups[i + length1] = group;
-  }
-
-  return groups;
-}
-
-export function getNewMonomer(rdkit: RDModule, mLib: IMonomerLib, rule: RuleReaction): [string, Monomer] {
-  const reacSmarts = rule.reaction;
-  const monomerName = rule.name;
-
-  const monomer1 = mLib.getMonomer('PEPTIDE', rule.firstMonomer);
-  if (!monomer1) throw new MonomerNotFoundError('PEPTIDE', rule.firstMonomer);
-  const monomer2 = mLib.getMonomer('PEPTIDE', rule.secondMonomer);
-  if (!monomer2) throw new MonomerNotFoundError('PEPTIDE', rule.secondMonomer);
-
-  const [mb1, mb2] = getMonomersMolBlocks(monomer1!, monomer2!);
-  const molBlock = getSyntheticMolBlock(rdkit, reacSmarts, mb1, mb2, monomerName);
-  const groups: RGroup[] = getNewGroups(monomer1!, monomer2!);
-
-  const resMonomer: Monomer = {
-    [REQ.SYMBOL]: monomerName,
-    [REQ.NAME]: monomerName,
-    [REQ.MOLFILE]: molBlock,
-    [REQ.AUTHOR]: '',
-    [REQ.ID]: 0,
-    [REQ.RGROUPS]: groups,
-    [REQ.SMILES]: '',
-    [REQ.POLYMER_TYPE]: 'PEPTIDE',
-    [REQ.MONOMER_TYPE]: 'Backbone',
-    [REQ.CREATE_DATE]: null,
-
-    // // @ts-ignore
-    // lib: {source: 'Reaction'},
-  };
-
-  resMonomer[OPT.META] = Object.assign(resMonomer[OPT.META] ?? {},
-    {'colors': {'default': {line: '#2083D5', text: '#2083D5', background: '#F2F2F5'}}});
-
-  return [monomerName, resMonomer];
-}
-
-export async function getOverriddenLibrary(rules: Rules): Promise<IMonomerLibBase> {
-  const monomerLibHelper = await getMonomerLibHelper();
-  const systemMonomerLib = monomerLibHelper.getMonomerLib();
-
-  const rdkit = await getRdKitModule();
-  const argLib: { [symbol: string]: Monomer } = {};
-
-  for (let i = 0; i < rules.reactionRules.length; i++) {
-    const [name, monomer] = getNewMonomer(rdkit, systemMonomerLib, rules.reactionRules[i]);
-    argLib[name] = monomer;
-  }
-
-  const overrideMonomerLibData: MonomerLibData = {[PolymerTypes.PEPTIDE]: argLib};
-  const overriddenMonomerLib = systemMonomerLib.override(overrideMonomerLibData,
-    'ST-PT-reactions.' + wu.repeat(1).map(() => Math.floor((Math.random() * 36)).toString(36)).take(4).toArray().join(''));
-  return overriddenMonomerLib;
-}
-
-/** Gets 0-based in-index (simple polymer) of out-index (continuous) {@link idx} */
-export function getInnerIdx(outIdx: number, monomers: string[][]): [number, number] {
-  // let prevSpCount = 0;
-  // for (let spI = 0; spI < monomers.length && idx >= (prevSpCount + monomers[spI].length); ++spI)
-  //   prevSpCount += monomers[spI].length;
-  // return idx - prevSpCount;
-  let inIdx = outIdx;
-  let spIdx: number;
-  for (spIdx = 0; spIdx < monomers.length && inIdx >= monomers[spIdx].length; ++spIdx)
-    inIdx -= monomers[spIdx].length;
-  return [inIdx, spIdx];
-}
-
-/** Gets 0-based out-index of 0-based in-index {@link inIdx} monomer of simple polymer {@link spIdx} */
-export function getOuterIdx(inIdx: number, spIdx: number, monomers: string[][]): number {
-  let outIdx = 0;
-  for (let i = 0; i < spIdx; ++i)
-    outIdx += monomers[i].length;
-  return outIdx + inIdx;
-}

package/src/polytool/conversion/pt-conversion.ts

@@ -0,0 +1,27 @@
+
+import {errInfo} from '@datagrok-libraries/bio/src/utils/err-info';
+import {IHelmHelper} from '@datagrok-libraries/bio/src/helm/helm-helper';
+
+import {Rules, RuleReaction} from '../pt-rules';
+import {InvalidReactionError, MonomerNotFoundError} from '../types';
+import {Chain} from './pt-chain';
+
+import {_package} from '../../package';
+
+/** The main PolyTool convert engine. Returns list of Helms. Covered with tests. */
+export function doPolyToolConvert(sequences: string[], rules: Rules, helmHelper: IHelmHelper): string[] {
+  const helms = new Array<string>(sequences.length);
+  for (let i = 0; i < sequences.length; i++) {
+    try {
+      if (sequences[i] == null) { helms[i] = ''; } else {
+        const chain = Chain.fromNotation(sequences[i], rules, helmHelper);
+        helms[i] = chain.getHelm();
+      }
+    } catch (err: any) {
+      const [errMsg, errStack] = errInfo(err);
+      _package.logger.error(errMsg, undefined, errStack);
+      helms[i] = '';
+    }
+  }
+  return helms;
+}