Cross your genomic regions against the ReMap catalogue of transcription factor binding sites. You can annotate any BED file and calculate statistical enrichments in TFs in your regions based on their binding sites.