Molecular modeling and molecular dynamics simulation of the human A2B adenosine receptor. The study of the possible binding modes of the A2B receptor antagonists

A molecular model of the human A(2B) adenosine receptor containing seven transmembrane alpha helices connected by three intracellular and three extracellular hydrophilic loops had been constructed. A molecular docking of seven structurally diverse xanthine antagonists of the A(2B) receptor was performed, and the differences in their binding modes were investigated. The 1 ns molecular dynamics (MD) simulations of several obtained ligand-receptor complexes inserted into the phospholipid bilayer were carried out. The conformational changes of the A(2B) receptor occurring during MD simulations were explored, and the stable binding modes of the studied antagonists were determined. According to the models presented in this work, the involvement of the His251, Asn282, Ser92, Thr89, and some aromatic residues in ligand recognition was determined. The obtained binding modes of the A(2B) antagonists demonstrate good agreement with the site-directed mutagenesis data.