Difference between revisions of "User:Remig/plico/modify"

From Jmol
Jump to navigation Jump to search
(Jmol script to modify amino acids and nucleotides by mouse actions)
 
m (Remove debug stuff)
Line 11: Line 11:
 
Copy and paste the following to a text editor and save to your scripts folder as modify.spt:
 
Copy and paste the following to a text editor and save to your scripts folder as modify.spt:
 
<pre>#  modify - Jmol script by Ron Mignery
 
<pre>#  modify - Jmol script by Ron Mignery
#  v1.0 beta    2/22/2014 for Jmol 14
+
#  v1.1 beta    2/22/2014 for Jmol 14
 
#
 
#
 
# See atom hover text for options:
 
# See atom hover text for options:
Line 30: Line 30:
 
# Return L tetrahedron point if i1<i2<i3, else R point
 
# Return L tetrahedron point if i1<i2<i3, else R point
 
function getTet(i1, i2, i3, dist) {
 
function getTet(i1, i2, i3, dist) {
print format("getTet(%d, %d, %d, dist", i1, i2, i3)
 
 
     var v1 = {atomIndex=i3}.xyz - {atomIndex=i2}.xyz
 
     var v1 = {atomIndex=i3}.xyz - {atomIndex=i2}.xyz
 
     var v2 = {atomIndex=i1}.xyz - {atomIndex=i2}.xyz
 
     var v2 = {atomIndex=i1}.xyz - {atomIndex=i2}.xyz
Line 48: Line 47:
 
# return vector parallel to previous bond
 
# return vector parallel to previous bond
 
function getTet1(pIdx, idx) {
 
function getTet1(pIdx, idx) {
print format("getTet1(%d, %d)", pIdx, idx)
 
 
     var cset = connected({atomIndex=pIdx}) and not {atomIndex=idx}
 
     var cset = connected({atomIndex=pIdx}) and not {atomIndex=idx}
 
     return {atomIndex=idx}.xyz + {atomIndex=pIdx}.xyz - cset[1].xyz
 
     return {atomIndex=idx}.xyz + {atomIndex=pIdx}.xyz - cset[1].xyz
Line 64: Line 62:
  
 
     return lpt + {atomIndex=i2}.xyz
 
     return lpt + {atomIndex=i2}.xyz
}
 
 
# Selected must include third parameter but not the first
 
function setAngleIdx (statorIdx, pivotIdx, rotorIdx, toangle) {
 
    try {
 
        var v1={atomIndex=statorIdx}.xyz - {atomIndex=pivotIdx}.xyz
 
        var v2={atomIndex=rotorIdx}.xyz - {atomIndex=pivotIdx}.xyz
 
        var axis = cross(v1, v2) + {atomIndex=pivotIdx}.xyz
 
        var curangle =  angle({atomIndex=statorIdx},
 
            {atomIndex=pivotIdx}, {atomIndex=rotorIdx})
 
        rotateselected @axis {atomIndex = pivotIdx} @{curangle-toangle}
 
    }
 
    catch {
 
    }
 
 
}
 
}
  
 
function addToIdx(idx, aElement, greek, addOs) {
 
function addToIdx(idx, aElement, greek, addOs) {
print format("addToIdx(%d, %s, %s, n)", idx, aElement, greek)
 
 
     var iGroup = {atomIndex=idx}.group
 
     var iGroup = {atomIndex=idx}.group
 
     var iResno = {atomIndex=idx}.resno
 
     var iResno = {atomIndex=idx}.resno
Line 217: Line 200:
  
 
function toDidx(idx) {
 
function toDidx(idx) {
print format("toDidx(%d)", idx)
 
 
     delay
 
     delay
 
     gIsPlanar = FALSE
 
     gIsPlanar = FALSE
Line 234: Line 216:
 
     var a1 = angle(n1, c6, c5, c4)\10
 
     var a1 = angle(n1, c6, c5, c4)\10
 
     var a2 = angle(c6, c5, c4, n3)\10
 
     var a2 = angle(c6, c5, c4, n3)\10
print format("a1=%d a2=%d", a1, a2)
+
   
 
     # N1-C6-C5-C4 = 0.5    C6-C5-C4-N3 = 0.8
 
     # N1-C6-C5-C4 = 0.5    C6-C5-C4-N3 = 0.8
 
     if ((a1 == 0) or (a1 == -1)) {
 
     if ((a1 == 0) or (a1 == -1)) {
Line 396: Line 378:
  
 
function modify1MB() {
 
function modify1MB() {
print "modify1MB"
 
 
     var idx = _atomPicked
 
     var idx = _atomPicked
 
     var isAmino = ({(atomIndex=idx) and {amino}}.size > 0)
 
     var isAmino = ({(atomIndex=idx) and {amino}}.size > 0)
Line 558: Line 539:
 
     setHoverLabels()
 
     setHoverLabels()
 
      
 
      
print "modify2MB"
 
 
}
 
}
  
 
function selectChainMB() {
 
function selectChainMB() {
print "selectChainMB"
 
 
     var cset = {atomIndex=_atomPicked}
 
     var cset = {atomIndex=_atomPicked}
 
     select cSet
 
     select cSet

Revision as of 18:02, 22 February 2014

Modify allows the user to make common amino acid and nucleotide modification by mouse actions. Hover text above each atom gives the available options including methylation on nucleosides, single-bond oxygens, and nitrogens, hydroxylation of PRO, formylation of MET, ribose/deoxyribose interconversion, phophorylation of SER, THR, TYR and HIS, and conversions to pseudouridine, dihyrouridine, inosine, and xanthine.

Some features unrelated to modification have been placed here as well to serve other Plico functions. Water molecules may be added and moved around to serve as pivot points for Tug when no other atoms are where the pivot point is wanted. Also the Jmol parameters allowMoveAtoms and allowRotateSelected are set true while modify is active to allow free rotation and translation of selected atoms.

Modify is a member of the Plico suite of protein folding tools described here. It may be installed and accessed as a macro with the file:

Title=PLICO Modify
Script=script <path to your scripts folder>/modify.spt;plicoModify

saved as plicoModify.macro in your .jmol/macros folder as described in Macro.

Copy and paste the following to a text editor and save to your scripts folder as modify.spt:

#   modify - Jmol script by Ron Mignery
#   v1.1 beta    2/22/2014 for Jmol 14
#
# See atom hover text for options:
# methyl on nucleosides, single-bond oxygens, and terminal N
# OH on ribose and proline, PO3 on SER, THR, TYR and HIS (2)
# pseudouridine and dihyrouridine, inosine, xanthine
# HOH singlet anywhere for phantom pivot and target atoms
# Enable jmol move for duration

var kCtolerance = 1.85
var gAllowMoveAtoms = allowMoveAtoms
var gAllowRotateSelected = allowRotateSelected
var gSelsave = ({})
var gScheme = defaultColorScheme
var gAltScheme = ((gScheme == "Jmol") ? "Rasmol" : "Jmol")
var gIsPlanar = FALSE

# Return L tetrahedron point if i1<i2<i3, else R point
function getTet(i1, i2, i3, dist) {
    var v1 = {atomIndex=i3}.xyz - {atomIndex=i2}.xyz
    var v2 = {atomIndex=i1}.xyz - {atomIndex=i2}.xyz
    var axis = cross(v1, v2)
    var pma = ({atomIndex=i1}.xyz + {atomIndex=i3}.xyz)/2
    var pmo = {atomIndex=i2}.xyz + {atomIndex=i2}.xyz - pma
    var pt = pmo + (axis/axis)

    var v = pt - {atomIndex=i2}.xyz
    var cdist = distance(pt, {atomIndex=i2})
    var factor = (dist/cdist)
    var lpt = v * factor

    return lpt + {atomIndex=i2}.xyz
}

# return vector parallel to previous bond
function getTet1(pIdx, idx) {
    var cset = connected({atomIndex=pIdx}) and not {atomIndex=idx}
    return {atomIndex=idx}.xyz + {atomIndex=pIdx}.xyz - cset[1].xyz
}

function getTrigonal(i1, i2, i3, dist) {
    var v1 = {atomIndex=i1}.xyz - {atomIndex=i2}.xyz
    var v2 = {atomIndex=i3}.xyz - {atomIndex=i2}.xyz
    var pt = {atomIndex=i2}.xyz - (v1 + v2)

    var v = pt - {atomIndex=i2}.xyz
    var cdist = distance(pt, {atomIndex=i2})
    var factor = (dist/cdist)
    var lpt = (v * factor)

    return lpt + {atomIndex=i2}.xyz
}

function addToIdx(idx, aElement, greek, addOs) {
    var iGroup = {atomIndex=idx}.group
    var iResno = {atomIndex=idx}.resno
    var iChain = {atomIndex=idx}.chain
    var iAtomName = {atomIndex=idx}.atomName
    var iAtomno = {atomIndex=idx}.atomno
    var aAtomno = {chain=iChain}.atomno.max + 1
    var cset = connected({atomIndex=idx})
    var pt = {0 0 0}
    if  ((iGroup = "HIS") or ((({atomIndex=idx} and {amino}).size == 0) 
        and (iAtomName[1][2] != "OP") and (iAtomName[3] != "\'"))) {
        pt = getTrigonal(cSet[1].atomIndex, idx, cSet[2].atomIndex, 1.5) 
    }
    else {
        if ((cSet.size == 2) and not greek.find("\'")) {
            pt = getTet(cSet[1].atomIndex, idx, cSet[2].atomIndex, 1.5) 
        }
        else if (cSet.size > 2) {
            pt = getTet(cSet[2].atomIndex, idx, cSet[1].atomIndex, 1.5) 
        }
        else {
            pt = getTet1(cSet[1].atomIndex, idx) 
        }
    }
    if (greek.size < 2) {
        greek += iAtomName[2]
    }
    var ls = "data \"append addto\"\n"
    ls += format("ATOM  %5d  %-4s%3s ", aAtomno, aElement + greek, iGroup)
    ls += format("%s%4d    %8.3f", iChain, iResno, pt.x)          
    ls += format("%8.3f%8.3f\n", pt.y, pt.z)
    ls += "end \"append addto\""
    script inline @{ls}
        
    var aIdx = {(atomno=aAtomno) and (chain=iChain)}.atomIndex
    connect {atomIndex=idx} {atomIndex=aIdx}

    if (addOs > 0) {    
        for (var i = 1; i <= addOs; i++) {
            aAtomno++
            var ccSet = connected({atomIndex=aIdx})
            if (addOs == 1) {
                pt = ({atomIndex=aIdx}.xyz +
                    ({atomIndex=aIdx}.xyz - {atomIndex=idx}.xyz))
            }
            else if (addOs == 2) {
                pt = getTrigonal(ccSet[1].atomIndex, aIdx, ccSet[2].atomIndex, 1.5) 
            }
            else {
                if (ccSet.size == 2) {
                    pt = getTet(ccSet[1].atomIndex, aIdx, ccSet[2].atomIndex, 1.5) 
                }
                else if (ccSet.size == 3) {
                    pt = getTet(ccSet[2].atomIndex, aIdx, ccSet[1].atomIndex, 1.5) 
                }
                else {
                    pt = getTet1(ccSet[1].atomIndex, aIdx) 
                }
            }
            ls = "data \"append addto\"\n"
            ls += format("ATOM  %5d  O%s%d ", aAtomno,
                {atomIndex=aIdx}.AtomName[1], i)
            ls += format("%3s %s%4d    ", iGroup, iChain, iResno)          
            ls += format("%8.3f%8.3f%8.3f\n",  pt.x, pt.y, pt.z)
            ls += "end \"append addto\""
            script inline @{ls}
            oIdx = {(atomno=aAtomno) and (chain=iChain)}.atomIndex
            connect {atomIndex=aIdx} {atomIndex=oIdx}
        }
    }
    setHoverLabels()
}

function removeFromIdx(idx) {
    var cSet = connected({atomIndex=idx})
    if ({atomIndex=idx}.atomName[2] != "M") {
        for (var i = 1; i <= cSet.size; i++) {
            if ((cSet[i].atomName[2] == "P") or (cSet[i].atomName[2] == "O")) {
                delete @{cSet[i]}
            }
        }
    }
}

function toXidx(idx, newElement, newGroup) {
    var cSet = {resno=@{{atomIndex=idx}.resno}}
    var ls = "data \"append toX\"\n"
    for (var i = 1; i <= cSet.size; i++) {
        var iName = cSet[i].atomName
        if (idx == cSet[i].atomIndex) {
            iName[1] = newElement
        }
        ls += format("ATOM  %5d  %-3s ", cSet[i].atomno, iName)
        ls += format("%3s %s%4d    ", newGroup,
            {atomIndex=idx}.chain, cSet[i].resno)          
        ls += format("%8.3f%8.3f%8.3f\n",  cSet[i].x, cSet[i].y, cSet[i].z)
    }
    ls += "end \"append toX\""
    delete cSet
    script inline @{ls}
    connect
}

function toPUidx(idx, undoIt) {
    var iResno = {atomIndex=idx}.resno
    var a1idx = {(resno=iResno) and (atomName="O2")}.atomIndex    
    var a2idx = {(resno=iResno) and (atomName="O4")}.atomIndex
    var xyz = {0 0 0} + {atomIndex=a1idx}.xyz
    {atomIndex=a1idx}.xyz = {0 0 0} + {atomIndex=a2idx}.xyz    
    {atomIndex=a2idx}.xyz = xyz
    
    a1idx = {(resno=iResno) and (atomName="C2")}.atomIndex
    a2idx = {(resno=iResno) and (atomName="C4")}.atomIndex
    xyz = {0 0 0} + {atomIndex=a1idx}.xyz
    {atomIndex=a1idx}.xyz = {0 0 0} + {atomIndex=a2idx}.xyz    
    {atomIndex=a2idx}.xyz = xyz
    
    a1idx = {(resno=iResno) and (atomName="N1")}.atomIndex
    a2idx = {(resno=iResno) and (atomName="C5")}.atomIndex
    xyz = {0 0 0} + {atomIndex=a1idx}.xyz
    {atomIndex=a1idx}.xyz = {0 0 0} + {atomIndex=a2idx}.xyz    
    {atomIndex=a2idx}.xyz = xyz
    
    var a3idx = {(resno=iResno) and (atomName="c1\'")}.atomIndex
    connect (atomIndex=a3idx) (atomindex=a1idx) DELETE
    connect (atomIndex=a3idx) (atomindex=a2idx)
    
    if (undoit) {
        toXidx(idx, "C", "U")
    }
    else {
        toXidx(idx, "N", "PU")
    }
}

function toDidx(idx) {
    delay
    gIsPlanar = FALSE
    var iResno = {atomIndex=idx}.resno
    var iChain = {atomIndex=idx}.chain
    var rSet = {(resno=iResno) and (chain=iChain)}
    var n1 = {rSet and (atomName="N1")}
    var c2 = {rSet and (atomName="C2")}   
    var o2 = {rSet and (atomName="O2")}   
    var n3 = {rSet and (atomName="N3")}
    var c4 = {rSet and (atomName="C4")}   
    var o4 = {rSet and (atomName="O4")}   
    var c5 = {rSet and (atomName="C5")}
    var c6 = {rSet and (atomName="C6")}

    var a1 = angle(n1, c6, c5, c4)\10
    var a2 = angle(c6, c5, c4, n3)\10
    
    # N1-C6-C5-C4 = 0.5     C6-C5-C4-N3 = 0.8
    if ((a1 == 0) or (a1 == -1)) {
        # Planar to Chair 1    
        axis = (c2.xyz - c6.xyz + n1.xyz)
        select c2 or o2 or n3 or c4 or o4 or c5 or c6
        rotateSelected @n1 @axis 40
        select o2 or n3 or c4 or o4 or c5
        rotateSelected @c2 @c6 40
        select c4 or o4
        rotateSelected @c5 @n3 40
    }
    
    # N1-C6-C5-C4 = -41.0     C6-C5-C4-N3 = -46.0
    else if ((a1 == -5) and ((a2 == 5) or (a2 == 4))) {
        # Chair 1 to boat 1 twist 1 `
        select c4 or o4
        rotateSelected @c5 @n3 -80
        select o2 or n3 or c4 or o4 or c5
        rotateSelected @n1 @c4 25
    }
    
    # N1-C6-C5-C4 = 41.6     C6-C5-C4-N3 = 23.3
    else if (((a1 == 3) or (a1 == 4)) and (a2 == -7)) {
        # Boat 1 twist 1 to boat 1 twist 2
        rotateSelected @n1 @c4 -50
    }
    
    # N1-C6-C5-C4 = -40.2     C6-C5-C4-N3 = 70.0
    else if (((a1 == -4) or (a1 == -5)) and ((a2 == -3) or (a2 == -2))) {
        # Boat1 twist 2 to chair 2    
        rotateSelected @n1 @c4 25
        axis = (c2.xyz - c6.xyz + n1.xyz)
        select c2 or o2 or n3 or c4 or o4 or c5 or c6
        rotateSelected @n1 @axis -80
        select o2 or n3 or c4 or o4 or c5
        rotateSelected @c2 @c6 -80
    }
    
    # N1-C6-C5-C4 = 41.8     C6-C5-C4-N3 = -44.10
    else if ((a1 == 4) and (a2 == -5)) {
        # Chair 2 to boat 2 twist 1` `
        select c4 or o4
        rotateSelected @c5 @n3 80
        select o2 or n3 or c4 or o4 or c5
        rotateSelected @n1 @c4 25
    }
    
    # N1-C6-C5-C4 = 41.2     C6-C5-C4-N3 = -68.5
    else if ((a1 == 4) and ((a2 == 1) or (a2 == 2))) {
        # Boat 2 twist 1 to boat 2 twist 2
        rotateSelected @n1 @c4 -50
    }
    
    # N1-C6-C5-C4 = -41.3     C6-C5-C4-N3 = 64.5
    else if ((a1 == -5) and (a2 == 6)) {
        # Boat 2 twist 2 to planar
        rotateSelected @n1 @c4 25
        select c4 or o4
        rotateSelected @c5 @n3 -40
        axis = (c2.xyz - c6.xyz + n1.xyz)
        select c2 or o2 or n3 or c4 or o4 or c5 or c6
#quit        
        rotateSelected @n1 @axis 40
        select o2 or n3 or c4 or o4 or c5
        rotateSelected @c2 @c6 40
        gIsPlanar = TRUE
    }
    
    else {
        print format("a1=%d a2=%d", a1, a2)
    }
    
}

function setHoverLabels() {
    select all
    set hoverLabel "%U"
    
    select ({(atomName="?O?")} or {(atomName="?P?") and {amino}})
    set hoverLabel "%U|1:  -X"
    
    select (oxygen or sulfur) and not *.O and (bondCount = 1)
    select {selected} or (nitrogen and (bondCount < 3))
    set hoverLabel "%U|1:  -CH3"
    select (not {amino} and not oxygen and (bondCount < 3))
    set hoverLabel "%U|1:  -CH3"
    
    select ({(atomName="?M?")} or {group="HOH"})
    set hoverLabel "%U|1:  -X"
    
    # MET N fMET
    select (MET.N and (bondCount == 1))
    set hoverLabel "%U|1:  -CH3|2: -CO"
    select (MET.N and (bondCount == 2))
    set hoverLabel "%U|1:  -CH3|2: -X"
        
    # ASP OD1 OD2 P
    select ASP.OD1 or ASP.OD2 or HIS.ND1 or SER.OG or THR.OG1 OR TYR.OH
    set hoverLabel "%U|1:  -CH3|2:  -PO3"
    
    # PRO CG HYP
    select PRO.CG
    set hoverLabel "%U|1: -CH3|2:  -OH"
    
    # A N6 I
    select [A].N6 or [I].O6
    set hoverLabel "%U|1:  -CH3|2: A <--> I"
    
    # G N2 X
    select [G].N2 or [X].O2
    set hoverLabel "%U|1:  -CH3|2: G <--> X"
    
    # U N1 PU
    select [U].N1 or [PU].C5
    set hoverLabel "%U|1:  -CH3|2: U <--> PU"
    
    # U C5-6 HU
    select [U].C6
    set hoverLabel "%U|1:  -CH3|2: U --> D"
    
    # U C5-6 HU
    select [D].C6
    set hoverLabel "%U|1:  ->Conformers|2: D --> U"
    
    # U C5 M T
    select [U].C5 OR [T].C5
    set hoverLabel "%U|1:  -CH3|2: U <--> T"
    
    # *.?2' -X
    select ??.?2'
    set hoverLabel "%U|1:  -CH3|2: X <--> DX"
    refresh
 }

function addWaterMB() {
    var idx = _atomPicked
    var iChain = {atomIndex=idx}.chain
    var pt = {0 0 0}
    var cSet = connected({atomIndex=idx})
    if (cSet.size == 2) {
        pt = getTet(cSet[1].atomIndex, idx, cSet[2].atomIndex, 2.3) 
    }
    else if (cSet.size == 3) {
        pt = getTet(cSet[2].atomIndex, idx, cSet[1].atomIndex, 2.3) 
    }
    else {
        var v = {atomIndex=idx}.xyz - cSet[1].xyz
        pt = {atomIndex=idx}.xyz + v + v
    }
    if (within(kCtolerance, FALSE, pt).size == 0) {
        var ls = "data \"append addto\"\n"
        ls += format("ATOM  %5d  O   HOH %s%4d    ",
             {*}.atomno.max + 1, iChain, {*}.resno.max + 1)          
        ls += format("%8.3f%8.3f%8.3f\n", pt.x, pt.y, pt.z)
        ls += "end \"append addto\""
        script inline @{ls}
        setHoverLabels()
    }
}

function modify1MB() {
    var idx = _atomPicked
    var isAmino = ({(atomIndex=idx) and {amino}}.size > 0)
    var iName = {atomIndex=idx}.atomName
    var iGroup = {atomIndex=idx}.group
    
    if (iGroup == "HOH") {
        delete {atomIndex=idx}
    }
    
    else if (iGroup == "D") {
        toDidx(idx)
        if(gIsPlanar) {
            toDidx(idx)
        }
    }
    
    else if ((iName[2]="M") or (iName[2]="O")
        or (isAmino and ((iName[1]="P") or (iName[2]=="P")))) {
        removeFromIdx(idx)
        delete {atomIndex=idx}
    }
    
    else {
        if ({atomIndex=idx}.element="P") {
            if (isAmino) {
                removeFromIdx(idx)
            }
        }
        else if (({atomIndex=idx}.element="O")
            and ({atomIndex=idx}.atomName!="O")
            and ({atomIndex=idx}.bondCount=1)) {
            addToIdx(idx, "C", "M", 0)
        }
        else if (({atomIndex=idx}.element!="O")
            and ({atomIndex=idx}.bondCount < 3)) {
            addToIdx(idx, "C", "M", 0)
        }
        else if (isAmino) {
            removeFromIdx(idx)
        }
    }
}

function modify2MB() {
    var idx = _atomPicked
    var iGroup = {atomIndex=idx}.group
    var iName = {atomIndex=idx}.atomName
    var iResno = {atomIndex=idx}.resno
    var iChain = {atomIndex=idx}.chain
    
    # MET.N fMET
    if ((iGroup="MET") and (iName="N")) {
        if ({atomIndex=idx}.bondCount == 1) {
            addToIdx(idx, "C", "F", 1)
        }
        else {
             delete {(resno=iResno) and (chain=iChain) and (atomName="CF")}
             delete {(resno=iResno) and (chain=iChain) and (atomName="OC1")}
        }
    }
    
    # ASP.OD1 or ASP.O or HIS.ND1 or SER.OG or THR.OG1 OR TYR.OH PO3
    else if (((iGroup="HIS") and (iName="ND1"))
        or ((iGroup="SER") and (iName="OG"))
        or ((iGroup="THR") and (iName="OG1"))
        or ((iGroup="TYR") and (iName="OH"))
        or ((iGroup="ASP") and (iName[1][2]="OD"))) {
        
        if (({atomIndex=idx}.bondCount == 1)
            or ((iName="ND1") and ({atomIndex=idx}.bondCount == 2))) {
            addToIdx(idx, "P", "", 3)
        }
        else {
            removeFromIdx(idx)
        }
    }
    
    # PRO.CG HYP
    else if ((iGroup="PRO") and (iName="CG")) {
        if ({atomIndex=idx}.bondCount == 2) {
            addToIdx(idx, "O", "H", 0)
        }
        else {
            removeFromIdx(idx)
        }
    }
    
    # A.N6 I
    else if ((iGroup="A") and (iName="N6")) {
        toXidx(idx, "O", "I")
    }
    else if ((iGroup="I") and (iName="O6")) {
        toXidx(idx, "N", "A")
    }
    
    # G.N2 X
    else if ((iGroup="G") and (iName="N2")) {
        toXidx(idx, "O", "X")
    }
    else if ((iGroup="X") and (iName="O2")) {
        toXidx(idx, "N", "G")
    }
    
    # U.N1 PU
    else if ((iGroup="U") and (iName="N1")) {
        toPUidx(idx, FALSE)
    }
    else if ((iGroup="PU") and (iName="C5")) {
        toPUidx(idx, TRUE)
    }
    
    # U.C5 T
    else if ((iGroup="U") and (iName="C5")) {
        addToIdx(idx, "C", "7 ", 0)
        toXidx(idx, "C", "T")
    }
    else if ((iGroup="T") and (iName="c5")) {
        delete {(resno=iResno) and (atomName="C7")}
        toXidx(idx, "C", "U")
    }
    
    # U.C6 HU
    else if ((iGroup="U") and ((iName="C5") or (iName="C6"))) {
        toXidx(idx, "C", "D")
        toDidx(idx)
    }
    else if ((iGroup="D") and ((iName="C5") or (iName="C6"))) {
        for (var i = 0; i < 10; i++) {
            toDidx(idx)
            if (gIsPlanar) {
                break
            }
        }
        toXidx(idx, "C", "U")
    }
    
    # O2'
    else if (iName="O2\'") {
        var cIdx = connected({atomIndex=idx})[1].atomIndex
        delete {atomIndex=idx}
        toXidx(cIdx, "C", "D"+iGroup)
    }
    
    # C2'
    else if (iName="C2\'") {
        var ccSet = connected({atomIndex=idx})
        if ((ccSet.size < 3) and (iGroup.size > 1)) {
            addToIdx(idx, "O", "2\'", 0)
            toXidx(idx, "C", iGroup[2])
        }
        else if (iGroup.size == 1) {
            for (var i = 1; i <= ccSet.size; i++) {
                if (ccSet[i].element == "O") {
                    delete {atomIndex=@{ccSet[i].atomIndex}}
                    toXidx(idx, "C", "D"+iGroup)
                }
            }
        }
    }
    setHoverLabels()
    
}

function selectChainMB() {
    var cset = {atomIndex=_atomPicked}
    select cSet
    while (cSet.size > 0) {
        cSet = connected({selected}) and not {selected}
        select {selected} or cSet
    }
    color {selected} @gAltScheme
}

function modifyExitMB() {
    unbind
    set echo TOP LEFT
    echo
    set allowRotateSelected gAllowRotateSelected
    set allowMoveAtoms gAllowMoveAtoms
    set hoverDelay gHoverDelay
    select all
    set hoverLabel = ""
    color {selected} @gScheme
    select {gSelsave}
}
    
# Top level of Modify
function plicoModify() {
    gSelsave = {selected}
    gAllowMoveAtoms = allowMoveAtoms
    gAllowRotateSelected = allowRotateSelected
    gHoverDelay = hoverDelay
    set allowMoveAtoms TRUE    
    set allowRotateSelected TRUE
    
    gPlico = "MODIFY"
    
    gScheme = defaultColorScheme
    gAltScheme = ((gScheme == "Jmol") ? "Rasmol" : "Jmol")
    set echo TOP LEFT
    background ECHO yellow
    
    gEcho = ("________MODIFY________|ALT-CLICK=change 1|SHIFT-CLICK=change 2|" +
        "ALT-SHIFT-CLICK=select chain|ALT-CTRL-CLICK=add water|ALT-DRAG=rotate selected|" +
        "ALT-SHIFT_DRAG=move selected|DOUBLE-CLICK=exit")
    echo @gEcho
    gChain = ""
    unbind
    
    bind "ALT-LEFT-CLICK" "_pickAtom";
    bind "ALT-LEFT-CLICK" "+:modify1MB";
    bind "SHIFT-LEFT-CLICK" "_pickAtom";
    bind "SHIFT-LEFT-CLICK" "+:modify2MB";
    bind "ALT-SHIFT-LEFT-CLICK" "_pickAtom";
    bind "ALT-SHIFT-LEFT-CLICK" "+:selectChainMB";
    bind "ALT-CTRL-LEFT-CLICK" "_pickAtom";
    bind "ALT-CTRL-LEFT-CLICK" "+:addWaterMB";
    bind "DOUBLE" "modifyExitMB";
    
    set hoverDelay 0.001
    hover on
    setHoverLabels()
}

# End of MODIFY.SPT

Contributors

Remig