Dear Elaine,<br><br>Thanks again for the good information. <br>I have done it the way you said and it is really a lot of work, but i was thinking was it the right way or not.<br>When I select aromatic ring, all aromatic-aromatic contacts between both residues and ligand aromatic atoms are displayed.<br>
<br>Best regards,<br><br>--<br>Shahid.<br><br><div class="gmail_quote">On Thu, Oct 29, 2009 at 6:11 PM, Elaine Meng <span dir="ltr"><<a href="mailto:meng@cgl.ucsf.edu">meng@cgl.ucsf.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi Shahid,<br>
For (d) there is not really any trick, it is just using the residue and atom names. For example, if you knew C1-C6 in LIG are an aromatic ring, for aromatic-aromatic contacts to PHE you could look for lines in the contacts output with PHE CG,CD1,CD2,CE1,CE2,CZ such as:<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
LIG 1 C1 PHE 168.A CD2 -0.163 3.488<br>
</blockquote>
<br>
As I mentioned this would be a lot of work since you would have to know all the atom/res names for a particular type of functional group such as aromatic ring. Also, I don't think this would work for your system because looking at what you sent below, your LIG atoms are not named uniquely (all the carbons are just named C, all the nitrogens named N, etc.).<br>
Best,<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<div class="im"><br>
<br>
On Oct 29, 2009, at 4:10 AM, M. Shahid wrote:<br>
<br>
</div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div class="im">
Dear Elaine,<br>
Dear Eric,<br>
<br>
Thanks a lot for suggestion of using an alternate way of classifying the contact interactions.<br>
I am currently using it with step by step selection of different types. And its working great!!! thanks again.<br>
<br>
Could you please also give me some hints on the other way you mentioned in (d) below in conversation?<br>
<br>
Best regards,<br>
<br>
--<br>
Shahid.<br>
<br>
On Mon, Oct 26, 2009 at 9:02 PM, Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu" target="_blank">meng@cgl.ucsf.edu</a>> wrote:<br>
<br></div><div class="im">
(d) yet another approach is to just get all the contacts of all types and filter them afterwards. However, that would require you to figure out the types of interactions based on the residue and atom names, which might be a lot of work. If using this approach, you might want to remove all the hydrogens first, since they would generally provide a large amount of redundant information.<br>
<br>
<br></div><div><div></div><div class="h5">
On Oct 24, 2009, at 8:00 AM, M. Shahid wrote:<br>
<br>
Dear All,<br>
I have a question regarding the contact interactions between a protein-ligand complex.<br>
<br>
I can retrieve the contacts by the findclash command in --nogui mode as below:<br>
chimera --nogui protligcomplex.pdb <a href="http://clash.com" target="_blank">clash.com</a><br>
<br>
and similarly I can find the hbonds.<br>
<br>
The output I am getting is as below:<br>
--------------------------------------------------------------------------------<br>
32 contacts<br>
atom1 atom2 overlap distance<br>
LIG 1 H ALA 265.A O 1.045 1.435<br>
LIG 1 O ALA 265.A O 0.550 2.430<br>
LIG 1 N GLU 169.A OE2 0.439 2.666<br>
LIG 1 H GLU 169.A CD 0.116 2.584<br>
LIG 1 O ASN 253.A 1HD2 0.094 2.386<br>
LIG 1 H ALA 265.A C 0.090 2.610<br>
LIG 1 C ALA 265.A O 0.059 3.121<br>
LIG 1 N ASN 253.A OD1 -0.037 3.142<br>
LIG 1 C ALA 265.A O -0.039 3.219<br>
LIG 1 N GLU 169.A HG3 -0.072 2.697<br>
LIG 1 C PHE 168.A HB2 -0.081 2.781<br>
LIG 1 N GLU 169.A CG -0.113 3.438<br>
LIG 1 N GLU 169.A CD -0.125 3.450<br>
LIG 1 N GLU 169.A CD -0.128 3.453<br>
LIG 1 C MET 270.A HG2 -0.140 2.840<br>
LIG 1 N PHE 168.A CD2 -0.163 3.488<br>
.......... .... ... .. .<br>
</div></div></blockquote>
<br>
</blockquote></div><br>