In a classic mode of estrogen action, estrogen receptors (ERs) bind to estrogen responsive element (ERE) to activate gene transcription. A perfect ERE contains a 13-base pair sequence of a palindromic repeat separated by a three-base spacer, 5′-GGTCAnnnTGACC-3′. In addition to the consensus or wild-type ERE (wtERE), naturally occurring EREs often have one or two base pairs’ alternation. Based on the newly constructed Thermodynamic Modeling of ChIP-seq (TherMos) model, binding energy between ER[Formula: see text] and a series of 34-bp mutant EREs (mutERE) was simulated to predict the binding affinity between ERs and EREs with single base pair deviation at different sites of the 13-bp inverted sequence. Experimentally, dual polarization interferometry (DPI) method was developed to measure ER[Formula: see text]–mutEREs binding affinity. On a biotin-NeutrAvidin (NA)-biotin treated DPI chip, wtERE is immobilized. In a direct binding assay, ER[Formula: see text]–wtERE binding affinity is determined. In a competition assay, ER[Formula: see text] was preincubated with mutant EREs before being added for competitive binding to the immobilized wtERE. This competition strategy provided a successful platform to evaluate the binding affinity variation among large number of ERE with different base mutations. The experimental result correlates well with the mathematically predicted binding energy with a Spearman correlation coefficient of 0.97.