Hi Anders,
With cmd.get_povray() (i.e. cmd.get_povray()[1]) you'll
get all vertex-coordinates describing the surface, but I
wouldn't know how you could use it to map something on
there neatly.
Cheers,
Tsjerk
On Wed, 22 Jun 2005 08:55:05 +0200
madsen <mad...@ccs.ki.ku.dk> wrote:
Dear Robert, Tsjerk, Jean-Philippe and Tim
thanks for your replies.
The colouring of atoms according to B-factors
or some other property is nice, but I really need
sub-atomic resolution of the properties that I would
like to map on the surface.
So, as Jean-Philippe and Tim are discussing, I would
like
to access the defining coordinates of the surface, or by
some
other means be able to hack the colouring of the
surface.
There is apparantly no standard trick for doing just
that.
I will have a look at povscript+.
Thanks for the hint!
regards
Anders
On Tue, Jun 21, 2005 at 02:23:24PM -0400, Robert
Campbell wrote:
Anders,
* Anders Madsen <mad...@ccs.ki.ku.dk> [2005-06-21 16:33]
wrote:
> Dear PyMol users,
> I would like to colour the beautiful
> PyMol molecular surfaces according to
> properties of the atoms involved.
> Has anyone experimented with this?
>
> E.g. by accessing the coordinates that define
> the surface and assign a colour to each triangle
> involved in the surface mesh?
If you only need to colour it at the resolution of the
atoms (as opposed
to individual surface vertices), then you can do so by
altering the
B-factor or occupancy values and colouring the structure
on those. I
have a script for reading a file of data and applying it
to the B-factor
or occupancy. Then you just need to colour using the
built-in spectrum
method or using my color_b.py script. See the
data2bfactor.py script in:
http://adelie.biochem.queensu.ca/~rlc/work/pymol
Within that script, four functions are defined:
data2b_atom
data2b_res
data2q_atom
data2q_res
The "2b" functions alter the B-factor and the "2q"
functions alter the
occupancy. The "_atom" functions are used to apply a
value to
individual atoms, while the "_res" functions are used to
apply a single
value to the whole residue.
You just need to create a data file that looks like:
A 1 ALA N data
A 1 ALA CA data
A 1 ALA C data
A 1 ALA O data
A 1 ALA CB data
etc. (or omit the chain if not defined in your
structure)
or
A 1 ALA data
A 2 SER data
A 3 PRO data
etc. (again omit the chain if not defined).
Cheers,
Rob
--
Robert L. Campbell, Ph.D.
<r...@post.queensu.ca>
Senior Research Associate
phone: 613-533-6821
Dept. of Biochemistry, Queen's University,
fax: 613-533-2497
Kingston, ON K7L 3N6 Canada
http://adelie.biochem.queensu.ca/~rlc
PGP Fingerprint: 9B49 3D3F A489 05DC B35C 8E33 F238
A8F5 F635 C0E2
-------------------------------------------------------
SF.Net email is sponsored by: Discover Easy Linux
Migration Strategies
from IBM. Find simple to follow Roadmaps,
straightforward articles,
informative Webcasts and more! Get everything you need
to get up to
speed, fast.
http://ads.osdn.com/?ad_id=7477&alloc_id=16492&op=click
_______________________________________________
PyMOL-users mailing list
PyMOL-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/pymol-users
--
Anders Østergaard Madsen
Centre for Crystallographic Studies
Department of Chemistry, University of Copenhagen
Universitetsparken 5 DK-2100 KBH Denmark
home: Røde Mellemvej 18, 1. DK-2300 Kbh. S
email: mad...@ccs.ki.ku.dk
-------------------------------------------------------
SF.Net email is sponsored by: Discover Easy Linux
Migration Strategies
from IBM. Find simple to follow Roadmaps,
straightforward articles,
informative Webcasts and more! Get everything you need
to get up to
speed, fast.
http://ads.osdn.com/?ad_idt77&alloc_id492&opÿick
_______________________________________________
PyMOL-users mailing list
PyMOL-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/pymol-users