<div dir="ltr"><br><br><div class="gmail_quote">---------- Forwarded message ----------<br>From: <b class="gmail_sendername">Francesco Pietra</b> <span dir="ltr"><<a href="mailto:chiendarret@gmail.com">chiendarret@gmail.com</a>></span><br>
Date: Tue, Aug 27, 2013 at 10:45 AM<br>Subject: Re: [Chimera-users] Fwd: Make bonds in a large polypeptide<br>To: UCSF Chimera Mailing List <<a href="mailto:chimera-users@cgl.ucsf.edu">chimera-users@cgl.ucsf.edu</a>><br>
<br><br><div dir="ltr"><div><div><div><div><div><div>Hi Elaine:<br></div><br></div>Matching was very good by using MultiSeq plugin of VMD.<br><br></div>The situation before "Join Models" is represented by the couple of strands on the left. Following "Join Models", by the couple of strands at the right. Highlighted in red Ala-Leu before bonding and after bonding. That was done with distances and dihedrals as in the integer subunit from which they come. In my hands, formation of a regular alpha-helix could improve a little on playing with distance/dihedral values, but only little.<br>
<br></div><div>Putting both fragments as a single model, command "make bond' between NH of Leu and CO of Ala refused to work. Therefore, I asked whether 450 residues are too much (for example, they are too much for Molefacture in VMD).<br>
</div><div><br></div>As "Join Models" is certainly a proved tool, I am not submitting any bug. There must be something wrong either with files or doing. Looking for alternatives.<br><br></div>Thanks a lot for your interest.<span class="HOEnZb"><font color="#888888"><br>
<br></font></span></div><span class="HOEnZb"><font color="#888888">francesco<br></font></span></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Aug 26, 2013 at 8:44 PM, 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 Francesco,<br>
Add bond is only for when the parts are already in the same model. Join Models is for when they are in two separate models.<br>
<br>
Join Models doesn't change the existing dihedrals. It just joins the models using the distance and angle values you give it, which includes a new dihedral around the new bond. I don't know what else to say other than to be careful to select the right atoms, as described, and enter the right values:<br>
<<a href="http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/editing/editing.html#join" target="_blank">http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/editing/editing.html#join</a>><br>
<br>
Also if you are joining peptides it should be easier to use the "C-N peptide bond" option, because then you don't have to add hydrogens and decide which to replace (as in the "other bond" option).<br>
<br>
If you are sure it is changing some dihedral that you already have, you could use "Help… Report a Bug" and attach your session before the join (with atoms already selected). In the report please describe all of the values you used in the join, and say which dihedral changes, and its value measured before and after the join.<br>
I hope this helps,<br>
Elaine<br>
----------<br>
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>
<div><div><br>
On Aug 26, 2013, at 11:32 AM, Francesco Pietra <<a href="mailto:chiendarret@gmail.com" target="_blank">chiendarret@gmail.com</a>> wrote:<br>
<br>
> Must add that I had also tried "Join Models" at no avail. Restraing concerned dihedrals, through tmsd-colvars applied to the CA atom, so that the side chains could move. With "Join Models" I was therefore unable to get the continuous alpha helix shape: formation of the new peptide bond resulted in a kinked alpha helix, no matter how I dealt with the "Join Models" entries. What I need is that the new peptide bond is formed without altering the dihedrals based on the CA atoms.<br>
><br>
> francesco pietra<br>
><br>
> ---------- Forwarded message ----------<br>
> From: Francesco Pietra <<a href="mailto:chiendarret@gmail.com" target="_blank">chiendarret@gmail.com</a>><br>
> Date: Mon, Aug 26, 2013 at 8:02 PM<br>
> Subject: Make bonds in a large polypeptide<br>
> To: chimera <<a href="mailto:chimera-users@cgl.ucsf.edu" target="_blank">chimera-users@cgl.ucsf.edu</a>><br>
><br>
><br>
> Hello;<br>
> Is there a limit in the number of atoms when using the function "make a bond"? I was trying that in a polypeptide of 450 standard residues (I made it a single model) and it did not work.<br>
><br>
> Just to not be cryptic, that need arises from having modeled by special techniques the initial portion of the polypeptide. Now I would like to replace that portion with the modeled one. I had restrained with rmsd-colvars the dihedrals of a short stretch of alpha helix, so that I could align the modeled portion with the polypeptide.<br>
><br>
> thanks<br>
> francesco pietra<br>
<br>
</div></div></blockquote></div><br></div>
</div></div></div><br></div>