<div dir="ltr">Thanks so much Elaine! I finally got this to work. It took me a while to figure out how to text edit pdb files.<div><br></div><div>If anyone on the forum is interested, here is what I did: </div><div>
<br></div><div>1. Opened the original pdb file in text edit containing the receptor and two ligands (3 chains total).</div><div>2. Inserted a linker sequence between the two ligands with x,y,z coordinates set to zero and occupation set to -1. </div>

<div>3. Used a python script to renumber all the atoms in the correct order give both ligands and the linker the same chain ID. At this point, I had ligand 1 - linker - ligand two, all in the same chain with atoms in sequential order. I also deleted all the secondary structure information (HELIX and SHEET lines at the beginning of the pdb file) and manually added these back in later. Pymol will add secondary structure back in automatically, but it seems Chimera has trouble with this sometimes. </div>

<div>4. Saved the file in .pdb format, and opened in Chimera.</div><div>5. Opened a fasta file with the amino acid sequence of my ligand1-linker-ligand2 construct. </div><div>6. Highlighted the linker sequence in the sequence file and opened Modeller (loops/refinement), selected "chimera selection region", "models to generate" = 5, "apply". </div>

<div><br></div><div>This was all done with the receptor in place. Of the 5 generated models, 3 contained linkers that did not clash with the receptor and had acceptable ramachandran plots. </div><div><br>
</div><div>Thanks again for all your help!</div><div><br></div><div>Best,</div><div><br></div><div>Gerald Maxwell Cherf</div><div>PhD candidate</div><div>Stanford University, Bioengineering Dept. </div>
<div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Aug 12, 2013 at 11:44 AM, Elaine Meng <span dir="ltr"><<a href="mailto:meng@cgl.ucsf.edu" target="_blank">meng@cgl.ucsf.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Gerald,<br>
You don't need to do any linker-building, you can use Modeller to fill it in for you, keeping the ends (ligands) in their input positions.  There are probably a few alternative routes, but what I would do is:<br>
<br>
(1) Save PDB containing the peptide ends (the two ligands) in the positions that you want.  What you want at the end is just these two ligands with the same chain ID, and no other atoms.  That can be achieved either by text-editing the PDB after writing it out, or by changing chain IDs if necessary (Tools… Structure Editing… Change Chain IDs) of the ligands before writing it out.  The PDB save dialog gives you the option of saving only selected atoms, so you could select the ligands only and use that option.<br>

<br>
 ** Here's probably the biggest caveat of the approach, which I realized after I sent the last message:  via the Chimera-Modeller interface it is not possible to model one chain while taking other additional chains into account.  In other words, it won't help to include the protein in the file because it will be ignored.  Modeller itself is capable of doing that but the Chimera-Modeller interface only accesses a subset of Modeller's many capabilities… you'd have to run Modeller yourself independent of Chimera, which is much more difficult.  However, all may not be lost: the ensemble of possible solutions (linker conformations) may include some that happen to avoid intersecting with the protein.<br>

<br>
(2) Create plain text sequence file (FASTA, named something.fa) which contains the sequence of entire peptide: the ligand ends and the missing linker portion.  I would just create it in a text editor.  FASTA format is very simple, basically one line starting with > followed by a comment or description (or nothing, whatever you want) followed by one or more lines containing the amino acid sequence.<br>

<br>
(3) In Chimera, open the new PDB file containing the two ligands, and open the FASTA sequence file.  Then from the sequence window menu, choose Structure… Modeller (loops/refinement).  In the resulting dialog, you will want to model "non-terminal missing structure" or "all missing structure" (equivalent in this case).  You will probably aso want to change "Allow this many residues adjacent to missing regions to move" to 0 and increase the number of models to generate … I'd try something much higher like 20 or 50 in hopes of getting something that doesn't intersect with the protein.<br>

<<a href="http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/multalignviewer/modeller.html#building" target="_blank">http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/multalignviewer/modeller.html#building</a>><br>

<div class="im"><br>
I hope this helps,<br>
Elaine<br>
</div>----------<br>
<div class="im HOEnZb">Elaine C. Meng, Ph.D.<br>
UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab<br>
Department of Pharmaceutical Chemistry<br>
University of California, San Francisco<br>
<br>
</div><div class="HOEnZb"><div class="h5">On Aug 9, 2013, at 7:19 PM, Maxwell Cherf <<a href="mailto:maxcherf@gmail.com">maxcherf@gmail.com</a>> wrote:<br>
<br>
> Hi Elaine,<br>
><br>
> Using modeller sounds like the exact solution I'm looking for. Both the ligands and the linker I will connect them with are peptides.<br>
><br>
> I tried using modeller, but I haven't gotten it to work yet. Here's what I've done:<br>
><br>
> 1. built peptide linker using the "build structure" function<br>
> 2. attached the linker to one ligand1 using the "join" function<br>
> 3. combined the ligand1-linker with the second ligand into one model using the "combine" function<br>
> 4. added an unrealistic ~40 angstrom bond between the linker and second ligand (now we have ligand1 - linker - 40A bond - ligand 2, all in one model).<br>
> 5. I then tried using the "model/refine loop" function to find acceptable orientations for the linker (including the 40A bond) I incorporated, but I was unable to do this because the two ligands were still separate chains, and the "model/refine loop" function requires them to be in the same chain. If I use the "join" command instead of the "combine" command to attach the ligand1-linker to ligand2, I get one chain and then the model/refine loop function works perfectly to remodel the linker. However, using the join function changes the orientation of the ligands relative to the receptor, and I need the ligands to remain docked properly to the receptor.<br>

><br>
> Is there a way to join both ligands into the same chain without changing their orientation? Or, do you know of an alternate approach for adding a linker between the two ligands without changing their orientation using modeller?<br>

><br>
> Thanks again,<br>
> Gerald<br>
><br>
><br>
<br>
</div></div></blockquote></div><br></div>