Homology modeling, molecular docking and molecular dynamics (MD) simulation methods were used to build a reliable model for A 1 AR (as one of the G protein-coupled receptors-GPCRs) and to explore the structural features and binding mechanism of ligands to this receptor. A model of A 1 AR was built and inserted in a hydrated lipid bilayer, and 20-ns MD simulation was performed to examine the stability of the best model. In this study, RG-14718 as the best A 1 AR agonist and bamifylline as a selective antagonist of A 1 AR have been docked into the active site of the A 1 AR. After docking, two 20-ns MD simulation was performed on the A 1 ARligand complex to explore effects of the presence of lipid membrane in the vicinity of the A 1 AR-ligand complex. At the end of the MD simulation, a change in the position and orientation of the ligand in the binding site was observed. This important observation indicated that the application of MD simulation after docking of ligands is useful. Thr270, His278 and Asn70 were crucial residues for hydrogen bonds with these ligands. Phe171, Glu172, Tyr271 and Ile274 were determined to be involved in ligand-receptor binding. The results obtained are in good agreement with most of the site-directed mutagenesis data reported by others. Our results show that molecular modeling and rational drug design for adenosine targets is a possible approach.