JG,
I predict you *will not* find a good way in the current state. When the brightness exceeds a certain value, no amount of color saturation will help.
A good rule of thumb is around 0.8 . Above this value it is almost impossible to manage color in a pleasing way. So the issue is how bright the stars are relative to the chrominance information. One technique, the one I typically employ, is to use strong non-linear brightening schemes to manage the brightness of things that are important to color (to make certain they do not exceed that 0.8 brightness value). Another technique is to dim the stars below the impossible bright threshold and then recolor them from the original RGB information that is stretched in a strongly non-linear way (to keep the color). Typically this is done using a star mask. This only works if the stars are not saturated from the get go. If the stars have the high values from the originally acquired images... well...then you really have to get creative. One method is to using convolution and purposely blur the stars a little to even out the flat-topped PSF. Finally, it is possible to use HDRMT to also moderate the star values- but this one is pretty messy and tough to manage for bright stars.
-adam
