Estimating genome-wide regulatory activity from multi-omics data sets using mathematical optimization