[Chimera-users] Analysis of the protein-protein binding interphaces in MARTINI models

[Elaine Meng]

Just try the Chimera findclash command and figure out what options you want.  When you have one command that works the way you want, you can use bash to make a plain text file containing many such commands (Chimera command script file, *.com).  The command would be something like 

findclash :.b test other

… but you would have to give additional options with the values you want to use (overlapCutoff, hbondAllowance etc.).  See the documentation:
<http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/findclash.html>
<http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/indexcommand.html#cmdfile>

We can’t do the exact script for you because everybody’s situation is slightly different.  You have actually try entering some commands.
Elaine


> On Dec 5, 2017, at 1:24 AM, James Starlight <jmsstarlight at gmail.com> wrote:
> 
> Thank you so much for the suggestions, Elaine!
> 
> In fact I have assigned chains (I guess as chain IDs) using Chimera's
> GUI interface, which I have mentioned in the another topic.
> For the simplicity I would use Chimera's scripting interface for the
> batch mode (will be grateful for the example of the script because I
> am also very bad in the programming!) e.g to check the contacts
> between chain A and "all other atoms" in only one pdb  file (not a
> trajectory).
> 
> BTW, via bash I can produce several such scripts (for each of the
> chains) and execute them in a loop manner, doesn't it?
> e.g to check contacts between chain B and "all other atoms".
>                                   between chain C and "all other atoms" etc ..
> 
> 
> James
> 
> 2017-12-04 18:57 GMT+01:00 Elaine Meng <meng at cgl.ucsf.edu>:
>> Hi James,
>> I don’t know how you have assigned these letters (as aliases? chain IDs? residue names?) but first figure out the “findclash” command that does the same thing you were doing with the Find Clashes/Contacts GUI.
>> <http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/findclash.html>
>> 
>> There are two main kinds of Chimera scripts, Chimera commands and Python.  Python can have Chimera commands embedded, and is generally needed if you want to substitute arguments into the commands or loop through many variations of the command or a large set of PDB files.  Python is also needed for actions for which no Chimera command is available.
>> 
>> I’m not a programmer myself so I can only point to information others have generated… there’s always the programmer’s guide:
>> <http://www.rbvi.ucsf.edu/chimera/docs/ProgrammersGuide/index.html>
>> 
>> In this case after figuring out the findclash command options you want, I guess you could write a python script to substitute in the letters into the command to loop through the monomers.  There is a page on how to loop through multiple different structures, which is different, but may still be helpful:
>> <https://www.cgl.ucsf.edu/chimera/docs/ProgrammersGuide/basicPrimer.html>
>> 
>> There are lot of example Chimera python scripts for miscellaneous tasks here:
>> <http://plato.cgl.ucsf.edu/trac/chimera/wiki/Scripts>
>> 
>> After you figured out a script for an individual PDB, you could adapt it to specify a per-frame script in MD Movie for your trajectory.
>> <http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/movie/framemovie.html>
>> <http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/movie/movie.html#per-frame>
>> 
>> I hope this helps,
>> Elaine
>> ----------
>> Elaine C. Meng, Ph.D.
>> UCSF Chimera(X) team
>> Department of Pharmaceutical Chemistry
>> University of California, San Francisco
>> 
>>> On Dec 4, 2017, at 5:37 AM, James Starlight <jmsstarlight at gmail.com> wrote:
>>> 
>>> One question regarding Find Clashes/Contacts.
>>> 
>>> assuming that I have assigned the unique letters for each of the monomer
>>> in my CG system. Now I am using this selection ("the individual
>>> chain" vs "all other systems") to calculate contacts using -2.2 cutoff
>>> for MARTINI (operating with the backbone atoms only). Is it possible
>>> to run the same routine in batch mode i) for the individual pdb as
>>> input as well as ii) for the md trajectory (to calculate occurence of
>>> the contacts along the trajectory)? I would be very grateful for some
>>> example script file for batch Chimera!
>>> 
>>> James
>> 
> 
> _______________________________________________
> Chimera-users mailing list: Chimera-users at cgl.ucsf.edu
> Manage subscription: http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users
>