
Dear Chimera Developers, First, thank you all for such a wonderful program! Also, thanks to the users' list for such quick, thorough, and always-helpful answers! Second, a little background: I have a gromacs MD simulation (well several, actually) of a small ligand diffusing through a protein. Because the ligand requires the entire simulation to traverse the protein cavity, I was thinking of using the Surfnet tool to accumulate or "build-up" the cavity interface between the ligand and protein using a per-frame python script a la Eric's suggestion in the Dec 2009 [Chimera-users] Surfnet command line thread on this mailing list. I was thinking I could then "prune" the accumulated SurfaceModel using the surfvol.py script available on the http://plato.cgl.ucsf.edu/trac/chimera/wiki/Scripts download page. Now, I can generate a Surfnet surface easily (as Eric points out) for each MD frame using the python commands from Surfnet import interface_surfnet interface_surfnet("LigandAtoms","ProteinAtoms") But how, if possible, can I combine all the individual SurfaceModels into one surface? Also, and not secondarily, does this seem like a reasonable approach to obtain the ligand-protein cavity surface? Thanks again for all your generous support. -- William J. McDonald Postdoctoral Scholar Department of Chemistry and Biochemistry University of California, Santa Cruz