icn3d 3.29.17 → 3.30.0

package/icn3d.js

@@ -5810,6 +5810,26 @@ class ClickMenu {
         $("#" + me.pre + "menulist").html(html);
     }
 
+    async setIgTemplate(template) { let me = this.icn3dui, ic = me.icn3d;
+      ic.bRunRefnumAgain = true;
+
+      // reset for the selection
+      let residueArray = ic.resid2specCls.atoms2residues(Object.keys(ic.hAtoms));
+      for(let i = 0, il = residueArray.length; i < il; ++i) {
+         let resid = residueArray[i];
+
+         if(ic.resid2refnum) delete ic.resid2refnum[resid];
+         if(ic.resid2refnum_ori) delete ic.resid2refnum_ori[resid];
+         if(ic.resid2domainid) delete ic.resid2domainid[resid];
+      }
+
+      let bSelection = true;
+      // await ic.refnumCls.showIgRefNum(template);
+      if(!ic.bAnnoShown) await ic.showAnnoCls.showAnnotations();
+      await ic.annotationCls.setAnnoTabIg(bSelection, template);
+
+      // ic.bRunRefnumAgain = false;
+    }
 
     clickMenu1() { let me = this.icn3dui; me.icn3d;
         if(me.bNode) return;
@@ -7375,6 +7395,8 @@ class ClickMenu {
         });
 
         me.myEventCls.onIds("#" + me.pre + "mn6_igrefYes", "click", async function(e) { let ic = me.icn3d; //e.preventDefault();
+            ic.bRunRefnumAgain = true;
+
             thisClass.setLogCmd('ig refnum on', true);
             // await ic.refnumCls.showIgRefNum();
             // thisClass.setLogCmd('set annotation ig', true);
@@ -7391,21 +7413,22 @@ class ClickMenu {
             //    ic.hlUpdateCls.updateHlAll();
             //    ic.drawCls.draw();
             // }
+
+            // ic.bRunRefnumAgain = false;
          });
 
         me.myEventCls.onIds("#" + me.pre + "mn6_igrefTpl", "click", async function(e) { me.icn3d; //e.preventDefault();
             me.htmlCls.dialogCls.openDlg('dl_igrefTpl', 'Choose an Ig template');
          });
 
-        me.myEventCls.onIds("#" + me.pre + "mn6_igrefTpl_apply", "click", async function(e) { let ic = me.icn3d; //e.preventDefault();
+        me.myEventCls.onIds("#" + me.pre + "mn6_igrefTpl_apply", "click", async function(e) { me.icn3d; //e.preventDefault();
             if(!me.cfg.notebook) dialog.dialog( "close" );
-                  
+
             let template = $("#" + me.pre + "refTpl").val();
+
+            await thisClass.setIgTemplate(template);
+
             thisClass.setLogCmd('ig template ' + template, true);
-            let bSelection = true;
-            // await ic.refnumCls.showIgRefNum(template);
-            if(!ic.bAnnoShown) await ic.showAnnoCls.showAnnotations();
-            await ic.annotationCls.setAnnoTabIg(bSelection, template);
          });
 
          me.myEventCls.onIds("#" + me.pre + "mn6_alignrefTpl", "click", async function(e) { me.icn3d; //e.preventDefault();
@@ -7414,18 +7437,26 @@ class ClickMenu {
 
          me.myEventCls.onIds("#" + me.pre + "mn6_alignrefTpl_apply", "click", async function(e) { let ic = me.icn3d; //e.preventDefault();
             if(!me.cfg.notebook) dialog.dialog( "close" );
-                  
+   
             let template = $("#" + me.pre + "refTpl2").val();
 
+            let selAtoms = me.hashUtilsCls.cloneHash(ic.hAtoms);
+
             // load the template
             let url = me.htmlCls.baseUrl + "icn3d/refpdb/" + template + ".pdb";
             await ic.pdbParserCls.downloadUrl(url, 'pdb', undefined, template);
             thisClass.setLogCmd('load url ' + url + ' | type pdb', true);
+   
+            let structure = template.replace(/_/g, '').substr(0,4);
+
+            let chainid = ic.structures[structure][0];
+
+            ic.hAtoms = me.hashUtilsCls.unionHash(selAtoms, ic.chains[chainid]);
 
             // align the template with the selection
             me.cfg.aligntool = 'tmalign';
             await ic.realignParserCls.realignOnStructAlign();
-            thisClass.setLogCmd('realign on tmalign', true);
+            thisClass.setLogCmd('realign on tmalign', true);   
          });
 
          me.myEventCls.onIds("#" + me.pre + "mn6_igrefNo", "click", async function(e) { let ic = me.icn3d; //e.preventDefault();
@@ -37609,9 +37640,10 @@ class AnnoIg {
 
     //Show the annotations of CDD domains and binding sites.
     async showIg(chnid, template) { let ic = this.icn3d; ic.icn3dui;
-        if(!ic.bRunRefnum || Object.keys(ic.atoms).length > Object.keys(ic.hAtoms).length) {
+        // if(!ic.bRunRefnum || Object.keys(ic.atoms).length > Object.keys(ic.hAtoms).length) {
+        if(ic.bRunRefnumAgain) {
             await ic.refnumCls.showIgRefNum(template);
-            ic.bRunRefnum = true;    
+            // ic.bRunRefnum = true;    
         }
 
         let type = 'ig';
@@ -37629,6 +37661,8 @@ class AnnoIg {
         $("#" + ic.pre + "dt_" + type + "_" + chnid).html(html);
         $("#" + ic.pre + "ov_" + type + "_" + chnid).html(html2);
         $("#" + ic.pre + "tt_" + type + "_" + chnid).html(html3);
+
+        ic.bRunRefnumAgain = false;
     }
 
     showAllRefNum(giSeq, chnid) {  let ic = this.icn3d; ic.icn3dui;
@@ -38131,6 +38165,7 @@ class AnnoDomain {
         let pdbid = pdbArray[index];
         //let url = me.htmlCls.baseUrl + "mmdb/mmdb_strview.cgi?v=2&program=icn3d&domain&molinfor&uid=" + pdbid;
 
+/*
         if(!ic.bResetAnno && index == 0 && ic.mmdb_data !== undefined) {      
             for(let chnid in ic.protein_chainid) {
                 if(chnid.indexOf(pdbid) !== -1) {
@@ -38146,6 +38181,7 @@ class AnnoDomain {
             }
         }
         else {
+*/
             // calculate 3D domains on-the-fly
             //ic.protein_chainid[chainArray[i]] 
             let data = {};
@@ -38187,7 +38223,7 @@ class AnnoDomain {
 
             ic.bAjax3ddomain = true;
             ic.bAjaxDoneArray[index] = true;          
-        }
+        // }
     }
 
     //Show the annotations of 3D domains.
@@ -39222,6 +39258,10 @@ class Domain3d {
     constructor(icn3d) {
 		this.icn3d = icn3d;
 
+		this.init3ddomain();
+	}
+
+	init3ddomain() { let ic = this.icn3d; ic.icn3dui;
         //this.dcut = 8; // threshold for C-alpha interactions
 
 		// It seemed the threshold 7 angstrom works better
@@ -39229,7 +39269,8 @@ class Domain3d {
 		this.dcut = 8; // threshold for C-alpha interactions
 
 		// added by Jiyao
-		this.min_contacts = 5; //3;			// minimum number of contacts to be considered as neighbors
+		// pdbid 1CD8 requires a min contact 4, not 5
+		this.min_contacts = 4; //5; //3;			// minimum number of contacts to be considered as neighbors
 
 		this.MAX_SSE = 512;
 
@@ -39241,6 +39282,7 @@ class Domain3d {
 
         //let this.elt_size[this.MAX_SSE];			// element sizes in residues
         this.elt_size = [];
+
         this.elt_size.length = this.MAX_SSE;
 
         //let this.group_num[this.MAX_SSE];			// indicates required element groupings
@@ -39850,6 +39892,8 @@ class Domain3d {
 	// x0, y0, z0: array of x,y,z coordinates of C-alpha atoms
 	//c2b_NewSplitChain(chnid, dcut) { let ic = this.icn3d, me = ic.icn3dui;
 	c2b_NewSplitChain(atoms, dcut) { let ic = this.icn3d; ic.icn3dui;
+		this.init3ddomain();
+
 		let x0 = [], y0 = [], z0 = [], resiArray = [];
 
 		//substruct: array of secondary structures, each of which has the keys: From (1-based), To (1-based), Sheet (0 or 1), also add these paras: x1, y1, z1, x2, y2, z2
@@ -42383,7 +42427,11 @@ class Annotation {
         this.updateSsbond();
         this.updateCrosslink();
         await this.updateTransmem();
+
+        ic.bRunRefnumAgain = true;
         await this.updateIg();
+        // ic.bRunRefnumAgain = false;
+
         this.updateInteraction();
     }
     hideAnnoTabAll() {  let ic = this.icn3d; ic.icn3dui;
@@ -42465,10 +42513,13 @@ class Annotation {
             $("[id^=" + ic.pre + "transmem]").show();
         }
         if($("#" + ic.pre + "anno_ig").length && $("#" + ic.pre + "anno_ig")[0].checked) {
-            ic.bRunRefnum = false;
+            // ic.bRunRefnumAgain = true;
+
             await this.updateIg();
 
             $("[id^=" + ic.pre + "ig]").show();
+
+            // ic.bRunRefnumAgain = false;
         }
     }
     setAnnoTabCustom() {  let ic = this.icn3d; ic.icn3dui;
@@ -42729,12 +42780,15 @@ class Annotation {
 
         me.myEventCls.onIds("#" + ic.pre + "anno_ig", "click", async function(e) {
             if($("#" + ic.pre + "anno_ig").length && $("#" + ic.pre + "anno_ig")[0].checked) {
-                if(Object.keys(ic.atoms).length > Object.keys(ic.hAtoms).length) {
-                    ic.bRunRefnum = false;
-                }
+                // if(Object.keys(ic.atoms).length > Object.keys(ic.hAtoms).length) {
+                //     ic.bRunRefnum = false;
+                // }
 
+                ic.bRunRefnumAgain = true;
                 await thisClass.setAnnoTabIg();
                 me.htmlCls.clickMenuCls.setLogCmd("set annotation ig", true);
+
+                // ic.bRunRefnumAgain = false;
             }
             else {
                 thisClass.hideAnnoTabIg();
@@ -59959,11 +60013,17 @@ class ApplyCommand {
           await ic.annotationCls.setAnnoTabTransmem();
       }
       else if(command == 'set annotation ig') {
+          ic.bRunRefnumAgain = true;
           await ic.annotationCls.setAnnoTabIg();
+          // ic.bRunRefnumAgain = false;
       }
       else if(command == 'ig refnum on') {
+        ic.bRunRefnumAgain = true;
+
         if(!ic.bAnnoShown) await ic.showAnnoCls.showAnnotations();
         await ic.annotationCls.setAnnoTabIg(true);
+
+        // ic.bRunRefnumAgain = false;
     }
       else if(command == 'highlight level up') {
           ic.resid2specCls.switchHighlightLevelUp();
@@ -62451,9 +62511,7 @@ class LoadScript {
           // }
           else if(command.indexOf('ig template') == 0 ) { 
             let template = command.substr(command.lastIndexOf(' ') + 1);
-            // await ic.refnumCls.showIgRefNum(template);
-            if(!ic.bAnnoShown) await ic.showAnnoCls.showAnnotations();
-            await ic.annotationCls.setAnnoTabIg(true, template);
+            await me.htmlCls.clickMenuCls.setIgTemplate(template);
           }
           else if(command.indexOf('set annotation 3ddomain') == 0) { // the command may have "|||{"factor"...
               if(Object.keys(ic.proteins).length > 0) {
@@ -62656,8 +62714,12 @@ class LoadScript {
                     }
                     else if(lastCommand.indexOf('ig refnum on') == 0) {
                         // await ic.refnumCls.showIgRefNum();
+                        ic.bRunRefnumAgain = true;
+
                         if(!ic.bAnnoShown) await ic.showAnnoCls.showAnnotations();
                         await ic.annotationCls.setAnnoTabIg(true);
+
+                        // ic.bRunRefnumAgain = false;
                     }
                     else if(lastCommand.indexOf('set annotation 3ddomain') == 0) {
                         thisClass.applyCommand3ddomain(lastCommand);
@@ -65654,7 +65716,7 @@ class Dssp {
 
     async hideIgRefNum() { let ic = this.icn3d; ic.icn3dui;
         ic.bShowRefnum = false;
-        ic.bRunRefnum = false;
+        // ic.bRunRefnum = false;
 
         // redo all ref numbers
         ic.resid2refnum = {};
@@ -66257,7 +66319,7 @@ class Dssp {
                 strandHash[strand] = 1;
             }
             let score = parseFloat(queryData[0].score);
-            //!!!
+
             // if the TM score difference is within 0.1 and more strands are found, use the template with more strands
             // if(!domainid2segs.hasOwnProperty(domainid) || 
             //     (score >= parseFloat(ic.domainid2score[domainid].split('_')[0]) + tmAdjust)
@@ -66299,7 +66361,7 @@ class Dssp {
                     refpdbnameList.sort(function(a, b) {
                         return refpdbname2score[b] - refpdbname2score[a]
                     });
-                    // top 3 templates
+                    // top 4 templates
                     ic.domainid2refpdbname[domainid] = refpdbnameList.slice(0,4);
                 }
             }
@@ -66313,7 +66375,7 @@ class Dssp {
 
         let domainid2segs = this.parseAlignData_part1(dataArray, domainidpairArray, bRound1);
 
-        // !!!no more Igs to detect
+        // no more Igs to detect
         if(Object.keys(domainid2segs).length == 0) {
             bNoMoreIg = true;
             return bNoMoreIg;
@@ -67062,7 +67124,7 @@ class Dssp {
                             }
                             
                             // in case A1550 is not found, but A1551 is found
-                            if(!bFoundAnchor && (lastTwo == 51 || lastTwo == 52 || lastTwo == 53 || lastTwo == 54) ) {
+                            if(!bFoundAnchor && (lastTwo >= 46 && lastTwo <= 54) ) {
                                 let offset = lastTwo - 50;
                                 strandArray[strandCnt].anchorRefnum = refnum - offset;
                                 strandArray[strandCnt].resCntBfAnchor = resCntBfAnchor - offset;
@@ -75593,7 +75655,7 @@ class iCn3DUI {
     //even when multiple iCn3D viewers are shown together.
     this.pre = this.cfg.divid + "_";
 
-    this.REVISION = '3.29.5';
+    this.REVISION = '3.30.0';
 
     // In nodejs, iCn3D defines "window = {navigator: {}}"
     this.bNode = (Object.keys(window).length < 2) ? true : false;