zap;
load li.pdb;
wireframe 5;
select all;
spacefill 30%;

connect 3.05 (selected)(selected) create; 
select all; color [225,225,191]; color bonds [180,180,180]; select none;

select atomno=17; color [225,0,0];

select atomno=2, atomno=19, atomno=1, atomno=12, atomno=16, atomno=4, atomno=7, atomno=13; color [100,100,100]; select none;

select atomno=18, atomno=15, atomno=25, atomno=11, atomno=6, atomno=21;
color [225,225,0]; select none;