@datagrok/bio 2.25.17 → 2.26.1

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Files changed (50) hide show
  1. package/CHANGELOG.md +4 -0
  2. package/dist/282.js +2 -0
  3. package/dist/282.js.map +1 -0
  4. package/dist/287.js +2 -0
  5. package/dist/287.js.map +1 -0
  6. package/dist/288.js +2 -0
  7. package/dist/288.js.map +1 -0
  8. package/dist/422.js +2 -0
  9. package/dist/422.js.map +1 -0
  10. package/dist/455.js +1 -1
  11. package/dist/455.js.map +1 -1
  12. package/dist/767.js +2 -0
  13. package/dist/767.js.map +1 -0
  14. package/dist/package-test.js +5 -5
  15. package/dist/package-test.js.map +1 -1
  16. package/dist/package.js +3 -3
  17. package/dist/package.js.map +1 -1
  18. package/files/samples/antibodies.csv +494 -0
  19. package/package.json +2 -2
  20. package/src/package-api.ts +28 -0
  21. package/src/package.g.ts +31 -1
  22. package/src/package.ts +40 -1
  23. package/src/tests/substructure-filters-tests.ts +1 -0
  24. package/src/utils/annotations/annotation-actions.ts +130 -0
  25. package/src/utils/annotations/annotation-manager-ui.ts +118 -0
  26. package/src/utils/annotations/annotation-manager.ts +163 -0
  27. package/src/utils/annotations/liability-scanner-ui.ts +88 -0
  28. package/src/utils/annotations/liability-scanner.ts +147 -0
  29. package/src/utils/annotations/numbering-ui.ts +472 -0
  30. package/src/utils/antibody-numbering (WIP)/alignment.ts +578 -0
  31. package/src/utils/antibody-numbering (WIP)/annotator.ts +120 -0
  32. package/src/utils/antibody-numbering (WIP)/data/blosum62.ts +55 -0
  33. package/src/utils/antibody-numbering (WIP)/data/consensus-aho.ts +155 -0
  34. package/src/utils/antibody-numbering (WIP)/data/consensus-imgt.ts +162 -0
  35. package/src/utils/antibody-numbering (WIP)/data/consensus-kabat.ts +157 -0
  36. package/src/utils/antibody-numbering (WIP)/data/consensus-martin.ts +152 -0
  37. package/src/utils/antibody-numbering (WIP)/data/consensus.ts +36 -0
  38. package/src/utils/antibody-numbering (WIP)/data/regions.ts +63 -0
  39. package/src/utils/antibody-numbering (WIP)/index.ts +31 -0
  40. package/src/utils/antibody-numbering (WIP)/testdata.ts +5356 -0
  41. package/src/utils/antibody-numbering (WIP)/types.ts +69 -0
  42. package/src/utils/context-menu.ts +42 -2
  43. package/src/utils/get-region-func-editor.ts +18 -2
  44. package/src/utils/get-region.ts +167 -17
  45. package/src/utils/sequence-column-input.ts +57 -0
  46. package/src/viewers/vd-regions-viewer.ts +2 -0
  47. package/src/widgets/representations.ts +53 -2
  48. package/src/widgets/sequence-scrolling-widget.ts +28 -18
  49. package/test-console-output-1.log +587 -551
  50. package/test-record-1.mp4 +0 -0
package/CHANGELOG.md CHANGED
@@ -1,5 +1,9 @@
1
1
  # Bio changelog
2
2
 
3
+ ## 2.26.0 (2026-03-04)
4
+
5
+ * Antibodies support: Numbering schemes, highlighting, allignment, liabilities detection, extracting regions, and more.
6
+
3
7
  ## 2.25.15 (2026-02-26)
4
8
 
5
9
  * Monomer collections
