Modeling Stochasticity and Variability in Gene Regulatory Networks