[Chimera-users] Surface not including all molecules?

[Elaine Meng]

Hi Gustaf,
Another possibility is to simulate a density map from the atoms and show an isosurface of that density map.  For example, command

molmap :pyr 2.5 gridSpacing 0.3

would simulate a density map from the atoms in your file and open the Volume Viewer tool for adjusting the density contour level using the vertical bar on the histogram.  The possible disadvantage is that the surface is defined differently than the molecular solvent-excluded surface, and there are adjustable parameters: the contour level, and the "resolution" (2.5 in the example commands, you could try different values) used to simulate the density map. Smaller gridSpacing gives a smoother contour surface.

The molmap command has several additional options:
<http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/molmap.html>

Volume Viewer docs:
<http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/volumeviewer/framevolumeviewer.html>

Best,
Elaine

On Dec 16, 2011, at 12:52 PM, Gustaf Olsson wrote:

> It work, just a little to well. When I try I only get the "outer" surface for all molecules. Either trying to find a way to make the criteria more strict to leave more blank space "within" the sphere. Or divide the PYR residues into "clusters", creating surfaces for "aggregated" groups of molecules!
> 
> But thanks a lot for the startingpoint!
> 
> Best regards
> // Gustaf
> 
> On Dec 16, 2011, at 19:37 , Elaine Meng wrote:
> 
>> Hi Gustaf,
>> Your "structure" appears to be many separate small molecules.  In that case, you would probably have to tell Chimera that they should be enclosed together, since usually people would want different surfaces for the different molecules.  That could be done with the command surfcat.  Assuming you want all the atoms to be lumped together, something like:
>> 
>> surfcat mystuff #
>> surface mystuff
>> 
>> However, that surface calculation is likely to be problematic (there may be numerical failures) because the structure is very open, i.e. loosely packed: there are many holes and cavities if all those atoms are put together in one surface category.
>> 
>> When I tried it, I did get a surface plus a warning message.  One big surface was calculated successfully, but the calculation for other completely enclosed bubbles failed numerically.  The big surface does show many holes and a large interior cavity that is open to the exterior.  If you got that and it shows what you want, you are done.
>> 
>> However, if you needed the additional small bubbles, or you didn't get any surface at all when you tried this, you would need to try some tweaks to avoid the numerical failure.
>> 
>> Unfortunately a numerical failure is specific to everything (molecule, computer, all surface parameters, possibly even molecule orientation) and there is no one recipe to fix it in all situations.  Some things to try are changing the atomic VDW radii slightly and/or the molecular surface vertex density or probe radius.
>> 
>> For example, choose Favorites... Preferences, in Preferences go to Category: New Surfaces, and change settings before trying to calculate a surface.  To change vdw radii, could be something like command:
>> 
>> vdwdefine +.01
>> 
>> or to go back to default radii,
>> 
>> ~vdwdefine
>> 
>> The "split" command isn't relevant here, as it splits chains into separate models, but all the atoms in this structure are apparently chain "X".  If you did have multiple chains, but you wanted all the atoms to be in one surface, you still wouldn't want to split, because that would force calculation of a separate surface for each chain.
>> 
>> I hope this helps,
>> Elaine
>> -----
>> Elaine C. Meng, Ph.D.                       
>> UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab
>> Department of Pharmaceutical Chemistry
>> University of California, San Francisco
>> 
>> On Dec 16, 2011, at 12:54 AM, Gustaf Olsson wrote:
>> 
>>> Hello
>>> I'm trying to create a "cavity" model for visualization, I've done it before (a long time ago) and it worked great. But this time I'm having some problems. When trying with all molecules, I got some error regarding the calculation, so I split the model up into pieces and now it almost works. 
>>> 
>>> The problem I'm having is when I try creating a surface of one of the molecular species, I only get three "blobs" and not one big surface for all molecules. Any suggestions are most welcome.
>>> 
>>> Don't know if it works, but I'll attach the PDB file for the molecules not being displayed as a surface properly.
>>> 
>>> <shell_pyr.pdb>
>>> ______________________________
>>> Gustaf Olsson
>>> p://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users
>> 
> 
> 
> _______________________________________________
> Chimera-users mailing list
> Chimera-users at cgl.ucsf.edu
> http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users