icn3d 3.26.7 → 3.26.8

package/icn3d.module.js

@@ -6932,6 +6932,21 @@ class ClickMenu {
         thisClass.setLogCmd("export pdb hydrogen", true);
        });
 
+       me.myEventCls.onIds("#" + me.pre + "m1_exportIgstrand", "click", async function(e) { let ic = me.icn3d; //e.preventDefault();
+       ic.refnumCls.exportRefnum('igstrand');
+       thisClass.setLogCmd("export refnum igstrand", true);
+      });
+
+      me.myEventCls.onIds("#" + me.pre + "m1_exportKabat", "click", async function(e) { let ic = me.icn3d; //e.preventDefault();
+         ic.refnumCls.exportRefnum('kabat');
+         thisClass.setLogCmd("export refnum kabat", true);
+      });
+
+      me.myEventCls.onIds("#" + me.pre + "m1_exportImgt", "click", async function(e) { let ic = me.icn3d; //e.preventDefault();
+      ic.refnumCls.exportRefnum('imgt');
+      thisClass.setLogCmd("export refnum imgt", true);
+      });
+
         me.myEventCls.onIds("#" + me.pre + "mn1_exportStl", "click", function(e) { let ic = me.icn3d; //e.preventDefault();
            thisClass.setLogCmd("export stl file", false);
            //ic.threeDPrintCls.hideStabilizer();
@@ -9504,6 +9519,16 @@ class SetMenu {
             html += this.getLink('mn1_exportSecondary', 'Secondary Structure', undefined, 2);
         }
 
+        //!!!
+        /*
+        html += this.getMenuText('m1_exportrefnum', 'Reference Numbers', undefined, undefined, 2);
+        html += "<ul>";
+        html += this.getLink('mn1_exportIgstrand', 'Ig Strand', undefined, 3);
+        html += this.getLink('mn1_exportKabat', 'Kabat', undefined, 3);
+        html += this.getLink('mn1_exportImgt', 'IMGT', undefined, 3);
+        html += "</ul>";
+        */
+
         html += "<li><br/></li>";
 
         html += "</ul>";
@@ -13760,7 +13785,7 @@ class Events {
             if(esmfold_fasta.indexOf('>') != -1) { //FASTA with header
                 let pos = esmfold_fasta.indexOf('\n');
                 ic.esmTitle = esmfold_fasta.substr(1, pos - 1).trim();
-                pdbid = ic.esmTitle.substr(0, ic.esmTitle.indexOf(' '));
+                pdbid = (ic.esmTitle.indexOf(' ') != -1) ? ic.esmTitle.substr(0, ic.esmTitle.indexOf(' ')) : ic.esmTitle;
                 if(pdbid.length < 6) pdbid = pdbid.padEnd(6, '-');
 
                 esmfold_fasta = esmfold_fasta.substr(pos + 1);
@@ -42993,6 +43018,7 @@ class ShowSeq {
                 break;
             }
         }
+
         if(kabat_or_imgt == 1 && !bKabatFound) {
             return {html: '', html3: ''};
         }
@@ -45336,9 +45362,14 @@ class LineGraph {
                     domainAtomsArray.push(currAtoms);
 
                     let residueArray = ic.resid2specCls.atoms2residues(Object.keys(currAtoms));
+
+                    let atomFirst = ic.firstAtomObjCls.getFirstAtomObj(currAtoms);
+                    let atomLast = ic.firstAtomObjCls.getLastAtomObj(currAtoms);
+                    let resiSum = parseInt(atomFirst.resi) + parseInt(atomLast.resi);
+
                     for(let n = 0, nl = residueArray.length; n < nl; ++n) {
                         let resid = residueArray[n];
-                        ic.resid2domainid[resid] = chainid + '-0' + '_' + Object.keys(currAtoms).length; 
+                        ic.resid2domainid[resid] = chainid + '-0' + '_' + resiSum; 
                     }
                 }
                 else {                 
@@ -45358,13 +45389,17 @@ class LineGraph {
 
                         domainAtomsArray.push(domainAtoms);
 
+                        let atomFirst = ic.firstAtomObjCls.getFirstAtomObj(domainAtoms);
+                        let atomLast = ic.firstAtomObjCls.getLastAtomObj(domainAtoms);
+                        let resiSum = parseInt(atomFirst.resi) + parseInt(atomLast.resi);
+
                         for(let m = 0, ml = segArray.length; m < ml; m += 2) {
                             let startResi = segArray[m];
                             let endResi = segArray[m+1];
                             for(let n = parseInt(startResi); n <= parseInt(endResi); ++n) {
                                 let resid = chainid + '_' + pos2resi[n];
                                 //domainAtoms = me.hashUtilsCls.unionHash(domainAtoms, ic.residues[resid]);
-                                ic.resid2domainid[resid] = chainid + '-' + k + '_' + Object.keys(domainAtoms).length; 
+                                ic.resid2domainid[resid] = chainid + '-' + k + '_' + resiSum; 
                             }
                         }
                     }
@@ -45404,7 +45439,7 @@ class LineGraph {
             // dataArray2 = await allPromise;
 
             //split arrays into chunks of 100 jobs
-            let n = (me.bNode) ? 16 : 96;
+            let n = (me.bNode) ? 96 : 96;
 
             for(let i = 0, il = parseInt((ajaxArray.length - 1) / n + 1); i < il; ++i) {
                 let currAjaxArray = [];
@@ -45627,6 +45662,8 @@ class LineGraph {
         if(!ic.resid2refnum) ic.resid2refnum = {};
         if(!ic.refnum2residArray) ic.refnum2residArray = {};
         if(!ic.chainsMapping) ic.chainsMapping = {};
+
+        if(!ic.refPdbList) ic.refPdbList = [];
         for(let chainid in chainid2segs) {
             let segArray = chainid2segs[chainid];
 
@@ -45635,7 +45672,8 @@ class LineGraph {
             let chainList = this.getTemplateList(chainid);
 
             //if(!me.bNode) console.log("The reference PDB(s) for chain " + chainid + " are " + chainList);
-            console.log("The reference PDB(s) for chain " + chainid + " are " + chainList);
+            if(!me.bNode) console.log("The reference PDB(s) for chain " + chainid + " are " + chainList);
+            ic.refPdbList.push("The reference PDB(s) for chain " + chainid + " are " + chainList);
 
             let prevStrand;
             let bCd19 = refpdbnameArray.length == 1 && refpdbnameArray[0] == 'CD19_6al5A_human_C2orV-n1';
@@ -45818,9 +45856,121 @@ class LineGraph {
         ic.annotationCls.setAnnoViewAndDisplay('detailed view');
     }
 
-    rmStrandFromRefnumlabel(refnumLabel) {
+    rmStrandFromRefnumlabel(refnumLabel) { let ic = this.icn3d; ic.icn3dui;
         return (!refnumLabel) ? refnumLabel : refnumLabel.replace(/'/g, '').replace(/\*/g, '').replace(/\^/g, '').replace(/\+/g, '').replace(/\-/g, '').substr(1); // C', C''
     }
+
+    exportRefnum(type) { let ic = this.icn3d, me = ic.icn3dui;
+        let refData = '';
+
+        // 1. show IgStrand ref numbers
+        if(type == 'igstrand' || type == 'IgStrand') {
+            // iGStrand reference numbers were adjusted when showing in sequences
+            if(me.bNode) {
+                for(let chnid in ic.chains) {
+                    let atom = ic.firstAtomObjCls.getFirstAtomObj(ic.chains[chnid]);
+                    if(ic.proteins.hasOwnProperty(atom.serial)) {
+                        let giSeq = [];
+                        for(let i = 0; i < ic.chainsSeq[chnid].length; ++i) {
+                            giSeq.push(ic.chainsSeq[chnid][i].name);
+                        }
+                        ic.showSeqCls.showRefNum(giSeq, chnid);
+                    }
+                }
+            }
+    
+            let resid2refnum = {};
+            for(let resid in ic.resid2refnum) {
+                let atom = ic.firstAtomObjCls.getFirstAtomObj(ic.residues[resid]);
+                if(!atom) continue;
+                
+                let resn = me.utilsCls.residueName2Abbr(atom.resn.substr(0, 3));
+        
+                let domainid = ic.resid2domainid[resid];
+                let refnumLabel = ic.resid2refnum[resid];
+        
+                if(refnumLabel) {
+                    let refnumStr_ori = ic.refnumCls.rmStrandFromRefnumlabel(refnumLabel);
+                    (ic.domainid2ig2kabat[domainid]) ? ic.domainid2ig2kabat[domainid][refnumStr_ori] : undefined;
+                }
+        
+                if(ic.resid2refnum[resid]) {
+                    if(ic.residIgLoop.hasOwnProperty(resid)) { // loop
+                    resid2refnum[resid + '_' + resn] = ic.resid2refnum[resid] + '_loop';
+                    }
+                    else {
+                    resid2refnum[resid + '_' + resn] = ic.resid2refnum[resid];
+                    }
+                }
+            }
+
+            refData += '{"ref PDB" : ' + JSON.stringify(ic.refPdbList) + ",\n";
+    
+            refData += '"data": {\n';
+            for(let chnid in ic.chains) {
+            let atom = ic.firstAtomObjCls.getFirstAtomObj(ic.chains[chnid]);
+            if(ic.proteins.hasOwnProperty(atom.serial)) {
+                for(let i = 0; i < ic.chainsSeq[chnid].length; ++i) {
+                const resid = chnid + '_' + ic.chainsSeq[chnid][i].resi + '_' + ic.chainsSeq[chnid][i].name;
+    
+                refData += "'" + resid + "': '" + resid2refnum[resid] + "',\n";
+                }
+            }
+            }
+            refData += '}\n';
+            refData += '}\n';
+        }
+        // 2. show Kabat ref numbers
+        else if(type == 'kabat' || type == 'Kabat') {
+            let resid2kabat = {};
+            for(let resid in ic.resid2refnum) {
+            let domainid = ic.resid2domainid[resid];
+            let refnumStr, refnumLabel = ic.resid2refnum[resid];
+    
+            let atom = ic.firstAtomObjCls.getFirstAtomObj(ic.residues[resid]);
+            let resn = me.utilsCls.residueName2Abbr(atom.resn.substr(0, 3));
+    
+            if(refnumLabel) {
+                let refnumStr_ori = ic.refnumCls.rmStrandFromRefnumlabel(refnumLabel);
+                refnumStr = (ic.domainid2ig2kabat[domainid]) ? ic.domainid2ig2kabat[domainid][refnumStr_ori] : undefined;
+            }
+    
+            resid2kabat[resid + '_' + resn] = refnumStr;
+            }
+            
+            refData += JSON.stringify(resid2kabat);
+        }
+        // 3. show IMGT ref numbers
+        else if(type == 'imgt'|| type == 'IMGT') {
+            let resid2imgt = {};
+            for(let resid in ic.resid2refnum) {
+            let domainid = ic.resid2domainid[resid];
+            let refnumStr, refnumLabel = ic.resid2refnum[resid];
+    
+            let atom = ic.firstAtomObjCls.getFirstAtomObj(ic.residues[resid]);
+            let resn = me.utilsCls.residueName2Abbr(atom.resn.substr(0, 3));
+    
+            if(refnumLabel) {
+                let refnumStr_ori = ic.refnumCls.rmStrandFromRefnumlabel(refnumLabel);
+                refnumStr = (ic.domainid2ig2imgt[domainid]) ? ic.domainid2ig2imgt[domainid][refnumStr_ori] : undefined;
+            }
+    
+            resid2imgt[resid + '_' + resn] = refnumStr;
+            }
+            
+            refData += JSON.stringify(resid2imgt);
+        }
+
+
+        if(!me.bNode) {
+            let file_pref = Object.keys(me.utilsCls.getHlStructures()).join(',');
+    
+            ic.saveFileCls.saveFile(file_pref + '_refnum_' + type + '.txt', 'text', [refData]);
+        }
+        else {
+            return refData;
+        }
+    }
  }
 
 /**
@@ -58355,6 +58505,11 @@ class ApplyCommand {
       else if(command == 'export pdb hydrogen') {
         await ic.scapCls.exportPdbProfix(true);
       }
+      else if(command.indexOf('export refnum ') != -1) {
+        let type = command.substr(14);
+        
+        ic.refnumCls.exportRefnum(type);
+      }
       else if(command == 'export secondary structure') {
          me.htmlCls.setHtmlCls.exportSecondary();
       }

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "icn3d",
-  "version": "3.26.7",
+  "version": "3.26.8",
   "main": "./icn3d.js",
   "exports": {
     ".": {