It's easy to do this in PyMOL. Just set ray_trace_mode=2 and it converts to B&W. See more tips and tricks on their wiki. Atom labeling is also very easy in PyMOL.

This looks like a rather old style of drawing. All modern software I am aware of uses by default a more elaborate style, e.g. with shading for a 3D-look. Older software, e.g. ATOMS 3, created this style. However, handling the data and creating the plot has become so much more comfortable in modern software, that it would be a really bad idea to go back to the old stuff.

I think it is possible to get this flat style from VESTA (http://jp-minerals.org/vesta/en/) by disenabling the "lighting"-option in menu point "overall appearance".    

Thanks for the answer Daniel. I will give it a go. 

I'm able to generate some pretty good Tachyon images using VMD for when I am working on MD simulation. But I find the labelling of atoms in VMD is very primitive.  

Here's my best effort , I think the setting that I used can be seen on the picture attached below. Few things to highlight are:

Have fun!

It's easy to do this in PyMOL. Just set ray_trace_mode=2 and it converts to B&W. See more tips and tricks on their wiki. Atom labeling is also very easy in PyMOL.

Two options:

Get a nice colorful VMD representationand then post-process to B&W using the GIMP.

Change options in VMD to get the right colors (Grey, black, white) for your atoms

Either way is straightforward procedure...

The slacker option:  send the pdb file and I will return you the B&W stuff...

:)

Thanks all, 

I tried with Avogadro (1.1.1) just now, was going well until I tried changing the hydrogen bond colour from yellow (default) to black. Unfortunately, on a 100% fresh Mac Yosemite 10.10.1  (two weeks old), it doesn't work! I fired an email off to the developers.

In the meanwhile I will try PyMol and perhaps go back to an attempt with VMD and see if I can actually make the labels look half decent. 

An very late answer but I think that Gabedit will be the best option