As crucial members of the G-protein coupled receptor (GPCR) superfamily, alpha (1)-adrenergic receptors (alpha(1)-ARs) are recognized to intervene the actions of endogenous catecholamines such as norepinephrine and epinephrine. So far three distinct alpha(1)-AR subtypes, alpha(1A), alpha(1B) and alpha(1D), have been characterized by functional analysis, radio-ligand binding and molecular biology studies. The alpha(1)-ARs are of therapeutic interest because of their distinct and critical roles in many physiological processes, containing hypertension, benign prostatic hyperplasia, smooth muscle contraction, myocardial inotropy and chronotropy, and hepatic glucose metabolism. Accordingly, designing subtype-selective antagonists for each of the three alpha(1)-AR subtypes has been an enthusiastic region of medicinal research. Even though a large number of studies on GPCRs have been conducted, understanding of how known antagonists bind to alpha(1)-ARs still remains sketchy and has been a serious impediment to search for potent and subtype-selective alpha(1)-AR antagonists because of the lack of detailed experimental structural knowledge. This review deliberates the simulation of alpha(1)-ARs and their interactions with antagonists by using ligand-based (pharmacophore identification and QSAR modeling) and structure-based (comparative modeling and molecular docking) approaches. Combined with experimental data, these computational attempts could improve our understanding of the structural basis of antagonist binding and the molecular basis of receptor activation, thus offering a more reasonable approach in the design of drugs targeting alpha(1)-ARs.