package/dist/282.js ADDED
@@ -0,0 +1,2 @@
1
+ "use strict";(self.webpackChunkbio=self.webpackChunkbio||[]).push([[282],{1282(e,t,a){a.d(t,{showNumberingSchemeDialog:()=>p});var n,i,o=a(4328),s=a(7389),l=a(6082),r=a(2003),g=a(4517);!function(e){e.IMGT="IMGT",e.Kabat="Kabat",e.Chothia="Chothia",e.AHo="AHo"}(n||(n={})),function(e){e.Heavy="Heavy",e.Light_Kappa="Light_Kappa",e.Light_Lambda="Light_Lambda"}(i||(i={}));i.Heavy,i.Heavy,i.Heavy,i.Heavy,i.Heavy,i.Heavy,i.Heavy,i.Heavy,i.Light_Kappa,i.Light_Kappa,i.Light_Kappa,i.Light_Kappa,i.Light_Kappa,i.Light_Kappa,i.Light_Kappa,i.Light_Kappa,i.Light_Lambda,i.Light_Lambda,i.Light_Lambda,i.Light_Lambda,i.Light_Lambda,i.Light_Lambda,i.Light_Lambda,i.Light_Lambda,i.Heavy,i.Heavy,i.Heavy,i.Heavy,i.Heavy,i.Heavy,i.Heavy,i.Heavy,i.Light_Kappa,i.Light_Kappa,i.Light_Kappa,i.Light_Kappa,i.Light_Kappa,i.Light_Kappa,i.Light_Kappa,i.Light_Kappa,i.Light_Lambda,i.Light_Lambda,i.Light_Lambda,i.Light_Lambda,i.Light_Lambda,i.Light_Lambda,i.Light_Lambda,i.Light_Lambda,i.Heavy,i.Heavy,i.Heavy,i.Heavy,i.Heavy,i.Heavy,i.Heavy,i.Heavy,i.Light_Kappa,i.Light_Kappa,i.Light_Kappa,i.Light_Kappa,i.Light_Kappa,i.Light_Kappa,i.Light_Kappa,i.Light_Kappa,i.Light_Lambda,i.Light_Lambda,i.Light_Lambda,i.Light_Lambda,i.Light_Lambda,i.Light_Lambda,i.Light_Lambda,i.Light_Lambda,i.Heavy,i.Heavy,i.Heavy,i.Heavy,i.Heavy,i.Heavy,i.Heavy,i.Heavy,i.Light_Kappa,i.Light_Kappa,i.Light_Kappa,i.Light_Kappa,i.Light_Kappa,i.Light_Kappa,i.Light_Kappa,i.Light_Kappa,i.Light_Lambda,i.Light_Lambda,i.Light_Lambda,i.Light_Lambda,i.Light_Lambda,i.Light_Lambda,i.Light_Lambda,i.Light_Lambda;n.IMGT,n.Kabat,n.Chothia,n.AHo;var h=a(3736);function p(){const e=o.shell.tv?.dataFrame;if(!e)return void o.shell.warning("No table open");const t=e.columns.bySemTypeAll(l.SEMTYPE.MACROMOLECULE);if(0===t.length)return void o.shell.warning("No macromolecule columns found");const i=function(){const e=[],t=l.Func.find({meta:{role:"antibodyNumbering"}});if(0===t.length)throw o.shell.error("No antibody numbering engines found. Make sure that Proteomics plugin is installed and up to date."),new Error("No external antibody numbering engines found. Make sure that Proteomics plugin is installed and up to date.");for(const a of t){const t=a.package?.name??"",n=a.friendlyName||a.name;e.push({label:n,key:t?`${t}:${a.name}`:a.name,func:a})}return e.push({label:"Built-in (TypeScript)",key:"__builtin__",func:null}),e}(),p=i.map(e=>e.label),m=Object.values(n),L=s.input.table("Table",{value:e}),d=s.input.column("Sequence",{table:e,value:t[0],filter:e=>e.semType===l.SEMTYPE.MACROMOLECULE}),f=s.input.choice("Scheme",{value:n.IMGT,items:m}),_=s.input.choice("Engine",{value:p[0],items:p}),b=s.dialog({title:"Apply Antibody Numbering"}).add(s.inputs([L,d,f,_])).onOK(async()=>{const t=d.value,n=f.value,s=_.value,p=i.find(e=>e.label===s)??i[i.length-1],m=l.TaskBarProgressIndicator.create(`Applying ${n} numbering...`);try{let i;i=p.func?await p.func.apply({df:e,seqCol:t,scheme:n.toLowerCase()}):await async function(e,t){const{numberSequences:n,extractSequence:i}=await a.e(767).then(a.bind(a,6767)),o=t.toLowerCase(),s=[];for(let t=0;t<e.length;t++){const a=e.get(t);s.push(i(a??""))}return function(e){const t=e.length,a=l.Column.fromType(l.COLUMN_TYPE.STRING,"position_names",t),n=l.Column.fromType(l.COLUMN_TYPE.STRING,"chain_type",t),i=l.Column.fromType(l.COLUMN_TYPE.STRING,"annotations_json",t),o=l.Column.fromType(l.COLUMN_TYPE.STRING,"numbering_detail",t),s=l.Column.fromType(l.COLUMN_TYPE.STRING,"numbering_map",t);for(let l=0;l<t;l++){const t=e[l];t.error&&t.percentIdentity<.3?(a.set(l,""),n.set(l,""),i.set(l,"[]"),o.set(l,""),s.set(l,"")):(a.set(l,t.positionNames),n.set(l,t.chainType),i.set(l,JSON.stringify(t.annotations)),o.set(l,JSON.stringify(t.numberingDetail)),s.set(l,JSON.stringify(t.numberingMap)))}return l.DataFrame.fromColumns([a,n,i,o,s])}(n(s,o))}(t,n),function(e,t,a,n,i,s){if(!a||0===a.rowCount)return void o.shell.warning("No numbering results returned");const p=a.getCol("position_names"),m=a.getCol("chain_type"),L=a.getCol("annotations_json"),d=a.col("numbering_map");if(d)for(let e=0;e<a.rowCount;e++){const a=d.get(e);if(a)try{const n=JSON.parse(a),i=u(t.get(e)??""),o={};for(const[e,t]of Object.entries(n))t<i.length&&(o[e]=i[t]);d.set(e,JSON.stringify(o))}catch{}}let f=-1,_=-1;for(let e=0;e<a.rowCount;e++){const t=p.get(e);if(!t||0===t.length)continue;const a=L.get(e);let n=0;if(a)try{n=JSON.parse(a).length}catch{}n>_&&(_=n,f=e)}const b=f>=0?m.get(f)??"":"",y=f>=0?L.get(f)??"[]":"[]";if(f<0)return void o.shell.warning(`${s} could not number the sequences. Check that they are valid antibody variable region sequences.`);t.setTag(r.gp.numberingScheme,n);try{const e=JSON.parse(y),a=(0,h.Ln)(t).filter(e=>e.category!==g.eI.Structure);(0,h.fh)(t,[...a,...e])}catch(e){console.warn("Failed to set annotation definitions on original column:",e)}if(d)try{const n=(0,h.Lz)(e,t);for(let e=0;e<a.rowCount;e++){const t=d.get(e),a=L.get(e);if(!t||!a)continue;const i=JSON.parse(t),o=JSON.parse(a),s=[];for(const e of o){if(null==e.start||null==e.end)continue;const t=i[e.start],a=i[e.end];null!=t&&null!=a&&s.push({annotationId:e.id,positionIndex:t,endPositionIndex:a,positionName:e.start,matchedMonomers:""})}const l=(0,h.JG)(n,e)??[];(0,h.z5)(n,e,(0,h.bq)(l,s,!0,!1))}}catch(e){console.warn("Failed to store per-row region data on original column:",e)}const v=function(e,t,a){const n=a.col("numbering_map");if(!n)return null;const i=new Set,o=[],s=[],g=[];for(let e=0;e<a.rowCount;e++){const a=n.get(e),l=t.get(e)??"";if(a)try{const e=JSON.parse(a);o.push(e);for(const t of Object.keys(e))i.add(t);const t=Object.values(e),n=Math.min(...t),r=Math.max(...t);s.push(n),g.push(Math.max(0,l.length-r-1))}catch{o.push(null),s.push(0),g.push(0)}else o.push(null),s.push(0),g.push(0)}if(0===i.size)return null;const h=Math.max(0,...s),p=Math.max(0,...g),u=Array.from(i).slice().sort((e,t)=>{const[a,n]=c(e),[i,o]=c(t);return a!==i?a-i:n.localeCompare(o)}),m=[];for(let e=h;e>0;e--)m.push(`N-${e}`);const L=[];for(let e=1;e<=p;e++)L.push(`C+${e}`);const d=[...m,...u,...L],f=d.length,_=h,b=new Map;for(let e=0;e<u.length;e++)b.set(u[e],_+e);const y=e.columns.getUnusedName(`${t.name} (aligned)`),v=l.Column.fromType(l.COLUMN_TYPE.STRING,y,e.rowCount);for(let a=0;a<e.rowCount;a++){const e=a<o.length?o[a]:null,n=t.get(a)??"";if(!e){v.set(a,"-".repeat(f));continue}const i=new Array(f).fill("-");for(const[t,a]of Object.entries(e)){const e=b.get(t);null!=e&&a<n.length&&(i[e]=n[a])}const l=s[a],r=Object.values(e),h=Math.min(...r);for(let e=0;e<l;e++)i[_-l+e]=n[h-l+e];const p=g[a],c=Math.max(...r),m=_+u.length;for(let e=0;e<p;e++)i[m+e]=n[c+1+e];v.set(a,i.join(""))}return v.semType=l.SEMTYPE.MACROMOLECULE,v.setTag(r.gp.aligned,"SEQ.MSA"),v.setTag(r.gp.alphabet,r.YI.PT),v.meta.units=r.Hi.FASTA,v.setTag(l.Tags.CellRenderer,"sequence"),v.setTag(r.gp.positionNames,d.join(", ")),{alignedCol:v,unifiedPositions:d,preOffset:_}}(e,t,a);if(v&&(e.columns.insert(v.alignedCol,e.columns.toList().indexOf(t)+1),o.shell.tv?.dataFrame===e&&o.shell.tv.grid.scrollToCell(t,0),v.alignedCol.setTag(r.gp.numberingScheme,n),i))try{const e=JSON.parse(y),t=(0,h.Ln)(v.alignedCol).filter(e=>e.category!==g.eI.Structure);(0,h.fh)(v.alignedCol,[...t,...e])}catch(e){console.warn("Failed to set column-level annotations on aligned column:",e)}e.fireValuesChanged(),o.shell.info(`Numbering applied: ${n}, chain type: ${b}`)}(e,t,i,n,!0,p.label)}catch(e){o.shell.error(`Numbering failed: ${e.message??e}`),console.error(e)}finally{m.close()}});b.show()}function u(e){const t=[];for(let a=0;a<e.length;a++)"-"!==e[a]&&"."!==e[a]&&t.push(a);return t}function c(e){const t=e.match(/^(\d+)([A-Z]?)$/);return t?[parseInt(t[1],10),t[2]]:[1/0,e]}}}]);
2
+ //# sourceMappingURL=282.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"282.js","mappings":"mIAMWA,EAOAC,E,mDANX,SAAWD,GACPA,EAAsB,KAAI,OAC1BA,EAAuB,MAAI,QAC3BA,EAAyB,QAAI,UAC7BA,EAAqB,IAAI,KAC5B,CALD,CAKGA,IAAoBA,EAAkB,CAAC,IAE1C,SAAWC,GACPA,EAAiB,MAAI,QACrBA,EAAuB,YAAI,cAC3BA,EAAwB,aAAI,cAC/B,CAJD,CAIGA,IAAcA,EAAY,CAAC,IAKzBA,EAAUC,MACsED,EAAUC,MACPD,EAAUC,MACZD,EAAUC,MACRD,EAAUC,MACXD,EAAUC,MACPD,EAAUC,MACZD,EAAUC,MAE7FD,EAAUE,YACsEF,EAAUE,YACPF,EAAUE,YACZF,EAAUE,YACRF,EAAUE,YACXF,EAAUE,YACPF,EAAUE,YACZF,EAAUE,YAE7FF,EAAUG,aACsEH,EAAUG,aACPH,EAAUG,aACZH,EAAUG,aACRH,EAAUG,aACXH,EAAUG,aACPH,EAAUG,aACZH,EAAUG,aAK7FH,EAAUC,MACsED,EAAUC,MACPD,EAAUC,MACZD,EAAUC,MACRD,EAAUC,MACZD,EAAUC,MACPD,EAAUC,MACXD,EAAUC,MAE7FD,EAAUE,YACsEF,EAAUE,YACPF,EAAUE,YACZF,EAAUE,YACRF,EAAUE,YACZF,EAAUE,YACRF,EAAUE,YACXF,EAAUE,YAE5FF,EAAUG,aACsEH,EAAUG,aACPH,EAAUG,aACZH,EAAUG,aACRH,EAAUG,aACZH,EAAUG,aACRH,EAAUG,aACXH,EAAUG,aAK5FH,EAAUC,MACsED,EAAUC,MACPD,EAAUC,MACZD,EAAUC,MACRD,EAAUC,MACZD,EAAUC,MACPD,EAAUC,MACXD,EAAUC,MAE7FD,EAAUE,YACsEF,EAAUE,YACPF,EAAUE,YACZF,EAAUE,YACRF,EAAUE,YACZF,EAAUE,YACRF,EAAUE,YACXF,EAAUE,YAE5FF,EAAUG,aACsEH,EAAUG,aACPH,EAAUG,aACZH,EAAUG,aACRH,EAAUG,aACZH,EAAUG,aACRH,EAAUG,aACXH,EAAUG,aAK5FH,EAAUC,MACsED,EAAUC,MACPD,EAAUC,MACZD,EAAUC,MACRD,EAAUC,MACXD,EAAUC,MACPD,EAAUC,MACZD,EAAUC,MAE7FD,EAAUE,YACsEF,EAAUE,YACPF,EAAUE,YACZF,EAAUE,YACRF,EAAUE,YACXF,EAAUE,YACPF,EAAUE,YACZF,EAAUE,YAE7FF,EAAUG,aACsEH,EAAUG,aACPH,EAAUG,aACZH,EAAUG,aACRH,EAAUG,aACXH,EAAUG,aACPH,EAAUG,aACZH,EAAUG,aAK7FJ,EAAgBK,KAChBL,EAAgBM,MAChBN,EAAgBO,QAChBP,EAAgBQ,I,cC1Cd,SAASC,IACd,MAAMC,EAAK,QAAWC,IAAIC,UAC1B,IAAKF,EAEH,YADA,QAAWG,QAAQ,iBAIrB,MAAMC,EAAUJ,EAAGK,QAAQC,aAAa,UAAWC,eACnD,GAAuB,IAAnBH,EAAQI,OAEV,YADA,QAAWL,QAAQ,kCAIrB,MAAMM,EAnFR,WACE,MAAMA,EAA6B,GAE7BC,EAAQ,OAAQC,KAAK,CAACC,KAAM,CAACC,KAAM,uBACzC,GAAqB,IAAjBH,EAAMF,OAER,MADA,QAAWM,MAAM,sGACX,IAAIC,MAAM,+GAElB,IAAK,MAAMC,KAAKN,EAAO,CACrB,MAAMO,EAAUD,EAAEE,SAASC,MAAQ,GAC7BC,EAAQJ,EAAEK,cAAgBL,EAAEG,KAClCV,EAAQa,KAAK,CACXF,MAAOA,EACPG,IAAKN,EAAU,GAAGA,KAAWD,EAAEG,OAASH,EAAEG,KAC1CK,KAAMR,GAEV,CAIA,OADAP,EAAQa,KAAK,CAACF,MAjCa,wBAiCgBG,IAlClB,cAkC2CC,KAAM,OACnEf,CACT,CA8DkBgB,GACVC,EAAejB,EAAQkB,IAAKC,GAAMA,EAAER,OACpCS,EAAgBC,OAAOC,OAAOzC,GAE9B0C,EAAa,QAASC,MAAM,QAAS,CAACC,MAAOlC,IAC7CmC,EAAW,QAASC,OAAO,WAAY,CAC3CH,MAAOjC,EAAIkC,MAAO9B,EAAQ,GAC1BiC,OAASC,GAAmBA,EAAIC,UAAY,UAAWhC,gBAEnDiC,EAAc,QAASC,OAAO,SAAU,CAACP,MAAO5C,EAAgBK,KAAM+C,MAAOb,IAC7Ec,EAAc,QAASF,OAAO,SAAU,CAC5CP,MAAOR,EAAa,GAAIgB,MAAOhB,IAK3BkB,EAAS,SAAU,CAACC,MAAO,6BAC9BC,IAAI,SAAU,CAACd,EAAYG,EAAUK,EAAaG,KAClDI,KAAKC,UACJ,MAAMC,EAASd,EAASD,MAClBgB,EAAaV,EAAYN,MACzBiB,EAAgBR,EAAYT,MAC5BkB,EAAS3C,EAAQE,KAAMiB,GAAMA,EAAER,QAAU+B,IAAkB1C,EAAQA,EAAQD,OAAS,GACpF6C,EAAK,2BAA4BC,OAAO,YAAYJ,kBAC1D,IACE,IAAIK,EAEFA,EADEH,EAAO5B,WACM4B,EAAO5B,KAAKgC,MAAM,CAACxD,GAAIA,EAAIiD,OAAQA,EAAQQ,OAAQP,EAAWQ,sBAxDvFV,eACEC,EAA2BC,GAE3B,MAAM,gBAACS,EAAe,gBAAEC,SAAyB,8BAC3CH,EAASP,EAAWQ,cAEpBG,EAAsB,GAC5B,IAAK,IAAIC,EAAI,EAAGA,EAAIb,EAAOzC,OAAQsD,IAAK,CACtC,MAAMC,EAAMd,EAAOe,IAAIF,GACvBD,EAAUvC,KAAKsC,EAAgBG,GAAO,IACxC,CAGA,OA1CK,SAAqCE,GAC1C,MAAMC,EAAID,EAAQzD,OACZ2D,EAAW,SAAUC,SAAS,cAAeC,OAAQ,iBAAkBH,GACvEI,EAAa,SAAUF,SAAS,cAAeC,OAAQ,aAAcH,GACrEK,EAAY,SAAUH,SAAS,cAAeC,OAAQ,mBAAoBH,GAC1EM,EAAY,SAAUJ,SAAS,cAAeC,OAAQ,mBAAoBH,GAC1EO,EAAS,SAAUL,SAAS,cAAeC,OAAQ,gBAAiBH,GAE1E,IAAK,IAAIJ,EAAI,EAAGA,EAAII,EAAGJ,IAAK,CAC1B,MAAMY,EAAIT,EAAQH,GACdY,EAAE5D,OAAS4D,EAAEC,gBAAkB,IACjCR,EAASS,IAAId,EAAG,IAChBQ,EAAWM,IAAId,EAAG,IAClBS,EAAUK,IAAId,EAAG,MACjBU,EAAUI,IAAId,EAAG,IACjBW,EAAOG,IAAId,EAAG,MAEdK,EAASS,IAAId,EAAGY,EAAEG,eAClBP,EAAWM,IAAId,EAAGY,EAAEI,WACpBP,EAAUK,IAAId,EAAGiB,KAAKC,UAAUN,EAAEO,cAClCT,EAAUI,IAAId,EAAGiB,KAAKC,UAAUN,EAAEQ,kBAClCT,EAAOG,IAAId,EAAGiB,KAAKC,UAAUN,EAAES,eAEnC,CAEA,OAAO,YAAaC,YAAY,CAACjB,EAAUG,EAAYC,EAAWC,EAAWC,GAC/E,CAgBSY,CADS1B,EAAgBE,EAAWJ,GAE7C,CA4CyB6B,CAAoBrC,EAAQC,GAiLrD,SACElD,EAAkBiD,EAA2BM,EAC7CL,EAAoBqC,EAA0BC,GAE9C,IAAKjC,GAA8B,IAApBA,EAAOkC,SAEpB,YADA,QAAWtF,QAAQ,iCAIrB,MAAMuF,EAAcnC,EAAOoC,OAAO,kBAC5BC,EAAerC,EAAOoC,OAAO,cAC7BE,EAAetC,EAAOoC,OAAO,oBAC7BG,EAAkBvC,EAAOjB,IAAI,iBAKnC,GAAiBwD,EACf,IAAK,IAAIhC,EAAI,EAAGA,EAAIP,EAAOkC,SAAU3B,IAAK,CACxC,MAAMiC,EAASD,EAAgB9B,IAAIF,GACnC,GAAKiC,EACL,IACE,MAAMC,EAAuCjB,KAAKkB,MAAMF,GAElDG,EAAaC,EADJlD,EAAOe,IAAIF,IAAM,IAE1BsC,EAAmC,CAAC,EAC1C,IAAK,MAAOC,EAASC,KAAgBxE,OAAOyE,QAAQP,GAC9CM,EAAcJ,EAAW1F,SAC3B4F,EAASC,GAAWH,EAAWI,IAEnCR,EAAgBlB,IAAId,EAAGiB,KAAKC,UAAUoB,GACxC,CAAE,MAAmB,CACvB,CAIF,IAAII,GAAc,EACdC,GAAkB,EAEtB,IAAK,IAAI3C,EAAI,EAAGA,EAAIP,EAAOkC,SAAU3B,IAAK,CACxC,MAAM4C,EAAKhB,EAAY1B,IAAIF,GAC3B,IAAK4C,GAAoB,IAAdA,EAAGlG,OAAc,SAC5B,MAAMmG,EAAKd,EAAa7B,IAAIF,GAC5B,IAAI8C,EAAQ,EACZ,GAAID,EACF,IAAMC,EAAQ7B,KAAKkB,MAAMU,GAAInG,MAAQ,CAAE,MAAmB,CACxDoG,EAAQH,IACVA,EAAiBG,EACjBJ,EAAa1C,EAEjB,CAEA,MAAMgB,EAAY0B,GAAc,EAAKZ,EAAa5B,IAAIwC,IAAe,GAAM,GACrEK,EAAkBL,GAAc,EAAKX,EAAa7B,IAAIwC,IAAe,KAAQ,KAEnF,GAAIA,EAAa,EAEf,YADA,QAAWrG,QAAQ,GAAGqF,mGAKxBvC,EAAO6D,OAAO,KAAQC,gBAAiB7D,GAKrC,IACE,MAAM8D,EAAqCjC,KAAKkB,MAAMY,GAChDI,GAAW,QAAqBhE,GAAQZ,OAAQ6E,GAAMA,EAAEC,WAAa,KAAmBC,YAC9F,QAAqBnE,EAAQ,IAAIgE,KAAaD,GAChD,CAAE,MAAOK,GACPC,QAAQC,KAAK,2DAA4DF,EAC3E,CAGF,GAAuBvB,EACrB,IACE,MAAM0B,GAAW,QAA4BxH,EAAIiD,GACjD,IAAK,IAAIa,EAAI,EAAGA,EAAIP,EAAOkC,SAAU3B,IAAK,CACxC,MAAMiC,EAASD,EAAgB9B,IAAIF,GAC7B2D,EAAe5B,EAAa7B,IAAIF,GACtC,IAAKiC,IAAW0B,EAAc,SAC9B,MAAMzB,EAAuCjB,KAAKkB,MAAMF,GAClD2B,EAA8B3C,KAAKkB,MAAMwB,GAEzCE,EAAiC,GACvC,IAAK,MAAMC,KAAUF,EAAY,CAC/B,GAAoB,MAAhBE,EAAOC,OAA+B,MAAdD,EAAOE,IAAa,SAChD,MAAMC,EAAe/B,EAAa4B,EAAOC,OACnCG,EAAahC,EAAa4B,EAAOE,KACnB,MAAhBC,GAAsC,MAAdC,GAC5BL,EAAWrG,KAAK,CACd2G,aAAcL,EAAOM,GACrBC,cAAeJ,EACfK,iBAAkBJ,EAClBK,aAAcT,EAAOC,MACrBS,gBAAiB,IAErB,CAEA,MAAMC,GAAe,QAAkBf,EAAU1D,IAAM,IACvD,QAAkB0D,EAAU1D,GAAG,QAAayE,EAAcZ,GAAY,GAAM,GAC9E,CACF,CAAE,MAAON,GACPC,QAAQC,KAAK,0DAA2DF,EAC1E,CAIF,MAAMmB,EArOR,SACExI,EAAkBiD,EAA2BM,GAE7C,MAAMuC,EAAkBvC,EAAOjB,IAAI,iBACnC,IAAKwD,EAAiB,OAAO,KAG7B,MAAM2C,EAAW,IAAIC,IACfC,EAA6C,GAC7CC,EAAuB,GACvBC,EAAwB,GAE9B,IAAK,IAAI/E,EAAI,EAAGA,EAAIP,EAAOkC,SAAU3B,IAAK,CACxC,MAAMiC,EAASD,EAAgB9B,IAAIF,GAC7BgF,EAAS7F,EAAOe,IAAIF,IAAM,GAChC,GAAKiC,EAML,IACE,MAAMC,EAAuCjB,KAAKkB,MAAMF,GACxD4C,EAAQrH,KAAK0E,GACb,IAAK,MAAM+C,KAAQjH,OAAOkH,KAAKhD,GAC7ByC,EAAS3F,IAAIiG,GAGf,MAAME,EAAcnH,OAAOC,OAAOiE,GAC5BkD,EAAUC,KAAKC,OAAOH,GACtBI,EAAUF,KAAKG,OAAOL,GAC5BL,EAAWtH,KAAK4H,GAChBL,EAAYvH,KAAK6H,KAAKG,IAAI,EAAGR,EAAOtI,OAAS6I,EAAU,GACzD,CAAE,MACAV,EAAQrH,KAAK,MACbsH,EAAWtH,KAAK,GAChBuH,EAAYvH,KAAK,EACnB,MArBEqH,EAAQrH,KAAK,MACbsH,EAAWtH,KAAK,GAChBuH,EAAYvH,KAAK,EAoBrB,CAEA,GAAsB,IAAlBmH,EAASc,KAAY,OAAO,KAEhC,MAAMC,EAAYL,KAAKG,IAAI,KAAMV,GAC3Ba,EAAaN,KAAKG,IAAI,KAAMT,GAG5Ba,EAAoCC,MAAMC,KAAKnB,GA1DxCoB,QAAQC,KAAK,CAAC5C,EAAG6C,KAC5B,MAAOC,EAAMC,GAAQC,EAAkBhD,IAChCiD,EAAMC,GAAQF,EAAkBH,GACvC,OAAIC,IAASG,EAAaH,EAAOG,EAC1BF,EAAKI,cAAcD,KAuDtBE,EAAqB,GAC3B,IAAK,IAAIC,EAAIf,EAAWe,EAAI,EAAGA,IAC7BD,EAAShJ,KAAK,KAAKiJ,KACrB,MAAMC,EAAsB,GAC5B,IAAK,IAAID,EAAI,EAAGA,GAAKd,EAAYc,IAC/BC,EAAUlJ,KAAK,KAAKiJ,KAEtB,MAAME,EAAmB,IAAIH,KAAaZ,KAAoBc,GACxDE,EAAWD,EAAiBjK,OAC5BmK,EAAYnB,EAGZoB,EAAkB,IAAIC,IAC5B,IAAK,IAAIC,EAAI,EAAGA,EAAIpB,EAAgBlJ,OAAQsK,IAC1CF,EAAgBhG,IAAI8E,EAAgBoB,GAAIH,EAAYG,GAGtD,MAAMC,EAAU/K,EAAGK,QAAQ2K,cAAc,GAAG/H,EAAO9B,kBAC7C8J,EAAa,SAAU7G,SAAS,cAAeC,OAAQ0G,EAAS/K,EAAGyF,UAEzE,IAAK,IAAI3B,EAAI,EAAGA,EAAI9D,EAAGyF,SAAU3B,IAAK,CACpC,MAAMnC,EAAMmC,EAAI6E,EAAQnI,OAASmI,EAAQ7E,GAAK,KACxCgF,EAAS7F,EAAOe,IAAIF,IAAM,GAEhC,IAAKnC,EAAK,CACRsJ,EAAWrG,IAAId,EAAG,IAAIoH,OAAOR,IAC7B,QACF,CAEA,MAAMS,EAAU,IAAIxB,MAAce,GAAUU,KAAK,KAGjD,IAAK,MAAO/E,EAASgF,KAAYvJ,OAAOyE,QAAQ5E,GAAM,CACpD,MAAM2J,EAAOV,EAAgB5G,IAAIqC,GACrB,MAARiF,GAAgBD,EAAUvC,EAAOtI,SACnC2K,EAAQG,GAAQxC,EAAOuC,GAC3B,CAGA,MAAME,EAAS3C,EAAW9E,GACpBmF,EAAcnH,OAAOC,OAAOJ,GAC5BuH,EAAUC,KAAKC,OAAOH,GAC5B,IAAK,IAAIsB,EAAI,EAAGA,EAAIgB,EAAQhB,IAC1BY,EAAQR,EAAYY,EAAShB,GAAKzB,EAAOI,EAAUqC,EAAShB,GAG9D,MAAMiB,EAAU3C,EAAY/E,GACtBuF,EAAUF,KAAKG,OAAOL,GACtBwC,EAAYd,EAAYjB,EAAgBlJ,OAC9C,IAAK,IAAI+J,EAAI,EAAGA,EAAIiB,EAASjB,IAC3BY,EAAQM,EAAYlB,GAAKzB,EAAOO,EAAU,EAAIkB,GAEhDU,EAAWrG,IAAId,EAAGqH,EAAQO,KAAK,IACjC,CAUA,OAPAT,EAAW1I,QAAU,UAAWhC,cAChC0K,EAAWnE,OAAO,KAAQqE,QAAS,WACnCF,EAAWnE,OAAO,KAAQ6E,SAAU,KAASC,IAC7CX,EAAWrK,KAAKiL,MAAQ,KAASC,MACjCb,EAAWnE,OAAO,OAAQiF,aAAc,YACxCd,EAAWnE,OAAO,KAAQjC,cAAe4F,EAAiBiB,KAAK,OAExD,CAACT,aAAYR,mBAAkBE,YACxC,CAsHoBqB,CAAoBhM,EAAIiD,EAAQM,GAClD,GAAIiF,IACFxI,EAAGK,QAAQ4L,OAAOzD,EAAUyC,WAAYjL,EAAGK,QAAQ6L,SAASC,QAAQlJ,GAAU,GAC1E,QAAWhD,IAAIC,YAAcF,GAC/B,QAAWC,GAAGmM,KAAKC,aAAapJ,EAAQ,GAE1CuF,EAAUyC,WAAWnE,OAAO,KAAQC,gBAAiB7D,GAEjDqC,GACF,IACE,MAAMyB,EAAqCjC,KAAKkB,MAAMY,GAChDI,GAAW,QAAqBuB,EAAUyC,YAAY5I,OAAQ6E,GAAMA,EAAEC,WAAa,KAAmBC,YAC5G,QAAqBoB,EAAUyC,WAAY,IAAIhE,KAAaD,GAkB9D,CAAE,MAAOK,GACPC,QAAQC,KAAK,4DAA6DF,EAC5E,CAIJrH,EAAGsM,oBAEH,QAAWC,KAAK,sBAAsBrJ,kBAA2B4B,IACnE,CAnUQ0H,CAAsBxM,EAAIiD,EAAQM,EAAQL,GAAY,EAAME,EAAOhC,MAUrE,CAAE,MAAOiG,GACP,QAAWvG,MAAM,qBAAqBuG,EAAIoF,SAAWpF,KACrDC,QAAQxG,MAAMuG,EAChB,C,QACEhE,EAAGqJ,OACL,IAGJ9J,EAAO+J,MACT,CAKA,SAASxG,EAAyByG,GAChC,MAAMjL,EAAgB,GACtB,IAAK,IAAIkL,EAAI,EAAGA,EAAID,EAAUpM,OAAQqM,IACf,MAAjBD,EAAUC,IAA+B,MAAjBD,EAAUC,IACpClL,EAAIL,KAAKuL,GAEb,OAAOlL,CACT,CAIA,SAASuI,EAAkBnB,GACzB,MAAM+D,EAAQ/D,EAAK+D,MAAM,mBACzB,OAAKA,EACE,CAACC,SAASD,EAAM,GAAI,IAAKA,EAAM,IADnB,CAACE,IAAUjE,EAEhC,C","sources":["webpack://bio/./node_modules/@datagrok-libraries/bio/src/utils/macromolecule/numbering-schemes.js","webpack://bio/./src/utils/annotations/numbering-ui.ts"],"sourcesContent":["/** Static antibody numbering scheme region definitions.\n *\n * Contains FR/CDR boundaries for IMGT, Kabat, Chothia, and AHo schemes.\n * These are the canonical position ranges used after a numbering tool\n * (e.g. ANARCI) has assigned scheme-specific position names.\n */\nexport var NumberingScheme;\n(function (NumberingScheme) {\n NumberingScheme[\"IMGT\"] = \"IMGT\";\n NumberingScheme[\"Kabat\"] = \"Kabat\";\n NumberingScheme[\"Chothia\"] = \"Chothia\";\n NumberingScheme[\"AHo\"] = \"AHo\";\n})(NumberingScheme || (NumberingScheme = {}));\nexport var ChainType;\n(function (ChainType) {\n ChainType[\"Heavy\"] = \"Heavy\";\n ChainType[\"Light_Kappa\"] = \"Light_Kappa\";\n ChainType[\"Light_Lambda\"] = \"Light_Lambda\";\n})(ChainType || (ChainType = {}));\n/** IMGT region boundaries — the reference numbering scheme.\n * @see https://www.imgt.org/IMGTScientificChart/Numbering/IMGTIGVLsuperfamily.html\n */\nexport const IMGT_REGIONS = {\n [ChainType.Heavy]: [\n { name: 'FR1', type: 'FR', startPosition: '1', endPosition: '26', chainType: ChainType.Heavy },\n { name: 'CDR1', type: 'CDR', startPosition: '27', endPosition: '38', chainType: ChainType.Heavy },\n { name: 'FR2', type: 'FR', startPosition: '39', endPosition: '55', chainType: ChainType.Heavy },\n { name: 'CDR2', type: 'CDR', startPosition: '56', endPosition: '65', chainType: ChainType.Heavy },\n { name: 'FR3', type: 'FR', startPosition: '66', endPosition: '104', chainType: ChainType.Heavy },\n { name: 'CDR3', type: 'CDR', startPosition: '105', endPosition: '117', chainType: ChainType.Heavy },\n { name: 'FR4', type: 'FR', startPosition: '118', endPosition: '128', chainType: ChainType.Heavy },\n ],\n [ChainType.Light_Kappa]: [\n { name: 'FR1', type: 'FR', startPosition: '1', endPosition: '26', chainType: ChainType.Light_Kappa },\n { name: 'CDR1', type: 'CDR', startPosition: '27', endPosition: '38', chainType: ChainType.Light_Kappa },\n { name: 'FR2', type: 'FR', startPosition: '39', endPosition: '55', chainType: ChainType.Light_Kappa },\n { name: 'CDR2', type: 'CDR', startPosition: '56', endPosition: '65', chainType: ChainType.Light_Kappa },\n { name: 'FR3', type: 'FR', startPosition: '66', endPosition: '104', chainType: ChainType.Light_Kappa },\n { name: 'CDR3', type: 'CDR', startPosition: '105', endPosition: '117', chainType: ChainType.Light_Kappa },\n { name: 'FR4', type: 'FR', startPosition: '118', endPosition: '127', chainType: ChainType.Light_Kappa },\n ],\n [ChainType.Light_Lambda]: [\n { name: 'FR1', type: 'FR', startPosition: '1', endPosition: '26', chainType: ChainType.Light_Lambda },\n { name: 'CDR1', type: 'CDR', startPosition: '27', endPosition: '38', chainType: ChainType.Light_Lambda },\n { name: 'FR2', type: 'FR', startPosition: '39', endPosition: '55', chainType: ChainType.Light_Lambda },\n { name: 'CDR2', type: 'CDR', startPosition: '56', endPosition: '65', chainType: ChainType.Light_Lambda },\n { name: 'FR3', type: 'FR', startPosition: '66', endPosition: '104', chainType: ChainType.Light_Lambda },\n { name: 'CDR3', type: 'CDR', startPosition: '105', endPosition: '117', chainType: ChainType.Light_Lambda },\n { name: 'FR4', type: 'FR', startPosition: '118', endPosition: '127', chainType: ChainType.Light_Lambda },\n ],\n};\n/** Kabat region boundaries (differ mainly in CDR definitions). */\nexport const KABAT_REGIONS = {\n [ChainType.Heavy]: [\n { name: 'FR1', type: 'FR', startPosition: '1', endPosition: '30', chainType: ChainType.Heavy },\n { name: 'CDR1', type: 'CDR', startPosition: '31', endPosition: '35', chainType: ChainType.Heavy },\n { name: 'FR2', type: 'FR', startPosition: '36', endPosition: '49', chainType: ChainType.Heavy },\n { name: 'CDR2', type: 'CDR', startPosition: '50', endPosition: '65', chainType: ChainType.Heavy },\n { name: 'FR3', type: 'FR', startPosition: '66', endPosition: '94', chainType: ChainType.Heavy },\n { name: 'CDR3', type: 'CDR', startPosition: '95', endPosition: '102', chainType: ChainType.Heavy },\n { name: 'FR4', type: 'FR', startPosition: '103', endPosition: '113', chainType: ChainType.Heavy },\n ],\n [ChainType.Light_Kappa]: [\n { name: 'FR1', type: 'FR', startPosition: '1', endPosition: '23', chainType: ChainType.Light_Kappa },\n { name: 'CDR1', type: 'CDR', startPosition: '24', endPosition: '34', chainType: ChainType.Light_Kappa },\n { name: 'FR2', type: 'FR', startPosition: '35', endPosition: '49', chainType: ChainType.Light_Kappa },\n { name: 'CDR2', type: 'CDR', startPosition: '50', endPosition: '56', chainType: ChainType.Light_Kappa },\n { name: 'FR3', type: 'FR', startPosition: '57', endPosition: '88', chainType: ChainType.Light_Kappa },\n { name: 'CDR3', type: 'CDR', startPosition: '89', endPosition: '97', chainType: ChainType.Light_Kappa },\n { name: 'FR4', type: 'FR', startPosition: '98', endPosition: '107', chainType: ChainType.Light_Kappa },\n ],\n [ChainType.Light_Lambda]: [\n { name: 'FR1', type: 'FR', startPosition: '1', endPosition: '23', chainType: ChainType.Light_Lambda },\n { name: 'CDR1', type: 'CDR', startPosition: '24', endPosition: '34', chainType: ChainType.Light_Lambda },\n { name: 'FR2', type: 'FR', startPosition: '35', endPosition: '49', chainType: ChainType.Light_Lambda },\n { name: 'CDR2', type: 'CDR', startPosition: '50', endPosition: '56', chainType: ChainType.Light_Lambda },\n { name: 'FR3', type: 'FR', startPosition: '57', endPosition: '88', chainType: ChainType.Light_Lambda },\n { name: 'CDR3', type: 'CDR', startPosition: '89', endPosition: '97', chainType: ChainType.Light_Lambda },\n { name: 'FR4', type: 'FR', startPosition: '98', endPosition: '107', chainType: ChainType.Light_Lambda },\n ],\n};\n/** Chothia region boundaries. */\nexport const CHOTHIA_REGIONS = {\n [ChainType.Heavy]: [\n { name: 'FR1', type: 'FR', startPosition: '1', endPosition: '25', chainType: ChainType.Heavy },\n { name: 'CDR1', type: 'CDR', startPosition: '26', endPosition: '32', chainType: ChainType.Heavy },\n { name: 'FR2', type: 'FR', startPosition: '33', endPosition: '51', chainType: ChainType.Heavy },\n { name: 'CDR2', type: 'CDR', startPosition: '52', endPosition: '56', chainType: ChainType.Heavy },\n { name: 'FR3', type: 'FR', startPosition: '57', endPosition: '94', chainType: ChainType.Heavy },\n { name: 'CDR3', type: 'CDR', startPosition: '95', endPosition: '102', chainType: ChainType.Heavy },\n { name: 'FR4', type: 'FR', startPosition: '103', endPosition: '113', chainType: ChainType.Heavy },\n ],\n [ChainType.Light_Kappa]: [\n { name: 'FR1', type: 'FR', startPosition: '1', endPosition: '25', chainType: ChainType.Light_Kappa },\n { name: 'CDR1', type: 'CDR', startPosition: '26', endPosition: '32', chainType: ChainType.Light_Kappa },\n { name: 'FR2', type: 'FR', startPosition: '33', endPosition: '49', chainType: ChainType.Light_Kappa },\n { name: 'CDR2', type: 'CDR', startPosition: '50', endPosition: '52', chainType: ChainType.Light_Kappa },\n { name: 'FR3', type: 'FR', startPosition: '53', endPosition: '90', chainType: ChainType.Light_Kappa },\n { name: 'CDR3', type: 'CDR', startPosition: '91', endPosition: '96', chainType: ChainType.Light_Kappa },\n { name: 'FR4', type: 'FR', startPosition: '97', endPosition: '107', chainType: ChainType.Light_Kappa },\n ],\n [ChainType.Light_Lambda]: [\n { name: 'FR1', type: 'FR', startPosition: '1', endPosition: '25', chainType: ChainType.Light_Lambda },\n { name: 'CDR1', type: 'CDR', startPosition: '26', endPosition: '32', chainType: ChainType.Light_Lambda },\n { name: 'FR2', type: 'FR', startPosition: '33', endPosition: '49', chainType: ChainType.Light_Lambda },\n { name: 'CDR2', type: 'CDR', startPosition: '50', endPosition: '52', chainType: ChainType.Light_Lambda },\n { name: 'FR3', type: 'FR', startPosition: '53', endPosition: '90', chainType: ChainType.Light_Lambda },\n { name: 'CDR3', type: 'CDR', startPosition: '91', endPosition: '96', chainType: ChainType.Light_Lambda },\n { name: 'FR4', type: 'FR', startPosition: '97', endPosition: '107', chainType: ChainType.Light_Lambda },\n ],\n};\n/** AHo region boundaries. */\nexport const AHO_REGIONS = {\n [ChainType.Heavy]: [\n { name: 'FR1', type: 'FR', startPosition: '1', endPosition: '24', chainType: ChainType.Heavy },\n { name: 'CDR1', type: 'CDR', startPosition: '25', endPosition: '40', chainType: ChainType.Heavy },\n { name: 'FR2', type: 'FR', startPosition: '41', endPosition: '55', chainType: ChainType.Heavy },\n { name: 'CDR2', type: 'CDR', startPosition: '56', endPosition: '78', chainType: ChainType.Heavy },\n { name: 'FR3', type: 'FR', startPosition: '79', endPosition: '108', chainType: ChainType.Heavy },\n { name: 'CDR3', type: 'CDR', startPosition: '109', endPosition: '138', chainType: ChainType.Heavy },\n { name: 'FR4', type: 'FR', startPosition: '139', endPosition: '149', chainType: ChainType.Heavy },\n ],\n [ChainType.Light_Kappa]: [\n { name: 'FR1', type: 'FR', startPosition: '1', endPosition: '24', chainType: ChainType.Light_Kappa },\n { name: 'CDR1', type: 'CDR', startPosition: '25', endPosition: '40', chainType: ChainType.Light_Kappa },\n { name: 'FR2', type: 'FR', startPosition: '41', endPosition: '55', chainType: ChainType.Light_Kappa },\n { name: 'CDR2', type: 'CDR', startPosition: '56', endPosition: '78', chainType: ChainType.Light_Kappa },\n { name: 'FR3', type: 'FR', startPosition: '79', endPosition: '108', chainType: ChainType.Light_Kappa },\n { name: 'CDR3', type: 'CDR', startPosition: '109', endPosition: '138', chainType: ChainType.Light_Kappa },\n { name: 'FR4', type: 'FR', startPosition: '139', endPosition: '149', chainType: ChainType.Light_Kappa },\n ],\n [ChainType.Light_Lambda]: [\n { name: 'FR1', type: 'FR', startPosition: '1', endPosition: '24', chainType: ChainType.Light_Lambda },\n { name: 'CDR1', type: 'CDR', startPosition: '25', endPosition: '40', chainType: ChainType.Light_Lambda },\n { name: 'FR2', type: 'FR', startPosition: '41', endPosition: '55', chainType: ChainType.Light_Lambda },\n { name: 'CDR2', type: 'CDR', startPosition: '56', endPosition: '78', chainType: ChainType.Light_Lambda },\n { name: 'FR3', type: 'FR', startPosition: '79', endPosition: '108', chainType: ChainType.Light_Lambda },\n { name: 'CDR3', type: 'CDR', startPosition: '109', endPosition: '138', chainType: ChainType.Light_Lambda },\n { name: 'FR4', type: 'FR', startPosition: '139', endPosition: '149', chainType: ChainType.Light_Lambda },\n ],\n};\n/** Lookup table: scheme → chain type → region defs */\nexport const SCHEME_REGIONS = {\n [NumberingScheme.IMGT]: IMGT_REGIONS,\n [NumberingScheme.Kabat]: KABAT_REGIONS,\n [NumberingScheme.Chothia]: CHOTHIA_REGIONS,\n [NumberingScheme.AHo]: AHO_REGIONS,\n};\n//# sourceMappingURL=numbering-schemes.js.map","/* eslint-disable max-len */\nimport * as grok from 'datagrok-api/grok';\nimport * as ui from 'datagrok-api/ui';\nimport * as DG from 'datagrok-api/dg';\n\nimport {TAGS as bioTAGS, ALIGNMENT, ALPHABET, NOTATION} from '@datagrok-libraries/bio/src/utils/macromolecule';\nimport {\n SeqAnnotation, SeqAnnotationHit, AnnotationCategory,\n} from '@datagrok-libraries/bio/src/utils/macromolecule/annotations';\nimport {NumberingScheme} from '@datagrok-libraries/bio/src/utils/macromolecule/numbering-schemes';\nimport {\n setColumnAnnotations, getColumnAnnotations,\n getOrCreateAnnotationColumn, getRowAnnotations, setRowAnnotations, mergeRowHits,\n} from './annotation-manager';\nimport {_package} from '../../package';\nimport type {NumberingResult, Scheme} from '../antibody-numbering (WIP)';\nimport {VdRegionsViewer} from '../../viewers/vd-regions-viewer';\nimport { VdRegion, VdRegionType } from '@datagrok-libraries/bio/src/viewers/vd-regions';\n\nconst BUILTIN_ENGINE_KEY = '__builtin__';\nconst BUILTIN_ENGINE_LABEL = 'Built-in (TypeScript)';\n\n/** An engine entry: either a dynamically discovered DG.Func or the built-in TS engine. */\ninterface NumberingEngine {\n /** Display label for the dropdown */\n label: string;\n /** Unique key — nqName for DG.Func engines, BUILTIN_ENGINE_KEY for built-in */\n key: string;\n /** The DG.Func to call, or null for the built-in engine */\n func: DG.Func | null;\n}\n\n/** Discovers all registered antibody numbering engines + the built-in TS engine.\n * Dynamic engines (meta.role = 'antibodyNumbering') come first; built-in is last. */\nfunction discoverEngines(): NumberingEngine[] {\n const engines: NumberingEngine[] = [];\n\n const funcs = DG.Func.find({meta: {role: 'antibodyNumbering'}});\n if (funcs.length === 0) {\n grok.shell.error('No antibody numbering engines found. Make sure that Proteomics plugin is installed and up to date.');\n throw new Error('No external antibody numbering engines found. Make sure that Proteomics plugin is installed and up to date.');\n }\n for (const f of funcs) {\n const pkgName = f.package?.name ?? '';\n const label = f.friendlyName || f.name;\n engines.push({\n label: label,\n key: pkgName ? `${pkgName}:${f.name}` : f.name,\n func: f,\n });\n }\n\n // Built-in TS engine is always last\n engines.push({label: BUILTIN_ENGINE_LABEL, key: BUILTIN_ENGINE_KEY, func: null});\n return engines;\n}\n\n/** Converts TS NumberingResult[] to a DG.DataFrame matching the expected output shape.\n * Columns: position_names, chain_type, annotations_json, numbering_detail, numbering_map. */\nexport function numberingResultsToDataFrame(results: NumberingResult[]): DG.DataFrame {\n const n = results.length;\n const posNames = DG.Column.fromType(DG.COLUMN_TYPE.STRING, 'position_names', n);\n const chainTypes = DG.Column.fromType(DG.COLUMN_TYPE.STRING, 'chain_type', n);\n const annotJson = DG.Column.fromType(DG.COLUMN_TYPE.STRING, 'annotations_json', n);\n const numDetail = DG.Column.fromType(DG.COLUMN_TYPE.STRING, 'numbering_detail', n);\n const numMap = DG.Column.fromType(DG.COLUMN_TYPE.STRING, 'numbering_map', n);\n\n for (let i = 0; i < n; i++) {\n const r = results[i];\n if (r.error && r.percentIdentity < 0.3) {\n posNames.set(i, '');\n chainTypes.set(i, '');\n annotJson.set(i, '[]');\n numDetail.set(i, '');\n numMap.set(i, '');\n } else {\n posNames.set(i, r.positionNames);\n chainTypes.set(i, r.chainType);\n annotJson.set(i, JSON.stringify(r.annotations));\n numDetail.set(i, JSON.stringify(r.numberingDetail));\n numMap.set(i, JSON.stringify(r.numberingMap));\n }\n }\n\n return DG.DataFrame.fromColumns([posNames, chainTypes, annotJson, numDetail, numMap]);\n}\n\n/** Runs the built-in TS numbering engine on all rows of a sequence column. */\nasync function runBuiltinNumbering(\n seqCol: DG.Column<string>, schemeName: string,\n): Promise<DG.DataFrame> {\n const {numberSequences, extractSequence} = await import('../antibody-numbering (WIP)');\n const scheme = schemeName.toLowerCase() as Scheme;\n\n const sequences: string[] = [];\n for (let i = 0; i < seqCol.length; i++) {\n const raw = seqCol.get(i);\n sequences.push(extractSequence(raw ?? ''));\n }\n\n const results = numberSequences(sequences, scheme);\n return numberingResultsToDataFrame(results);\n}\n\nexport function showNumberingSchemeDialog(): void {\n const df = grok.shell.tv?.dataFrame;\n if (!df) {\n grok.shell.warning('No table open');\n return;\n }\n\n const seqCols = df.columns.bySemTypeAll(DG.SEMTYPE.MACROMOLECULE);\n if (seqCols.length === 0) {\n grok.shell.warning('No macromolecule columns found');\n return;\n }\n\n const engines = discoverEngines();\n const engineLabels = engines.map((e) => e.label);\n const schemeChoices = Object.values(NumberingScheme);\n\n const tableInput = ui.input.table('Table', {value: df});\n const seqInput = ui.input.column('Sequence', {\n table: df, value: seqCols[0],\n filter: (col: DG.Column) => col.semType === DG.SEMTYPE.MACROMOLECULE,\n });\n const schemeInput = ui.input.choice('Scheme', {value: NumberingScheme.IMGT, items: schemeChoices});\n const engineInput = ui.input.choice('Engine', {\n value: engineLabels[0], items: engineLabels,\n });\n // const populateRegions = ui.input.bool('Populate FR/CDR regions', {value: true});\n // const openVdRegions = ui.input.bool('Open VD Regions viewer', {value: true});\n\n const dialog = ui.dialog({title: 'Apply Antibody Numbering'})\n .add(ui.inputs([tableInput, seqInput, schemeInput, engineInput]))\n .onOK(async () => {\n const seqCol = seqInput.value!;\n const schemeName = schemeInput.value!;\n const selectedLabel = engineInput.value!;\n const engine = engines.find((e) => e.label === selectedLabel) ?? engines[engines.length - 1];\n const pi = DG.TaskBarProgressIndicator.create(`Applying ${schemeName} numbering...`);\n try {\n let result: DG.DataFrame;\n if (engine.func)\n result = await engine.func.apply({df: df, seqCol: seqCol, scheme: schemeName.toLowerCase()});\n else\n result = await runBuiltinNumbering(seqCol, schemeName);\n\n applyNumberingResults(df, seqCol, result, schemeName, true, engine.label);\n\n // // Open VD Regions viewer\n // if (openVdRegions.value && grok.shell.tv) {\n // try {\n // await grok.shell.tv.dataFrame.plot.fromType('VdRegions', {});\n // } catch (err) {\n // console.warn('Could not open VD Regions viewer:', err);\n // }\n // }\n } catch (err: any) {\n grok.shell.error(`Numbering failed: ${err.message ?? err}`);\n console.error(err);\n } finally {\n pi.close();\n }\n });\n\n dialog.show();\n}\n\n/** Builds a map from ungapped character index to gapped character index.\n * Used when the source column is already aligned (MSA) — numbering engines strip gaps,\n * so their output indices refer to the ungapped sequence, not the gapped original. */\nfunction buildUngappedToGappedMap(gappedSeq: string): number[] {\n const map: number[] = [];\n for (let g = 0; g < gappedSeq.length; g++) {\n if (gappedSeq[g] !== '-' && gappedSeq[g] !== '.')\n map.push(g);\n }\n return map;\n}\n\n/** Parses a position code into [numericPart, insertionLetter] for sorting.\n * E.g. \"27\" → [27, \"\"], \"111A\" → [111, \"A\"], \"27B\" → [27, \"B\"]. */\nfunction parsePositionCode(code: string): [number, string] {\n const match = code.match(/^(\\d+)([A-Z]?)$/);\n if (!match) return [Infinity, code];\n return [parseInt(match[1], 10), match[2]];\n}\n\n/** Sorts position codes in scheme order: numeric ascending, then insertion letter. */\nfunction sortPositionCodes(codes: string[]): string[] {\n return codes.slice().sort((a, b) => {\n const [numA, insA] = parsePositionCode(a);\n const [numB, insB] = parsePositionCode(b);\n if (numA !== numB) return numA - numB;\n return insA.localeCompare(insB);\n });\n}\n\n/** Builds unified position list from all rows and creates an aligned sequence column.\n * Includes flanking residues (before/after the numbered region) padded with gaps.\n * Layout: [pre-region gaps+residues] [scheme-aligned region] [post-region residues+gaps]\n * @returns aligned column, full position list (including flanking), and the pre-region offset. */\nfunction createAlignedColumn(\n df: DG.DataFrame, seqCol: DG.Column<string>, result: DG.DataFrame,\n): {alignedCol: DG.Column<string>; unifiedPositions: string[]; preOffset: number} | null {\n const numberingMapCol = result.col('numbering_map');\n if (!numberingMapCol) return null;\n\n // Pass 1: collect all scheme position codes and per-row flanking lengths\n const allCodes = new Set<string>();\n const rowMaps: (Record<string, number> | null)[] = [];\n const rowPreLens: number[] = []; // chars before first numbered position\n const rowPostLens: number[] = []; // chars after last numbered position\n\n for (let i = 0; i < result.rowCount; i++) {\n const mapStr = numberingMapCol.get(i);\n const rawSeq = seqCol.get(i) ?? '';\n if (!mapStr) {\n rowMaps.push(null);\n rowPreLens.push(0);\n rowPostLens.push(0);\n continue;\n }\n try {\n const posToCharIdx: Record<string, number> = JSON.parse(mapStr);\n rowMaps.push(posToCharIdx);\n for (const code of Object.keys(posToCharIdx))\n allCodes.add(code);\n\n // Find min/max char indices that were numbered\n const charIndices = Object.values(posToCharIdx);\n const minChar = Math.min(...charIndices);\n const maxChar = Math.max(...charIndices);\n rowPreLens.push(minChar); // chars before first numbered\n rowPostLens.push(Math.max(0, rawSeq.length - maxChar - 1)); // chars after last numbered\n } catch {\n rowMaps.push(null);\n rowPreLens.push(0);\n rowPostLens.push(0);\n }\n }\n\n if (allCodes.size === 0) return null;\n\n const maxPreLen = Math.max(0, ...rowPreLens);\n const maxPostLen = Math.max(0, ...rowPostLens);\n\n // Build position names: [pre-flanking] + [scheme positions] + [post-flanking]\n const schemePositions = sortPositionCodes(Array.from(allCodes));\n const preNames: string[] = [];\n for (let p = maxPreLen; p > 0; p--)\n preNames.push(`N-${p}`);\n const postNames: string[] = [];\n for (let p = 1; p <= maxPostLen; p++)\n postNames.push(`C+${p}`);\n\n const unifiedPositions = [...preNames, ...schemePositions, ...postNames];\n const totalLen = unifiedPositions.length;\n const preOffset = maxPreLen; // scheme region starts at this index in the aligned string\n\n // Map scheme position codes → index in the full unified list\n const posToUnifiedIdx = new Map<string, number>();\n for (let s = 0; s < schemePositions.length; s++)\n posToUnifiedIdx.set(schemePositions[s], preOffset + s);\n\n // Pass 2: build aligned sequences\n const colName = df.columns.getUnusedName(`${seqCol.name} (aligned)`);\n const alignedCol = DG.Column.fromType(DG.COLUMN_TYPE.STRING, colName, df.rowCount);\n\n for (let i = 0; i < df.rowCount; i++) {\n const map = i < rowMaps.length ? rowMaps[i] : null;\n const rawSeq = seqCol.get(i) ?? '';\n\n if (!map) {\n alignedCol.set(i, '-'.repeat(totalLen));\n continue;\n }\n\n const aligned = new Array<string>(totalLen).fill('-');\n\n // Place scheme-numbered residues\n for (const [posCode, charIdx] of Object.entries(map)) {\n const uIdx = posToUnifiedIdx.get(posCode);\n if (uIdx != null && charIdx < rawSeq.length)\n aligned[uIdx] = rawSeq[charIdx];\n }\n\n // Place pre-region flanking residues (right-aligned within the pre-region block)\n const preLen = rowPreLens[i];\n const charIndices = Object.values(map);\n const minChar = Math.min(...charIndices);\n for (let p = 0; p < preLen; p++)\n aligned[preOffset - preLen + p] = rawSeq[minChar - preLen + p];\n\n // Place post-region flanking residues (left-aligned within the post-region block)\n const postLen = rowPostLens[i];\n const maxChar = Math.max(...charIndices);\n const postStart = preOffset + schemePositions.length;\n for (let p = 0; p < postLen; p++)\n aligned[postStart + p] = rawSeq[maxChar + 1 + p];\n\n alignedCol.set(i, aligned.join(''));\n }\n\n // Set macromolecule tags on the aligned column\n alignedCol.semType = DG.SEMTYPE.MACROMOLECULE;\n alignedCol.setTag(bioTAGS.aligned, ALIGNMENT.SEQ_MSA);\n alignedCol.setTag(bioTAGS.alphabet, ALPHABET.PT);\n alignedCol.meta.units = NOTATION.FASTA;\n alignedCol.setTag(DG.Tags.CellRenderer, 'sequence');\n alignedCol.setTag(bioTAGS.positionNames, unifiedPositions.join(', '));\n\n return {alignedCol, unifiedPositions, preOffset};\n}\n\n/** Applies numbering results (from either engine) to the sequence column and dataframe.\n *\n * Annotation strategy:\n * - Original column: row-level region spans only (char indices from numbering_map).\n * No column-level position names (they differ per row).\n * - Aligned column: column-level annotations only (all rows share unified positions).\n * Position names tag set to the unified list. */\nfunction applyNumberingResults(\n df: DG.DataFrame, seqCol: DG.Column<string>, result: DG.DataFrame,\n schemeName: string, populateRegions: boolean, engineLabel: string,\n): void {\n if (!result || result.rowCount === 0) {\n grok.shell.warning('No numbering results returned');\n return;\n }\n\n const posNamesCol = result.getCol('position_names');\n const chainTypeCol = result.getCol('chain_type');\n const annotJsonCol = result.getCol('annotations_json');\n const numberingMapCol = result.col('numbering_map');\n\n // If the source column is already aligned (MSA), remap numbering_map indices\n // from ungapped to gapped, since numbering engines strip gaps before processing.\n const isAligned = true; // always treat as aligned to handle remapping\n if (isAligned && numberingMapCol) {\n for (let i = 0; i < result.rowCount; i++) {\n const mapStr = numberingMapCol.get(i);\n if (!mapStr) continue;\n try {\n const posToCharIdx: Record<string, number> = JSON.parse(mapStr);\n const rawSeq = seqCol.get(i) ?? '';\n const ungapToGap = buildUngappedToGappedMap(rawSeq);\n const remapped: Record<string, number> = {};\n for (const [posCode, ungappedIdx] of Object.entries(posToCharIdx)) {\n if (ungappedIdx < ungapToGap.length)\n remapped[posCode] = ungapToGap[ungappedIdx];\n }\n numberingMapCol.set(i, JSON.stringify(remapped));\n } catch { /* skip */ }\n }\n }\n\n // Pick the row with the most annotations for column-level representative data\n let bestRowIdx = -1;\n let bestAnnotCount = -1;\n\n for (let i = 0; i < result.rowCount; i++) {\n const pn = posNamesCol.get(i);\n if (!pn || pn.length === 0) continue;\n const aj = annotJsonCol.get(i);\n let count = 0;\n if (aj)\n try { count = JSON.parse(aj).length; } catch { /* skip */ }\n if (count > bestAnnotCount) {\n bestAnnotCount = count;\n bestRowIdx = i;\n }\n }\n\n const chainType = bestRowIdx >= 0 ? (chainTypeCol.get(bestRowIdx) ?? '') : '';\n const annotationsJson = bestRowIdx >= 0 ? (annotJsonCol.get(bestRowIdx) ?? '[]') : '[]';\n\n if (bestRowIdx < 0) {\n grok.shell.warning(`${engineLabel} could not number the sequences. Check that they are valid antibody variable region sequences.`);\n return;\n }\n\n // Mark scheme on original column (no position names — they differ per row)\n seqCol.setTag(bioTAGS.numberingScheme, schemeName);\n\n // --- Original column: column-level annotation definitions (needed for renderer\n // to resolve annotation IDs → colors/names) + row-level region spans ---\n if (populateRegions) {\n try {\n const regionAnnotations: SeqAnnotation[] = JSON.parse(annotationsJson);\n const existing = getColumnAnnotations(seqCol).filter((a) => a.category !== AnnotationCategory.Structure);\n setColumnAnnotations(seqCol, [...existing, ...regionAnnotations]);\n } catch (err) {\n console.warn('Failed to set annotation definitions on original column:', err);\n }\n }\n\n if (populateRegions && numberingMapCol) {\n try {\n const annotCol = getOrCreateAnnotationColumn(df, seqCol);\n for (let i = 0; i < result.rowCount; i++) {\n const mapStr = numberingMapCol.get(i);\n const rowAnnotJson = annotJsonCol.get(i);\n if (!mapStr || !rowAnnotJson) continue;\n const posToCharIdx: Record<string, number> = JSON.parse(mapStr);\n const rowRegions: SeqAnnotation[] = JSON.parse(rowAnnotJson);\n\n const regionHits: SeqAnnotationHit[] = [];\n for (const region of rowRegions) {\n if (region.start == null || region.end == null) continue;\n const startCharIdx = posToCharIdx[region.start];\n const endCharIdx = posToCharIdx[region.end];\n if (startCharIdx == null || endCharIdx == null) continue;\n regionHits.push({\n annotationId: region.id,\n positionIndex: startCharIdx,\n endPositionIndex: endCharIdx,\n positionName: region.start,\n matchedMonomers: '',\n });\n }\n\n const existingHits = getRowAnnotations(annotCol, i) ?? [];\n setRowAnnotations(annotCol, i, mergeRowHits(existingHits, regionHits, true, false));\n }\n } catch (err) {\n console.warn('Failed to store per-row region data on original column:', err);\n }\n }\n\n // --- Aligned column: column-level annotations only ---\n const alignment = createAlignedColumn(df, seqCol, result);\n if (alignment) {\n df.columns.insert(alignment.alignedCol, df.columns.toList().indexOf(seqCol) + 1);\n if (grok.shell.tv?.dataFrame === df)\n grok.shell.tv.grid.scrollToCell(seqCol, 0);\n\n alignment.alignedCol.setTag(bioTAGS.numberingScheme, schemeName);\n\n if (populateRegions) {\n try {\n const regionAnnotations: SeqAnnotation[] = JSON.parse(annotationsJson);\n const existing = getColumnAnnotations(alignment.alignedCol).filter((a) => a.category !== AnnotationCategory.Structure);\n setColumnAnnotations(alignment.alignedCol, [...existing, ...regionAnnotations]);\n // chunk for vd regions viewer if that becomes a desired feature in the future\n // if (grok.shell.tv?.dataFrame === df) {\n // (async () => {\n // const vdRegionsViewer: VdRegionsViewer = await grok.shell.tv.dataFrame.plot.fromType('VdRegions',) as VdRegionsViewer;\n // vdRegionsViewer.chains = [chainType];\n // vdRegionsViewer.setData(regionAnnotations.map((a, i) => ({\n // type: a.name?.toLowerCase().startsWith('cdr') ? VdRegionType.CDR : VdRegionType.FR,\n // name: a.name,\n // chain: chainType,\n // order: i,\n // sequenceColumnName: alignment.alignedCol.name,\n // positionStartName: a.start ?? '',\n // positionEndName: a.end ?? '',\n // } satisfies VdRegion)));\n // grok.shell.tv.addViewer(vdRegionsViewer);\n // })();\n // }\n } catch (err) {\n console.warn('Failed to set column-level annotations on aligned column:', err);\n }\n }\n }\n\n df.fireValuesChanged();\n\n grok.shell.info(`Numbering applied: ${schemeName}, chain type: ${chainType}`);\n}\n"],"names":["NumberingScheme","ChainType","Heavy","Light_Kappa","Light_Lambda","IMGT","Kabat","Chothia","AHo","showNumberingSchemeDialog","df","tv","dataFrame","warning","seqCols","columns","bySemTypeAll","MACROMOLECULE","length","engines","funcs","find","meta","role","error","Error","f","pkgName","package","name","label","friendlyName","push","key","func","discoverEngines","engineLabels","map","e","schemeChoices","Object","values","tableInput","table","value","seqInput","column","filter","col","semType","schemeInput","choice","items","engineInput","dialog","title","add","onOK","async","seqCol","schemeName","selectedLabel","engine","pi","create","result","apply","scheme","toLowerCase","numberSequences","extractSequence","sequences","i","raw","get","results","n","posNames","fromType","STRING","chainTypes","annotJson","numDetail","numMap","r","percentIdentity","set","positionNames","chainType","JSON","stringify","annotations","numberingDetail","numberingMap","fromColumns","numberingResultsToDataFrame","runBuiltinNumbering","populateRegions","engineLabel","rowCount","posNamesCol","getCol","chainTypeCol","annotJsonCol","numberingMapCol","mapStr","posToCharIdx","parse","ungapToGap","buildUngappedToGappedMap","remapped","posCode","ungappedIdx","entries","bestRowIdx","bestAnnotCount","pn","aj","count","annotationsJson","setTag","numberingScheme","regionAnnotations","existing","a","category","Structure","err","console","warn","annotCol","rowAnnotJson","rowRegions","regionHits","region","start","end","startCharIdx","endCharIdx","annotationId","id","positionIndex","endPositionIndex","positionName","matchedMonomers","existingHits","alignment","allCodes","Set","rowMaps","rowPreLens","rowPostLens","rawSeq","code","keys","charIndices","minChar","Math","min","maxChar","max","size","maxPreLen","maxPostLen","schemePositions","Array","from","slice","sort","b","numA","insA","parsePositionCode","numB","insB","localeCompare","preNames","p","postNames","unifiedPositions","totalLen","preOffset","posToUnifiedIdx","Map","s","colName","getUnusedName","alignedCol","repeat","aligned","fill","charIdx","uIdx","preLen","postLen","postStart","join","alphabet","PT","units","FASTA","CellRenderer","createAlignedColumn","insert","toList","indexOf","grid","scrollToCell","fireValuesChanged","info","applyNumberingResults","message","close","show","gappedSeq","g","match","parseInt","Infinity"],"sourceRoot":""}
package/dist/287.js ADDED
@@ -0,0 +1,2 @@
1
+ "use strict";(self.webpackChunkbio=self.webpackChunkbio||[]).push([[287],{9287(e,o,t){t.d(o,{showAnnotationManagerDialog:()=>d});var n=t(4328),i=t(7389),l=t(6082),a=t(4517),r=t(3736);const s={[a.eI.Structure]:"Structure (FR/CDR)",[a.eI.Liability]:"Liability",[a.eI.PTM]:"Post-translational Modification",[a.eI.Custom]:"Custom"},c={[a.Hq.High]:"High",[a.Hq.Medium]:"Medium",[a.Hq.Low]:"Low",[a.Hq.Info]:"Info"};function d(){const e=n.shell.tv?.dataFrame;if(!e)return void n.shell.warning("No table open");const o=e.columns.bySemTypeAll(l.SEMTYPE.MACROMOLECULE);if(0===o.length)return void n.shell.warning("No macromolecule columns found");let t=o[0];const a=i.input.column("Sequence Column",{table:e,value:t,filter:e=>e.semType===l.SEMTYPE.MACROMOLECULE,onValueChanged:e=>{t=e,u()}}),d=i.divV([],{style:{maxHeight:"380px",overflowY:"auto",paddingRight:"8px"}});function u(){d.innerHTML="";const o=(0,r.Ln)(t);if(0!==o.length)for(const n of o){const o=s[n.category]??n.category,a=n.severity?` [${c[n.severity]??n.severity}]`:"",p=n.start&&n.end?` (${n.start}-${n.end})`:"",g=n.sourceScheme?` ${n.sourceScheme}`:"",h=i.iconFA("trash",()=>{const o=(0,r.Ln)(t).filter(e=>e.id!==n.id);(0,r.fh)(t,o),e.fireValuesChanged(),u()});h.style.cursor="pointer",h.style.color="#999",h.style.marginLeft="8px";const f=n.color??"#ccc";let m=f;const y=i.div([],{style:{width:"12px",height:"12px",borderRadius:"2px",backgroundColor:m,display:"inline-block",marginRight:"6px",flexShrink:"0",cursor:"pointer"}});i.colorPicker(l.Color.fromHtml(n.color??"#ccc"),e=>{m=l.Color.toHtml(e)},y,()=>{const o=(0,r.Ln)(t).map(e=>e.id===n.id?{...e,color:m}:e);(0,r.fh)(t,o),e.fireValuesChanged(),u()},()=>{m=f,y.style.backgroundColor=m});const x=i.divH([y,i.divText(`${n.name}${p}${g}${a}`,{style:{flex:"1",fontSize:"12px",padding:"4px"}}),i.divText(o,{style:{color:"#888",fontSize:"11px",marginRight:"8px"}}),h],{style:{alignItems:"center",padding:"4px 0",borderBottom:"1px solid #eee"}});d.append(x)}else d.append(i.divText("No annotations on this column.",{style:{color:"#888",padding:"8px"}}))}u();const p=i.button("Clear All",()=>{(0,r.OW)(e,t),e.fireValuesChanged(),u(),n.shell.info("All annotations cleared")});i.dialog({title:"Manage Annotations"}).add(i.inputs([a])).add(i.h3("Annotations")).add(d).add(i.divH([p],{style:{marginTop:"8px"}})).onOK(()=>{}).show()}}}]);
2
+ //# sourceMappingURL=287.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"287.js","mappings":"uLAUA,MAAMA,EAAyC,CAC7C,CAAC,KAAmBC,WAAY,qBAChC,CAAC,KAAmBC,WAAY,YAChC,CAAC,KAAmBC,KAAM,kCAC1B,CAAC,KAAmBC,QAAS,UAGzBC,EAAyC,CAC7C,CAAC,KAAkBC,MAAO,OAC1B,CAAC,KAAkBC,QAAS,SAC5B,CAAC,KAAkBC,KAAM,MACzB,CAAC,KAAkBC,MAAO,QAGrB,SAASC,IACd,MAAMC,EAAK,QAAWC,IAAIC,UAC1B,IAAKF,EAEH,YADA,QAAWG,QAAQ,iBAIrB,MAAMC,EAAUJ,EAAGK,QAAQC,aAAa,UAAWC,eACnD,GAAuB,IAAnBH,EAAQI,OAEV,YADA,QAAWL,QAAQ,kCAIrB,IAAIM,EAAcL,EAAQ,GAC1B,MAAMM,EAAW,QAASC,OAAO,kBAAmB,CAClDC,MAAOZ,EAAIa,MAAOJ,EAClBK,OAASC,GAAmBA,EAAIC,UAAY,UAAWT,cACvDU,eAAiBF,IAAUN,EAAcM,EAAMG,OAG3CC,EAAU,OAAQ,GAAI,CAACC,MAAO,CAACC,UAAW,QAASC,UAAW,OAAQC,aAAc,SAE1F,SAASL,IACPC,EAAQK,UAAY,GACpB,MAAMC,GAAc,QAAqBhB,GACzC,GAA2B,IAAvBgB,EAAYjB,OAKhB,IAAK,MAAMkB,KAASD,EAAa,CAC/B,MAAME,EAAWtC,EAAeqC,EAAME,WAAaF,EAAME,SACnDC,EAAWH,EAAMI,SAAW,KAAKpC,EAAegC,EAAMI,WAAaJ,EAAMI,YAAc,GACvFC,EAAaL,EAAMM,OAASN,EAAMO,IAAM,KAAKP,EAAMM,SAASN,EAAMO,OAAS,GAC3EC,EAAcR,EAAMS,aAAe,IAAIT,EAAMS,eAAiB,GAE9DC,EAAY,SAAU,QAAS,KACnC,MAAMC,GAAU,QAAqB5B,GAAaK,OAAQwB,GAAMA,EAAEC,KAAOb,EAAMa,KAC/E,QAAqB9B,EAAa4B,GAClCrC,EAAGwC,oBACHtB,MAEFkB,EAAUhB,MAAMqB,OAAS,UACzBL,EAAUhB,MAAMsB,MAAQ,OACxBN,EAAUhB,MAAMuB,WAAa,MAC7B,MAAMC,EAAgBlB,EAAMgB,OAAS,OACrC,IAAIG,EAAeD,EACnB,MAAME,EAAc,MAAO,GAAI,CAAC1B,MAAO,CACrC2B,MAAO,OAAQC,OAAQ,OAAQC,aAAc,MAC7CC,gBAAiBL,EAAcM,QAAS,eAAgBC,YAAa,MACrEC,WAAY,IAAKZ,OAAQ,aAG3B,cAAe,QAASa,SAAS5B,EAAMgB,OAAS,QAAUa,IACxDV,EAAe,QAASW,OAAOD,IAC9BT,EAAa,KACd,MAAMT,GAAU,QAAqB5B,GAAagD,IAAKnB,GAAMA,EAAEC,KAAOb,EAAMa,GAAK,IAAID,EAAGI,MAAOG,GAAgBP,IAC/G,QAAqB7B,EAAa4B,GAClCrC,EAAGwC,oBACHtB,KACC,KACD2B,EAAeD,EACfE,EAAY1B,MAAM8B,gBAAkBL,IAGtC,MAAMa,EAAM,OAAQ,CAClBZ,EACA,UAAW,GAAGpB,EAAMiC,OAAO5B,IAAaG,IAAcL,IAAY,CAACT,MAAO,CAACwC,KAAM,IAAKC,SAAU,OAAQC,QAAS,SACjH,UAAWnC,EAAU,CAACP,MAAO,CAACsB,MAAO,OAAQmB,SAAU,OAAQT,YAAa,SAC5EhB,GACC,CAAChB,MAAO,CAAC2C,WAAY,SAAUD,QAAS,QAASE,aAAc,oBAElE7C,EAAQ8C,OAAOP,EACjB,MA/CEvC,EAAQ8C,OAAO,UAAW,iCAAkC,CAAC7C,MAAO,CAACsB,MAAO,OAAQoB,QAAS,SAgDjG,CAEA5C,IAEA,MAAMgD,EAAW,SAAU,YAAa,MACtC,QAAiBlE,EAAIS,GACrBT,EAAGwC,oBACHtB,IACA,QAAWiD,KAAK,6BAGH,SAAU,CAACC,MAAO,uBAC9BC,IAAI,SAAU,CAAC3D,KACf2D,IAAI,KAAM,gBACVA,IAAIlD,GACJkD,IAAI,OAAQ,CAACH,GAAW,CAAC9C,MAAO,CAACkD,UAAW,UAC5CC,KAAK,QAEDC,MACT,C","sources":["webpack://bio/./src/utils/annotations/annotation-manager-ui.ts"],"sourcesContent":["/* eslint-disable max-len */\nimport * as grok from 'datagrok-api/grok';\nimport * as ui from 'datagrok-api/ui';\nimport * as DG from 'datagrok-api/dg';\n\nimport {\n SeqAnnotation, AnnotationCategory, LiabilitySeverity,\n} from '@datagrok-libraries/bio/src/utils/macromolecule/annotations';\nimport {getColumnAnnotations, setColumnAnnotations, clearAnnotations} from './annotation-manager';\n\nconst categoryLabels: Record<string, string> = {\n [AnnotationCategory.Structure]: 'Structure (FR/CDR)',\n [AnnotationCategory.Liability]: 'Liability',\n [AnnotationCategory.PTM]: 'Post-translational Modification',\n [AnnotationCategory.Custom]: 'Custom',\n};\n\nconst severityLabels: Record<string, string> = {\n [LiabilitySeverity.High]: 'High',\n [LiabilitySeverity.Medium]: 'Medium',\n [LiabilitySeverity.Low]: 'Low',\n [LiabilitySeverity.Info]: 'Info',\n};\n\nexport function showAnnotationManagerDialog(): void {\n const df = grok.shell.tv?.dataFrame;\n if (!df) {\n grok.shell.warning('No table open');\n return;\n }\n\n const seqCols = df.columns.bySemTypeAll(DG.SEMTYPE.MACROMOLECULE);\n if (seqCols.length === 0) {\n grok.shell.warning('No macromolecule columns found');\n return;\n }\n\n let selectedCol = seqCols[0];\n const colInput = ui.input.column('Sequence Column', {\n table: df, value: selectedCol,\n filter: (col: DG.Column) => col.semType === DG.SEMTYPE.MACROMOLECULE,\n onValueChanged: (col) => { selectedCol = col!; refreshList(); },\n });\n\n const listDiv = ui.divV([], {style: {maxHeight: '380px', overflowY: 'auto', paddingRight: '8px'}});\n\n function refreshList(): void {\n listDiv.innerHTML = '';\n const annotations = getColumnAnnotations(selectedCol);\n if (annotations.length === 0) {\n listDiv.append(ui.divText('No annotations on this column.', {style: {color: '#888', padding: '8px'}}));\n return;\n }\n\n for (const annot of annotations) {\n const catLabel = categoryLabels[annot.category] ?? annot.category;\n const sevLabel = annot.severity ? ` [${severityLabels[annot.severity] ?? annot.severity}]` : '';\n const rangeLabel = annot.start && annot.end ? ` (${annot.start}-${annot.end})` : '';\n const schemeLabel = annot.sourceScheme ? ` ${annot.sourceScheme}` : '';\n\n const removeBtn = ui.iconFA('trash', () => {\n const updated = getColumnAnnotations(selectedCol).filter((a) => a.id !== annot.id);\n setColumnAnnotations(selectedCol, updated);\n df.fireValuesChanged();\n refreshList();\n });\n removeBtn.style.cursor = 'pointer';\n removeBtn.style.color = '#999';\n removeBtn.style.marginLeft = '8px';\n const originalColor = annot.color ?? '#ccc';\n let currentColor = originalColor;\n const colorSwatch = ui.div([], {style: {\n width: '12px', height: '12px', borderRadius: '2px',\n backgroundColor: currentColor, display: 'inline-block', marginRight: '6px',\n flexShrink: '0', cursor: 'pointer',\n }});\n\n ui.colorPicker(DG.Color.fromHtml(annot.color ?? '#ccc'), (newColor) => {\n currentColor = DG.Color.toHtml(newColor);\n }, colorSwatch, () => {\n const updated = getColumnAnnotations(selectedCol).map((a) => a.id === annot.id ? {...a, color: currentColor} : a);\n setColumnAnnotations(selectedCol, updated);\n df.fireValuesChanged();\n refreshList();\n }, () => {\n currentColor = originalColor;\n colorSwatch.style.backgroundColor = currentColor;\n });\n\n const row = ui.divH([\n colorSwatch,\n ui.divText(`${annot.name}${rangeLabel}${schemeLabel}${sevLabel}`, {style: {flex: '1', fontSize: '12px', padding: '4px'}}),\n ui.divText(catLabel, {style: {color: '#888', fontSize: '11px', marginRight: '8px'}}),\n removeBtn,\n ], {style: {alignItems: 'center', padding: '4px 0', borderBottom: '1px solid #eee'}});\n\n listDiv.append(row);\n }\n }\n\n refreshList();\n\n const clearBtn = ui.button('Clear All', () => {\n clearAnnotations(df, selectedCol);\n df.fireValuesChanged();\n refreshList();\n grok.shell.info('All annotations cleared');\n });\n\n const dialog = ui.dialog({title: 'Manage Annotations'})\n .add(ui.inputs([colInput]))\n .add(ui.h3('Annotations'))\n .add(listDiv)\n .add(ui.divH([clearBtn], {style: {marginTop: '8px'}}))\n .onOK(() => {});\n\n dialog.show();\n}\n"],"names":["categoryLabels","Structure","Liability","PTM","Custom","severityLabels","High","Medium","Low","Info","showAnnotationManagerDialog","df","tv","dataFrame","warning","seqCols","columns","bySemTypeAll","MACROMOLECULE","length","selectedCol","colInput","column","table","value","filter","col","semType","onValueChanged","refreshList","listDiv","style","maxHeight","overflowY","paddingRight","innerHTML","annotations","annot","catLabel","category","sevLabel","severity","rangeLabel","start","end","schemeLabel","sourceScheme","removeBtn","updated","a","id","fireValuesChanged","cursor","color","marginLeft","originalColor","currentColor","colorSwatch","width","height","borderRadius","backgroundColor","display","marginRight","flexShrink","fromHtml","newColor","toHtml","map","row","name","flex","fontSize","padding","alignItems","borderBottom","append","clearBtn","info","title","add","marginTop","onOK","show"],"sourceRoot":""}
package/dist/288.js ADDED
@@ -0,0 +1,2 @@
1
+ "use strict";(self.webpackChunkbio=self.webpackChunkbio||[]).push([[288],{9288(n,t,e){e.d(t,{extractAnnotatedRegion:()=>r,filterByLiabilityHits:()=>a});var o=e(4328),i=e(6082),l=e(4517),s=e(3736);function a(n,t){const e=(0,s.VT)(t.name);let l=null;try{l=n.columns.byName(e)}catch{}if(!l)return void o.shell.warning("No annotation data found. Run liability scanning first.");const a=(0,s.wi)(l),r=i.BitSet.create(n.rowCount);for(let t=0;t<n.rowCount;t++)a[t]&&a[t].length>0&&r.set(t,!0);n.filter.copyFrom(r),o.shell.info(`Filtered to ${r.trueCount} rows with liability hits`)}function r(n,t,e,a){const r=(0,s.Ln)(t).find(n=>n.name===e&&n.category===l.eI.Structure);if(!r)return o.shell.warning(`Region annotation "${e}" not found.`),null;const u=a.getSeqHandler(t),d=`${t.name}(${e})`,c=(0,s.VT)(t.name);let f=null;try{f=n.columns.byName(c)}catch{}if(f){const t=(0,s.wi)(f);if(t.some(n=>n?.some(n=>n.annotationId===r.id&&null!=n.endPositionIndex))){const l=i.Column.fromType(i.COLUMN_TYPE.STRING,d,n.rowCount);for(let e=0;e<n.rowCount;e++){const n=t[e],o=n?.find(n=>n.annotationId===r.id&&null!=n.endPositionIndex);if(o){const n=u.getSplitted(e),t=[];for(let e=o.positionIndex;e<=o.endPositionIndex;e++)e<n.length&&t.push(n.getOriginal(e));l.set(e,t.join(u.separator||""))}else l.set(e,"")}return n.columns.add(l),o.data.detectSemanticTypes(n),o.shell.info(`Extracted region ${e} as column "${d}"`),l}}if(null==r.start||null==r.end)return o.shell.warning(`Region annotation "${e}" has no position range.`),null;const g=u.posList.indexOf(r.start),m=u.posList.indexOf(r.end);if(g<0||m<0)return o.shell.warning(`Position names "${r.start}" or "${r.end}" not found in position list.`),null;const h=u.getRegion(g,m,d);return n.columns.add(h),o.data.detectSemanticTypes(n),o.shell.info(`Extracted region ${e} as column "${d}"`),h}}}]);
2
+ //# sourceMappingURL=288.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"288.js","mappings":"oMAQO,SAASA,EAAsBC,EAAkBC,GACtD,MAAMC,GAAe,QAAwBD,EAAOE,MACpD,IAAIC,EAAqC,KACzC,IACEA,EAAWJ,EAAGK,QAAQC,OAAOJ,EAC/B,CAAE,MAAwB,CAE1B,IAAKE,EAEH,YADA,QAAWG,QAAQ,2DAIrB,MAAMC,GAAU,QAAuBJ,GACjCK,EAAK,SAAUC,OAAOV,EAAGW,UAC/B,IAAK,IAAIC,EAAI,EAAGA,EAAIZ,EAAGW,SAAUC,IAC3BJ,EAAQI,IAAMJ,EAAQI,GAAIC,OAAS,GACrCJ,EAAGK,IAAIF,GAAG,GAGdZ,EAAGe,OAAOC,SAASP,GACnB,QAAWQ,KAAK,eAAeR,EAAGS,qCACpC,CA6BO,SAASC,EACdnB,EACAC,EACAmB,EACAC,GAEA,MACMC,GADc,QAAqBrB,GACfsB,KAAMC,GAC9BA,EAAErB,OAASiB,GAAkBI,EAAEC,WAAa,KAAmBC,WAEjE,IAAKJ,EAEH,OADA,QAAWf,QAAQ,sBAAsBa,iBAClC,KAGT,MAAMO,EAAKN,EAAUO,cAAc3B,GAC7B4B,EAAU,GAAG5B,EAAOE,QAAQiB,KAG5BlB,GAAe,QAAwBD,EAAOE,MACpD,IAAIC,EAAqC,KACzC,IAAMA,EAAWJ,EAAGK,QAAQC,OAAOJ,EAAoC,CAAE,MAAwB,CAEjG,GAAIE,EAAU,CACZ,MAAM0B,GAAa,QAAuB1B,GAI1C,GAHyB0B,EAAWC,KAAMC,GACxCA,GAAID,KAAME,GAAwBA,EAAEC,eAAiBZ,EAAMa,IAA4B,MAAtBF,EAAEG,mBAE/C,CACpB,MAAMC,EAAS,SAAUC,SAAS,cAAeC,OAAQV,EAAS7B,EAAGW,UACrE,IAAK,IAAIC,EAAI,EAAGA,EAAIZ,EAAGW,SAAUC,IAAK,CACpC,MAAM4B,EAAUV,EAAWlB,GACrB6B,EAAYD,GAASjB,KAAMU,GAC/BA,EAAEC,eAAiBZ,EAAMa,IAA4B,MAAtBF,EAAEG,kBACnC,GAAIK,EAAW,CACb,MAAMC,EAAWf,EAAGgB,YAAY/B,GAC1BgC,EAAkB,GACxB,IAAK,IAAIC,EAAIJ,EAAUK,cAAeD,GAAKJ,EAAUL,iBAAmBS,IAClEA,EAAIH,EAAS7B,QACf+B,EAAMG,KAAKL,EAASM,YAAYH,IAEpCR,EAAOvB,IAAIF,EAAGgC,EAAMK,KAAKtB,EAAGuB,WAAa,IAC3C,MACEb,EAAOvB,IAAIF,EAAG,GAClB,CAIA,OAHAZ,EAAGK,QAAQ8C,IAAId,GACf,OAAUe,oBAAoBpD,GAC9B,QAAWiB,KAAK,oBAAoBG,gBAA6BS,MAC1DQ,CACT,CACF,CAGA,GAAmB,MAAff,EAAM+B,OAA8B,MAAb/B,EAAMgC,IAE/B,OADA,QAAW/C,QAAQ,sBAAsBa,6BAClC,KAGT,MAAMmC,EAAW5B,EAAG6B,QAAQC,QAAQnC,EAAM+B,OACpCK,EAAS/B,EAAG6B,QAAQC,QAAQnC,EAAMgC,KAExC,GAAIC,EAAW,GAAKG,EAAS,EAE3B,OADA,QAAWnD,QAAQ,mBAAmBe,EAAM+B,cAAc/B,EAAMgC,oCACzD,KAGT,MAAMjB,EAASV,EAAGgC,UAAUJ,EAAUG,EAAQ7B,GAI9C,OAHA7B,EAAGK,QAAQ8C,IAAId,GACf,OAAUe,oBAAoBpD,GAC9B,QAAWiB,KAAK,oBAAoBG,gBAA6BS,MAC1DQ,CACT,C","sources":["webpack://bio/./src/utils/annotations/annotation-actions.ts"],"sourcesContent":["import * as grok from 'datagrok-api/grok';\nimport * as DG from 'datagrok-api/dg';\n\nimport {AnnotationCategory, SeqAnnotationHit} from '@datagrok-libraries/bio/src/utils/macromolecule/annotations';\nimport {ISeqHelper} from '@datagrok-libraries/bio/src/utils/seq-helper';\nimport {getAnnotationColumnName, getColumnAnnotations, cacheAllRowAnnotations} from './annotation-manager';\n\n/** Filters the DataFrame to show only rows that have at least one liability hit. */\nexport function filterByLiabilityHits(df: DG.DataFrame, seqCol: DG.Column<string>): void {\n const annotColName = getAnnotationColumnName(seqCol.name);\n let annotCol: DG.Column<string> | null = null;\n try {\n annotCol = df.columns.byName(annotColName) as DG.Column<string>;\n } catch { /* not found */ }\n\n if (!annotCol) {\n grok.shell.warning('No annotation data found. Run liability scanning first.');\n return;\n }\n\n const rowData = cacheAllRowAnnotations(annotCol);\n const bs = DG.BitSet.create(df.rowCount);\n for (let i = 0; i < df.rowCount; i++) {\n if (rowData[i] && rowData[i]!.length > 0)\n bs.set(i, true);\n }\n\n df.filter.copyFrom(bs);\n grok.shell.info(`Filtered to ${bs.trueCount} rows with liability hits`);\n}\n\n/** Selects all rows that contain a specific annotation hit. */\nexport function selectRowsWithAnnotation(df: DG.DataFrame, seqCol: DG.Column<string>, annotationId: string): void {\n const annotColName = getAnnotationColumnName(seqCol.name);\n let annotCol: DG.Column<string> | null = null;\n try {\n annotCol = df.columns.byName(annotColName) as DG.Column<string>;\n } catch { /* not found */ }\n\n if (!annotCol) {\n grok.shell.warning('No annotation data found.');\n return;\n }\n\n const rowData = cacheAllRowAnnotations(annotCol);\n const bs = DG.BitSet.create(df.rowCount);\n for (let i = 0; i < df.rowCount; i++) {\n if (rowData[i]?.some((h) => h.annotationId === annotationId))\n bs.set(i, true);\n }\n\n df.selection.copyFrom(bs);\n grok.shell.info(`Selected ${bs.trueCount} rows with ${annotationId} hits`);\n}\n\n/** Extracts a named region annotation as a new column.\n * Uses per-row region spans from the companion column when available (unaligned data),\n * falls back to column-level position names (aligned/MSA data). */\nexport function extractAnnotatedRegion(\n df: DG.DataFrame,\n seqCol: DG.Column<string>,\n annotationName: string,\n seqHelper: ISeqHelper,\n): DG.Column<string> | null {\n const annotations = getColumnAnnotations(seqCol);\n const annot = annotations.find((a) =>\n a.name === annotationName && a.category === AnnotationCategory.Structure);\n\n if (!annot) {\n grok.shell.warning(`Region annotation \"${annotationName}\" not found.`);\n return null;\n }\n\n const sh = seqHelper.getSeqHandler(seqCol);\n const colName = `${seqCol.name}(${annotationName})`;\n\n // Try per-row extraction using companion column region spans\n const annotColName = getAnnotationColumnName(seqCol.name);\n let annotCol: DG.Column<string> | null = null;\n try { annotCol = df.columns.byName(annotColName) as DG.Column<string>; } catch { /* not found */ }\n\n if (annotCol) {\n const allRowData = cacheAllRowAnnotations(annotCol);\n const hasPerRowRegions = allRowData.some((rd) =>\n rd?.some((h: SeqAnnotationHit) => h.annotationId === annot.id && h.endPositionIndex != null));\n\n if (hasPerRowRegions) {\n const regCol = DG.Column.fromType(DG.COLUMN_TYPE.STRING, colName, df.rowCount);\n for (let i = 0; i < df.rowCount; i++) {\n const rowHits = allRowData[i];\n const regionHit = rowHits?.find((h: SeqAnnotationHit) =>\n h.annotationId === annot.id && h.endPositionIndex != null);\n if (regionHit) {\n const splitted = sh.getSplitted(i);\n const parts: string[] = [];\n for (let p = regionHit.positionIndex; p <= regionHit.endPositionIndex!; p++) {\n if (p < splitted.length)\n parts.push(splitted.getOriginal(p));\n }\n regCol.set(i, parts.join(sh.separator || ''));\n } else\n regCol.set(i, '');\n }\n df.columns.add(regCol);\n grok.data.detectSemanticTypes(df);\n grok.shell.info(`Extracted region ${annotationName} as column \"${colName}\"`);\n return regCol;\n }\n }\n\n // Fall back to column-level position names (aligned/MSA data)\n if (annot.start == null || annot.end == null) {\n grok.shell.warning(`Region annotation \"${annotationName}\" has no position range.`);\n return null;\n }\n\n const startIdx = sh.posList.indexOf(annot.start);\n const endIdx = sh.posList.indexOf(annot.end);\n\n if (startIdx < 0 || endIdx < 0) {\n grok.shell.warning(`Position names \"${annot.start}\" or \"${annot.end}\" not found in position list.`);\n return null;\n }\n\n const regCol = sh.getRegion(startIdx, endIdx, colName);\n df.columns.add(regCol);\n grok.data.detectSemanticTypes(df);\n grok.shell.info(`Extracted region ${annotationName} as column \"${colName}\"`);\n return regCol;\n}\n"],"names":["filterByLiabilityHits","df","seqCol","annotColName","name","annotCol","columns","byName","warning","rowData","bs","create","rowCount","i","length","set","filter","copyFrom","info","trueCount","extractAnnotatedRegion","annotationName","seqHelper","annot","find","a","category","Structure","sh","getSeqHandler","colName","allRowData","some","rd","h","annotationId","id","endPositionIndex","regCol","fromType","STRING","rowHits","regionHit","splitted","getSplitted","parts","p","positionIndex","push","getOriginal","join","separator","add","detectSemanticTypes","start","end","startIdx","posList","indexOf","endIdx","getRegion"],"sourceRoot":""}
package/dist/422.js ADDED
@@ -0,0 +1,2 @@
1
+ "use strict";(self.webpackChunkbio=self.webpackChunkbio||[]).push([[422],{9422(e,t,i){i.d(t,{showLiabilityScannerDialog:()=>c});var n=i(4328),a=i(7389),o=i(6082),l=i(3785),r=i(4517),d=i(3736);const s=[{id:"deamid-ng",name:"Deamidation (NG)",pattern:/NG/g,length:2,severity:r.Hq.High,ruleCategory:"deamidation",color:r.D5.liability.deamidation,enabled:!0},{id:"deamid-ns",name:"Deamidation (NS)",pattern:/NS/g,length:2,severity:r.Hq.Medium,ruleCategory:"deamidation",color:r.D5.liability.deamidation,enabled:!0},{id:"deamid-na",name:"Deamidation (NA)",pattern:/NA/g,length:2,severity:r.Hq.Low,ruleCategory:"deamidation",color:r.D5.liability.deamidation,enabled:!0},{id:"deamid-nd",name:"Deamidation (ND)",pattern:/ND/g,length:2,severity:r.Hq.Low,ruleCategory:"deamidation",color:r.D5.liability.deamidation,enabled:!0},{id:"deamid-nt",name:"Deamidation (NT)",pattern:/NT/g,length:2,severity:r.Hq.Low,ruleCategory:"deamidation",color:r.D5.liability.deamidation,enabled:!0},{id:"isom-dg",name:"Isomerization (DG)",pattern:/DG/g,length:2,severity:r.Hq.High,ruleCategory:"isomerization",color:r.D5.liability.isomerization,enabled:!0},{id:"isom-ds",name:"Isomerization (DS)",pattern:/DS/g,length:2,severity:r.Hq.Medium,ruleCategory:"isomerization",color:r.D5.liability.isomerization,enabled:!0},{id:"oxid-m",name:"Oxidation (Met)",pattern:/M/g,length:1,severity:r.Hq.Medium,ruleCategory:"oxidation",color:r.D5.liability.oxidation,enabled:!0},{id:"oxid-w",name:"Oxidation (Trp)",pattern:/W/g,length:1,severity:r.Hq.Low,ruleCategory:"oxidation",color:r.D5.liability.oxidation,enabled:!0},{id:"glyco-nxst",name:"N-glycosylation",pattern:/N[^P][ST]/g,length:3,severity:r.Hq.High,ruleCategory:"glycosylation",color:r.D5.liability.glycosylation,enabled:!0},{id:"free-cys",name:"Free Cysteine",pattern:/C/g,length:1,severity:r.Hq.Info,ruleCategory:"freeCysteine",color:r.D5.liability.freeCysteine,enabled:!1}];function u(e,t){const i=e.getSplitted(t),n=new Array(i.length);for(let e=0;e<i.length;e++)n[e]=i.getOriginal(e);return n.join("")}const g={[r.Hq.High]:"High",[r.Hq.Medium]:"Medium",[r.Hq.Low]:"Low",[r.Hq.Info]:"Info"};function c(){const e=n.shell.tv?.dataFrame;if(!e)return void n.shell.warning("No table open");const t=e.columns.bySemTypeAll(o.SEMTYPE.MACROMOLECULE);if(0===t.length)return void n.shell.warning("No macromolecule columns found");const i=s.map(e=>({...e,pattern:new RegExp(e.pattern.source,"g")})),c=a.input.table("Table",{value:e}),m=a.input.column("Sequence",{table:e,value:t[0],filter:e=>e.semType===o.SEMTYPE.MACROMOLECULE}),y=[],p=a.divV([]);for(const e of i){const t=a.input.bool(e.name,{value:e.enabled,tooltipText:`Severity: ${g[e.severity]??e.severity}`});y.push({rule:e,input:t}),p.append(t.root)}const h=a.input.bool("Highlight in cell renderer",{value:!0}),b=a.input.bool("Create annotation column",{value:!0}),f=a.input.bool("Create summary count column",{value:!1}),v=a.dialog({title:"Scan Sequence Liabilities"}).add(a.inputs([c,m])).add(a.h3("Rules")).add(p).add(a.h3("Output")).add(a.inputs([h,b,f])).onOK(()=>{try{const t=m.value,a=l._package.seqHelper.getSeqHandler(t);for(const{rule:e,input:t}of y)e.enabled=t.value;const o=function(e,t,i){const n=i.filter(e=>e.enabled),a=t.posList,o=new Map,l=new Array(e.length);let d=0;for(let i=0;i<e.length;i++){const e=u(t,i),r=[];for(const t of n){let i;for(t.pattern.lastIndex=0;null!==(i=t.pattern.exec(e));)r.push({annotationId:t.id,positionIndex:i.index,positionName:i.index<a.length?a[i.index]:void 0,matchedMonomers:i[0]}),o.set(t.id,(o.get(t.id)??0)+1),d++}l[i]=r}return{annotations:n.filter(e=>o.has(e.id)).map(e=>({id:e.id,name:e.name,description:`${e.ruleCategory} liability pattern (${o.get(e.id)} hits)`,start:null,end:null,visualType:1===e.length?r.ao.Point:r.ao.Motif,category:r.eI.Liability,color:e.color,severity:e.severity,motifPattern:e.pattern.source,autoGenerated:!0})),rowData:l,totalHits:d}}(t,a,i);(b.value||h.value)&&function(e,t,i){const n=(0,d.Ln)(t).filter(e=>e.category!==r.eI.Liability);(0,d.fh)(t,[...n,...i.annotations]);const a=(0,d.Lz)(e,t);for(let e=0;e<i.rowData.length;e++){const t=(0,d.JG)(a,e)??[];(0,d.z5)(a,e,(0,d.bq)(t,i.rowData[e],!1,!0))}}(e,t,o),f.value&&function(e,t,i){const n=`${t.name}_liability_count`,a=i.rowData.map(e=>e.length),o=e.columns.addNewInt(n);for(let e=0;e<a.length;e++)o.set(e,a[e])}(e,t,o),n.shell.info(`Liability scan: ${o.totalHits} hits found across ${o.annotations.length} rules`),e.fireValuesChanged()}catch(e){n.shell.error(`Liability scan failed: ${e.message??e}`),console.error(e)}});v.show()}}}]);
2
+ //# sourceMappingURL=422.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"422.js","mappings":"gMA4BO,MAAMA,EAA2C,CACtD,CAACC,GAAI,YAAaC,KAAM,mBAAoBC,QAAS,MAAOC,OAAQ,EAAGC,SAAU,KAAkBC,KAAMC,aAAc,cAAeC,MAAO,KAAkBC,UAAUC,YAAaC,SAAS,GAC/L,CAACV,GAAI,YAAaC,KAAM,mBAAoBC,QAAS,MAAOC,OAAQ,EAAGC,SAAU,KAAkBO,OAAQL,aAAc,cAAeC,MAAO,KAAkBC,UAAUC,YAAaC,SAAS,GACjM,CAACV,GAAI,YAAaC,KAAM,mBAAoBC,QAAS,MAAOC,OAAQ,EAAGC,SAAU,KAAkBQ,IAAKN,aAAc,cAAeC,MAAO,KAAkBC,UAAUC,YAAaC,SAAS,GAC9L,CAACV,GAAI,YAAaC,KAAM,mBAAoBC,QAAS,MAAOC,OAAQ,EAAGC,SAAU,KAAkBQ,IAAKN,aAAc,cAAeC,MAAO,KAAkBC,UAAUC,YAAaC,SAAS,GAC9L,CAACV,GAAI,YAAaC,KAAM,mBAAoBC,QAAS,MAAOC,OAAQ,EAAGC,SAAU,KAAkBQ,IAAKN,aAAc,cAAeC,MAAO,KAAkBC,UAAUC,YAAaC,SAAS,GAC9L,CAACV,GAAI,UAAWC,KAAM,qBAAsBC,QAAS,MAAOC,OAAQ,EAAGC,SAAU,KAAkBC,KAAMC,aAAc,gBAAiBC,MAAO,KAAkBC,UAAUK,cAAeH,SAAS,GACnM,CAACV,GAAI,UAAWC,KAAM,qBAAsBC,QAAS,MAAOC,OAAQ,EAAGC,SAAU,KAAkBO,OAAQL,aAAc,gBAAiBC,MAAO,KAAkBC,UAAUK,cAAeH,SAAS,GACrM,CAACV,GAAI,SAAUC,KAAM,kBAAmBC,QAAS,KAAMC,OAAQ,EAAGC,SAAU,KAAkBO,OAAQL,aAAc,YAAaC,MAAO,KAAkBC,UAAUM,UAAWJ,SAAS,GACxL,CAACV,GAAI,SAAUC,KAAM,kBAAmBC,QAAS,KAAMC,OAAQ,EAAGC,SAAU,KAAkBQ,IAAKN,aAAc,YAAaC,MAAO,KAAkBC,UAAUM,UAAWJ,SAAS,GACrL,CAACV,GAAI,aAAcC,KAAM,kBAAmBC,QAAS,aAAcC,OAAQ,EAAGC,SAAU,KAAkBC,KAAMC,aAAc,gBAAiBC,MAAO,KAAkBC,UAAUO,cAAeL,SAAS,GAC1M,CAACV,GAAI,WAAYC,KAAM,gBAAiBC,QAAS,KAAMC,OAAQ,EAAGC,SAAU,KAAkBY,KAAMV,aAAc,eAAgBC,MAAO,KAAkBC,UAAUS,aAAcP,SAAS,IAI9L,SAASQ,EAAmBC,EAAiBC,GAC3C,MAAMC,EAAWF,EAAGG,YAAYF,GAC1BG,EAAkB,IAAIC,MAAMH,EAASlB,QAC3C,IAAK,IAAIsB,EAAI,EAAGA,EAAIJ,EAASlB,OAAQsB,IACnCF,EAAME,GAAKJ,EAASK,YAAYD,GAClC,OAAOF,EAAMI,KAAK,GACpB,CCtCA,MAAMC,EAAyC,CAC7C,CAAC,KAAkBvB,MAAO,OAC1B,CAAC,KAAkBM,QAAS,SAC5B,CAAC,KAAkBC,KAAM,MACzB,CAAC,KAAkBI,MAAO,QAGrB,SAASa,IACd,MAAMC,EAAK,QAAWC,IAAIC,UAC1B,IAAKF,EAEH,YADA,QAAWG,QAAQ,iBAIrB,MAAMC,EAAUJ,EAAGK,QAAQC,aAAa,UAAWC,eACnD,GAAuB,IAAnBH,EAAQ/B,OAEV,YADA,QAAW8B,QAAQ,kCAIrB,MAAMK,EAAQvC,EAAwBwC,IAAKC,IAAM,IAAKA,EAAGtC,QAAS,IAAIuC,OAAOD,EAAEtC,QAAQwC,OAAQ,QAEzFC,EAAa,QAASC,MAAM,QAAS,CAACC,MAAOf,IAC7CgB,EAAW,QAASC,OAAO,WAAY,CAC3CH,MAAOd,EAAIe,MAAOX,EAAQ,GAC1Bc,OAASC,GAAmBA,EAAIC,UAAY,UAAWb,gBAInDc,EAAoE,GACpEC,EAAW,OAAQ,IACzB,IAAK,MAAMC,KAAQf,EAAO,CACxB,MAAMgB,EAAQ,QAASC,KAAKF,EAAKpD,KAAM,CACrC4C,MAAOQ,EAAK3C,QACZ8C,YAAa,aAAa5B,EAAeyB,EAAKjD,WAAaiD,EAAKjD,aAElE+C,EAAWM,KAAK,CAACJ,OAAMK,MAAOJ,IAC9BF,EAASO,OAAOL,EAAMM,KACxB,CAEA,MAAMC,EAAiB,QAASN,KAAK,6BAA8B,CAACV,OAAO,IACrEiB,EAAgB,QAASP,KAAK,2BAA4B,CAACV,OAAO,IAClEkB,EAAe,QAASR,KAAK,8BAA+B,CAACV,OAAO,IAEpEmB,EAAS,SAAU,CAACC,MAAO,8BAC9BC,IAAI,SAAU,CAACvB,EAAYG,KAC3BoB,IAAI,KAAM,UACVA,IAAId,GACJc,IAAI,KAAM,WACVA,IAAI,SAAU,CAACL,EAAgBC,EAAeC,KAC9CI,KAAK,KACJ,IACE,MAAMC,EAAStB,EAASD,MAClB1B,EAAK,EAAAkD,SAASC,UAAUC,cAAcH,GAG5C,IAAK,MAAM,KAACf,EAAI,MAAEK,KAAUP,EAC1BE,EAAK3C,QAAUgD,EAAMb,MAEvB,MAAM2B,EDXP,SACLvB,EACA9B,EACAmB,GAEA,MAAMmC,EAAenC,EAAMU,OAAQR,GAAMA,EAAE9B,SACrCgE,EAAUvD,EAAGuD,QAGbC,EAAgB,IAAIC,IAEpBC,EAA+B,IAAIrD,MAAMyB,EAAI9C,QACnD,IAAI2E,EAAY,EAEhB,IAAK,IAAI1D,EAAS,EAAGA,EAAS6B,EAAI9C,OAAQiB,IAAU,CAClD,MAAM2D,EAAM7D,EAAmBC,EAAIC,GAC7B4D,EAA2B,GAEjC,IAAK,MAAM3B,KAAQoB,EAAc,CAG/B,IAAIQ,EACJ,IAFA5B,EAAKnD,QAAQgF,UAAY,EAEmB,QAApCD,EAAQ5B,EAAKnD,QAAQiF,KAAKJ,KAChCC,EAAKvB,KAAK,CACR2B,aAAc/B,EAAKrD,GACnBqF,cAAeJ,EAAMK,MACrBC,aAAcN,EAAMK,MAAQZ,EAAQvE,OAASuE,EAAQO,EAAMK,YAASE,EACpEC,gBAAiBR,EAAM,KAEzBN,EAAce,IAAIrC,EAAKrD,IAAK2E,EAAcgB,IAAItC,EAAKrD,KAAO,GAAK,GAC/D8E,GAEJ,CACAD,EAAQzD,GAAU4D,CACpB,CAmBA,MAAO,CAACY,YAhB6BnB,EAClCzB,OAAQR,GAAMmC,EAAckB,IAAIrD,EAAExC,KAClCuC,IAAKC,IAAM,CACVxC,GAAIwC,EAAExC,GACNC,KAAMuC,EAAEvC,KACR6F,YAAa,GAAGtD,EAAElC,mCAAmCqE,EAAcgB,IAAInD,EAAExC,YACzE+F,MAAO,KACPC,IAAK,KACLC,WAAyB,IAAbzD,EAAErC,OAAe,KAAqB+F,MAAQ,KAAqBC,MAC/EC,SAAU,KAAmBC,UAC7B9F,MAAOiC,EAAEjC,MACTH,SAAUoC,EAAEpC,SACZkG,aAAc9D,EAAEtC,QAAQwC,OACxB6D,eAAe,KAGE1B,UAASC,YAChC,CC3CuB0B,CAAgBpC,EAAQjD,EAAImB,IAEvCwB,EAAcjB,OAASgB,EAAehB,QD4C3C,SACLf,EACAsC,EACAI,GAGA,MAAMiC,GAAW,QAAqBrC,GACnCpB,OAAQ0D,GAAMA,EAAEN,WAAa,KAAmBC,YACnD,QAAqBjC,EAAQ,IAAIqC,KAAajC,EAAOoB,cAGrD,MAAMe,GAAW,QAA4B7E,EAAIsC,GACjD,IAAK,IAAI3C,EAAI,EAAGA,EAAI+C,EAAOK,QAAQ1E,OAAQsB,IAAK,CAC9C,MAAMmF,GAAe,QAAkBD,EAAUlF,IAAM,IACvD,QAAkBkF,EAAUlF,GAAG,QAAamF,EAAcpC,EAAOK,QAAQpD,IAAI,GAAO,GACtF,CACF,CC3DUoF,CAA0B/E,EAAIsC,EAAQI,GAEpCT,EAAalB,OD4DlB,SACLf,EACAsC,EACAI,GAEA,MAAMsC,EAAU,GAAG1C,EAAOnE,uBACpB8G,EAASvC,EAAOK,QAAQtC,IAAKyC,GAASA,EAAK7E,QAC3C8C,EAAMnB,EAAGK,QAAQ6E,UAAUF,GACjC,IAAK,IAAIrF,EAAI,EAAGA,EAAIsF,EAAO5G,OAAQsB,IACjCwB,EAAIyC,IAAIjE,EAAGsF,EAAOtF,GAEtB,CCtEUwF,CAA6BnF,EAAIsC,EAAQI,GAE3C,QAAW0C,KAAK,mBAAmB1C,EAAOM,+BAA+BN,EAAOoB,YAAYzF,gBAC5F2B,EAAGqF,mBACL,CAAE,MAAOC,GACP,QAAWC,MAAM,0BAA0BD,EAAIE,SAAWF,KAC1DG,QAAQF,MAAMD,EAChB,IAGJpD,EAAOwD,MACT,C","sources":["webpack://bio/./src/utils/annotations/liability-scanner.ts","webpack://bio/./src/utils/annotations/liability-scanner-ui.ts"],"sourcesContent":["/* eslint-disable max-len */\nimport * as DG from 'datagrok-api/dg';\n\nimport {ISeqHandler} from '@datagrok-libraries/bio/src/utils/macromolecule/seq-handler';\nimport {\n SeqAnnotation, SeqAnnotationHit, RowAnnotationData,\n AnnotationVisualType, AnnotationCategory, LiabilitySeverity,\n ANNOTATION_COLORS,\n} from '@datagrok-libraries/bio/src/utils/macromolecule/annotations';\nimport {\n getOrCreateAnnotationColumn, setColumnAnnotations, setRowAnnotations,\n getColumnAnnotations, getRowAnnotations, mergeRowHits,\n} from './annotation-manager';\n\n/** A single liability scanning rule. */\nexport interface LiabilityRule {\n id: string;\n name: string;\n pattern: RegExp;\n length: number;\n severity: LiabilitySeverity;\n /** Sub-category for grouping (e.g. \"deamidation\", \"oxidation\") */\n ruleCategory: string;\n color: string;\n enabled: boolean;\n}\n\n/** Built-in liability rules for antibody engineering. */\nexport const BUILTIN_LIABILITY_RULES: LiabilityRule[] = [\n {id: 'deamid-ng', name: 'Deamidation (NG)', pattern: /NG/g, length: 2, severity: LiabilitySeverity.High, ruleCategory: 'deamidation', color: ANNOTATION_COLORS.liability.deamidation, enabled: true},\n {id: 'deamid-ns', name: 'Deamidation (NS)', pattern: /NS/g, length: 2, severity: LiabilitySeverity.Medium, ruleCategory: 'deamidation', color: ANNOTATION_COLORS.liability.deamidation, enabled: true},\n {id: 'deamid-na', name: 'Deamidation (NA)', pattern: /NA/g, length: 2, severity: LiabilitySeverity.Low, ruleCategory: 'deamidation', color: ANNOTATION_COLORS.liability.deamidation, enabled: true},\n {id: 'deamid-nd', name: 'Deamidation (ND)', pattern: /ND/g, length: 2, severity: LiabilitySeverity.Low, ruleCategory: 'deamidation', color: ANNOTATION_COLORS.liability.deamidation, enabled: true},\n {id: 'deamid-nt', name: 'Deamidation (NT)', pattern: /NT/g, length: 2, severity: LiabilitySeverity.Low, ruleCategory: 'deamidation', color: ANNOTATION_COLORS.liability.deamidation, enabled: true},\n {id: 'isom-dg', name: 'Isomerization (DG)', pattern: /DG/g, length: 2, severity: LiabilitySeverity.High, ruleCategory: 'isomerization', color: ANNOTATION_COLORS.liability.isomerization, enabled: true},\n {id: 'isom-ds', name: 'Isomerization (DS)', pattern: /DS/g, length: 2, severity: LiabilitySeverity.Medium, ruleCategory: 'isomerization', color: ANNOTATION_COLORS.liability.isomerization, enabled: true},\n {id: 'oxid-m', name: 'Oxidation (Met)', pattern: /M/g, length: 1, severity: LiabilitySeverity.Medium, ruleCategory: 'oxidation', color: ANNOTATION_COLORS.liability.oxidation, enabled: true},\n {id: 'oxid-w', name: 'Oxidation (Trp)', pattern: /W/g, length: 1, severity: LiabilitySeverity.Low, ruleCategory: 'oxidation', color: ANNOTATION_COLORS.liability.oxidation, enabled: true},\n {id: 'glyco-nxst', name: 'N-glycosylation', pattern: /N[^P][ST]/g, length: 3, severity: LiabilitySeverity.High, ruleCategory: 'glycosylation', color: ANNOTATION_COLORS.liability.glycosylation, enabled: true},\n {id: 'free-cys', name: 'Free Cysteine', pattern: /C/g, length: 1, severity: LiabilitySeverity.Info, ruleCategory: 'freeCysteine', color: ANNOTATION_COLORS.liability.freeCysteine, enabled: false},\n];\n\n/** Extracts a canonical single-letter string from a sequence handler for a given row. */\nfunction getCanonicalString(sh: ISeqHandler, rowIdx: number): string {\n const splitted = sh.getSplitted(rowIdx);\n const chars: string[] = new Array(splitted.length);\n for (let i = 0; i < splitted.length; i++)\n chars[i] = splitted.getOriginal(i);\n return chars.join('');\n}\n\nexport interface ScanLiabilitiesResult {\n annotations: SeqAnnotation[];\n rowData: RowAnnotationData[];\n totalHits: number;\n}\n\n/** Scans all rows of a macromolecule column for liability motifs.\n * Returns column-level SeqAnnotation entries + per-row SeqAnnotationHit arrays. */\nexport function scanLiabilities(\n col: DG.Column<string>,\n sh: ISeqHandler,\n rules: LiabilityRule[],\n): ScanLiabilitiesResult {\n const enabledRules = rules.filter((r) => r.enabled);\n const posList = sh.posList;\n\n // Track which rules had hits\n const ruleHitCounts = new Map<string, number>();\n\n const rowData: RowAnnotationData[] = new Array(col.length);\n let totalHits = 0;\n\n for (let rowIdx = 0; rowIdx < col.length; rowIdx++) {\n const seq = getCanonicalString(sh, rowIdx);\n const hits: SeqAnnotationHit[] = [];\n\n for (const rule of enabledRules) {\n // Reset regex lastIndex for global patterns\n rule.pattern.lastIndex = 0;\n let match: RegExpExecArray | null;\n while ((match = rule.pattern.exec(seq)) !== null) {\n hits.push({\n annotationId: rule.id,\n positionIndex: match.index,\n positionName: match.index < posList.length ? posList[match.index] : undefined,\n matchedMonomers: match[0],\n });\n ruleHitCounts.set(rule.id, (ruleHitCounts.get(rule.id) ?? 0) + 1);\n totalHits++;\n }\n }\n rowData[rowIdx] = hits;\n }\n\n // Build column-level annotations only for rules that had hits\n const annotations: SeqAnnotation[] = enabledRules\n .filter((r) => ruleHitCounts.has(r.id))\n .map((r) => ({\n id: r.id,\n name: r.name,\n description: `${r.ruleCategory} liability pattern (${ruleHitCounts.get(r.id)} hits)`,\n start: null,\n end: null,\n visualType: r.length === 1 ? AnnotationVisualType.Point : AnnotationVisualType.Motif,\n category: AnnotationCategory.Liability,\n color: r.color,\n severity: r.severity,\n motifPattern: r.pattern.source,\n autoGenerated: true,\n }));\n\n return {annotations, rowData, totalHits};\n}\n\n/** Applies liability scan results to the DataFrame (writes tags + companion column). */\nexport function applyLiabilityScanResults(\n df: DG.DataFrame,\n seqCol: DG.Column<string>,\n result: ScanLiabilitiesResult,\n): void {\n // Merge with existing annotations, removing old liability entries\n const existing = getColumnAnnotations(seqCol)\n .filter((a) => a.category !== AnnotationCategory.Liability);\n setColumnAnnotations(seqCol, [...existing, ...result.annotations]);\n\n // Write per-row data to hidden companion column, preserving region hits from numbering\n const annotCol = getOrCreateAnnotationColumn(df, seqCol);\n for (let i = 0; i < result.rowData.length; i++) {\n const existingHits = getRowAnnotations(annotCol, i) ?? [];\n setRowAnnotations(annotCol, i, mergeRowHits(existingHits, result.rowData[i], false, true));\n }\n}\n\n/** Creates a liability summary count column (total hits per row). */\nexport function createLiabilitySummaryColumn(\n df: DG.DataFrame,\n seqCol: DG.Column<string>,\n result: ScanLiabilitiesResult,\n): DG.Column<number> {\n const colName = `${seqCol.name}_liability_count`;\n const counts = result.rowData.map((hits) => hits.length);\n const col = df.columns.addNewInt(colName);\n for (let i = 0; i < counts.length; i++)\n col.set(i, counts[i]);\n return col;\n}\n","import * as grok from 'datagrok-api/grok';\nimport * as ui from 'datagrok-api/ui';\nimport * as DG from 'datagrok-api/dg';\n\nimport {_package} from '../../package';\nimport {\n BUILTIN_LIABILITY_RULES, LiabilityRule, scanLiabilities,\n applyLiabilityScanResults, createLiabilitySummaryColumn,\n} from './liability-scanner';\nimport {LiabilitySeverity} from '@datagrok-libraries/bio/src/utils/macromolecule/annotations';\n\nconst severityLabels: Record<string, string> = {\n [LiabilitySeverity.High]: 'High',\n [LiabilitySeverity.Medium]: 'Medium',\n [LiabilitySeverity.Low]: 'Low',\n [LiabilitySeverity.Info]: 'Info',\n};\n\nexport function showLiabilityScannerDialog(): void {\n const df = grok.shell.tv?.dataFrame;\n if (!df) {\n grok.shell.warning('No table open');\n return;\n }\n\n const seqCols = df.columns.bySemTypeAll(DG.SEMTYPE.MACROMOLECULE);\n if (seqCols.length === 0) {\n grok.shell.warning('No macromolecule columns found');\n return;\n }\n\n const rules = BUILTIN_LIABILITY_RULES.map((r) => ({...r, pattern: new RegExp(r.pattern.source, 'g')}));\n\n const tableInput = ui.input.table('Table', {value: df});\n const seqInput = ui.input.column('Sequence', {\n table: df, value: seqCols[0],\n filter: (col: DG.Column) => col.semType === DG.SEMTYPE.MACROMOLECULE,\n });\n\n // Rule checkboxes\n const ruleChecks: {rule: LiabilityRule; input: DG.InputBase<boolean>}[] = [];\n const rulesDiv = ui.divV([]);\n for (const rule of rules) {\n const check = ui.input.bool(rule.name, {\n value: rule.enabled,\n tooltipText: `Severity: ${severityLabels[rule.severity] ?? rule.severity}`,\n });\n ruleChecks.push({rule, input: check});\n rulesDiv.append(check.root);\n }\n\n const highlightInput = ui.input.bool('Highlight in cell renderer', {value: true});\n const annotColInput = ui.input.bool('Create annotation column', {value: true});\n const summaryInput = ui.input.bool('Create summary count column', {value: false});\n\n const dialog = ui.dialog({title: 'Scan Sequence Liabilities'})\n .add(ui.inputs([tableInput, seqInput]))\n .add(ui.h3('Rules'))\n .add(rulesDiv)\n .add(ui.h3('Output'))\n .add(ui.inputs([highlightInput, annotColInput, summaryInput]))\n .onOK(() => {\n try {\n const seqCol = seqInput.value!;\n const sh = _package.seqHelper.getSeqHandler(seqCol);\n\n // Apply checkbox state\n for (const {rule, input} of ruleChecks)\n rule.enabled = input.value;\n\n const result = scanLiabilities(seqCol, sh, rules);\n\n if (annotColInput.value || highlightInput.value)\n applyLiabilityScanResults(df, seqCol, result);\n\n if (summaryInput.value)\n createLiabilitySummaryColumn(df, seqCol, result);\n\n grok.shell.info(`Liability scan: ${result.totalHits} hits found across ${result.annotations.length} rules`);\n df.fireValuesChanged();\n } catch (err: any) {\n grok.shell.error(`Liability scan failed: ${err.message ?? err}`);\n console.error(err);\n }\n });\n\n dialog.show();\n}\n"],"names":["BUILTIN_LIABILITY_RULES","id","name","pattern","length","severity","High","ruleCategory","color","liability","deamidation","enabled","Medium","Low","isomerization","oxidation","glycosylation","Info","freeCysteine","getCanonicalString","sh","rowIdx","splitted","getSplitted","chars","Array","i","getOriginal","join","severityLabels","showLiabilityScannerDialog","df","tv","dataFrame","warning","seqCols","columns","bySemTypeAll","MACROMOLECULE","rules","map","r","RegExp","source","tableInput","table","value","seqInput","column","filter","col","semType","ruleChecks","rulesDiv","rule","check","bool","tooltipText","push","input","append","root","highlightInput","annotColInput","summaryInput","dialog","title","add","onOK","seqCol","_package","seqHelper","getSeqHandler","result","enabledRules","posList","ruleHitCounts","Map","rowData","totalHits","seq","hits","match","lastIndex","exec","annotationId","positionIndex","index","positionName","undefined","matchedMonomers","set","get","annotations","has","description","start","end","visualType","Point","Motif","category","Liability","motifPattern","autoGenerated","scanLiabilities","existing","a","annotCol","existingHits","applyLiabilityScanResults","colName","counts","addNewInt","createLiabilitySummaryColumn","info","fireValuesChanged","err","error","message","console","show"],"sourceRoot":""}
package/dist/455.js CHANGED
@@ -1,2 +1,2 @@
1
- var bio;(()=>{"use strict";const t={V2K_RGP_SHIFT:8,V2K_RGP_LINE:"M RGP",V2K_A_LINE:"A ",V3K_COUNTS_SHIFT:14,V3K_IDX_SHIFT:7,V3K_HEADER_FIRST_LINE:"\nDatagrok macromolecule handler\n\n",V3K_HEADER_SECOND_LINE:" 0 0 0 0 0 0 999 V3000\n",V3K_BEGIN_CTAB_BLOCK:"M V30 BEGIN CTAB\n",V3K_END_CTAB_BLOCK:"M V30 END CTAB\n",V3K_BEGIN_COUNTS_LINE:"M V30 COUNTS ",V3K_COUNTS_LINE_ENDING:" 0 0 0\n",V3K_BEGIN_ATOM_BLOCK:"M V30 BEGIN ATOM\n",V3K_END_ATOM_BLOCK:"M V30 END ATOM\n",V3K_BEGIN_BOND_BLOCK:"M V30 BEGIN BOND\n",V3K_END_BOND_BLOCK:"M V30 END BOND\n",V3K_BOND_CONFIG:" CFG=",V3K_BEGIN_DATA_LINE:"M V30 ",V3K_END:"M END",PRECISION_FACTOR:1e4,DEOXYRIBOSE:{polymerType:"RNA",symbol:"d"},RIBOSE:{polymerType:"RNA",symbol:"r"},PHOSPHATE:{polymerType:"RNA",symbol:"p"},OXYGEN:"O",HYDROGEN:"H"};function o(t,o){return t[o.polymerType]?.[o.symbol]}class n extends Map{constructor(t){super(t)}}class e{constructor(t,o){this.molfile=t,this.monomers=o}static createEmpty(){return new e("",new n(null))}}var a,s,i;new RegExp("[rd]\\((\\w)\\)p?","g");class r{constructor(t,o,n){this.name=t,this.alphabet=o,this.cutoff=n}}!function(t){t.FASTA="fasta",t.SEPARATOR="separator",t.HELM="helm",t.CUSTOM="custom",t.BILN="biln"}(a||(a={})),function(t){t.DNA="DNA",t.RNA="RNA",t.PT="PT",t.UN="UN"}(s||(s={})),function(t){t.aligned="aligned",t.alphabet="alphabet",t.alphabetSize=".alphabetSize",t.alphabetIsMultichar=".alphabetIsMultichar",t.separator="separator",t.isHelmCompatible=".isHelmCompatible",t.positionNames=".positionNames",t.positionLabels=".positionLabels",t.regions=".regions",t.positionShift=".positionShift",t.selectedPosition=".selectedPosition",t.polymerTypeColumnName=".polymerTypeColumnName"}(i||(i={}));const l=new class{constructor(){this.fasta={peptide:new Set(["G","L","Y","S","E","Q","D","N","F","A","K","R","H","C","V","P","W","I","M","T"]),dna:new Set(["A","C","G","T"]),rna:new Set(["A","C","G","U"])}}};var c;function E(t){let o;switch(t){case c.BASE:case c.SUGAR:case c.LINKER:case c.NUCLEOTIDE:case"nucleotide":o="RNA";break;case c.AA:o="PEPTIDE";break;case c.CHEM:o="CHEM";break;case c.BLOB:o="BLOB";break;default:o="PEPTIDE",console.warn(`Unexpected HelmType '${t}'`)}return o}function h(a,i,r,l){if(0===a.length)return e.createEmpty();const c=m,{atomCount:h,bondCount:f,needsCapping:_}=c(a,i,r,l),S=new Array(h),d=new Array(f);let p,T=null,I=null;"PEPTIDE"===l?p=N:(p=b,T=r===s.DNA?o(i,t.DEOXYRIBOSE):o(i,t.RIBOSE),I=o(i,t.PHOSPHATE));const O={i:0,nodeShift:0,bondShift:0,backbonePositionShift:new Array(2).fill(0),branchPositionShift:new Array(2).fill(0),backboneAttachNode:0,branchAttachNode:0,flipFactor:1},u={sugar:T,phosphate:I,seqLength:a.length,atomCount:h,bondCount:f},C=new n,B=[];let L,y=0;for(O.i=0;O.i<u.seqLength;++O.i){const t=a[O.i];if(""===t.symbol)continue;const n=o(i,{symbol:t.symbol,polymerType:E(t.biotype)});L=n.terminalR2Atom;const e=O.nodeShift,s=O.bondShift;p(n,S,d,O,u),n.stereoAtoms?.forEach(t=>B.push(t+y)),y+=n.atoms.x.length;const r=O.nodeShift-e,l=new Array(r);for(let t=0;t<r;++t)l[t]=e+t;const c=O.bondShift-s,h=new Array(c);for(let t=0;t<c;++t)h[t]=s+t;C.set(O.i,{biotype:t.biotype,symbol:t.symbol,atoms:l,bonds:h})}_&&function(o,n,e,a,s=t.OXYGEN){const i=e.nodeShift+1;o[a.atomCount]=t.V3K_BEGIN_DATA_LINE+i+" "+(s??t.OXYGEN)+" "+A(e.backbonePositionShift[0])+" "+e.flipFactor*A(e.backbonePositionShift[1])+" 0.000000 0\n";const r=e.backboneAttachNode,l=i;n[a.bondCount]=t.V3K_BEGIN_DATA_LINE+e.bondShift+" 1 "+r+" "+l+"\n"}(S,d,O,u,L??t.OXYGEN);const D=t.V3K_BEGIN_COUNTS_LINE+h+" "+f+t.V3K_COUNTS_LINE_ENDING;let g="";return g+=t.V3K_HEADER_FIRST_LINE,g+=t.V3K_HEADER_SECOND_LINE,g+=t.V3K_BEGIN_CTAB_BLOCK,g+=D,g+=t.V3K_BEGIN_ATOM_BLOCK,g+=S.join(""),g+=t.V3K_END_ATOM_BLOCK,g+=t.V3K_BEGIN_BOND_BLOCK,g+=d.join(""),g+=t.V3K_END_BOND_BLOCK,B.length>0&&(g+=function(t){const o=[];let n=`M V30 MDLV30/STEABS ATOMS=(${t.length}`;for(let e=0;e<t.length;e++){const a=`${n} ${t[e]}`;a.length>76?(o.push(`${n} -\n`),n=`M V30 ${t[e]}`):n=a,e===t.length-1&&o.push(`${n})\n`)}return`M V30 BEGIN COLLECTION\n${o.join("")}M V30 END COLLECTION\n`}(B)),g+=t.V3K_END_CTAB_BLOCK,g+=t.V3K_END,{molfile:g,monomers:C}}function N(t,o,n,e){e.flipFactor=(-1)**(e.i%2),f(t,o,n,e)}function f(o,n,e,a){!function(o,n,e){for(let a=0;a<o.atoms.atomTypes.length;++a){const s=e.nodeShift+a+1;n[e.nodeShift+a]=t.V3K_BEGIN_DATA_LINE+s+" "+o.atoms.atomTypes[a]+" "+A(e.backbonePositionShift[0]+o.atoms.x[a])+" "+A(e.backbonePositionShift[1]+e.flipFactor*o.atoms.y[a])+" "+o.atoms.kwargs[a]}}(o,n,a),_(o,e,a),function(o,n,e){if(0!==e.backboneAttachNode){const a=e.bondShift,s=e.backboneAttachNode,i=o.meta.terminalNodes[0]+e.nodeShift;n[e.bondShift-1]=t.V3K_BEGIN_DATA_LINE+a+" 1 "+s+" "+i+"\n"}}(o,e,a),null!==o.meta.branchShift&&o.meta.terminalNodes.length>2&&function(t,o){o.branchAttachNode=o.nodeShift+t.meta.terminalNodes[2];for(let n=0;n<2;++n)o.branchPositionShift[n]=o.backbonePositionShift[n]+t.meta.branchShift[n]}(o,a),function(t,o){o.backboneAttachNode=o.nodeShift+t.meta.terminalNodes[1],o.bondShift+=t.bonds.atomPairs.length+1,o.nodeShift+=t.atoms.atomTypes.length,o.backbonePositionShift[0]+=t.meta.backboneShift?.[0]??0,o.backbonePositionShift[1]+=o.flipFactor*(t.meta.backboneShift?.[1]??0)}(o,a)}function b(o,n,e,a,s){if(0===a.i)f(s.sugar,n,e,a);else for(const t of[s.phosphate,s.sugar])f(t,n,e,a);!function(o,n,e,a){(function(o,n,e){for(let a=0;a<o.atoms.atomTypes.length;++a){const s=e.nodeShift+a+1;n[e.nodeShift+a]=t.V3K_BEGIN_DATA_LINE+s+" "+o.atoms.atomTypes[a]+" "+A(e.branchPositionShift[0]+o.atoms.x[a])+" "+A(e.branchPositionShift[1]+e.flipFactor*o.atoms.y[a])+" "+o.atoms.kwargs[a]}})(o,n,a),_(o,e,a),function(o,n,e){const a=e.bondShift,s=e.branchAttachNode,i=o.meta.terminalNodes[0]+e.nodeShift;n[a-1]=t.V3K_BEGIN_DATA_LINE+a+" 1 "+s+" "+i+"\n"}(o,e,a);const s=a.bondShift,i=a.branchAttachNode,r=o.meta.terminalNodes[0]+a.nodeShift;e[s-1]=t.V3K_BEGIN_DATA_LINE+s+" 1 "+i+" "+r+"\n",a.bondShift+=o.bonds.atomPairs.length+1,a.nodeShift+=o.atoms.atomTypes.length}(o,n,e,a)}function _(o,n,e){for(let a=0;a<o.bonds.atomPairs.length;++a){const s=e.bondShift+a+1,i=o.bonds.atomPairs[a][0]+e.nodeShift,r=o.bonds.atomPairs[a][1]+e.nodeShift;let l="";if(o.bonds.bondConfiguration.has(a)){let t=o.bonds.bondConfiguration.get(a);e.flipFactor<0&&(t=1===t?3:1),l=" CFG="+t}const c=o.bonds.kwargs.has(a)?" "+o.bonds.kwargs.get(a):"";n[e.bondShift+a]=t.V3K_BEGIN_DATA_LINE+s+" "+o.bonds.bondTypes[a]+" "+i+" "+r+l+c+"\n"}}function m(n,e,a,i){let r=0,l=0,c=0,h=!0,N=null;for(const t of n)if(""!==t.symbol){if("*"==t.symbol)throw new Error("Gap canonical symbol is '', not '*");N=o(e,{symbol:t.symbol,polymerType:E(t.biotype)}),r+=N.atoms.x.length,l+=N.bonds.bondTypes.length,c++}if("PEPTIDE"===i)r+=1,l+=c,c>0&&((N?.meta?.rNodes?.length??0)<2||N?.terminalR2Atom?.toLowerCase()===t.HYDROGEN.toLowerCase())&&(h=!1,r-=1,l-=1);else{const i=a===s.DNA?o(e,t.DEOXYRIBOSE):o(e,t.RIBOSE),c=o(e,t.PHOSPHATE);r+=(n.length-1)*c.atoms.x.length,r+=n.length*i.atoms.x.length,r+=1,l+=(n.length-1)*c.bonds.bondTypes.length,l+=n.length*i.bonds.bondTypes.length,l-=1,l+=3*n.length}return{atomCount:r,bondCount:l,needsCapping:h}}function A(o){return Math.round(t.PRECISION_FACTOR*o)/t.PRECISION_FACTOR}new r(s.PT,l.fasta.peptide,.5),new r(s.DNA,l.fasta.dna,.55),new r(s.RNA,l.fasta.rna,.55),a.FASTA,a.SEPARATOR,a.HELM,a.BILN,function(t){t.BASE="HELM_BASE",t.SUGAR="HELM_SUGAR",t.LINKER="HELM_LINKER",t.AA="HELM_AA",t.CHEM="HELM_CHEM",t.BLOB="HELM_BLOB",t.NUCLEOTIDE="HELM_NUCLETIDE"}(c||(c={})),onmessage=t=>{const{seqList:o,monomersDict:n,alphabet:e,polymerType:a,start:s,end:i}=t.data,r=new Array(i-s),l=new Array(0);for(let t=s;t<i;++t)try{const i=o[t];r[t-s]=h(i,n,e,a)}catch(o){const n=`Cannot get molfile of row #${t}: ${o instanceof Error?o.message:o.toString()}.`;l.push(n)}postMessage({molfiles:r,warnings:l})},bio={}})();
1
+ var bio;(()=>{"use strict";const t={V2K_RGP_SHIFT:8,V2K_RGP_LINE:"M RGP",V2K_A_LINE:"A ",V3K_COUNTS_SHIFT:14,V3K_IDX_SHIFT:7,V3K_HEADER_FIRST_LINE:"\nDatagrok macromolecule handler\n\n",V3K_HEADER_SECOND_LINE:" 0 0 0 0 0 0 999 V3000\n",V3K_BEGIN_CTAB_BLOCK:"M V30 BEGIN CTAB\n",V3K_END_CTAB_BLOCK:"M V30 END CTAB\n",V3K_BEGIN_COUNTS_LINE:"M V30 COUNTS ",V3K_COUNTS_LINE_ENDING:" 0 0 0\n",V3K_BEGIN_ATOM_BLOCK:"M V30 BEGIN ATOM\n",V3K_END_ATOM_BLOCK:"M V30 END ATOM\n",V3K_BEGIN_BOND_BLOCK:"M V30 BEGIN BOND\n",V3K_END_BOND_BLOCK:"M V30 END BOND\n",V3K_BOND_CONFIG:" CFG=",V3K_BEGIN_DATA_LINE:"M V30 ",V3K_END:"M END",PRECISION_FACTOR:1e4,DEOXYRIBOSE:{polymerType:"RNA",symbol:"d"},RIBOSE:{polymerType:"RNA",symbol:"r"},PHOSPHATE:{polymerType:"RNA",symbol:"p"},OXYGEN:"O",HYDROGEN:"H"};function o(t,o){return t[o.polymerType]?.[o.symbol]}class n extends Map{constructor(t){super(t)}}class e{constructor(t,o){this.molfile=t,this.monomers=o}static createEmpty(){return new e("",new n(null))}}var a,s,i;new RegExp("[rd]\\((\\w)\\)p?","g");class r{constructor(t,o,n){this.name=t,this.alphabet=o,this.cutoff=n}}!function(t){t.FASTA="fasta",t.SEPARATOR="separator",t.HELM="helm",t.CUSTOM="custom",t.BILN="biln"}(a||(a={})),function(t){t.DNA="DNA",t.RNA="RNA",t.PT="PT",t.UN="UN"}(s||(s={})),function(t){t.aligned="aligned",t.alphabet="alphabet",t.alphabetSize=".alphabetSize",t.alphabetIsMultichar=".alphabetIsMultichar",t.separator="separator",t.isHelmCompatible=".isHelmCompatible",t.positionNames=".positionNames",t.positionLabels=".positionLabels",t.regions=".regions",t.positionShift=".positionShift",t.selectedPosition=".selectedPosition",t.polymerTypeColumnName=".polymerTypeColumnName",t.annotations=".annotations",t.numberingScheme=".numberingScheme",t.annotationColumnName=".annotationColumnName"}(i||(i={}));const l=new class{constructor(){this.fasta={peptide:new Set(["G","L","Y","S","E","Q","D","N","F","A","K","R","H","C","V","P","W","I","M","T"]),dna:new Set(["A","C","G","T"]),rna:new Set(["A","C","G","U"])}}};var c;function h(t){let o;switch(t){case c.BASE:case c.SUGAR:case c.LINKER:case c.NUCLEOTIDE:case"nucleotide":o="RNA";break;case c.AA:o="PEPTIDE";break;case c.CHEM:o="CHEM";break;case c.BLOB:o="BLOB";break;default:o="PEPTIDE",console.warn(`Unexpected HelmType '${t}'`)}return o}function E(a,i,r,l){if(0===a.length)return e.createEmpty();const c=_,{atomCount:E,bondCount:f,needsCapping:m}=c(a,i,r,l),S=new Array(E),d=new Array(f);let p,T=null,I=null;"PEPTIDE"===l?p=N:(p=b,T=r===s.DNA?o(i,t.DEOXYRIBOSE):o(i,t.RIBOSE),I=o(i,t.PHOSPHATE));const u={i:0,nodeShift:0,bondShift:0,backbonePositionShift:new Array(2).fill(0),branchPositionShift:new Array(2).fill(0),backboneAttachNode:0,branchAttachNode:0,flipFactor:1},O={sugar:T,phosphate:I,seqLength:a.length,atomCount:E,bondCount:f},C=new n,B=[];let L,y=0;for(u.i=0;u.i<O.seqLength;++u.i){const t=a[u.i];if(""===t.symbol)continue;const n=o(i,{symbol:t.symbol,polymerType:h(t.biotype)});L=n.terminalR2Atom;const e=u.nodeShift,s=u.bondShift;p(n,S,d,u,O),n.stereoAtoms?.forEach(t=>B.push(t+y)),y+=n.atoms.x.length;const r=u.nodeShift-e,l=new Array(r);for(let t=0;t<r;++t)l[t]=e+t;const c=u.bondShift-s,E=new Array(c);for(let t=0;t<c;++t)E[t]=s+t;C.set(u.i,{biotype:t.biotype,symbol:t.symbol,atoms:l,bonds:E})}m&&function(o,n,e,a,s=t.OXYGEN){const i=e.nodeShift+1;o[a.atomCount]=t.V3K_BEGIN_DATA_LINE+i+" "+(s??t.OXYGEN)+" "+A(e.backbonePositionShift[0])+" "+e.flipFactor*A(e.backbonePositionShift[1])+" 0.000000 0\n";const r=e.backboneAttachNode,l=i;n[a.bondCount]=t.V3K_BEGIN_DATA_LINE+e.bondShift+" 1 "+r+" "+l+"\n"}(S,d,u,O,L??t.OXYGEN);const D=t.V3K_BEGIN_COUNTS_LINE+E+" "+f+t.V3K_COUNTS_LINE_ENDING;let g="";return g+=t.V3K_HEADER_FIRST_LINE,g+=t.V3K_HEADER_SECOND_LINE,g+=t.V3K_BEGIN_CTAB_BLOCK,g+=D,g+=t.V3K_BEGIN_ATOM_BLOCK,g+=S.join(""),g+=t.V3K_END_ATOM_BLOCK,g+=t.V3K_BEGIN_BOND_BLOCK,g+=d.join(""),g+=t.V3K_END_BOND_BLOCK,B.length>0&&(g+=function(t){const o=[];let n=`M V30 MDLV30/STEABS ATOMS=(${t.length}`;for(let e=0;e<t.length;e++){const a=`${n} ${t[e]}`;a.length>76?(o.push(`${n} -\n`),n=`M V30 ${t[e]}`):n=a,e===t.length-1&&o.push(`${n})\n`)}return`M V30 BEGIN COLLECTION\n${o.join("")}M V30 END COLLECTION\n`}(B)),g+=t.V3K_END_CTAB_BLOCK,g+=t.V3K_END,{molfile:g,monomers:C}}function N(t,o,n,e){e.flipFactor=(-1)**(e.i%2),f(t,o,n,e)}function f(o,n,e,a){!function(o,n,e){for(let a=0;a<o.atoms.atomTypes.length;++a){const s=e.nodeShift+a+1;n[e.nodeShift+a]=t.V3K_BEGIN_DATA_LINE+s+" "+o.atoms.atomTypes[a]+" "+A(e.backbonePositionShift[0]+o.atoms.x[a])+" "+A(e.backbonePositionShift[1]+e.flipFactor*o.atoms.y[a])+" "+o.atoms.kwargs[a]}}(o,n,a),m(o,e,a),function(o,n,e){if(0!==e.backboneAttachNode){const a=e.bondShift,s=e.backboneAttachNode,i=o.meta.terminalNodes[0]+e.nodeShift;n[e.bondShift-1]=t.V3K_BEGIN_DATA_LINE+a+" 1 "+s+" "+i+"\n"}}(o,e,a),null!==o.meta.branchShift&&o.meta.terminalNodes.length>2&&function(t,o){o.branchAttachNode=o.nodeShift+t.meta.terminalNodes[2];for(let n=0;n<2;++n)o.branchPositionShift[n]=o.backbonePositionShift[n]+t.meta.branchShift[n]}(o,a),function(t,o){o.backboneAttachNode=o.nodeShift+t.meta.terminalNodes[1],o.bondShift+=t.bonds.atomPairs.length+1,o.nodeShift+=t.atoms.atomTypes.length,o.backbonePositionShift[0]+=t.meta.backboneShift?.[0]??0,o.backbonePositionShift[1]+=o.flipFactor*(t.meta.backboneShift?.[1]??0)}(o,a)}function b(o,n,e,a,s){if(0===a.i)f(s.sugar,n,e,a);else for(const t of[s.phosphate,s.sugar])f(t,n,e,a);!function(o,n,e,a){(function(o,n,e){for(let a=0;a<o.atoms.atomTypes.length;++a){const s=e.nodeShift+a+1;n[e.nodeShift+a]=t.V3K_BEGIN_DATA_LINE+s+" "+o.atoms.atomTypes[a]+" "+A(e.branchPositionShift[0]+o.atoms.x[a])+" "+A(e.branchPositionShift[1]+e.flipFactor*o.atoms.y[a])+" "+o.atoms.kwargs[a]}})(o,n,a),m(o,e,a),function(o,n,e){const a=e.bondShift,s=e.branchAttachNode,i=o.meta.terminalNodes[0]+e.nodeShift;n[a-1]=t.V3K_BEGIN_DATA_LINE+a+" 1 "+s+" "+i+"\n"}(o,e,a);const s=a.bondShift,i=a.branchAttachNode,r=o.meta.terminalNodes[0]+a.nodeShift;e[s-1]=t.V3K_BEGIN_DATA_LINE+s+" 1 "+i+" "+r+"\n",a.bondShift+=o.bonds.atomPairs.length+1,a.nodeShift+=o.atoms.atomTypes.length}(o,n,e,a)}function m(o,n,e){for(let a=0;a<o.bonds.atomPairs.length;++a){const s=e.bondShift+a+1,i=o.bonds.atomPairs[a][0]+e.nodeShift,r=o.bonds.atomPairs[a][1]+e.nodeShift;let l="";if(o.bonds.bondConfiguration.has(a)){let t=o.bonds.bondConfiguration.get(a);e.flipFactor<0&&(t=1===t?3:1),l=" CFG="+t}const c=o.bonds.kwargs.has(a)?" "+o.bonds.kwargs.get(a):"";n[e.bondShift+a]=t.V3K_BEGIN_DATA_LINE+s+" "+o.bonds.bondTypes[a]+" "+i+" "+r+l+c+"\n"}}function _(n,e,a,i){let r=0,l=0,c=0,E=!0,N=null;for(const t of n)if(""!==t.symbol){if("*"==t.symbol)throw new Error("Gap canonical symbol is '', not '*");N=o(e,{symbol:t.symbol,polymerType:h(t.biotype)}),r+=N.atoms.x.length,l+=N.bonds.bondTypes.length,c++}if("PEPTIDE"===i)r+=1,l+=c,c>0&&((N?.meta?.rNodes?.length??0)<2||N?.terminalR2Atom?.toLowerCase()===t.HYDROGEN.toLowerCase())&&(E=!1,r-=1,l-=1);else{const i=a===s.DNA?o(e,t.DEOXYRIBOSE):o(e,t.RIBOSE),c=o(e,t.PHOSPHATE);r+=(n.length-1)*c.atoms.x.length,r+=n.length*i.atoms.x.length,r+=1,l+=(n.length-1)*c.bonds.bondTypes.length,l+=n.length*i.bonds.bondTypes.length,l-=1,l+=3*n.length}return{atomCount:r,bondCount:l,needsCapping:E}}function A(o){return Math.round(t.PRECISION_FACTOR*o)/t.PRECISION_FACTOR}new r(s.PT,l.fasta.peptide,.5),new r(s.DNA,l.fasta.dna,.55),new r(s.RNA,l.fasta.rna,.55),a.FASTA,a.SEPARATOR,a.HELM,a.BILN,function(t){t.BASE="HELM_BASE",t.SUGAR="HELM_SUGAR",t.LINKER="HELM_LINKER",t.AA="HELM_AA",t.CHEM="HELM_CHEM",t.BLOB="HELM_BLOB",t.NUCLEOTIDE="HELM_NUCLETIDE"}(c||(c={})),onmessage=t=>{const{seqList:o,monomersDict:n,alphabet:e,polymerType:a,start:s,end:i}=t.data,r=new Array(i-s),l=new Array(0);for(let t=s;t<i;++t)try{const i=o[t];r[t-s]=E(i,n,e,a)}catch(o){const n=`Cannot get molfile of row #${t}: ${o instanceof Error?o.message:o.toString()}.`;l.push(n)}postMessage({molfiles:r,warnings:l})},bio={}})();
2
2
  //# sourceMappingURL=455.js.map
