<html><head><meta http-equiv="Content-Type" content="text/html charset=windows-1252"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hi Elaine,<div class=""><br class=""></div><div class="">As always “quick on the ball” and extremely helpful.</div><div class=""><br class=""></div><div class="">I followed your suggestions for the first snapshot of each trajectory. </div><div class=""><br class=""></div><div class="">I got the following Chimera Warning <font color="#ff2600" class="">“Calculation of some surface components failed. Using only single surface component. This may give inaccurate areas if surfaces of either set of atoms or the combined set are disconnected”.</font> </div><div class=""><br class=""></div><div class="">How important is this?</div><div class=""><br class=""></div><div class="">Other than that, from the reply log, it seems that the protein and ligand are recognised. Here’s a summary of the log:</div><div class=""><br class=""></div><div class="">1. Monomer. Command: measure buriedArea :126 protein</div><div class=""><div class=""><br class=""></div><div class="">Assigning chain ID A to 125 residues, e.g. ASP <font color="#0433ff" class=""># correct. The monomer has 125 residues</font></div><div class="">Model 0 (prod_100ns.nc) appears to be a protein without secondary structure assignments.</div></div><div class=""><br class=""></div><div class="">Several WARNINGS</div><div class=""><br class=""></div><div class=""><div class="">Calculation of some surface components failed. Using only single surface component. This may give inaccurate areas if surfaces of either set of atoms or the combined set are disconnected.</div><div class=""><br class=""></div><div class="">Buried solvent accessible surface area</div><div class=""> B1SAS = 324.016, B2SAS = 223.111, BaveSAS = 273.564</div><div class=""> (A1 = 401.975, A2 = 8118.7, A12 = 7973.55 = 77.959 + 7895.59)</div><div class="">Buried solvent excluded surface area</div><div class=""> B1SES = 96.114, B2SES = 177.797, BaveSES = 136.955</div><div class=""> (A1 = 207.413, A2 = 7117.15, A12 = 7050.65 = 111.299 + 6939.35)</div></div><div class=""><br class=""></div><div class="">2. Dimer. Command: measure buriedArea :126 protein</div><div class=""><br class=""></div><div class=""><div class="">Assigning chain ID A to 125 residues, e.g. ASP</div><div class="">Assigning chain ID B to 125 residues, e.g. ASP</div></div><div class=""><font color="#0433ff" class=""># above assignment of chains is correct</font></div><div class=""><br class=""></div><div class="">Several WARNINGS</div><div class=""><br class=""></div><div class=""><div class="">Buried solvent accessible surface area</div><div class=""> B1SAS = 395.601, B2SAS = 202.981, BaveSAS = 299.291</div><div class=""> (A1 = 405.778, A2 = 14917.3, A12 = 14724.5 = 10.1763 + 14714.3)</div><div class="">Buried solvent excluded surface area</div><div class=""> B1SES = 169.315, B2SES = 300.475, BaveSES = 234.895</div><div class=""> (A1 = 208.956, A2 = 13673.7, A12 = 13412.9 = 39.6405 + 13373.3)</div></div><div class=""><br class=""></div><div class="">If I read the above correctly, the buried solvent accessible area for the ligand in the monomer complex is ~ 78 angstrom and for the ligand in the dimer complex ~ 10 angstrom.</div><div class=""><br class=""></div><div class="">Many many thanks for your kind help</div><div class=""><br class=""></div><div class="">George</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On 9 Nov 2016, at 22:40, Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu" class="">meng@cgl.ucsf.edu</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div class="">Hi George,<br class="">That earlier post is probably not that relevant to your situation. I’ll try to outline what you might try. First you would need to figure out command-line specifications that work on your trajectories, for:<br class=""><br class="">(1) the ligand - sometimes “ligand” will work, but not always, as that might also include other stuff you don’t want like sulfate ions. Just check with “sel ligand” and see if that is the only thing that gets selected. Otherwise you might need to use its residue name or number, e.g. :HEM or :128<br class=""><br class="">(2) the protein monomer in one trajectory - “protein” might work but you should check to make sure that it selects everything you want and nothing else. If not, you may need to use a residue number range instead.<br class=""><br class="">(3) the protein dimer in the other trajectory - again “protein” might work<br class=""><br class="">Once you’ve figured out those things, let’s call them spec1, spec2, and spec3, you can try out the buriedArea command on some single frame of the trajectory. For the first trajectory, command:<br class=""><br class="">measure buriedArea spec1 spec2<br class="">… for example:<br class="">measure buriedArea ligand protein<br class=""><br class=""><<a href="http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/measure.html#buriedArea" class="">http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/measure.html#buriedArea</a>><br class=""><br class="">… and for the second trajectory, measure buriedArea spec1 spec3. If the command works (you see plausible results in Reply Log), then you could put the command in a per-frame script in MD Movie to have it run at every frame. The results in the Reply Log will include the area of the two things together and each individually, as well as the resulting buried area, and two sets of results, one for SES (solvent-excluded surface) and one for SAS (solvent-accessible surface). You may have used a per-frame script before for H-bond calculations. <br class=""><br class=""><<a href="http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/movie/movie.html#per-frame" class="">http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/movie/movie.html#per-frame</a>><br class=""><br class="">One problem is that sometimes the surface calculation may fail, especially during dynamics. You may not be able to get results for every frame, so it is probably good to have the script echo frame number into the Reply Log so you can be sure about which results you are getting. In other words, the per-frame command script could include a command like:<br class=""><br class="">echo frame # <FRAME><br class=""><br class="">I hope this helps,<br class="">Elaine<br class="">-----<br class="">Elaine C. Meng, Ph.D. <br class="">UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab<br class="">Department of Pharmaceutical Chemistry<br class="">University of California, San Francisco<br class=""><br class=""><br class=""><br class="">On Nov 9, 2016, at 12:10 PM, George Tzotzos <<a href="mailto:gtzotzos@me.com" class="">gtzotzos@me.com</a>> wrote:<br class=""><br class=""><blockquote type="cite" class="">I’m dealing with two separate trajectories involving:<br class=""><br class="">1. Ligand A in complex with the monomeric subunit of a given protein A<br class="">2. the same ligand in complex with the dimer of protein A<br class=""><br class="">I would like to measure the buried (or solvent accessible) surface area of ligand A in both cases.<br class=""><br class="">I referred to a relevant posting in the list (<a href="http://www.cgl.ucsf.edu/pipermail/chimera-users/2016-June/012472.html" class="">http://www.cgl.ucsf.edu/pipermail/chimera-users/2016-June/012472.html</a>) but I’m not sure how this could apply in the above case. <br class=""><br class="">I’d be grateful for any suggestion on how to deal with this type of measurement.<br class=""><br class="">Many thanks in advance for any suggestions<br class=""><br class="">George<br class=""><br class=""><br class="">_______________________________________________<br class="">Chimera-users mailing list: <a href="mailto:Chimera-users@cgl.ucsf.edu" class="">Chimera-users@cgl.ucsf.edu</a><br class="">Manage subscription: <a href="http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users" class="">http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a><br class=""><br class=""></blockquote><br class=""></div></div></blockquote></div><br class=""></div></body></html>