A forward modeling approach to AGN variability-method description and early applications