package/dist/455.js.map CHANGED
@@ -1 +1 @@
1
- {"version":3,"file":"455.js","mappings":"2BAAO,MAAMA,EAAqB,CAE9BC,cAAe,EACfC,aAAc,SACdC,WAAY,MAEZC,iBAAkB,GAClBC,cAAe,EACfC,sBAAuB,uCACvBC,uBAAwB,4CACxBC,qBAAsB,sBACtBC,mBAAoB,oBACpBC,sBAAuB,iBACvBC,uBAAwB,WACxBC,qBAAsB,sBACtBC,mBAAoB,oBACpBC,qBAAsB,sBACtBC,mBAAoB,oBACpBC,gBAAiB,QACjBC,oBAAqB,UACrBC,QAAS,SACTC,iBAAkB,IAElBC,YAAa,CAAEC,YAAa,MAA8BC,OAAQ,KAClEC,OAAQ,CAAEF,YAAa,MAA8BC,OAAQ,KAC7DE,UAAW,CAAEH,YAAa,MAA8BC,OAAQ,KAChEG,OAAQ,IACRC,SAAU,KC3BP,SAASC,EAAYC,EAAMC,GAC9B,OAAOD,EAAKC,EAAOR,eAAeQ,EAAOP,OAC7C,CAUO,MAAMQ,UAAmBC,IAC5B,WAAAC,CAAYC,GACRC,MAAMD,EACV,EAGG,MAAME,EACT,WAAAH,CAAYI,EAASC,GACjBC,KAAKF,QAAUA,EACfE,KAAKD,SAAWA,CACpB,CACA,kBAAOE,GAAgB,OAAO,IAAIJ,EAAe,GAAI,IAAIL,EAAW,MAAQ,E,ICrBrEU,EAQAC,EAQAC,ECyByB,IAAIC,OAAO,oBAA0E,KCdlH,MAAMC,EACT,WAAAZ,CAAYa,EAAMC,EAAUC,GACxBT,KAAKO,KAAOA,EACZP,KAAKQ,SAAWA,EAChBR,KAAKS,OAASA,CAClB,GF/BJ,SAAWP,GACPA,EAAgB,MAAI,QACpBA,EAAoB,UAAI,YACxBA,EAAe,KAAI,OACaA,EAAiB,OAAI,SACrBA,EAAe,KAAI,MACtD,CAND,CAMGA,IAAaA,EAAW,CAAC,IAE5B,SAAWC,GACPA,EAAc,IAAI,MAClBA,EAAc,IAAI,MAClBA,EAAa,GAAI,KAEjBA,EAAa,GAAI,IACpB,CAND,CAMGA,IAAaA,EAAW,CAAC,IAE5B,SAAWC,GACPA,EAAc,QAAI,UAClBA,EAAe,SAAI,WACnBA,EAAmB,aAAI,gBACvBA,EAA0B,oBAAI,uBAC9BA,EAAgB,UAAI,YACpBA,EAAuB,iBAAI,oBAC3BA,EAAoB,cAAI,iBACxBA,EAAqB,eAAI,kBACzBA,EAAc,QAAI,WAClBA,EAAoB,cAAI,iBACxBA,EAAuB,iBAAI,oBAC3BA,EAA4B,sBAAI,wBACnC,CAbD,CAaGA,IAASA,EAAO,CAAC,IAEb,MAIMM,EAAY,IAAI,MACzB,WAAAhB,GACIM,KAAKW,MAAQ,CACTC,QAAS,IAAIC,IAAI,CACb,IAAK,IAAK,IAAK,IAAK,IAAK,IAAK,IAAK,IAAK,IAAK,IAC7C,IAAK,IAAK,IAAK,IAAK,IAAK,IAAK,IAAK,IAAK,IAAK,MAEjDC,IAAK,IAAID,IAAI,CAAC,IAAK,IAAK,IAAK,MAC7BE,IAAK,IAAIF,IAAI,CAAC,IAAK,IAAK,IAAK,MAErC,GGhDG,IAAIG,ECgBJ,SAASC,EAAsBC,GAClC,IAAInC,EACJ,OAAQmC,GACJ,KAAKF,EAAUG,KACf,KAAKH,EAAUI,MACf,KAAKJ,EAAUK,OACf,KAAKL,EAAUM,WAEf,IAAK,aACDvC,EAAc,MACd,MACJ,KAAKiC,EAAUO,GACXxC,EAAc,UACd,MACJ,KAAKiC,EAAUQ,KACXzC,EAAc,OACd,MACJ,KAAKiC,EAAUS,KACX1C,EAAc,OACd,MACJ,QACIA,EAAc,UACd2C,QAAQC,KAAK,wBAAwBT,MAE7C,OAAOnC,CACX,CCHO,SAAS6C,EAAoBC,EAAYC,EAActB,EAAUzB,GACpE,GAA0B,IAAtB8C,EAAWE,OAEX,OAAOlC,EAAeI,cAG1B,MAAM+B,EAAuBC,GACvB,UAAEC,EAAS,UAAEC,EAAS,aAAEC,GAAiBJ,EAAqBH,EAAYC,EAActB,EAAUzB,GAElGsD,EAAmB,IAAIC,MAAMJ,GAC7BK,EAAmB,IAAID,MAAMH,GACnC,IAAIK,EACAC,EAAQ,KACRC,EAAY,KACI,YAAhB3D,EACAyD,EAAuBG,GAEvBH,EAAuBI,EACvBH,EAASjC,IAAaL,EAAS0C,IAAOxD,EAAYyC,EAAc,EAAEhD,aAAeO,EAAYyC,EAAc,EAAE7C,QAC7GyD,EAAYrD,EAAYyC,EAAc,EAAE5C,YAE5C,MAAM4D,EAAI,CACNC,EAAG,EACHC,UAAW,EACXC,UAAW,EACXC,sBAAuB,IAAIZ,MAAM,GAAGa,KAAK,GACzCC,oBAAqB,IAAId,MAAM,GAAGa,KAAK,GACvCE,mBAAoB,EACpBC,iBAAkB,EAClBC,WAAY,GAEVC,EAAK,CACPf,MAAOA,EACPC,UAAWA,EACXe,UAAW5B,EAAWE,OACtBG,UAAWA,EACXC,UAAWA,GAETpC,EAAW,IAAIP,EACfkE,EAAmB,GACzB,IACIC,EADAC,EAAS,EAEb,IAAKd,EAAEC,EAAI,EAAGD,EAAEC,EAAIS,EAAGC,YAAaX,EAAEC,EAAG,CACrC,MAAMc,EAAahC,EAAWiB,EAAEC,GAChC,GL1BkB,KK0Bdc,EAAW7E,OACX,SACJ,MAAM8E,EAAUzE,EAAYyC,EAAc,CAAE9C,OAAQ6E,EAAW7E,OAAQD,YAAakC,EAAsB4C,EAAWE,WACrHJ,EAAyBG,EAAQE,eACjC,MAAMC,EAAanB,EAAEE,UACfkB,EAAapB,EAAEG,UACrBT,EAAqBsB,EAASzB,EAAkBE,EAAkBO,EAAGU,GAErEM,EAAQK,aAAaC,QAASrB,GAAMW,EAAiBW,KAAKtB,EAAIa,IAC9DA,GAAUE,EAAQQ,MAAMC,EAAExC,OAC1B,MAAMyC,EAAa1B,EAAEE,UAAYiB,EAC3BQ,EAAY,IAAInC,MAAMkC,GAC5B,IAAK,IAAIE,EAAM,EAAGA,EAAMF,IAAcE,EAClCD,EAAUC,GAAOT,EAAaS,EAClC,MAAMC,EAAa7B,EAAEG,UAAYiB,EAC3BU,EAAY,IAAItC,MAAMqC,GAC5B,IAAK,IAAIE,EAAM,EAAGA,EAAMF,IAAcE,EAClCD,EAAUC,GAAOX,EAAaW,EAClC9E,EAAS+E,IAAIhC,EAAEC,EAAG,CACdgB,QAASF,EAAWE,QACpB/E,OAAQ6E,EAAW7E,OACnBsF,MAAOG,EAAWM,MAAOH,GAEjC,CAEIxC,GAiDR,SAA8BC,EAAkBE,EAAkBO,EAAGU,EAAIwB,EAAkB,EAAE7F,QAEzF,MAAM8F,EAAUnC,EAAEE,UAAY,EAC9BX,EAAiBmB,EAAGtB,WAAa,EAAEvD,oBAAsBsG,EAAU,KAC9DD,GAAmB,EAAE7F,QAAU,IAAM+F,EAAcpC,EAAEI,sBAAsB,IAAM,IAClFJ,EAAES,WAAa2B,EAAcpC,EAAEI,sBAAsB,IAFxB,gBAIjC,MAAMiC,EAAYrC,EAAEO,mBACd+B,EAAaH,EACnB1C,EAAiBiB,EAAGrB,WAAa,EAAExD,oBAAsBmE,EAAEG,UAA1B,MACnBkC,EAAY,IAAMC,EAAa,IACjD,CA3DQC,CAAqBhD,EAAkBE,EAAkBO,EAAGU,EAAIG,GAA0B,EAAExE,QAChG,MAAMmG,EAAoB,EAAElH,sBAAwB8D,EAAY,IAAMC,EAAY,EAAE9D,uBAIpF,IAAIkH,EAAS,GAgBb,OAfAA,GAAU,EAAEvH,sBACZuH,GAAU,EAAEtH,uBACZsH,GAAU,EAAErH,qBACZqH,GAAUD,EACVC,GAAU,EAAEjH,qBACZiH,GAAUlD,EAAiBmD,KAAK,IAChCD,GAAU,EAAEhH,mBACZgH,GAAU,EAAE/G,qBACZ+G,GAAUhD,EAAiBiD,KAAK,IAChCD,GAAU,EAAE9G,mBACRiF,EAAiB3B,OAAS,IAC1BwD,GAMR,SAA4BE,GAGxB,MACMC,EAAY,GAClB,IAAIC,EAAmB,+BAA+BF,EAAW1D,SACjE,IAAK,IAAIgB,EAAI,EAAGA,EAAI0C,EAAW1D,OAAQgB,IAAK,CACxC,MAAM6C,EAAa,GAAGD,KAAoBF,EAAW1C,KACjD6C,EAAW7D,OALA,IAMX2D,EAAUrB,KAAK,GAAGsB,SAClBA,EAAmB,UAAUF,EAAW1C,MAGxC4C,EAAmBC,EAEnB7C,IAAM0C,EAAW1D,OAAS,GAC1B2D,EAAUrB,KAAK,GAAGsB,OAC1B,CACA,MAAO,4BAA4BD,EAAUF,KAAK,4BACtD,CAzBkBK,CAAmBnC,IACjC6B,GAAU,EAAEpH,mBACZoH,GAAU,EAAE3G,QAEL,CAAEkB,QAASyF,EAAQxF,SAAUA,EACxC,CAsCA,SAAS4C,EAAuBmB,EAASzB,EAAkBE,EAAkBO,GACzEA,EAAES,aAAe,KAAOT,EAAEC,EAAI,GAC9B+C,EAA6BhC,EAASzB,EAAkBE,EAAkBO,EAC9E,CACA,SAASgD,EAA6BhC,EAASzB,EAAkBE,EAAkBO,IAmDnF,SAAuBgB,EAASzB,EAAkBS,GAC9C,IAAK,IAAIiD,EAAI,EAAGA,EAAIjC,EAAQQ,MAAM0B,UAAUjE,SAAUgE,EAAG,CACrD,MAAMd,EAAUnC,EAAEE,UAAY+C,EAAI,EAClC1D,EAAiBS,EAAEE,UAAY+C,GAAK,EAAEpH,oBAAsBsG,EAAU,IAClEnB,EAAQQ,MAAM0B,UAAUD,GAAK,IAC7Bb,EAAcpC,EAAEI,sBAAsB,GAAKY,EAAQQ,MAAMC,EAAEwB,IAAM,IACjEb,EAAcpC,EAAEI,sBAAsB,GAAKJ,EAAES,WAAaO,EAAQQ,MAAM2B,EAAEF,IAC1E,IAAMjC,EAAQQ,MAAM4B,OAAOH,EACnC,CACJ,CAzDII,CAAcrC,EAASzB,EAAkBS,GAEzCsD,EAActC,EAASvB,EAAkBO,GAwF7C,SAAgCgB,EAASvB,EAAkBO,GACvD,GAA6B,IAAzBA,EAAEO,mBAA0B,CAC5B,MAAMgD,EAAUvD,EAAEG,UACZkC,EAAYrC,EAAEO,mBACd+B,EAAatB,EAAQwC,KAAKC,cAAc,GAAKzD,EAAEE,UACrDT,EAAiBO,EAAEG,UAAY,GAAK,EAAEtE,oBAAsB0H,EAAxB,MACtBlB,EAAY,IAAMC,EAAa,IACjD,CACJ,CA9FIoB,CAAuB1C,EAASvB,EAAkBO,GAEjB,OAA7BgB,EAAQwC,KAAKG,aAAwB3C,EAAQwC,KAAKC,cAAcxE,OAAS,GAqCjF,SAA+B+B,EAAShB,GACpCA,EAAEQ,iBAAmBR,EAAEE,UAAYc,EAAQwC,KAAKC,cAAc,GAC9D,IAAK,IAAIxD,EAAI,EAAGA,EAAI,IAAKA,EACrBD,EAAEM,oBAAoBL,GAAKD,EAAEI,sBAAsBH,GAAKe,EAAQwC,KAAKG,YAAY1D,EACzF,CAxCQ2D,CAAsB5C,EAAShB,GA6BvC,SAAuCgB,EAAShB,GAC5CA,EAAEO,mBAAqBP,EAAEE,UAAYc,EAAQwC,KAAKC,cAAc,GAChEzD,EAAEG,WAAaa,EAAQiB,MAAM4B,UAAU5E,OAAS,EAChDe,EAAEE,WAAac,EAAQQ,MAAM0B,UAAUjE,OACvCe,EAAEI,sBAAsB,IAAMY,EAAQwC,KAAKM,gBAAgB,IAAM,EACjE9D,EAAEI,sBAAsB,IAAMJ,EAAES,YAAcO,EAAQwC,KAAKM,gBAAgB,IAAM,EACrF,CAjCIC,CAA8B/C,EAAShB,EAC3C,CACA,SAASF,EAAwBkE,EAAYzE,EAAkBE,EAAkBO,EAAGU,GAGhF,GAAY,IAARV,EAAEC,EACF+C,EAA6BtC,EAAGf,MAAOJ,EAAkBE,EAAkBO,QAE3E,IAAK,MAAMgB,IAAW,CAACN,EAAGd,UAAWc,EAAGf,OACpCqD,EAA6BhC,EAASzB,EAAkBE,EAAkBO,IAItF,SAAoCgB,EAASzB,EAAkBE,EAAkBO,IAqCjF,SAA6BgB,EAASzB,EAAkBS,GACpD,IAAK,IAAIiD,EAAI,EAAGA,EAAIjC,EAAQQ,MAAM0B,UAAUjE,SAAUgE,EAAG,CACrD,MAAMd,EAAUnC,EAAEE,UAAY+C,EAAI,EAClC1D,EAAiBS,EAAEE,UAAY+C,GAAK,EAAEpH,oBAAsBsG,EAAU,IAClEnB,EAAQQ,MAAM0B,UAAUD,GAAK,IAC7Bb,EAAcpC,EAAEM,oBAAoB,GAAKU,EAAQQ,MAAMC,EAAEwB,IAAM,IAC/Db,EAAcpC,EAAEM,oBAAoB,GAAKN,EAAES,WAAaO,EAAQQ,MAAM2B,EAAEF,IACxE,IAAMjC,EAAQQ,MAAM4B,OAAOH,EACnC,CACJ,EA7CIgB,CAAoBjD,EAASzB,EAAkBS,GAC/CsD,EAActC,EAASvB,EAAkBO,GA4E7C,SAAkCkE,EAAezE,EAAkBO,GAC/D,MAAMuD,EAAUvD,EAAEG,UACZkC,EAAYrC,EAAEQ,iBACd8B,EAAa4B,EAAcV,KAAKC,cAAc,GAAKzD,EAAEE,UAC3DT,EAAiB8D,EAAU,GAAK,EAAE1H,oBAAsB0H,EAAxB,MAClBlB,EAAY,IAAMC,EAAa,IACjD,CAjFI6B,CAAyBnD,EAASvB,EAAkBO,GAEpD,MAAMuD,EAAUvD,EAAEG,UACZkC,EAAYrC,EAAEQ,iBACd8B,EAAatB,EAAQwC,KAAKC,cAAc,GAAKzD,EAAEE,UACrDT,EAAiB8D,EAAU,GAAK,EAAE1H,oBAAsB0H,EAAxB,MAClBlB,EAAY,IAAMC,EAAa,KAE7CtC,EAAEG,WAAaa,EAAQiB,MAAM4B,UAAU5E,OAAS,EAChDe,EAAEE,WAAac,EAAQQ,MAAM0B,UAAUjE,MAC3C,CAfImF,CAA2BJ,EAAYzE,EAAkBE,EAAkBO,EAC/E,CAgDA,SAASsD,EAActC,EAASvB,EAAkBO,GAE9C,IAAK,IAAIiD,EAAI,EAAGA,EAAIjC,EAAQiB,MAAM4B,UAAU5E,SAAUgE,EAAG,CACrD,MAAMM,EAAUvD,EAAEG,UAAY8C,EAAI,EAC5BZ,EAAYrB,EAAQiB,MAAM4B,UAAUZ,GAAG,GAAKjD,EAAEE,UAC9CoC,EAAatB,EAAQiB,MAAM4B,UAAUZ,GAAG,GAAKjD,EAAEE,UACrD,IAAImE,EAAU,GACd,GAAIrD,EAAQiB,MAAMqC,kBAAkBC,IAAItB,GAAI,CAExC,IAAIuB,EAAcxD,EAAQiB,MAAMqC,kBAAkBG,IAAIxB,GAClDjD,EAAES,WAAa,IACf+D,EAA+B,IAAhBA,EAAqB,EAAI,GAC5CH,EAAU,QAAUG,CACxB,CACA,MAAMpB,EAASpC,EAAQiB,MAAMmB,OAAOmB,IAAItB,GACpC,IAAMjC,EAAQiB,MAAMmB,OAAOqB,IAAIxB,GAAK,GACxCxD,EAAiBO,EAAEG,UAAY8C,GAAK,EAAEpH,oBAAsB0H,EAAU,IAClEvC,EAAQiB,MAAMyC,UAAUzB,GAAK,IAC7BZ,EAAY,IAAMC,EAAa+B,EAAUjB,EAAS,IAC1D,CACJ,CAyBA,SAASjE,EAA2BJ,EAAYC,EAActB,EAAUzB,GACpE,IAAImD,EAAY,EACZC,EAAY,EACZsF,EAAe,EACfrF,GAAe,EACfsF,EAAmB,KAEvB,IAAK,MAAM7D,KAAchC,EACrB,GLjPkB,KKiPdgC,EAAW7E,OAAf,CAEA,GAAyB,KAArB6E,EAAW7E,OACX,MAAM,IAAI2I,MAAM,sCACpBD,EAAmBrI,EAAYyC,EAAc,CAAE9C,OAAQ6E,EAAW7E,OAAQD,YAAakC,EAAsB4C,EAAWE,WACxH7B,GAAawF,EAAiBpD,MAAMC,EAAExC,OACtCI,GAAauF,EAAiB3C,MAAMyC,UAAUzF,OAC9C0F,GANY,CAShB,GAAoB,YAAhB1I,EAEAmD,GAAa,EAEbC,GAAasF,EAITA,EAAe,KACVC,GAAkBpB,MAAMsB,QAAQ7F,QAAU,GAAK,GAAK2F,GAAkB1D,gBAAgB6D,gBAAkB,EAAEzI,SAASyI,iBACpHzF,GAAe,EACfF,GAAa,EACbC,GAAa,OAIpB,CACD,MAAMM,EAASjC,IAAaL,EAAS0C,IACjCxD,EAAYyC,EAAc,EAAEhD,aAAeO,EAAYyC,EAAc,EAAE7C,QACrEyD,EAAYrD,EAAYyC,EAAc,EAAE5C,WAE9CgD,IAAcL,EAAWE,OAAS,GAAKW,EAAU4B,MAAMC,EAAExC,OAEzDG,GAAaL,EAAWE,OAASU,EAAM6B,MAAMC,EAAExC,OAE/CG,GAAa,EAEbC,IAAcN,EAAWE,OAAS,GAAKW,EAAUqC,MAAMyC,UAAUzF,OAEjEI,GAAaN,EAAWE,OAASU,EAAMsC,MAAMyC,UAAUzF,OAEvDI,GAAa,EAEbA,GAAiC,EAApBN,EAAWE,MAC5B,CACA,MAAO,CAAEG,YAAWC,YAAWC,eACnC,CAKO,SAAS8C,EAAcX,GAC1B,OAAOuD,KAAKC,MAAM,EAAElJ,iBAAmB0F,GAAK,EAAE1F,gBAClD,CL3SI,IAAIyB,EAAcH,EAAS6H,GAAItH,EAAUC,MAAMC,QAAS,IACxD,IAAIN,EAAcH,EAAS0C,IAAKnC,EAAUC,MAAMG,IAAK,KACrD,IAAIR,EAAcH,EAAS8H,IAAKvH,EAAUC,MAAMI,IAAK,KAKpDb,EAASgI,MACThI,EAASiI,UACTjI,EAASkI,KACTlI,EAASmI,KG5Dd,SAAWrH,GACPA,EAAgB,KAAI,YACpBA,EAAiB,MAAI,aACrBA,EAAkB,OAAI,cACtBA,EAAc,GAAI,UAClBA,EAAgB,KAAI,YACpBA,EAAgB,KAAI,YACpBA,EAAsB,WAAI,gBAC7B,CARD,CAQGA,IAAcA,EAAY,CAAC,IGR9BsH,UAAaC,IACT,MAAM,QAAEC,EAAO,aAAE1G,EAAY,SAAEtB,EAAQ,YAAEzB,EAAW,MAAE0J,EAAK,IAAEC,GAAQH,EAAMI,KACrEC,EAAa,IAAItG,MAAMoG,EAAMD,GAC7BI,EAAqB,IAAIvG,MAAM,GACrC,IAAK,IAAIwG,EAAOL,EAAOK,EAAOJ,IAAOI,EACjC,IACI,MAAMC,EAAMP,EAAQM,GACpBF,EAAWE,EAAOL,GAAS7G,EAAoBmH,EAAKjH,EAActB,EAAUzB,EAChF,CACA,MAAOiK,GACH,MACMC,EAAM,8BAA8BH,MAD3BE,aAAerB,MAAQqB,EAAIE,QAAUF,EAAIG,cAExDN,EAAmBxE,KAAK4E,EAC5B,CAEJG,YAAY,CAAEC,SAAUT,EAAYU,SAAUT,K","sources":["webpack://bio/./node_modules/@datagrok-libraries/bio/src/monomer-works/consts.js","webpack://bio/./node_modules/@datagrok-libraries/bio/src/monomer-works/types.js","webpack://bio/./node_modules/@datagrok-libraries/bio/src/utils/macromolecule/consts.js","webpack://bio/./node_modules/@datagrok-libraries/bio/src/utils/const.js","webpack://bio/./node_modules/@datagrok-libraries/bio/src/utils/macromolecule/types.js","webpack://bio/./node_modules/@datagrok-libraries/js-draw-lite/src/types/org.js","webpack://bio/./node_modules/@datagrok-libraries/bio/src/monomer-works/monomer-works.js","webpack://bio/./node_modules/@datagrok-libraries/bio/src/monomer-works/to-atomic-level-utils.js","webpack://bio/./node_modules/@datagrok-libraries/bio/src/monomer-works/seq-to-molfile-worker.js"],"sourcesContent":["export const monomerWorksConsts = {\n // constants for parsing molfile V2000\n V2K_RGP_SHIFT: 8,\n V2K_RGP_LINE: 'M RGP',\n V2K_A_LINE: 'A ',\n // constants for parsing/reconstruction of molfile V3000\n V3K_COUNTS_SHIFT: 14,\n V3K_IDX_SHIFT: 7,\n V3K_HEADER_FIRST_LINE: '\\nDatagrok macromolecule handler\\n\\n',\n V3K_HEADER_SECOND_LINE: ' 0 0 0 0 0 0 999 V3000\\n',\n V3K_BEGIN_CTAB_BLOCK: 'M V30 BEGIN CTAB\\n',\n V3K_END_CTAB_BLOCK: 'M V30 END CTAB\\n',\n V3K_BEGIN_COUNTS_LINE: 'M V30 COUNTS ',\n V3K_COUNTS_LINE_ENDING: ' 0 0 0\\n',\n V3K_BEGIN_ATOM_BLOCK: 'M V30 BEGIN ATOM\\n',\n V3K_END_ATOM_BLOCK: 'M V30 END ATOM\\n',\n V3K_BEGIN_BOND_BLOCK: 'M V30 BEGIN BOND\\n',\n V3K_END_BOND_BLOCK: 'M V30 END BOND\\n',\n V3K_BOND_CONFIG: ' CFG=',\n V3K_BEGIN_DATA_LINE: 'M V30 ',\n V3K_END: 'M END',\n PRECISION_FACTOR: 10000, // HELMCoreLibrary has 4 significant digits after decimal point in atom coordinates\n // symbols for the corresponding monomers in HELM library\n DEOXYRIBOSE: { polymerType: \"RNA\" /* PolymerTypes.RNA */, symbol: 'd' },\n RIBOSE: { polymerType: \"RNA\" /* PolymerTypes.RNA */, symbol: 'r' },\n PHOSPHATE: { polymerType: \"RNA\" /* PolymerTypes.RNA */, symbol: 'p' },\n OXYGEN: 'O',\n HYDROGEN: 'H',\n};\n//# sourceMappingURL=consts.js.map","export function getMolGraph(dict, libKey) {\n return dict[libKey.polymerType]?.[libKey.symbol];\n}\nexport function hasMolGraph(dict, libKey) {\n return !!dict[libKey.polymerType]?.[libKey.symbol];\n}\nexport function setMolGraph(dict, libKey, value) {\n let pt = dict[libKey.polymerType];\n if (!pt)\n pt = dict[libKey.polymerType] = {};\n pt[libKey.symbol] = value;\n}\nexport class MonomerMap extends Map {\n constructor(entries) {\n super(entries);\n }\n}\n/** @property monomers key - helm seq position, */\nexport class MolfileWithMap {\n constructor(molfile, monomers) {\n this.molfile = molfile;\n this.monomers = monomers;\n }\n static createEmpty() { return new MolfileWithMap('', new MonomerMap(null)); }\n}\n//# sourceMappingURL=types.js.map","import { CandidateType } from './types';\n/** enum type to simplify setting \"user-friendly\" notation if necessary */\nexport var NOTATION;\n(function (NOTATION) {\n NOTATION[\"FASTA\"] = \"fasta\";\n NOTATION[\"SEPARATOR\"] = \"separator\";\n NOTATION[\"HELM\"] = \"helm\";\n /* Requires notation handler */ NOTATION[\"CUSTOM\"] = \"custom\";\n /* Requires notation handler */ NOTATION[\"BILN\"] = \"biln\";\n})(NOTATION || (NOTATION = {}));\nexport var ALPHABET;\n(function (ALPHABET) {\n ALPHABET[\"DNA\"] = \"DNA\";\n ALPHABET[\"RNA\"] = \"RNA\";\n ALPHABET[\"PT\"] = \"PT\";\n /** Unknown */\n ALPHABET[\"UN\"] = \"UN\";\n})(ALPHABET || (ALPHABET = {}));\nexport var TAGS;\n(function (TAGS) {\n TAGS[\"aligned\"] = \"aligned\";\n TAGS[\"alphabet\"] = \"alphabet\";\n TAGS[\"alphabetSize\"] = \".alphabetSize\";\n TAGS[\"alphabetIsMultichar\"] = \".alphabetIsMultichar\";\n TAGS[\"separator\"] = \"separator\";\n TAGS[\"isHelmCompatible\"] = \".isHelmCompatible\";\n TAGS[\"positionNames\"] = \".positionNames\";\n TAGS[\"positionLabels\"] = \".positionLabels\";\n TAGS[\"regions\"] = \".regions\";\n TAGS[\"positionShift\"] = \".positionShift\";\n TAGS[\"selectedPosition\"] = \".selectedPosition\";\n TAGS[\"polymerTypeColumnName\"] = \".polymerTypeColumnName\";\n})(TAGS || (TAGS = {}));\nexport { TAGS as BioTags };\nexport const positionSeparator = ', ';\nexport const monomerRe = /(?:\\[([A-Za-z0-9_\\-,()]+)\\])|([A-Za-z\\-])/g;\nexport const helmRe = /(PEPTIDE1|DNA1|RNA1)\\{([^}]+)}/g;\nexport const helmPp1Re = /\\[([^\\[\\]]+)]/g;\nexport const Alphabets = new class {\n constructor() {\n this.fasta = {\n peptide: new Set([\n 'G', 'L', 'Y', 'S', 'E', 'Q', 'D', 'N', 'F', 'A',\n 'K', 'R', 'H', 'C', 'V', 'P', 'W', 'I', 'M', 'T',\n ]),\n dna: new Set(['A', 'C', 'G', 'T']),\n rna: new Set(['A', 'C', 'G', 'U']),\n };\n }\n}();\nexport const candidateAlphabets = [\n new CandidateType(ALPHABET.PT, Alphabets.fasta.peptide, 0.50),\n new CandidateType(ALPHABET.DNA, Alphabets.fasta.dna, 0.55),\n new CandidateType(ALPHABET.RNA, Alphabets.fasta.rna, 0.55),\n];\n/** Canonical gap symbol */\nexport const GAP_SYMBOL = '';\nexport const GapOriginals = {\n [NOTATION.FASTA]: '-',\n [NOTATION.SEPARATOR]: '',\n [NOTATION.HELM]: '*',\n [NOTATION.BILN]: '',\n};\nexport const MONOMER_MOTIF_SPLITTER = ' , ';\nexport const NOTATION_PROVIDER_CONSTRUCTOR_ROLE = 'notationProviderConstructor';\n//# sourceMappingURL=consts.js.map","// core fields of HELM library object used in toAtomicLevel function\nexport const HELM_CORE_FIELDS = [\n \"symbol\" /* HELM_FIELDS.SYMBOL */,\n \"molfile\" /* HELM_FIELDS.MOLFILE */,\n \"rgroups\" /* HELM_FIELDS.RGROUPS */,\n \"name\" /* HELM_FIELDS.NAME */,\n // HELM_FIELDS.MONOMER_TYPE, // add if terminal monomers for PEPTIDEs to be\n // supported\n];\nexport const SDF_MONOMER_NAME = 'MonomerName';\n// todo: ideally, keys should be expressed via constants\nexport const jsonSdfMonomerLibDict = {\n 'monomerType': null, // -> Backbone\n 'smiles': null,\n 'name': 'Name',\n 'author': null,\n 'molfile': 'molecule',\n 'naturalAnalog': 'MonomerNaturalAnalogCode',\n 'rgroups': 'MonomerCaps',\n 'createDate': null,\n 'id': null,\n 'polymerType': 'MonomerType',\n 'symbol': 'MonomerName'\n};\nexport const DUMMY_MONOMER = {\n 'monomerType': 'Backbone',\n 'smiles': '',\n 'name': '',\n 'author': 'Datagrok',\n 'molfile': '',\n 'naturalAnalog': '',\n 'rgroups': [],\n 'createDate': null,\n 'id': 0,\n 'polymerType': 'PEPTIDE',\n 'symbol': ''\n};\n// range of hex nubers used in PepSea library to endode monomers\nexport const MONOMER_ENCODE_MIN = 0x100;\nexport const MONOMER_ENCODE_MAX = 0x40A;\nexport const RIBOSE_SYMBOL = 'r';\nexport const DEOXYRIBOSE_SYMBOL = 'd';\nexport const PHOSPHATE_SYMBOL = 'p';\nexport const HELM_WRAPPERS_REGEXP = new RegExp(`[${RIBOSE_SYMBOL}${DEOXYRIBOSE_SYMBOL}]\\\\((\\\\w)\\\\)${PHOSPHATE_SYMBOL}?`, 'g');\n//# sourceMappingURL=const.js.map","import { NOTATION_PROVIDER_CONSTRUCTOR_ROLE } from './consts';\nexport class NotationProviderBase {\n /** Name of the custom notation */\n static get notationName() {\n return 'Custom';\n }\n ;\n /** flag to let bio know if this provider implements method for converting helm to it */\n static get implementsFromHelm() {\n return false;\n }\n ;\n /** Method for converting HELM to this notation */\n static convertFromHelm(helm, options) {\n throw new Error(`Method convertFromHelm not implemented for this notation provider`);\n }\n ;\n static async getProviderConstructors() {\n // this is terrible, I know, but otherwise this gets put in webworkers and fails due to DG resolution)))\n // @ts-ignore\n if (window?.DG) {\n // @ts-ignore\n const constFuncs = window.DG.Func.find({ meta: { role: NOTATION_PROVIDER_CONSTRUCTOR_ROLE } });\n return Promise.all(constFuncs.map((f) => f.apply({})));\n }\n return [];\n }\n}\n/** Alphabet candidate type */\nexport class CandidateType {\n constructor(name, alphabet, cutoff) {\n this.name = name;\n this.alphabet = alphabet;\n this.cutoff = cutoff;\n }\n}\n/** Alphabet candidate similarity type */\nexport class CandidateSimType extends CandidateType {\n constructor(candidate, freq, similarity) {\n super(candidate.name, candidate.alphabet, candidate.cutoff);\n this.freq = freq;\n this.similarity = similarity;\n }\n}\n//# sourceMappingURL=types.js.map","export var HelmTypes;\n(function (HelmTypes) {\n HelmTypes[\"BASE\"] = \"HELM_BASE\";\n HelmTypes[\"SUGAR\"] = \"HELM_SUGAR\";\n HelmTypes[\"LINKER\"] = \"HELM_LINKER\";\n HelmTypes[\"AA\"] = \"HELM_AA\";\n HelmTypes[\"CHEM\"] = \"HELM_CHEM\";\n HelmTypes[\"BLOB\"] = \"HELM_BLOB\";\n HelmTypes[\"NUCLEOTIDE\"] = \"HELM_NUCLETIDE\";\n})(HelmTypes || (HelmTypes = {}));\n//# sourceMappingURL=org.js.map","import { HelmTypes } from '../helm/consts';\nexport class MonomerWorks {\n //private molfileConverter: IMolfileConverter;\n constructor(monomerLib /*, molfileConverter: IMolfileConverter*/) {\n this.monomerLib = monomerLib;\n //this.molfileConverter = molfileConverter;\n }\n //types according to Monomer possible\n getCappedRotatedMonomer(monomerType, monomerName) {\n // TODO: Check type of monomerType arg\n const monomer = this.monomerLib.getMonomer(monomerType, monomerName);\n if (monomer)\n return monomer.molfile; //TODO cap\n return null;\n }\n}\nexport function helmTypeToPolymerType(helmType) {\n let polymerType = undefined;\n switch (helmType) {\n case HelmTypes.BASE:\n case HelmTypes.SUGAR: // r - ribose, d - deoxyribose\n case HelmTypes.LINKER: // p - phosphate\n case HelmTypes.NUCLEOTIDE:\n // @ts-ignore\n case 'nucleotide':\n polymerType = \"RNA\" /* PolymerTypes.RNA */;\n break;\n case HelmTypes.AA:\n polymerType = \"PEPTIDE\" /* PolymerTypes.PEPTIDE */;\n break;\n case HelmTypes.CHEM:\n polymerType = \"CHEM\" /* PolymerTypes.CHEM */;\n break;\n case HelmTypes.BLOB:\n polymerType = \"BLOB\" /* PolymerTypes.BLOB */;\n break;\n default:\n polymerType = \"PEPTIDE\" /* PolymerTypes.PEPTIDE */;\n console.warn(`Unexpected HelmType '${helmType}'`);\n }\n return polymerType;\n}\n//# sourceMappingURL=monomer-works.js.map","/* eslint-disable max-len */\nimport { monomerWorksConsts as C } from './consts';\nimport { getMolGraph, MolfileWithMap, MonomerMap } from './types';\nimport { HELM_CORE_FIELDS, } from '../utils/const';\nimport { ALPHABET, GAP_SYMBOL } from '../utils/macromolecule/consts';\nimport { helmTypeToPolymerType } from './monomer-works';\n/** Get a mapping of peptide symbols to HELM monomer library objects with selected fields.\n * @param {IMonomerLib} monomerLib - Monomer library\n * @param {HELM_POLYMER_TYPE} polymerType - Polymer type\n * @param {ALPHABET} alphabet - Alphabet of the column\n * @return {Map<string, any>} - Mapping of peptide symbols to HELM monomer library objects with selected fields*/\nexport function getFormattedMonomerLib(monomerLib, polymerType, alphabet) {\n const map = new Map();\n for (const monomerSymbol of monomerLib.getMonomerSymbolsByType(polymerType)) {\n const it = monomerLib.getMonomer(polymerType, monomerSymbol);\n if (polymerType === \"RNA\" /* HELM_POLYMER_TYPE.RNA */ &&\n (it[\"monomerType\" /* HELM_FIELDS.MONOMER_TYPE */] === \"Branch\" /* HELM_MONOMER_TYPE.BRANCH */ ||\n alphabet === ALPHABET.DNA && it[\"symbol\" /* HELM_FIELDS.SYMBOL */] === C.DEOXYRIBOSE.symbol ||\n alphabet === ALPHABET.RNA && it[\"symbol\" /* HELM_FIELDS.SYMBOL */] === C.RIBOSE.symbol ||\n it[\"symbol\" /* HELM_FIELDS.SYMBOL */] === C.PHOSPHATE.symbol) ||\n polymerType === \"PEPTIDE\" /* HELM_POLYMER_TYPE.PEPTIDE */ &&\n it[\"monomerType\" /* HELM_FIELDS.MONOMER_TYPE */] !== \"Branch\" /* HELM_MONOMER_TYPE.BRANCH */) {\n const monomerObject = {};\n HELM_CORE_FIELDS.forEach((field) => {\n //@ts-ignore\n monomerObject[field] = it[field];\n });\n map.set(monomerSymbol, monomerObject);\n }\n }\n return map;\n}\n/** Translate a sequence of monomer symbols into Molfile V3000\n * @param {string[]} monomerSeq - Sequence of monomer symbols (canonical)\n * @param {Map<string, MolGraph>} monomersDict - Mapping of monomer symbols to MolGraph objects\n * @param {ALPHABET} alphabet - Alphabet of the column\n * @param {PolymerType} polymerType - Polymer type\n * @return {string} - Molfile V3000*/\nexport function monomerSeqToMolfile(monomerSeq, monomersDict, alphabet, polymerType) {\n if (monomerSeq.length === 0) {\n // throw new Error('monomerSeq is empty');\n return MolfileWithMap.createEmpty();\n }\n // define atom and bond counts, taking into account the bond type\n const getAtomAndBondCounts = getResultingAtomBondCounts;\n const { atomCount, bondCount, needsCapping } = getAtomAndBondCounts(monomerSeq, monomersDict, alphabet, polymerType);\n // create arrays to store lines of the resulting molfile\n const molfileAtomBlock = new Array(atomCount);\n const molfileBondBlock = new Array(bondCount);\n let addMonomerToMolblock;\n let sugar = null;\n let phosphate = null;\n if (polymerType === \"PEPTIDE\" /* HELM_POLYMER_TYPE.PEPTIDE */)\n addMonomerToMolblock = addAminoAcidToMolblock;\n else { // nucleotides\n addMonomerToMolblock = addNucleotideToMolblock;\n sugar = (alphabet === ALPHABET.DNA) ? getMolGraph(monomersDict, C.DEOXYRIBOSE) : getMolGraph(monomersDict, C.RIBOSE);\n phosphate = getMolGraph(monomersDict, C.PHOSPHATE);\n }\n const v = {\n i: 0,\n nodeShift: 0,\n bondShift: 0,\n backbonePositionShift: new Array(2).fill(0),\n branchPositionShift: new Array(2).fill(0),\n backboneAttachNode: 0,\n branchAttachNode: 0,\n flipFactor: 1,\n };\n const LC = {\n sugar: sugar,\n phosphate: phosphate,\n seqLength: monomerSeq.length,\n atomCount: atomCount,\n bondCount: bondCount,\n };\n const monomers = new MonomerMap();\n const steabsCollection = [];\n let nAtoms = 0;\n let lastMonomerCappingAtom = undefined;\n for (v.i = 0; v.i < LC.seqLength; ++v.i) {\n const seqMonomer = monomerSeq[v.i];\n if (seqMonomer.symbol === GAP_SYMBOL)\n continue;\n const monomer = getMolGraph(monomersDict, { symbol: seqMonomer.symbol, polymerType: helmTypeToPolymerType(seqMonomer.biotype) });\n lastMonomerCappingAtom = monomer.terminalR2Atom;\n const mAtomFirst = v.nodeShift;\n const mBondFirst = v.bondShift;\n addMonomerToMolblock(monomer, molfileAtomBlock, molfileBondBlock, v, LC);\n //adding stereo atoms to array for further STEABS block generation\n monomer.stereoAtoms?.forEach((i) => steabsCollection.push(i + nAtoms));\n nAtoms += monomer.atoms.x.length;\n const mAtomCount = v.nodeShift - mAtomFirst;\n const mAtomList = new Array(mAtomCount);\n for (let maI = 0; maI < mAtomCount; ++maI)\n mAtomList[maI] = mAtomFirst + maI;\n const mBondCount = v.bondShift - mBondFirst;\n const mBondList = new Array(mBondCount);\n for (let mbI = 0; mbI < mBondCount; ++mbI)\n mBondList[mbI] = mBondFirst + mbI;\n monomers.set(v.i, {\n biotype: seqMonomer.biotype,\n symbol: seqMonomer.symbol,\n atoms: mAtomList, bonds: mBondList\n });\n }\n // if the last monomer needs to be capped, add the terminal OH to the resulting molfile\n if (needsCapping)\n capResultingMolblock(molfileAtomBlock, molfileBondBlock, v, LC, lastMonomerCappingAtom ?? C.OXYGEN);\n const molfileCountsLine = C.V3K_BEGIN_COUNTS_LINE + atomCount + ' ' + bondCount + C.V3K_COUNTS_LINE_ENDING;\n // todo: possible optimization may be achieved by replacing .join('') with +=\n // since counterintuitively joining an array into a new string is reportedly\n // slower than using += as below\n let result = '';\n result += C.V3K_HEADER_FIRST_LINE;\n result += C.V3K_HEADER_SECOND_LINE;\n result += C.V3K_BEGIN_CTAB_BLOCK;\n result += molfileCountsLine;\n result += C.V3K_BEGIN_ATOM_BLOCK;\n result += molfileAtomBlock.join('');\n result += C.V3K_END_ATOM_BLOCK;\n result += C.V3K_BEGIN_BOND_BLOCK;\n result += molfileBondBlock.join('');\n result += C.V3K_END_BOND_BLOCK;\n if (steabsCollection.length > 0)\n result += getCollectionBlock(steabsCollection);\n result += C.V3K_END_CTAB_BLOCK;\n result += C.V3K_END;\n // return molfileParts.join('');\n return { molfile: result, monomers: monomers };\n}\nfunction getCollectionBlock(collection) {\n //one row in STEABS block can be no longer than 80 symbols\n //maxSymbols = 80 symbols minus ' -\\n' (4 symbols)\n const maxSymbols = 76;\n const rowsArray = [];\n let newCollectionRow = `M V30 MDLV30/STEABS ATOMS=(${collection.length}`;\n for (let i = 0; i < collection.length; i++) {\n const updatedRow = `${newCollectionRow} ${collection[i]}`;\n if (updatedRow.length > maxSymbols) {\n rowsArray.push(`${newCollectionRow} -\\n`);\n newCollectionRow = `M V30 ${collection[i]}`;\n }\n else\n newCollectionRow = updatedRow;\n //in case last atom was added - close the block\n if (i === collection.length - 1)\n rowsArray.push(`${newCollectionRow})\\n`);\n }\n return `M V30 BEGIN COLLECTION\\n${rowsArray.join('')}M V30 END COLLECTION\\n`;\n}\n/** Cap the resulting (after sewing up all the monomers) molfile with 'O'\n * @param {string[]} molfileAtomBlock - Array of lines of the resulting molfile atom block\n * @param {string[]} molfileBondBlock - Array of lines of the resulting molfile bond block\n * @param {LoopVariables} v - Loop variables\n * @param {LoopConstants} LC - Loop constants*/\nfunction capResultingMolblock(molfileAtomBlock, molfileBondBlock, v, LC, cappingAtomType = C.OXYGEN) {\n // add terminal oxygen\n const atomIdx = v.nodeShift + 1;\n molfileAtomBlock[LC.atomCount] = C.V3K_BEGIN_DATA_LINE + atomIdx + ' ' +\n (cappingAtomType ?? C.OXYGEN) + ' ' + keepPrecision(v.backbonePositionShift[0]) + ' ' +\n v.flipFactor * keepPrecision(v.backbonePositionShift[1]) + ' ' + '0.000000 0' + '\\n';\n // add terminal bond\n const firstAtom = v.backboneAttachNode;\n const secondAtom = atomIdx;\n molfileBondBlock[LC.bondCount] = C.V3K_BEGIN_DATA_LINE + v.bondShift + ' ' +\n 1 + ' ' + firstAtom + ' ' + secondAtom + '\\n';\n}\nfunction addAminoAcidToMolblock(monomer, molfileAtomBlock, molfileBondBlock, v) {\n v.flipFactor = (-1) ** (v.i % 2); // to flip every even monomer over OX\n addBackboneMonomerToMolblock(monomer, molfileAtomBlock, molfileBondBlock, v);\n}\nfunction addBackboneMonomerToMolblock(monomer, molfileAtomBlock, molfileBondBlock, v) {\n // todo: remove these comments to the docstrings of the corr. functions\n // construnct the lines of V3K molfile atom block\n fillAtomLines(monomer, molfileAtomBlock, v);\n // construct the lines of V3K molfile bond block\n fillBondLines(monomer, molfileBondBlock, v);\n // peptide bond\n fillChainExtendingBond(monomer, molfileBondBlock, v);\n // update branch variables if necessary\n if (monomer.meta.branchShift !== null && monomer.meta.terminalNodes.length > 2)\n updateBranchVariables(monomer, v);\n // update loop variables\n updateChainExtendingVariables(monomer, v);\n}\nfunction addNucleotideToMolblock(nucleobase, molfileAtomBlock, molfileBondBlock, v, LC) {\n // construnct the lines of V3K molfile atom block corresponding to phosphate\n // and sugar\n if (v.i === 0)\n addBackboneMonomerToMolblock(LC.sugar, molfileAtomBlock, molfileBondBlock, v);\n else {\n for (const monomer of [LC.phosphate, LC.sugar])\n addBackboneMonomerToMolblock(monomer, molfileAtomBlock, molfileBondBlock, v);\n }\n addBranchMonomerToMolblock(nucleobase, molfileAtomBlock, molfileBondBlock, v);\n}\nfunction addBranchMonomerToMolblock(monomer, molfileAtomBlock, molfileBondBlock, v) {\n fillBranchAtomLines(monomer, molfileAtomBlock, v);\n fillBondLines(monomer, molfileBondBlock, v);\n fillBackboneToBranchBond(monomer, molfileBondBlock, v);\n // C-N bond\n const bondIdx = v.bondShift;\n const firstAtom = v.branchAttachNode;\n const secondAtom = monomer.meta.terminalNodes[0] + v.nodeShift;\n molfileBondBlock[bondIdx - 1] = C.V3K_BEGIN_DATA_LINE + bondIdx + ' ' +\n 1 + ' ' + firstAtom + ' ' + secondAtom + '\\n';\n // update loop variables\n v.bondShift += monomer.bonds.atomPairs.length + 1;\n v.nodeShift += monomer.atoms.atomTypes.length;\n}\nfunction updateChainExtendingVariables(monomer, v) {\n v.backboneAttachNode = v.nodeShift + monomer.meta.terminalNodes[1];\n v.bondShift += monomer.bonds.atomPairs.length + 1;\n v.nodeShift += monomer.atoms.atomTypes.length;\n v.backbonePositionShift[0] += monomer.meta.backboneShift?.[0] ?? 0; // todo: non-null check\n v.backbonePositionShift[1] += v.flipFactor * (monomer.meta.backboneShift?.[1] ?? 0);\n}\nfunction updateBranchVariables(monomer, v) {\n v.branchAttachNode = v.nodeShift + monomer.meta.terminalNodes[2];\n for (let i = 0; i < 2; ++i)\n v.branchPositionShift[i] = v.backbonePositionShift[i] + monomer.meta.branchShift[i];\n}\nfunction fillAtomLines(monomer, molfileAtomBlock, v) {\n for (let j = 0; j < monomer.atoms.atomTypes.length; ++j) {\n const atomIdx = v.nodeShift + j + 1;\n molfileAtomBlock[v.nodeShift + j] = C.V3K_BEGIN_DATA_LINE + atomIdx + ' ' +\n monomer.atoms.atomTypes[j] + ' ' +\n keepPrecision(v.backbonePositionShift[0] + monomer.atoms.x[j]) + ' ' +\n keepPrecision(v.backbonePositionShift[1] + v.flipFactor * monomer.atoms.y[j]) +\n ' ' + monomer.atoms.kwargs[j];\n }\n}\n// todo: remove as quickfix\nfunction fillBranchAtomLines(monomer, molfileAtomBlock, v) {\n for (let j = 0; j < monomer.atoms.atomTypes.length; ++j) {\n const atomIdx = v.nodeShift + j + 1;\n molfileAtomBlock[v.nodeShift + j] = C.V3K_BEGIN_DATA_LINE + atomIdx + ' ' +\n monomer.atoms.atomTypes[j] + ' ' +\n keepPrecision(v.branchPositionShift[0] + monomer.atoms.x[j]) + ' ' +\n keepPrecision(v.branchPositionShift[1] + v.flipFactor * monomer.atoms.y[j]) +\n ' ' + monomer.atoms.kwargs[j];\n }\n}\nfunction fillBondLines(monomer, molfileBondBlock, v) {\n // construct the lines of V3K molfile bond block\n for (let j = 0; j < monomer.bonds.atomPairs.length; ++j) {\n const bondIdx = v.bondShift + j + 1;\n const firstAtom = monomer.bonds.atomPairs[j][0] + v.nodeShift;\n const secondAtom = monomer.bonds.atomPairs[j][1] + v.nodeShift;\n let bondCfg = '';\n if (monomer.bonds.bondConfiguration.has(j)) {\n // flip orientation when necessary\n let orientation = monomer.bonds.bondConfiguration.get(j);\n if (v.flipFactor < 0)\n orientation = (orientation === 1) ? 3 : 1;\n bondCfg = ' CFG=' + orientation;\n }\n const kwargs = monomer.bonds.kwargs.has(j) ?\n ' ' + monomer.bonds.kwargs.get(j) : '';\n molfileBondBlock[v.bondShift + j] = C.V3K_BEGIN_DATA_LINE + bondIdx + ' ' +\n monomer.bonds.bondTypes[j] + ' ' +\n firstAtom + ' ' + secondAtom + bondCfg + kwargs + '\\n';\n }\n}\nfunction fillChainExtendingBond(monomer, molfileBondBlock, v) {\n if (v.backboneAttachNode !== 0) {\n const bondIdx = v.bondShift;\n const firstAtom = v.backboneAttachNode;\n const secondAtom = monomer.meta.terminalNodes[0] + v.nodeShift;\n molfileBondBlock[v.bondShift - 1] = C.V3K_BEGIN_DATA_LINE + bondIdx + ' ' +\n 1 + ' ' + firstAtom + ' ' + secondAtom + '\\n';\n }\n}\n// todo: remove\nfunction fillBackboneToBranchBond(branchMonomer, molfileBondBlock, v) {\n const bondIdx = v.bondShift;\n const firstAtom = v.branchAttachNode;\n const secondAtom = branchMonomer.meta.terminalNodes[0] + v.nodeShift;\n molfileBondBlock[bondIdx - 1] = C.V3K_BEGIN_DATA_LINE + bondIdx + ' ' +\n 1 + ' ' + firstAtom + ' ' + secondAtom + '\\n';\n}\n/** Compute the atom/bond counts for the resulting molfile, depending on the\n * type of polymer (peptide/nucleotide)\n * @param {string[]}monomerSeq - the sequence of monomers\n * @param {Map<string, MolGraph>}monomersDict - the dictionary of monomers\n * @param {ALPHABET}alphabet - the alphabet of the monomers\n * @param {HELM_POLYMER_TYPE}polymerType - the type of polymer\n * @return {{atomCount: number, bondCount: number}} - the atom/bond counts*/\nfunction getResultingAtomBondCounts(monomerSeq, monomersDict, alphabet, polymerType) {\n let atomCount = 0;\n let bondCount = 0;\n let monomerCount = 0;\n let needsCapping = true;\n let lastMonomerGraph = null;\n // sum up all the atoms/nodes provided by the sequence\n for (const seqMonomer of monomerSeq) {\n if (seqMonomer.symbol === GAP_SYMBOL)\n continue; // Skip for gap/empty monomer in MSA\n if (seqMonomer.symbol == '*')\n throw new Error(`Gap canonical symbol is '', not '*`);\n lastMonomerGraph = getMolGraph(monomersDict, { symbol: seqMonomer.symbol, polymerType: helmTypeToPolymerType(seqMonomer.biotype) });\n atomCount += lastMonomerGraph.atoms.x.length;\n bondCount += lastMonomerGraph.bonds.bondTypes.length;\n monomerCount++;\n }\n // add extra values depending on the polymer type\n if (polymerType === \"PEPTIDE\" /* HELM_POLYMER_TYPE.PEPTIDE */) {\n // add the rightmost/terminating cap group 'OH' (i.e. 'O')\n atomCount += 1;\n // add chain-extending bonds (C-NH per each monomer pair and terminal C-OH)\n bondCount += monomerCount;\n // if the last monomer is something like NH2, which only has R1, there is no need to cap it\n // although, this should never happen, but hey... in other bits of code, there is a chunk that adds pseudo-R2 as hydrogen\n // we should also check, if the R2 of the last monomer is not hydrogen, that case should also be omitted\n if (monomerCount > 0) {\n if ((lastMonomerGraph?.meta?.rNodes?.length ?? 0) < 2 || lastMonomerGraph?.terminalR2Atom?.toLowerCase() === C.HYDROGEN.toLowerCase()) {\n needsCapping = false;\n atomCount -= 1; // remove the last atom (the terminal 'O')\n bondCount -= 1; // remove the last bond (the terminal C-OH)\n }\n }\n }\n else { // nucleotides\n const sugar = (alphabet === ALPHABET.DNA) ?\n getMolGraph(monomersDict, C.DEOXYRIBOSE) : getMolGraph(monomersDict, C.RIBOSE);\n const phosphate = getMolGraph(monomersDict, C.PHOSPHATE);\n // add phosphate per each pair of nucleobase symbols\n atomCount += (monomerSeq.length - 1) * phosphate.atoms.x.length;\n // add sugar per each nucleobase symbol\n atomCount += monomerSeq.length * sugar.atoms.x.length;\n // add the leftmost cap group 'OH' (i.e. 'O')\n atomCount += 1;\n // add bonds from phosphate monomers\n bondCount += (monomerSeq.length - 1) * phosphate.bonds.bondTypes.length;\n // add bonds from sugar monomers\n bondCount += monomerSeq.length * sugar.bonds.bondTypes.length;\n // exclude the first chain-extending bond O-P (absent, no 'leftmost' phosphate)\n bondCount -= 1;\n // add chain-extending and branch bonds (O-P, C-O and C-N per each nucleotide)\n bondCount += monomerSeq.length * 3;\n }\n return { atomCount, bondCount, needsCapping };\n}\n/** Keep precision upon floating point operations over atom coordinates\n * @param {number}x - the floating point number\n * @return {number} - the floating point number with the same precision\n */\nexport function keepPrecision(x) {\n return Math.round(C.PRECISION_FACTOR * x) / C.PRECISION_FACTOR;\n}\n//# sourceMappingURL=to-atomic-level-utils.js.map","import { monomerSeqToMolfile } from './to-atomic-level-utils';\nonmessage = (event) => {\n const { seqList, monomersDict, alphabet, polymerType, start, end } = event.data;\n const resMolList = new Array(end - start);\n const molfileWarningList = new Array(0);\n for (let rowI = start; rowI < end; ++rowI) {\n try {\n const seq = seqList[rowI];\n resMolList[rowI - start] = monomerSeqToMolfile(seq, monomersDict, alphabet, polymerType);\n }\n catch (err) {\n const errMsg = err instanceof Error ? err.message : err.toString();\n const msg = `Cannot get molfile of row #${rowI}: ${errMsg}.`;\n molfileWarningList.push(msg);\n }\n }\n postMessage({ molfiles: resMolList, warnings: molfileWarningList });\n};\n//# sourceMappingURL=seq-to-molfile-worker.js.map"],"names":["monomerWorksConsts","V2K_RGP_SHIFT","V2K_RGP_LINE","V2K_A_LINE","V3K_COUNTS_SHIFT","V3K_IDX_SHIFT","V3K_HEADER_FIRST_LINE","V3K_HEADER_SECOND_LINE","V3K_BEGIN_CTAB_BLOCK","V3K_END_CTAB_BLOCK","V3K_BEGIN_COUNTS_LINE","V3K_COUNTS_LINE_ENDING","V3K_BEGIN_ATOM_BLOCK","V3K_END_ATOM_BLOCK","V3K_BEGIN_BOND_BLOCK","V3K_END_BOND_BLOCK","V3K_BOND_CONFIG","V3K_BEGIN_DATA_LINE","V3K_END","PRECISION_FACTOR","DEOXYRIBOSE","polymerType","symbol","RIBOSE","PHOSPHATE","OXYGEN","HYDROGEN","getMolGraph","dict","libKey","MonomerMap","Map","constructor","entries","super","MolfileWithMap","molfile","monomers","this","createEmpty","NOTATION","ALPHABET","TAGS","RegExp","CandidateType","name","alphabet","cutoff","Alphabets","fasta","peptide","Set","dna","rna","HelmTypes","helmTypeToPolymerType","helmType","BASE","SUGAR","LINKER","NUCLEOTIDE","AA","CHEM","BLOB","console","warn","monomerSeqToMolfile","monomerSeq","monomersDict","length","getAtomAndBondCounts","getResultingAtomBondCounts","atomCount","bondCount","needsCapping","molfileAtomBlock","Array","molfileBondBlock","addMonomerToMolblock","sugar","phosphate","addAminoAcidToMolblock","addNucleotideToMolblock","DNA","v","i","nodeShift","bondShift","backbonePositionShift","fill","branchPositionShift","backboneAttachNode","branchAttachNode","flipFactor","LC","seqLength","steabsCollection","lastMonomerCappingAtom","nAtoms","seqMonomer","monomer","biotype","terminalR2Atom","mAtomFirst","mBondFirst","stereoAtoms","forEach","push","atoms","x","mAtomCount","mAtomList","maI","mBondCount","mBondList","mbI","set","bonds","cappingAtomType","atomIdx","keepPrecision","firstAtom","secondAtom","capResultingMolblock","molfileCountsLine","result","join","collection","rowsArray","newCollectionRow","updatedRow","getCollectionBlock","addBackboneMonomerToMolblock","j","atomTypes","y","kwargs","fillAtomLines","fillBondLines","bondIdx","meta","terminalNodes","fillChainExtendingBond","branchShift","updateBranchVariables","atomPairs","backboneShift","updateChainExtendingVariables","nucleobase","fillBranchAtomLines","branchMonomer","fillBackboneToBranchBond","addBranchMonomerToMolblock","bondCfg","bondConfiguration","has","orientation","get","bondTypes","monomerCount","lastMonomerGraph","Error","rNodes","toLowerCase","Math","round","PT","RNA","FASTA","SEPARATOR","HELM","BILN","onmessage","event","seqList","start","end","data","resMolList","molfileWarningList","rowI","seq","err","msg","message","toString","postMessage","molfiles","warnings"],"sourceRoot":""}
1
+ {"version":3,"file":"455.js","mappings":"2BAAO,MAAMA,EAAqB,CAE9BC,cAAe,EACfC,aAAc,SACdC,WAAY,MAEZC,iBAAkB,GAClBC,cAAe,EACfC,sBAAuB,uCACvBC,uBAAwB,4CACxBC,qBAAsB,sBACtBC,mBAAoB,oBACpBC,sBAAuB,iBACvBC,uBAAwB,WACxBC,qBAAsB,sBACtBC,mBAAoB,oBACpBC,qBAAsB,sBACtBC,mBAAoB,oBACpBC,gBAAiB,QACjBC,oBAAqB,UACrBC,QAAS,SACTC,iBAAkB,IAElBC,YAAa,CAAEC,YAAa,MAA8BC,OAAQ,KAClEC,OAAQ,CAAEF,YAAa,MAA8BC,OAAQ,KAC7DE,UAAW,CAAEH,YAAa,MAA8BC,OAAQ,KAChEG,OAAQ,IACRC,SAAU,KC3BP,SAASC,EAAYC,EAAMC,GAC9B,OAAOD,EAAKC,EAAOR,eAAeQ,EAAOP,OAC7C,CAUO,MAAMQ,UAAmBC,IAC5B,WAAAC,CAAYC,GACRC,MAAMD,EACV,EAGG,MAAME,EACT,WAAAH,CAAYI,EAASC,GACjBC,KAAKF,QAAUA,EACfE,KAAKD,SAAWA,CACpB,CACA,kBAAOE,GAAgB,OAAO,IAAIJ,EAAe,GAAI,IAAIL,EAAW,MAAQ,E,ICrBrEU,EAQAC,EAQAC,ECyByB,IAAIC,OAAO,oBAA0E,KCdlH,MAAMC,EACT,WAAAZ,CAAYa,EAAMC,EAAUC,GACxBT,KAAKO,KAAOA,EACZP,KAAKQ,SAAWA,EAChBR,KAAKS,OAASA,CAClB,GF/BJ,SAAWP,GACPA,EAAgB,MAAI,QACpBA,EAAoB,UAAI,YACxBA,EAAe,KAAI,OACaA,EAAiB,OAAI,SACrBA,EAAe,KAAI,MACtD,CAND,CAMGA,IAAaA,EAAW,CAAC,IAE5B,SAAWC,GACPA,EAAc,IAAI,MAClBA,EAAc,IAAI,MAClBA,EAAa,GAAI,KAEjBA,EAAa,GAAI,IACpB,CAND,CAMGA,IAAaA,EAAW,CAAC,IAE5B,SAAWC,GACPA,EAAc,QAAI,UAClBA,EAAe,SAAI,WACnBA,EAAmB,aAAI,gBACvBA,EAA0B,oBAAI,uBAC9BA,EAAgB,UAAI,YACpBA,EAAuB,iBAAI,oBAC3BA,EAAoB,cAAI,iBACxBA,EAAqB,eAAI,kBACzBA,EAAc,QAAI,WAClBA,EAAoB,cAAI,iBACxBA,EAAuB,iBAAI,oBAC3BA,EAA4B,sBAAI,yBAChCA,EAAkB,YAAI,eACtBA,EAAsB,gBAAI,mBAC1BA,EAA2B,qBAAI,uBAClC,CAhBD,CAgBGA,IAASA,EAAO,CAAC,IAEb,MAIMM,EAAY,IAAI,MACzB,WAAAhB,GACIM,KAAKW,MAAQ,CACTC,QAAS,IAAIC,IAAI,CACb,IAAK,IAAK,IAAK,IAAK,IAAK,IAAK,IAAK,IAAK,IAAK,IAC7C,IAAK,IAAK,IAAK,IAAK,IAAK,IAAK,IAAK,IAAK,IAAK,MAEjDC,IAAK,IAAID,IAAI,CAAC,IAAK,IAAK,IAAK,MAC7BE,IAAK,IAAIF,IAAI,CAAC,IAAK,IAAK,IAAK,MAErC,GGnDG,IAAIG,ECgBJ,SAASC,EAAsBC,GAClC,IAAInC,EACJ,OAAQmC,GACJ,KAAKF,EAAUG,KACf,KAAKH,EAAUI,MACf,KAAKJ,EAAUK,OACf,KAAKL,EAAUM,WAEf,IAAK,aACDvC,EAAc,MACd,MACJ,KAAKiC,EAAUO,GACXxC,EAAc,UACd,MACJ,KAAKiC,EAAUQ,KACXzC,EAAc,OACd,MACJ,KAAKiC,EAAUS,KACX1C,EAAc,OACd,MACJ,QACIA,EAAc,UACd2C,QAAQC,KAAK,wBAAwBT,MAE7C,OAAOnC,CACX,CCHO,SAAS6C,EAAoBC,EAAYC,EAActB,EAAUzB,GACpE,GAA0B,IAAtB8C,EAAWE,OAEX,OAAOlC,EAAeI,cAG1B,MAAM+B,EAAuBC,GACvB,UAAEC,EAAS,UAAEC,EAAS,aAAEC,GAAiBJ,EAAqBH,EAAYC,EAActB,EAAUzB,GAElGsD,EAAmB,IAAIC,MAAMJ,GAC7BK,EAAmB,IAAID,MAAMH,GACnC,IAAIK,EACAC,EAAQ,KACRC,EAAY,KACI,YAAhB3D,EACAyD,EAAuBG,GAEvBH,EAAuBI,EACvBH,EAASjC,IAAaL,EAAS0C,IAAOxD,EAAYyC,EAAc,EAAEhD,aAAeO,EAAYyC,EAAc,EAAE7C,QAC7GyD,EAAYrD,EAAYyC,EAAc,EAAE5C,YAE5C,MAAM4D,EAAI,CACNC,EAAG,EACHC,UAAW,EACXC,UAAW,EACXC,sBAAuB,IAAIZ,MAAM,GAAGa,KAAK,GACzCC,oBAAqB,IAAId,MAAM,GAAGa,KAAK,GACvCE,mBAAoB,EACpBC,iBAAkB,EAClBC,WAAY,GAEVC,EAAK,CACPf,MAAOA,EACPC,UAAWA,EACXe,UAAW5B,EAAWE,OACtBG,UAAWA,EACXC,UAAWA,GAETpC,EAAW,IAAIP,EACfkE,EAAmB,GACzB,IACIC,EADAC,EAAS,EAEb,IAAKd,EAAEC,EAAI,EAAGD,EAAEC,EAAIS,EAAGC,YAAaX,EAAEC,EAAG,CACrC,MAAMc,EAAahC,EAAWiB,EAAEC,GAChC,GLvBkB,KKuBdc,EAAW7E,OACX,SACJ,MAAM8E,EAAUzE,EAAYyC,EAAc,CAAE9C,OAAQ6E,EAAW7E,OAAQD,YAAakC,EAAsB4C,EAAWE,WACrHJ,EAAyBG,EAAQE,eACjC,MAAMC,EAAanB,EAAEE,UACfkB,EAAapB,EAAEG,UACrBT,EAAqBsB,EAASzB,EAAkBE,EAAkBO,EAAGU,GAErEM,EAAQK,aAAaC,QAASrB,GAAMW,EAAiBW,KAAKtB,EAAIa,IAC9DA,GAAUE,EAAQQ,MAAMC,EAAExC,OAC1B,MAAMyC,EAAa1B,EAAEE,UAAYiB,EAC3BQ,EAAY,IAAInC,MAAMkC,GAC5B,IAAK,IAAIE,EAAM,EAAGA,EAAMF,IAAcE,EAClCD,EAAUC,GAAOT,EAAaS,EAClC,MAAMC,EAAa7B,EAAEG,UAAYiB,EAC3BU,EAAY,IAAItC,MAAMqC,GAC5B,IAAK,IAAIE,EAAM,EAAGA,EAAMF,IAAcE,EAClCD,EAAUC,GAAOX,EAAaW,EAClC9E,EAAS+E,IAAIhC,EAAEC,EAAG,CACdgB,QAASF,EAAWE,QACpB/E,OAAQ6E,EAAW7E,OACnBsF,MAAOG,EAAWM,MAAOH,GAEjC,CAEIxC,GAiDR,SAA8BC,EAAkBE,EAAkBO,EAAGU,EAAIwB,EAAkB,EAAE7F,QAEzF,MAAM8F,EAAUnC,EAAEE,UAAY,EAC9BX,EAAiBmB,EAAGtB,WAAa,EAAEvD,oBAAsBsG,EAAU,KAC9DD,GAAmB,EAAE7F,QAAU,IAAM+F,EAAcpC,EAAEI,sBAAsB,IAAM,IAClFJ,EAAES,WAAa2B,EAAcpC,EAAEI,sBAAsB,IAFxB,gBAIjC,MAAMiC,EAAYrC,EAAEO,mBACd+B,EAAaH,EACnB1C,EAAiBiB,EAAGrB,WAAa,EAAExD,oBAAsBmE,EAAEG,UAA1B,MACnBkC,EAAY,IAAMC,EAAa,IACjD,CA3DQC,CAAqBhD,EAAkBE,EAAkBO,EAAGU,EAAIG,GAA0B,EAAExE,QAChG,MAAMmG,EAAoB,EAAElH,sBAAwB8D,EAAY,IAAMC,EAAY,EAAE9D,uBAIpF,IAAIkH,EAAS,GAgBb,OAfAA,GAAU,EAAEvH,sBACZuH,GAAU,EAAEtH,uBACZsH,GAAU,EAAErH,qBACZqH,GAAUD,EACVC,GAAU,EAAEjH,qBACZiH,GAAUlD,EAAiBmD,KAAK,IAChCD,GAAU,EAAEhH,mBACZgH,GAAU,EAAE/G,qBACZ+G,GAAUhD,EAAiBiD,KAAK,IAChCD,GAAU,EAAE9G,mBACRiF,EAAiB3B,OAAS,IAC1BwD,GAMR,SAA4BE,GAGxB,MACMC,EAAY,GAClB,IAAIC,EAAmB,+BAA+BF,EAAW1D,SACjE,IAAK,IAAIgB,EAAI,EAAGA,EAAI0C,EAAW1D,OAAQgB,IAAK,CACxC,MAAM6C,EAAa,GAAGD,KAAoBF,EAAW1C,KACjD6C,EAAW7D,OALA,IAMX2D,EAAUrB,KAAK,GAAGsB,SAClBA,EAAmB,UAAUF,EAAW1C,MAGxC4C,EAAmBC,EAEnB7C,IAAM0C,EAAW1D,OAAS,GAC1B2D,EAAUrB,KAAK,GAAGsB,OAC1B,CACA,MAAO,4BAA4BD,EAAUF,KAAK,4BACtD,CAzBkBK,CAAmBnC,IACjC6B,GAAU,EAAEpH,mBACZoH,GAAU,EAAE3G,QAEL,CAAEkB,QAASyF,EAAQxF,SAAUA,EACxC,CAsCA,SAAS4C,EAAuBmB,EAASzB,EAAkBE,EAAkBO,GACzEA,EAAES,aAAe,KAAOT,EAAEC,EAAI,GAC9B+C,EAA6BhC,EAASzB,EAAkBE,EAAkBO,EAC9E,CACA,SAASgD,EAA6BhC,EAASzB,EAAkBE,EAAkBO,IAmDnF,SAAuBgB,EAASzB,EAAkBS,GAC9C,IAAK,IAAIiD,EAAI,EAAGA,EAAIjC,EAAQQ,MAAM0B,UAAUjE,SAAUgE,EAAG,CACrD,MAAMd,EAAUnC,EAAEE,UAAY+C,EAAI,EAClC1D,EAAiBS,EAAEE,UAAY+C,GAAK,EAAEpH,oBAAsBsG,EAAU,IAClEnB,EAAQQ,MAAM0B,UAAUD,GAAK,IAC7Bb,EAAcpC,EAAEI,sBAAsB,GAAKY,EAAQQ,MAAMC,EAAEwB,IAAM,IACjEb,EAAcpC,EAAEI,sBAAsB,GAAKJ,EAAES,WAAaO,EAAQQ,MAAM2B,EAAEF,IAC1E,IAAMjC,EAAQQ,MAAM4B,OAAOH,EACnC,CACJ,CAzDII,CAAcrC,EAASzB,EAAkBS,GAEzCsD,EAActC,EAASvB,EAAkBO,GAwF7C,SAAgCgB,EAASvB,EAAkBO,GACvD,GAA6B,IAAzBA,EAAEO,mBAA0B,CAC5B,MAAMgD,EAAUvD,EAAEG,UACZkC,EAAYrC,EAAEO,mBACd+B,EAAatB,EAAQwC,KAAKC,cAAc,GAAKzD,EAAEE,UACrDT,EAAiBO,EAAEG,UAAY,GAAK,EAAEtE,oBAAsB0H,EAAxB,MACtBlB,EAAY,IAAMC,EAAa,IACjD,CACJ,CA9FIoB,CAAuB1C,EAASvB,EAAkBO,GAEjB,OAA7BgB,EAAQwC,KAAKG,aAAwB3C,EAAQwC,KAAKC,cAAcxE,OAAS,GAqCjF,SAA+B+B,EAAShB,GACpCA,EAAEQ,iBAAmBR,EAAEE,UAAYc,EAAQwC,KAAKC,cAAc,GAC9D,IAAK,IAAIxD,EAAI,EAAGA,EAAI,IAAKA,EACrBD,EAAEM,oBAAoBL,GAAKD,EAAEI,sBAAsBH,GAAKe,EAAQwC,KAAKG,YAAY1D,EACzF,CAxCQ2D,CAAsB5C,EAAShB,GA6BvC,SAAuCgB,EAAShB,GAC5CA,EAAEO,mBAAqBP,EAAEE,UAAYc,EAAQwC,KAAKC,cAAc,GAChEzD,EAAEG,WAAaa,EAAQiB,MAAM4B,UAAU5E,OAAS,EAChDe,EAAEE,WAAac,EAAQQ,MAAM0B,UAAUjE,OACvCe,EAAEI,sBAAsB,IAAMY,EAAQwC,KAAKM,gBAAgB,IAAM,EACjE9D,EAAEI,sBAAsB,IAAMJ,EAAES,YAAcO,EAAQwC,KAAKM,gBAAgB,IAAM,EACrF,CAjCIC,CAA8B/C,EAAShB,EAC3C,CACA,SAASF,EAAwBkE,EAAYzE,EAAkBE,EAAkBO,EAAGU,GAGhF,GAAY,IAARV,EAAEC,EACF+C,EAA6BtC,EAAGf,MAAOJ,EAAkBE,EAAkBO,QAE3E,IAAK,MAAMgB,IAAW,CAACN,EAAGd,UAAWc,EAAGf,OACpCqD,EAA6BhC,EAASzB,EAAkBE,EAAkBO,IAItF,SAAoCgB,EAASzB,EAAkBE,EAAkBO,IAqCjF,SAA6BgB,EAASzB,EAAkBS,GACpD,IAAK,IAAIiD,EAAI,EAAGA,EAAIjC,EAAQQ,MAAM0B,UAAUjE,SAAUgE,EAAG,CACrD,MAAMd,EAAUnC,EAAEE,UAAY+C,EAAI,EAClC1D,EAAiBS,EAAEE,UAAY+C,GAAK,EAAEpH,oBAAsBsG,EAAU,IAClEnB,EAAQQ,MAAM0B,UAAUD,GAAK,IAC7Bb,EAAcpC,EAAEM,oBAAoB,GAAKU,EAAQQ,MAAMC,EAAEwB,IAAM,IAC/Db,EAAcpC,EAAEM,oBAAoB,GAAKN,EAAES,WAAaO,EAAQQ,MAAM2B,EAAEF,IACxE,IAAMjC,EAAQQ,MAAM4B,OAAOH,EACnC,CACJ,EA7CIgB,CAAoBjD,EAASzB,EAAkBS,GAC/CsD,EAActC,EAASvB,EAAkBO,GA4E7C,SAAkCkE,EAAezE,EAAkBO,GAC/D,MAAMuD,EAAUvD,EAAEG,UACZkC,EAAYrC,EAAEQ,iBACd8B,EAAa4B,EAAcV,KAAKC,cAAc,GAAKzD,EAAEE,UAC3DT,EAAiB8D,EAAU,GAAK,EAAE1H,oBAAsB0H,EAAxB,MAClBlB,EAAY,IAAMC,EAAa,IACjD,CAjFI6B,CAAyBnD,EAASvB,EAAkBO,GAEpD,MAAMuD,EAAUvD,EAAEG,UACZkC,EAAYrC,EAAEQ,iBACd8B,EAAatB,EAAQwC,KAAKC,cAAc,GAAKzD,EAAEE,UACrDT,EAAiB8D,EAAU,GAAK,EAAE1H,oBAAsB0H,EAAxB,MAClBlB,EAAY,IAAMC,EAAa,KAE7CtC,EAAEG,WAAaa,EAAQiB,MAAM4B,UAAU5E,OAAS,EAChDe,EAAEE,WAAac,EAAQQ,MAAM0B,UAAUjE,MAC3C,CAfImF,CAA2BJ,EAAYzE,EAAkBE,EAAkBO,EAC/E,CAgDA,SAASsD,EAActC,EAASvB,EAAkBO,GAE9C,IAAK,IAAIiD,EAAI,EAAGA,EAAIjC,EAAQiB,MAAM4B,UAAU5E,SAAUgE,EAAG,CACrD,MAAMM,EAAUvD,EAAEG,UAAY8C,EAAI,EAC5BZ,EAAYrB,EAAQiB,MAAM4B,UAAUZ,GAAG,GAAKjD,EAAEE,UAC9CoC,EAAatB,EAAQiB,MAAM4B,UAAUZ,GAAG,GAAKjD,EAAEE,UACrD,IAAImE,EAAU,GACd,GAAIrD,EAAQiB,MAAMqC,kBAAkBC,IAAItB,GAAI,CAExC,IAAIuB,EAAcxD,EAAQiB,MAAMqC,kBAAkBG,IAAIxB,GAClDjD,EAAES,WAAa,IACf+D,EAA+B,IAAhBA,EAAqB,EAAI,GAC5CH,EAAU,QAAUG,CACxB,CACA,MAAMpB,EAASpC,EAAQiB,MAAMmB,OAAOmB,IAAItB,GACpC,IAAMjC,EAAQiB,MAAMmB,OAAOqB,IAAIxB,GAAK,GACxCxD,EAAiBO,EAAEG,UAAY8C,GAAK,EAAEpH,oBAAsB0H,EAAU,IAClEvC,EAAQiB,MAAMyC,UAAUzB,GAAK,IAC7BZ,EAAY,IAAMC,EAAa+B,EAAUjB,EAAS,IAC1D,CACJ,CAyBA,SAASjE,EAA2BJ,EAAYC,EAActB,EAAUzB,GACpE,IAAImD,EAAY,EACZC,EAAY,EACZsF,EAAe,EACfrF,GAAe,EACfsF,EAAmB,KAEvB,IAAK,MAAM7D,KAAchC,EACrB,GL9OkB,KK8OdgC,EAAW7E,OAAf,CAEA,GAAyB,KAArB6E,EAAW7E,OACX,MAAM,IAAI2I,MAAM,sCACpBD,EAAmBrI,EAAYyC,EAAc,CAAE9C,OAAQ6E,EAAW7E,OAAQD,YAAakC,EAAsB4C,EAAWE,WACxH7B,GAAawF,EAAiBpD,MAAMC,EAAExC,OACtCI,GAAauF,EAAiB3C,MAAMyC,UAAUzF,OAC9C0F,GANY,CAShB,GAAoB,YAAhB1I,EAEAmD,GAAa,EAEbC,GAAasF,EAITA,EAAe,KACVC,GAAkBpB,MAAMsB,QAAQ7F,QAAU,GAAK,GAAK2F,GAAkB1D,gBAAgB6D,gBAAkB,EAAEzI,SAASyI,iBACpHzF,GAAe,EACfF,GAAa,EACbC,GAAa,OAIpB,CACD,MAAMM,EAASjC,IAAaL,EAAS0C,IACjCxD,EAAYyC,EAAc,EAAEhD,aAAeO,EAAYyC,EAAc,EAAE7C,QACrEyD,EAAYrD,EAAYyC,EAAc,EAAE5C,WAE9CgD,IAAcL,EAAWE,OAAS,GAAKW,EAAU4B,MAAMC,EAAExC,OAEzDG,GAAaL,EAAWE,OAASU,EAAM6B,MAAMC,EAAExC,OAE/CG,GAAa,EAEbC,IAAcN,EAAWE,OAAS,GAAKW,EAAUqC,MAAMyC,UAAUzF,OAEjEI,GAAaN,EAAWE,OAASU,EAAMsC,MAAMyC,UAAUzF,OAEvDI,GAAa,EAEbA,GAAiC,EAApBN,EAAWE,MAC5B,CACA,MAAO,CAAEG,YAAWC,YAAWC,eACnC,CAKO,SAAS8C,EAAcX,GAC1B,OAAOuD,KAAKC,MAAM,EAAElJ,iBAAmB0F,GAAK,EAAE1F,gBAClD,CLxSI,IAAIyB,EAAcH,EAAS6H,GAAItH,EAAUC,MAAMC,QAAS,IACxD,IAAIN,EAAcH,EAAS0C,IAAKnC,EAAUC,MAAMG,IAAK,KACrD,IAAIR,EAAcH,EAAS8H,IAAKvH,EAAUC,MAAMI,IAAK,KAKpDb,EAASgI,MACThI,EAASiI,UACTjI,EAASkI,KACTlI,EAASmI,KG/Dd,SAAWrH,GACPA,EAAgB,KAAI,YACpBA,EAAiB,MAAI,aACrBA,EAAkB,OAAI,cACtBA,EAAc,GAAI,UAClBA,EAAgB,KAAI,YACpBA,EAAgB,KAAI,YACpBA,EAAsB,WAAI,gBAC7B,CARD,CAQGA,IAAcA,EAAY,CAAC,IGR9BsH,UAAaC,IACT,MAAM,QAAEC,EAAO,aAAE1G,EAAY,SAAEtB,EAAQ,YAAEzB,EAAW,MAAE0J,EAAK,IAAEC,GAAQH,EAAMI,KACrEC,EAAa,IAAItG,MAAMoG,EAAMD,GAC7BI,EAAqB,IAAIvG,MAAM,GACrC,IAAK,IAAIwG,EAAOL,EAAOK,EAAOJ,IAAOI,EACjC,IACI,MAAMC,EAAMP,EAAQM,GACpBF,EAAWE,EAAOL,GAAS7G,EAAoBmH,EAAKjH,EAActB,EAAUzB,EAChF,CACA,MAAOiK,GACH,MACMC,EAAM,8BAA8BH,MAD3BE,aAAerB,MAAQqB,EAAIE,QAAUF,EAAIG,cAExDN,EAAmBxE,KAAK4E,EAC5B,CAEJG,YAAY,CAAEC,SAAUT,EAAYU,SAAUT,K","sources":["webpack://bio/./node_modules/@datagrok-libraries/bio/src/monomer-works/consts.js","webpack://bio/./node_modules/@datagrok-libraries/bio/src/monomer-works/types.js","webpack://bio/./node_modules/@datagrok-libraries/bio/src/utils/macromolecule/consts.js","webpack://bio/./node_modules/@datagrok-libraries/bio/src/utils/const.js","webpack://bio/./node_modules/@datagrok-libraries/bio/src/utils/macromolecule/types.js","webpack://bio/./node_modules/@datagrok-libraries/js-draw-lite/src/types/org.js","webpack://bio/./node_modules/@datagrok-libraries/bio/src/monomer-works/monomer-works.js","webpack://bio/./node_modules/@datagrok-libraries/bio/src/monomer-works/to-atomic-level-utils.js","webpack://bio/./node_modules/@datagrok-libraries/bio/src/monomer-works/seq-to-molfile-worker.js"],"sourcesContent":["export const monomerWorksConsts = {\n // constants for parsing molfile V2000\n V2K_RGP_SHIFT: 8,\n V2K_RGP_LINE: 'M RGP',\n V2K_A_LINE: 'A ',\n // constants for parsing/reconstruction of molfile V3000\n V3K_COUNTS_SHIFT: 14,\n V3K_IDX_SHIFT: 7,\n V3K_HEADER_FIRST_LINE: '\\nDatagrok macromolecule handler\\n\\n',\n V3K_HEADER_SECOND_LINE: ' 0 0 0 0 0 0 999 V3000\\n',\n V3K_BEGIN_CTAB_BLOCK: 'M V30 BEGIN CTAB\\n',\n V3K_END_CTAB_BLOCK: 'M V30 END CTAB\\n',\n V3K_BEGIN_COUNTS_LINE: 'M V30 COUNTS ',\n V3K_COUNTS_LINE_ENDING: ' 0 0 0\\n',\n V3K_BEGIN_ATOM_BLOCK: 'M V30 BEGIN ATOM\\n',\n V3K_END_ATOM_BLOCK: 'M V30 END ATOM\\n',\n V3K_BEGIN_BOND_BLOCK: 'M V30 BEGIN BOND\\n',\n V3K_END_BOND_BLOCK: 'M V30 END BOND\\n',\n V3K_BOND_CONFIG: ' CFG=',\n V3K_BEGIN_DATA_LINE: 'M V30 ',\n V3K_END: 'M END',\n PRECISION_FACTOR: 10000, // HELMCoreLibrary has 4 significant digits after decimal point in atom coordinates\n // symbols for the corresponding monomers in HELM library\n DEOXYRIBOSE: { polymerType: \"RNA\" /* PolymerTypes.RNA */, symbol: 'd' },\n RIBOSE: { polymerType: \"RNA\" /* PolymerTypes.RNA */, symbol: 'r' },\n PHOSPHATE: { polymerType: \"RNA\" /* PolymerTypes.RNA */, symbol: 'p' },\n OXYGEN: 'O',\n HYDROGEN: 'H',\n};\n//# sourceMappingURL=consts.js.map","export function getMolGraph(dict, libKey) {\n return dict[libKey.polymerType]?.[libKey.symbol];\n}\nexport function hasMolGraph(dict, libKey) {\n return !!dict[libKey.polymerType]?.[libKey.symbol];\n}\nexport function setMolGraph(dict, libKey, value) {\n let pt = dict[libKey.polymerType];\n if (!pt)\n pt = dict[libKey.polymerType] = {};\n pt[libKey.symbol] = value;\n}\nexport class MonomerMap extends Map {\n constructor(entries) {\n super(entries);\n }\n}\n/** @property monomers key - helm seq position, */\nexport class MolfileWithMap {\n constructor(molfile, monomers) {\n this.molfile = molfile;\n this.monomers = monomers;\n }\n static createEmpty() { return new MolfileWithMap('', new MonomerMap(null)); }\n}\n//# sourceMappingURL=types.js.map","import { CandidateType } from './types';\n/** enum type to simplify setting \"user-friendly\" notation if necessary */\nexport var NOTATION;\n(function (NOTATION) {\n NOTATION[\"FASTA\"] = \"fasta\";\n NOTATION[\"SEPARATOR\"] = \"separator\";\n NOTATION[\"HELM\"] = \"helm\";\n /* Requires notation handler */ NOTATION[\"CUSTOM\"] = \"custom\";\n /* Requires notation handler */ NOTATION[\"BILN\"] = \"biln\";\n})(NOTATION || (NOTATION = {}));\nexport var ALPHABET;\n(function (ALPHABET) {\n ALPHABET[\"DNA\"] = \"DNA\";\n ALPHABET[\"RNA\"] = \"RNA\";\n ALPHABET[\"PT\"] = \"PT\";\n /** Unknown */\n ALPHABET[\"UN\"] = \"UN\";\n})(ALPHABET || (ALPHABET = {}));\nexport var TAGS;\n(function (TAGS) {\n TAGS[\"aligned\"] = \"aligned\";\n TAGS[\"alphabet\"] = \"alphabet\";\n TAGS[\"alphabetSize\"] = \".alphabetSize\";\n TAGS[\"alphabetIsMultichar\"] = \".alphabetIsMultichar\";\n TAGS[\"separator\"] = \"separator\";\n TAGS[\"isHelmCompatible\"] = \".isHelmCompatible\";\n TAGS[\"positionNames\"] = \".positionNames\";\n TAGS[\"positionLabels\"] = \".positionLabels\";\n TAGS[\"regions\"] = \".regions\";\n TAGS[\"positionShift\"] = \".positionShift\";\n TAGS[\"selectedPosition\"] = \".selectedPosition\";\n TAGS[\"polymerTypeColumnName\"] = \".polymerTypeColumnName\";\n TAGS[\"annotations\"] = \".annotations\";\n TAGS[\"numberingScheme\"] = \".numberingScheme\";\n TAGS[\"annotationColumnName\"] = \".annotationColumnName\";\n})(TAGS || (TAGS = {}));\nexport { TAGS as BioTags };\nexport const positionSeparator = ', ';\nexport const monomerRe = /(?:\\[([A-Za-z0-9_\\-,()]+)\\])|([A-Za-z\\-])/g;\nexport const helmRe = /(PEPTIDE1|DNA1|RNA1)\\{([^}]+)}/g;\nexport const helmPp1Re = /\\[([^\\[\\]]+)]/g;\nexport const Alphabets = new class {\n constructor() {\n this.fasta = {\n peptide: new Set([\n 'G', 'L', 'Y', 'S', 'E', 'Q', 'D', 'N', 'F', 'A',\n 'K', 'R', 'H', 'C', 'V', 'P', 'W', 'I', 'M', 'T',\n ]),\n dna: new Set(['A', 'C', 'G', 'T']),\n rna: new Set(['A', 'C', 'G', 'U']),\n };\n }\n}();\nexport const candidateAlphabets = [\n new CandidateType(ALPHABET.PT, Alphabets.fasta.peptide, 0.50),\n new CandidateType(ALPHABET.DNA, Alphabets.fasta.dna, 0.55),\n new CandidateType(ALPHABET.RNA, Alphabets.fasta.rna, 0.55),\n];\n/** Canonical gap symbol */\nexport const GAP_SYMBOL = '';\nexport const GapOriginals = {\n [NOTATION.FASTA]: '-',\n [NOTATION.SEPARATOR]: '',\n [NOTATION.HELM]: '*',\n [NOTATION.BILN]: '',\n};\nexport const MONOMER_MOTIF_SPLITTER = ' , ';\nexport const NOTATION_PROVIDER_CONSTRUCTOR_ROLE = 'notationProviderConstructor';\n//# sourceMappingURL=consts.js.map","// core fields of HELM library object used in toAtomicLevel function\nexport const HELM_CORE_FIELDS = [\n \"symbol\" /* HELM_FIELDS.SYMBOL */,\n \"molfile\" /* HELM_FIELDS.MOLFILE */,\n \"rgroups\" /* HELM_FIELDS.RGROUPS */,\n \"name\" /* HELM_FIELDS.NAME */,\n // HELM_FIELDS.MONOMER_TYPE, // add if terminal monomers for PEPTIDEs to be\n // supported\n];\nexport const SDF_MONOMER_NAME = 'MonomerName';\n// todo: ideally, keys should be expressed via constants\nexport const jsonSdfMonomerLibDict = {\n 'monomerType': null, // -> Backbone\n 'smiles': null,\n 'name': 'Name',\n 'author': null,\n 'molfile': 'molecule',\n 'naturalAnalog': 'MonomerNaturalAnalogCode',\n 'rgroups': 'MonomerCaps',\n 'createDate': null,\n 'id': null,\n 'polymerType': 'MonomerType',\n 'symbol': 'MonomerName'\n};\nexport const DUMMY_MONOMER = {\n 'monomerType': 'Backbone',\n 'smiles': '',\n 'name': '',\n 'author': 'Datagrok',\n 'molfile': '',\n 'naturalAnalog': '',\n 'rgroups': [],\n 'createDate': null,\n 'id': 0,\n 'polymerType': 'PEPTIDE',\n 'symbol': ''\n};\n// range of hex nubers used in PepSea library to endode monomers\nexport const MONOMER_ENCODE_MIN = 0x100;\nexport const MONOMER_ENCODE_MAX = 0x40A;\nexport const RIBOSE_SYMBOL = 'r';\nexport const DEOXYRIBOSE_SYMBOL = 'd';\nexport const PHOSPHATE_SYMBOL = 'p';\nexport const HELM_WRAPPERS_REGEXP = new RegExp(`[${RIBOSE_SYMBOL}${DEOXYRIBOSE_SYMBOL}]\\\\((\\\\w)\\\\)${PHOSPHATE_SYMBOL}?`, 'g');\n//# sourceMappingURL=const.js.map","import { NOTATION_PROVIDER_CONSTRUCTOR_ROLE } from './consts';\nexport class NotationProviderBase {\n /** Name of the custom notation */\n static get notationName() {\n return 'Custom';\n }\n ;\n /** flag to let bio know if this provider implements method for converting helm to it */\n static get implementsFromHelm() {\n return false;\n }\n ;\n /** Method for converting HELM to this notation */\n static convertFromHelm(helm, options) {\n throw new Error(`Method convertFromHelm not implemented for this notation provider`);\n }\n ;\n static async getProviderConstructors() {\n // this is terrible, I know, but otherwise this gets put in webworkers and fails due to DG resolution)))\n // @ts-ignore\n if (window?.DG) {\n // @ts-ignore\n const constFuncs = window.DG.Func.find({ meta: { role: NOTATION_PROVIDER_CONSTRUCTOR_ROLE } });\n return Promise.all(constFuncs.map((f) => f.apply({})));\n }\n return [];\n }\n}\n/** Alphabet candidate type */\nexport class CandidateType {\n constructor(name, alphabet, cutoff) {\n this.name = name;\n this.alphabet = alphabet;\n this.cutoff = cutoff;\n }\n}\n/** Alphabet candidate similarity type */\nexport class CandidateSimType extends CandidateType {\n constructor(candidate, freq, similarity) {\n super(candidate.name, candidate.alphabet, candidate.cutoff);\n this.freq = freq;\n this.similarity = similarity;\n }\n}\n//# sourceMappingURL=types.js.map","export var HelmTypes;\n(function (HelmTypes) {\n HelmTypes[\"BASE\"] = \"HELM_BASE\";\n HelmTypes[\"SUGAR\"] = \"HELM_SUGAR\";\n HelmTypes[\"LINKER\"] = \"HELM_LINKER\";\n HelmTypes[\"AA\"] = \"HELM_AA\";\n HelmTypes[\"CHEM\"] = \"HELM_CHEM\";\n HelmTypes[\"BLOB\"] = \"HELM_BLOB\";\n HelmTypes[\"NUCLEOTIDE\"] = \"HELM_NUCLETIDE\";\n})(HelmTypes || (HelmTypes = {}));\n//# sourceMappingURL=org.js.map","import { HelmTypes } from '../helm/consts';\nexport class MonomerWorks {\n //private molfileConverter: IMolfileConverter;\n constructor(monomerLib /*, molfileConverter: IMolfileConverter*/) {\n this.monomerLib = monomerLib;\n //this.molfileConverter = molfileConverter;\n }\n //types according to Monomer possible\n getCappedRotatedMonomer(monomerType, monomerName) {\n // TODO: Check type of monomerType arg\n const monomer = this.monomerLib.getMonomer(monomerType, monomerName);\n if (monomer)\n return monomer.molfile; //TODO cap\n return null;\n }\n}\nexport function helmTypeToPolymerType(helmType) {\n let polymerType = undefined;\n switch (helmType) {\n case HelmTypes.BASE:\n case HelmTypes.SUGAR: // r - ribose, d - deoxyribose\n case HelmTypes.LINKER: // p - phosphate\n case HelmTypes.NUCLEOTIDE:\n // @ts-ignore\n case 'nucleotide':\n polymerType = \"RNA\" /* PolymerTypes.RNA */;\n break;\n case HelmTypes.AA:\n polymerType = \"PEPTIDE\" /* PolymerTypes.PEPTIDE */;\n break;\n case HelmTypes.CHEM:\n polymerType = \"CHEM\" /* PolymerTypes.CHEM */;\n break;\n case HelmTypes.BLOB:\n polymerType = \"BLOB\" /* PolymerTypes.BLOB */;\n break;\n default:\n polymerType = \"PEPTIDE\" /* PolymerTypes.PEPTIDE */;\n console.warn(`Unexpected HelmType '${helmType}'`);\n }\n return polymerType;\n}\n//# sourceMappingURL=monomer-works.js.map","/* eslint-disable max-len */\nimport { monomerWorksConsts as C } from './consts';\nimport { getMolGraph, MolfileWithMap, MonomerMap } from './types';\nimport { HELM_CORE_FIELDS, } from '../utils/const';\nimport { ALPHABET, GAP_SYMBOL } from '../utils/macromolecule/consts';\nimport { helmTypeToPolymerType } from './monomer-works';\n/** Get a mapping of peptide symbols to HELM monomer library objects with selected fields.\n * @param {IMonomerLib} monomerLib - Monomer library\n * @param {HELM_POLYMER_TYPE} polymerType - Polymer type\n * @param {ALPHABET} alphabet - Alphabet of the column\n * @return {Map<string, any>} - Mapping of peptide symbols to HELM monomer library objects with selected fields*/\nexport function getFormattedMonomerLib(monomerLib, polymerType, alphabet) {\n const map = new Map();\n for (const monomerSymbol of monomerLib.getMonomerSymbolsByType(polymerType)) {\n const it = monomerLib.getMonomer(polymerType, monomerSymbol);\n if (polymerType === \"RNA\" /* HELM_POLYMER_TYPE.RNA */ &&\n (it[\"monomerType\" /* HELM_FIELDS.MONOMER_TYPE */] === \"Branch\" /* HELM_MONOMER_TYPE.BRANCH */ ||\n alphabet === ALPHABET.DNA && it[\"symbol\" /* HELM_FIELDS.SYMBOL */] === C.DEOXYRIBOSE.symbol ||\n alphabet === ALPHABET.RNA && it[\"symbol\" /* HELM_FIELDS.SYMBOL */] === C.RIBOSE.symbol ||\n it[\"symbol\" /* HELM_FIELDS.SYMBOL */] === C.PHOSPHATE.symbol) ||\n polymerType === \"PEPTIDE\" /* HELM_POLYMER_TYPE.PEPTIDE */ &&\n it[\"monomerType\" /* HELM_FIELDS.MONOMER_TYPE */] !== \"Branch\" /* HELM_MONOMER_TYPE.BRANCH */) {\n const monomerObject = {};\n HELM_CORE_FIELDS.forEach((field) => {\n //@ts-ignore\n monomerObject[field] = it[field];\n });\n map.set(monomerSymbol, monomerObject);\n }\n }\n return map;\n}\n/** Translate a sequence of monomer symbols into Molfile V3000\n * @param {string[]} monomerSeq - Sequence of monomer symbols (canonical)\n * @param {Map<string, MolGraph>} monomersDict - Mapping of monomer symbols to MolGraph objects\n * @param {ALPHABET} alphabet - Alphabet of the column\n * @param {PolymerType} polymerType - Polymer type\n * @return {string} - Molfile V3000*/\nexport function monomerSeqToMolfile(monomerSeq, monomersDict, alphabet, polymerType) {\n if (monomerSeq.length === 0) {\n // throw new Error('monomerSeq is empty');\n return MolfileWithMap.createEmpty();\n }\n // define atom and bond counts, taking into account the bond type\n const getAtomAndBondCounts = getResultingAtomBondCounts;\n const { atomCount, bondCount, needsCapping } = getAtomAndBondCounts(monomerSeq, monomersDict, alphabet, polymerType);\n // create arrays to store lines of the resulting molfile\n const molfileAtomBlock = new Array(atomCount);\n const molfileBondBlock = new Array(bondCount);\n let addMonomerToMolblock;\n let sugar = null;\n let phosphate = null;\n if (polymerType === \"PEPTIDE\" /* HELM_POLYMER_TYPE.PEPTIDE */)\n addMonomerToMolblock = addAminoAcidToMolblock;\n else { // nucleotides\n addMonomerToMolblock = addNucleotideToMolblock;\n sugar = (alphabet === ALPHABET.DNA) ? getMolGraph(monomersDict, C.DEOXYRIBOSE) : getMolGraph(monomersDict, C.RIBOSE);\n phosphate = getMolGraph(monomersDict, C.PHOSPHATE);\n }\n const v = {\n i: 0,\n nodeShift: 0,\n bondShift: 0,\n backbonePositionShift: new Array(2).fill(0),\n branchPositionShift: new Array(2).fill(0),\n backboneAttachNode: 0,\n branchAttachNode: 0,\n flipFactor: 1,\n };\n const LC = {\n sugar: sugar,\n phosphate: phosphate,\n seqLength: monomerSeq.length,\n atomCount: atomCount,\n bondCount: bondCount,\n };\n const monomers = new MonomerMap();\n const steabsCollection = [];\n let nAtoms = 0;\n let lastMonomerCappingAtom = undefined;\n for (v.i = 0; v.i < LC.seqLength; ++v.i) {\n const seqMonomer = monomerSeq[v.i];\n if (seqMonomer.symbol === GAP_SYMBOL)\n continue;\n const monomer = getMolGraph(monomersDict, { symbol: seqMonomer.symbol, polymerType: helmTypeToPolymerType(seqMonomer.biotype) });\n lastMonomerCappingAtom = monomer.terminalR2Atom;\n const mAtomFirst = v.nodeShift;\n const mBondFirst = v.bondShift;\n addMonomerToMolblock(monomer, molfileAtomBlock, molfileBondBlock, v, LC);\n //adding stereo atoms to array for further STEABS block generation\n monomer.stereoAtoms?.forEach((i) => steabsCollection.push(i + nAtoms));\n nAtoms += monomer.atoms.x.length;\n const mAtomCount = v.nodeShift - mAtomFirst;\n const mAtomList = new Array(mAtomCount);\n for (let maI = 0; maI < mAtomCount; ++maI)\n mAtomList[maI] = mAtomFirst + maI;\n const mBondCount = v.bondShift - mBondFirst;\n const mBondList = new Array(mBondCount);\n for (let mbI = 0; mbI < mBondCount; ++mbI)\n mBondList[mbI] = mBondFirst + mbI;\n monomers.set(v.i, {\n biotype: seqMonomer.biotype,\n symbol: seqMonomer.symbol,\n atoms: mAtomList, bonds: mBondList\n });\n }\n // if the last monomer needs to be capped, add the terminal OH to the resulting molfile\n if (needsCapping)\n capResultingMolblock(molfileAtomBlock, molfileBondBlock, v, LC, lastMonomerCappingAtom ?? C.OXYGEN);\n const molfileCountsLine = C.V3K_BEGIN_COUNTS_LINE + atomCount + ' ' + bondCount + C.V3K_COUNTS_LINE_ENDING;\n // todo: possible optimization may be achieved by replacing .join('') with +=\n // since counterintuitively joining an array into a new string is reportedly\n // slower than using += as below\n let result = '';\n result += C.V3K_HEADER_FIRST_LINE;\n result += C.V3K_HEADER_SECOND_LINE;\n result += C.V3K_BEGIN_CTAB_BLOCK;\n result += molfileCountsLine;\n result += C.V3K_BEGIN_ATOM_BLOCK;\n result += molfileAtomBlock.join('');\n result += C.V3K_END_ATOM_BLOCK;\n result += C.V3K_BEGIN_BOND_BLOCK;\n result += molfileBondBlock.join('');\n result += C.V3K_END_BOND_BLOCK;\n if (steabsCollection.length > 0)\n result += getCollectionBlock(steabsCollection);\n result += C.V3K_END_CTAB_BLOCK;\n result += C.V3K_END;\n // return molfileParts.join('');\n return { molfile: result, monomers: monomers };\n}\nfunction getCollectionBlock(collection) {\n //one row in STEABS block can be no longer than 80 symbols\n //maxSymbols = 80 symbols minus ' -\\n' (4 symbols)\n const maxSymbols = 76;\n const rowsArray = [];\n let newCollectionRow = `M V30 MDLV30/STEABS ATOMS=(${collection.length}`;\n for (let i = 0; i < collection.length; i++) {\n const updatedRow = `${newCollectionRow} ${collection[i]}`;\n if (updatedRow.length > maxSymbols) {\n rowsArray.push(`${newCollectionRow} -\\n`);\n newCollectionRow = `M V30 ${collection[i]}`;\n }\n else\n newCollectionRow = updatedRow;\n //in case last atom was added - close the block\n if (i === collection.length - 1)\n rowsArray.push(`${newCollectionRow})\\n`);\n }\n return `M V30 BEGIN COLLECTION\\n${rowsArray.join('')}M V30 END COLLECTION\\n`;\n}\n/** Cap the resulting (after sewing up all the monomers) molfile with 'O'\n * @param {string[]} molfileAtomBlock - Array of lines of the resulting molfile atom block\n * @param {string[]} molfileBondBlock - Array of lines of the resulting molfile bond block\n * @param {LoopVariables} v - Loop variables\n * @param {LoopConstants} LC - Loop constants*/\nfunction capResultingMolblock(molfileAtomBlock, molfileBondBlock, v, LC, cappingAtomType = C.OXYGEN) {\n // add terminal oxygen\n const atomIdx = v.nodeShift + 1;\n molfileAtomBlock[LC.atomCount] = C.V3K_BEGIN_DATA_LINE + atomIdx + ' ' +\n (cappingAtomType ?? C.OXYGEN) + ' ' + keepPrecision(v.backbonePositionShift[0]) + ' ' +\n v.flipFactor * keepPrecision(v.backbonePositionShift[1]) + ' ' + '0.000000 0' + '\\n';\n // add terminal bond\n const firstAtom = v.backboneAttachNode;\n const secondAtom = atomIdx;\n molfileBondBlock[LC.bondCount] = C.V3K_BEGIN_DATA_LINE + v.bondShift + ' ' +\n 1 + ' ' + firstAtom + ' ' + secondAtom + '\\n';\n}\nfunction addAminoAcidToMolblock(monomer, molfileAtomBlock, molfileBondBlock, v) {\n v.flipFactor = (-1) ** (v.i % 2); // to flip every even monomer over OX\n addBackboneMonomerToMolblock(monomer, molfileAtomBlock, molfileBondBlock, v);\n}\nfunction addBackboneMonomerToMolblock(monomer, molfileAtomBlock, molfileBondBlock, v) {\n // todo: remove these comments to the docstrings of the corr. functions\n // construnct the lines of V3K molfile atom block\n fillAtomLines(monomer, molfileAtomBlock, v);\n // construct the lines of V3K molfile bond block\n fillBondLines(monomer, molfileBondBlock, v);\n // peptide bond\n fillChainExtendingBond(monomer, molfileBondBlock, v);\n // update branch variables if necessary\n if (monomer.meta.branchShift !== null && monomer.meta.terminalNodes.length > 2)\n updateBranchVariables(monomer, v);\n // update loop variables\n updateChainExtendingVariables(monomer, v);\n}\nfunction addNucleotideToMolblock(nucleobase, molfileAtomBlock, molfileBondBlock, v, LC) {\n // construnct the lines of V3K molfile atom block corresponding to phosphate\n // and sugar\n if (v.i === 0)\n addBackboneMonomerToMolblock(LC.sugar, molfileAtomBlock, molfileBondBlock, v);\n else {\n for (const monomer of [LC.phosphate, LC.sugar])\n addBackboneMonomerToMolblock(monomer, molfileAtomBlock, molfileBondBlock, v);\n }\n addBranchMonomerToMolblock(nucleobase, molfileAtomBlock, molfileBondBlock, v);\n}\nfunction addBranchMonomerToMolblock(monomer, molfileAtomBlock, molfileBondBlock, v) {\n fillBranchAtomLines(monomer, molfileAtomBlock, v);\n fillBondLines(monomer, molfileBondBlock, v);\n fillBackboneToBranchBond(monomer, molfileBondBlock, v);\n // C-N bond\n const bondIdx = v.bondShift;\n const firstAtom = v.branchAttachNode;\n const secondAtom = monomer.meta.terminalNodes[0] + v.nodeShift;\n molfileBondBlock[bondIdx - 1] = C.V3K_BEGIN_DATA_LINE + bondIdx + ' ' +\n 1 + ' ' + firstAtom + ' ' + secondAtom + '\\n';\n // update loop variables\n v.bondShift += monomer.bonds.atomPairs.length + 1;\n v.nodeShift += monomer.atoms.atomTypes.length;\n}\nfunction updateChainExtendingVariables(monomer, v) {\n v.backboneAttachNode = v.nodeShift + monomer.meta.terminalNodes[1];\n v.bondShift += monomer.bonds.atomPairs.length + 1;\n v.nodeShift += monomer.atoms.atomTypes.length;\n v.backbonePositionShift[0] += monomer.meta.backboneShift?.[0] ?? 0; // todo: non-null check\n v.backbonePositionShift[1] += v.flipFactor * (monomer.meta.backboneShift?.[1] ?? 0);\n}\nfunction updateBranchVariables(monomer, v) {\n v.branchAttachNode = v.nodeShift + monomer.meta.terminalNodes[2];\n for (let i = 0; i < 2; ++i)\n v.branchPositionShift[i] = v.backbonePositionShift[i] + monomer.meta.branchShift[i];\n}\nfunction fillAtomLines(monomer, molfileAtomBlock, v) {\n for (let j = 0; j < monomer.atoms.atomTypes.length; ++j) {\n const atomIdx = v.nodeShift + j + 1;\n molfileAtomBlock[v.nodeShift + j] = C.V3K_BEGIN_DATA_LINE + atomIdx + ' ' +\n monomer.atoms.atomTypes[j] + ' ' +\n keepPrecision(v.backbonePositionShift[0] + monomer.atoms.x[j]) + ' ' +\n keepPrecision(v.backbonePositionShift[1] + v.flipFactor * monomer.atoms.y[j]) +\n ' ' + monomer.atoms.kwargs[j];\n }\n}\n// todo: remove as quickfix\nfunction fillBranchAtomLines(monomer, molfileAtomBlock, v) {\n for (let j = 0; j < monomer.atoms.atomTypes.length; ++j) {\n const atomIdx = v.nodeShift + j + 1;\n molfileAtomBlock[v.nodeShift + j] = C.V3K_BEGIN_DATA_LINE + atomIdx + ' ' +\n monomer.atoms.atomTypes[j] + ' ' +\n keepPrecision(v.branchPositionShift[0] + monomer.atoms.x[j]) + ' ' +\n keepPrecision(v.branchPositionShift[1] + v.flipFactor * monomer.atoms.y[j]) +\n ' ' + monomer.atoms.kwargs[j];\n }\n}\nfunction fillBondLines(monomer, molfileBondBlock, v) {\n // construct the lines of V3K molfile bond block\n for (let j = 0; j < monomer.bonds.atomPairs.length; ++j) {\n const bondIdx = v.bondShift + j + 1;\n const firstAtom = monomer.bonds.atomPairs[j][0] + v.nodeShift;\n const secondAtom = monomer.bonds.atomPairs[j][1] + v.nodeShift;\n let bondCfg = '';\n if (monomer.bonds.bondConfiguration.has(j)) {\n // flip orientation when necessary\n let orientation = monomer.bonds.bondConfiguration.get(j);\n if (v.flipFactor < 0)\n orientation = (orientation === 1) ? 3 : 1;\n bondCfg = ' CFG=' + orientation;\n }\n const kwargs = monomer.bonds.kwargs.has(j) ?\n ' ' + monomer.bonds.kwargs.get(j) : '';\n molfileBondBlock[v.bondShift + j] = C.V3K_BEGIN_DATA_LINE + bondIdx + ' ' +\n monomer.bonds.bondTypes[j] + ' ' +\n firstAtom + ' ' + secondAtom + bondCfg + kwargs + '\\n';\n }\n}\nfunction fillChainExtendingBond(monomer, molfileBondBlock, v) {\n if (v.backboneAttachNode !== 0) {\n const bondIdx = v.bondShift;\n const firstAtom = v.backboneAttachNode;\n const secondAtom = monomer.meta.terminalNodes[0] + v.nodeShift;\n molfileBondBlock[v.bondShift - 1] = C.V3K_BEGIN_DATA_LINE + bondIdx + ' ' +\n 1 + ' ' + firstAtom + ' ' + secondAtom + '\\n';\n }\n}\n// todo: remove\nfunction fillBackboneToBranchBond(branchMonomer, molfileBondBlock, v) {\n const bondIdx = v.bondShift;\n const firstAtom = v.branchAttachNode;\n const secondAtom = branchMonomer.meta.terminalNodes[0] + v.nodeShift;\n molfileBondBlock[bondIdx - 1] = C.V3K_BEGIN_DATA_LINE + bondIdx + ' ' +\n 1 + ' ' + firstAtom + ' ' + secondAtom + '\\n';\n}\n/** Compute the atom/bond counts for the resulting molfile, depending on the\n * type of polymer (peptide/nucleotide)\n * @param {string[]}monomerSeq - the sequence of monomers\n * @param {Map<string, MolGraph>}monomersDict - the dictionary of monomers\n * @param {ALPHABET}alphabet - the alphabet of the monomers\n * @param {HELM_POLYMER_TYPE}polymerType - the type of polymer\n * @return {{atomCount: number, bondCount: number}} - the atom/bond counts*/\nfunction getResultingAtomBondCounts(monomerSeq, monomersDict, alphabet, polymerType) {\n let atomCount = 0;\n let bondCount = 0;\n let monomerCount = 0;\n let needsCapping = true;\n let lastMonomerGraph = null;\n // sum up all the atoms/nodes provided by the sequence\n for (const seqMonomer of monomerSeq) {\n if (seqMonomer.symbol === GAP_SYMBOL)\n continue; // Skip for gap/empty monomer in MSA\n if (seqMonomer.symbol == '*')\n throw new Error(`Gap canonical symbol is '', not '*`);\n lastMonomerGraph = getMolGraph(monomersDict, { symbol: seqMonomer.symbol, polymerType: helmTypeToPolymerType(seqMonomer.biotype) });\n atomCount += lastMonomerGraph.atoms.x.length;\n bondCount += lastMonomerGraph.bonds.bondTypes.length;\n monomerCount++;\n }\n // add extra values depending on the polymer type\n if (polymerType === \"PEPTIDE\" /* HELM_POLYMER_TYPE.PEPTIDE */) {\n // add the rightmost/terminating cap group 'OH' (i.e. 'O')\n atomCount += 1;\n // add chain-extending bonds (C-NH per each monomer pair and terminal C-OH)\n bondCount += monomerCount;\n // if the last monomer is something like NH2, which only has R1, there is no need to cap it\n // although, this should never happen, but hey... in other bits of code, there is a chunk that adds pseudo-R2 as hydrogen\n // we should also check, if the R2 of the last monomer is not hydrogen, that case should also be omitted\n if (monomerCount > 0) {\n if ((lastMonomerGraph?.meta?.rNodes?.length ?? 0) < 2 || lastMonomerGraph?.terminalR2Atom?.toLowerCase() === C.HYDROGEN.toLowerCase()) {\n needsCapping = false;\n atomCount -= 1; // remove the last atom (the terminal 'O')\n bondCount -= 1; // remove the last bond (the terminal C-OH)\n }\n }\n }\n else { // nucleotides\n const sugar = (alphabet === ALPHABET.DNA) ?\n getMolGraph(monomersDict, C.DEOXYRIBOSE) : getMolGraph(monomersDict, C.RIBOSE);\n const phosphate = getMolGraph(monomersDict, C.PHOSPHATE);\n // add phosphate per each pair of nucleobase symbols\n atomCount += (monomerSeq.length - 1) * phosphate.atoms.x.length;\n // add sugar per each nucleobase symbol\n atomCount += monomerSeq.length * sugar.atoms.x.length;\n // add the leftmost cap group 'OH' (i.e. 'O')\n atomCount += 1;\n // add bonds from phosphate monomers\n bondCount += (monomerSeq.length - 1) * phosphate.bonds.bondTypes.length;\n // add bonds from sugar monomers\n bondCount += monomerSeq.length * sugar.bonds.bondTypes.length;\n // exclude the first chain-extending bond O-P (absent, no 'leftmost' phosphate)\n bondCount -= 1;\n // add chain-extending and branch bonds (O-P, C-O and C-N per each nucleotide)\n bondCount += monomerSeq.length * 3;\n }\n return { atomCount, bondCount, needsCapping };\n}\n/** Keep precision upon floating point operations over atom coordinates\n * @param {number}x - the floating point number\n * @return {number} - the floating point number with the same precision\n */\nexport function keepPrecision(x) {\n return Math.round(C.PRECISION_FACTOR * x) / C.PRECISION_FACTOR;\n}\n//# sourceMappingURL=to-atomic-level-utils.js.map","import { monomerSeqToMolfile } from './to-atomic-level-utils';\nonmessage = (event) => {\n const { seqList, monomersDict, alphabet, polymerType, start, end } = event.data;\n const resMolList = new Array(end - start);\n const molfileWarningList = new Array(0);\n for (let rowI = start; rowI < end; ++rowI) {\n try {\n const seq = seqList[rowI];\n resMolList[rowI - start] = monomerSeqToMolfile(seq, monomersDict, alphabet, polymerType);\n }\n catch (err) {\n const errMsg = err instanceof Error ? err.message : err.toString();\n const msg = `Cannot get molfile of row #${rowI}: ${errMsg}.`;\n molfileWarningList.push(msg);\n }\n }\n postMessage({ molfiles: resMolList, warnings: molfileWarningList });\n};\n//# sourceMappingURL=seq-to-molfile-worker.js.map"],"names":["monomerWorksConsts","V2K_RGP_SHIFT","V2K_RGP_LINE","V2K_A_LINE","V3K_COUNTS_SHIFT","V3K_IDX_SHIFT","V3K_HEADER_FIRST_LINE","V3K_HEADER_SECOND_LINE","V3K_BEGIN_CTAB_BLOCK","V3K_END_CTAB_BLOCK","V3K_BEGIN_COUNTS_LINE","V3K_COUNTS_LINE_ENDING","V3K_BEGIN_ATOM_BLOCK","V3K_END_ATOM_BLOCK","V3K_BEGIN_BOND_BLOCK","V3K_END_BOND_BLOCK","V3K_BOND_CONFIG","V3K_BEGIN_DATA_LINE","V3K_END","PRECISION_FACTOR","DEOXYRIBOSE","polymerType","symbol","RIBOSE","PHOSPHATE","OXYGEN","HYDROGEN","getMolGraph","dict","libKey","MonomerMap","Map","constructor","entries","super","MolfileWithMap","molfile","monomers","this","createEmpty","NOTATION","ALPHABET","TAGS","RegExp","CandidateType","name","alphabet","cutoff","Alphabets","fasta","peptide","Set","dna","rna","HelmTypes","helmTypeToPolymerType","helmType","BASE","SUGAR","LINKER","NUCLEOTIDE","AA","CHEM","BLOB","console","warn","monomerSeqToMolfile","monomerSeq","monomersDict","length","getAtomAndBondCounts","getResultingAtomBondCounts","atomCount","bondCount","needsCapping","molfileAtomBlock","Array","molfileBondBlock","addMonomerToMolblock","sugar","phosphate","addAminoAcidToMolblock","addNucleotideToMolblock","DNA","v","i","nodeShift","bondShift","backbonePositionShift","fill","branchPositionShift","backboneAttachNode","branchAttachNode","flipFactor","LC","seqLength","steabsCollection","lastMonomerCappingAtom","nAtoms","seqMonomer","monomer","biotype","terminalR2Atom","mAtomFirst","mBondFirst","stereoAtoms","forEach","push","atoms","x","mAtomCount","mAtomList","maI","mBondCount","mBondList","mbI","set","bonds","cappingAtomType","atomIdx","keepPrecision","firstAtom","secondAtom","capResultingMolblock","molfileCountsLine","result","join","collection","rowsArray","newCollectionRow","updatedRow","getCollectionBlock","addBackboneMonomerToMolblock","j","atomTypes","y","kwargs","fillAtomLines","fillBondLines","bondIdx","meta","terminalNodes","fillChainExtendingBond","branchShift","updateBranchVariables","atomPairs","backboneShift","updateChainExtendingVariables","nucleobase","fillBranchAtomLines","branchMonomer","fillBackboneToBranchBond","addBranchMonomerToMolblock","bondCfg","bondConfiguration","has","orientation","get","bondTypes","monomerCount","lastMonomerGraph","Error","rNodes","toLowerCase","Math","round","PT","RNA","FASTA","SEPARATOR","HELM","BILN","onmessage","event","seqList","start","end","data","resMolList","molfileWarningList","rowI","seq","err","msg","message","toString","postMessage","molfiles","warnings"],"sourceRoot":""}