<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
Dear Phenix Community,</div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
<br>
</div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
I'm at beginner level of using Phenix. I have difficulty to refine a protein-RNA complex in Phenix.</div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
<br>
</div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
I refined the protein model without RNA first, then I built the RNA model into the difference density map in Coot through Calculate&gt;Other Modelling Tools&gt;Ideal DNA/RNA, where I input the sequence and fit nucleotides into the density. At the 5' end it's guanosine
 triphosphate (GTP) but the modeler tool in Coot only generates nucleoside monophosphate.&nbsp;<span style="margin: 0px">So I built the RNA model from the second nucleotide at 5' end downstream until where RNA density was too poor to model. I did the real space
 refine zone in Coot to fit this main part of the RNA into its density the best possibly.</span></div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
<br>
</div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
For the GTP at 5' end, I realized there is GTP.cif restraint file in Phenix monomer library. So I obtained an arbitrary pdb file for GTP using eLBOW and corrected atom names matching the definition in GTP.cif file. Then I can also fit GTP into the RNA density
 at 5' end using&nbsp;<span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important">real space refine zone in Coot. Now I have a broken RNA model with GTP unlinked to the&nbsp;<span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important">ribose-</span>phosphate
 backbone.&nbsp;</span></div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
<span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important"><br>
</span></div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
For Phenix refinement, I supplied a&nbsp;cif_link_GTP-G.params file to define the linkage for the two RNA chains in the complex:</div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
<br>
</div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
<span style="margin: 0px">refinement.pdb_interpretation {<br>
</span>
<div style="margin: 0px">&nbsp; apply_cif_link {<br>
</div>
<div style="margin: 0px">&nbsp; &nbsp; data_link = p<br>
</div>
<div style="margin: 0px">&nbsp; &nbsp; residue_selection_1 = chain C and resname GTP and resid 1<br>
</div>
<div style="margin: 0px">&nbsp; &nbsp; residue_selection_2 = chain C and resname G and resid 2<br>
</div>
<div style="margin: 0px">&nbsp; }<br>
</div>
<div style="margin: 0px">&nbsp; apply_cif_link {<br>
</div>
<div style="margin: 0px">&nbsp; &nbsp; data_link = p<br>
</div>
<div style="margin: 0px">&nbsp; &nbsp; residue_selection_1 = chain D and resname GTP and resid 1<br>
</div>
<div style="margin: 0px">&nbsp; &nbsp; residue_selection_2 = chain D and resname G and resid 2<br>
</div>
<div style="margin: 0px">&nbsp; }<br>
</div>
<span style="margin: 0px">}</span><br>
</div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
<span style="margin: 0px"><br>
</span></div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
<span style="margin: 0px">After refinement, I noticed O3' of GTP appears connected to P of the following nucleotide G in the refined model in Coot. However, if I do the real space refine zone in Coot, the linkage becomes broken, i.e., I cannot refine the RNA
 model as a whole: Coot only allows separate real space refine for GTP and the rest. It seems Coot still cannot recognize GTP and the rest as one body. In the refinement settings there is a button for Automatic linking options.<span>&nbsp;</span><span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important"><span style="margin: 0px">&nbsp;</span>Interestingly,<span style="margin: 0px">&nbsp;</span></span>I
 noticed if I leave the automatic covalent linking as default, Phenix writes Link records in the output pdb after refinement. But it also links one of the protein residues to the GTP in chain D. In order to avoid the unwanted intermolecular linkage, I clicked
 Link none for<span>&nbsp;</span><span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important">Automatic linking options</span>, then there is no any link records written in the output pdb despite of the defined linkage in&nbsp;<span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important">cif_link_GTP-G.params
 file. So it seems for me Phenix ignores the input&nbsp;<span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important">cif_link_GTP-G.params file, it writes Link records only if the Automatic linking option is activated. I worried whether
 it could be due to the syntax of the&nbsp;<span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important">cif_link_GTP-G.params file. So I tried the following syntax in two separate files, the same issue was observed.</span></span></span></span></div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
<span style="margin: 0px"><span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important"><span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important"><span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important"><br>
</span></span></span></span></div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
<span style="margin: 0px"><span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important"><span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important"><span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important">cif_link_GTP-G_chC.params<br>
</span></span></span></span></div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
<span style="margin: 0px"><span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important"><span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important"><span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important"><span style="margin: 0px">refinement.pdb_interpretation.apply_cif_link
 {<br>
</span>
<div style="margin: 0px">&nbsp; &nbsp; data_link = p<br>
</div>
<div style="margin: 0px">&nbsp; &nbsp; residue_selection_1 = chain C and resname GTP and resid 1<br>
</div>
<div style="margin: 0px">&nbsp; &nbsp; residue_selection_2 = chain C and resname G and resid 2<br>
</div>
<span style="margin: 0px">}</span><br>
</span></span></span></span></div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
<span style="margin: 0px"><span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important"><span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important"><span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important"><span style="margin: 0px"><br>
</span></span></span></span></span></div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
<span style="margin: 0px"><span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important"><span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important"><span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important"><span style="margin: 0px">
<div style="margin: 0px"><span style="margin: 0px"><span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important"><span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important"><span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important">cif_link_GTP-G_chD.params<br>
</span></span></span></span></div>
<div style="margin: 0px"><span style="margin: 0px"><span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important"><span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important"><span style="margin: 0px; background-color: rgb(255, 255, 255); display: inline !important"><span style="margin: 0px">refinement.pdb_interpretation.apply_cif_link
 {<br>
</span>
<div style="margin: 0px">&nbsp; &nbsp; data_link = p<br>
</div>
<div style="margin: 0px">&nbsp; &nbsp; residue_selection_1 = chain D and resname GTP and resid 1<br>
</div>
<div style="margin: 0px">&nbsp; &nbsp; residue_selection_2 = chain D and resname G and resid 2<br>
</div>
<span style="margin: 0px">}</span></span></span></span></span></div>
<br>
</span></span></span></span></span></div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
I also need to point out that p linkage is pre-defined in the list of links (mon_lib_list.cif) in Phenix monomer library.<br>
</div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
<br>
</div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
I apologize for this long explanation of the problem. I hope it is clear.</div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
<br>
</div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
I really appreciate your advice in advance.</div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
<br>
</div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
Sincerely,</div>
<div style="margin: 0px; font-size: 12pt; font-family: Calibri, Arial, Helvetica, sans-serif; background-color: rgb(255, 255, 255)">
Zixian Li</div>
</div>
</body>
</html>