Estrogenic chemicals in the aquatic environment have been shown to cause a variety of reproductive anomalies in fish including full sex reversal, intersex, and altered population sex ratios. Two estrogens found in the aquatic environment, 17α-ethinylestradiol (EE2) and 17β-estradiol (E2), have been measured in wastewater treatment effluents and have been shown to cause adverse effects in fish. To further our understanding of how estrogen exposure affects reproductive endpoints in the male fathead minnow (FHM, Pimephales promelas), a physiologically based computational model was developed of the hypothalamic-pituitary-gonadal (HPG) axis. Apical reproductive endpoints in the model include plasma steroid hormone and vitellogenin concentrations. Using Markov chain Monte Carlo simulation, the model was calibrated with data from unexposed FHM, and FHM exposed to EE2 and E2. Independent experimental data sets were used to evaluate model predictions. We found good agreement between our model predictions and a variety of measured reproductive endpoints, although the model underpredicts unexposed FHM reproductive endpoint variances, and overpredicts variances in estrogen-exposed FHM. We conclude that this model provides a robust representation of the HPG axis in male FHM.