QuteMol allows production of nice realistic molecular pictures. PyMOL could do similar things albeit results are not strictly equivalent. Here is an example taken from the PyMOL wiki that produced the picture above.

fetch 1BTA
set_color oxygen, [1.0,0.4,0.4]
set_color nitrogen, [0.5,0.5,1 …