<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hi Shana,<div><span class="Apple-tab-span" style="white-space:pre">        </span>If you were starting from scratch, the procedure to add the CN would probably be simpler by using chimera.molEdit.addDihedralAtom() to add the C and then the N at the appropriate bond length, angle and dihedral and then rename the residue with: r.type = "MSCN". &nbsp;Nonetheless, you have come a long way down the road with your script and probably would not prefer to start over from scratch! &nbsp;So let's suppose you have your two relevant Residue instances in variables named 'cys' and 'cn'. &nbsp;The code to do what you want would be:</div><div><br></div><div>for a in cn.atoms:</div><div><span class="Apple-tab-span" style="white-space:pre">        </span>cn.removeAtom(a)</div><div><span class="Apple-tab-span" style="white-space:pre">        </span>cys.addAtom(a)</div><div>cys.type = "MSCN"</div><div># empty residues are problematic in various parts of Chimera, so...</div><div>cn.molecule.deleteResidue(cn)</div><div><br></div><div><span class="Apple-tab-span" style="white-space:pre">        </span>If you don't know how to get the 'cys' and 'cn' variables but you know an atom specifier that selects each, you can use this code to set the variables:</div><div><br></div><div>runCommand("sel <i>atom-spec-for-cys</i>")</div><div>cys = chimera.selection.currentResidues()[0]</div><div>runCommand("sel <i>atom-spec-for-cn</i>")</div><div>cn = chimera.selection.currentResidues()[0]</div><div><br></div><div><span class="Apple-tab-span" style="white-space:pre">        </span>There's a couple of further relevant things here. &nbsp;Strictly speaking, 4-letter residue names are not legal in PDB format (while there is a column where a 4th character could go, according to the format definition that column is supposed to be blank). &nbsp;Therefore if you export your modified structure to other software programs, they may or may not be able to handle the 4-character residue name without problems. &nbsp;Therefore I recommend you use the same residue name that the PDB uses for this moiety, 4CY (which can be found in the 3or0 entry).</div><div><span class="Apple-tab-span" style="white-space:pre">        </span>Along the same lines, 4CY can actually be found in the SwissSidechain rotamer library plugin for Chimera (<a href="http://www.swisssidechain.ch/visualization/chimera.html">SwissSidechain - A database of non-natural side chains</a>) and therefore if you install that plugin then the entire substitution process becomes a single command: &nbsp;swapnaa 4CY <i>cysteine-atom-spec</i>.</div><div><br></div><div>--Eric</div><div><br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;&nbsp;Eric Pettersen<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;&nbsp;UCSF Computer Graphics Lab<br>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;&nbsp;<a href="http://www.cgl.ucsf.edu">http://www.cgl.ucsf.edu</a><br><br></div><div><div><div>On Jun 18, 2014, at 6:50 AM, Shana Burstein &lt;<a href="mailto:sburstei@haverford.edu">sburstei@haverford.edu</a>&gt; wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div dir="ltr">Hello,<div>I am an undergraduate student at Haverford College and I am attempting to attach thiocyanate probes onto different sites along a protein by first replacing the existing amino acid with a cysteine and then by adding a C and N onto the S- of the Cys residue. I know it is possible to carry this out manually using Tools&gt;Structure Editing&gt;Modify Structure, but I am hoping to automate this process by calling an external .py script that will primarily contain a list of preexisting Chimera commands.&nbsp;</div>

<div><br></div><div>I have been successful at converting the targeted amino acid to Cys using "swapaa," building a CN residue called MSCN (using the createMolecule.py as a guide), combining the two models into one ("combine selected close True") and forming a bond between the Cys S and CN C atoms ("bond"). </div>

<div><br></div><div>At this point, I am seeking help in automating the renaming of the modified residue to MSCN. Currently, Chimera sees this residue as two different separate residues called CYS and MSCN. Is there a way to combine these two structures so that the label corresponds to MSCN and the residue number corresponds to CYS?&nbsp;</div>

<div><br></div><div>Sincerely,</div><div><br></div><div>Shana Burstein</div><div>Haverford College '15</div></div>
_______________________________________________<br>Chimera-users mailing list<br><a href="mailto:Chimera-users@cgl.ucsf.edu">Chimera-users@cgl.ucsf.edu</a><br>http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users<br></blockquote></div><br></div><br><br></body></html>