Alt click the stopwatch for the opacity property for each star layer. Add this expression:
wiggle(4,100)
Adjust the nubers to taste. first number is frequency per second, second number is amplitude.
--- EDIT ---
I think I misunderstood - reading your message again I now think your 8 large stars are on a single layer. Yes?
In that case - a different approach. Set up layers as follows:
1. Solid layer. Apply turbulent noise effect - adjust settings (contrast, complexity, scale etc) so you have random patches of white and black. Animate the evolution property so the pattern changes over time.
2. Your star layer
Create a Luma track matte between these layers.
Your stars will now vary in opacity over time.