Analysis of 300 ns (ns) molecular dynamics (MD) simulations of an adenosine A2a receptor (A2a AR) model, conducted in triplicate, in 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) and 1-palmitoyl-2-oleoylphosphatidylethanolamine (POPE) bilayers reveals significantly different protein dynamical behavior. Principal component analysis (PCA) shows that the dissimilarities stem from interhelical rather than intrahelical motions. The difference in the hydrophobic thicknesses of these simulated lipid bilayers is potentially a significant reason for the observed difference in results. The distinct lipid headgroups might also lead to different molecular interactions and hence different protein loop motions. Overall, the A2a AR shows higher mobility and flexibility in POPC as compared to POPE.
Molecular dynamics (MD) simulations of membrane-embedded G-protein coupled receptors (GPCRs) have rapidly gained popularity among the molecular simulation community in recent years, a trend which has an obvious link to the tremendous pharmaceutical importance of this group of receptors and the increasing availability of crystal structures. In view of the widespread use of this technique, it is of fundamental importance to ensure the reliability and robustness of the methodologies so they yield valid results and enable sufficiently accurate predictions to be made. In this work, 200 ns simulations of the A2a adenosine receptor (A2a AR) have been produced and evaluated in the light of these requirements. The conformational dynamics of the target protein, as obtained from replicate simulations in both the presence and absence of an inverse agonist ligand (ZM241385), have been investigated and compared using principal component analysis (PCA). Results show that, on this time scale, convergence of the replicates is not readily evident and dependent on the types of the protein motions considered. Thus rates of inter- as opposed to intrahelical relaxation and sampling can be different. When studied individually, we find that helices III and IV have noticeably greater stability than helices I, II, V, VI, and VII in the apo form. The addition of the inverse agonist ligand greatly improves the stability of all helices.
Human A3 adenosine receptor hA3AR has been implicated in gastrointestinal cancer, where its cellular expression has been found increased, thus suggesting its potential as a molecular target for novel anticancer compounds. Observation made in our previous work indicated the importance of the carbonyl group of amide in the indolylpyrimidylpiperazine (IPP) for its human A2A adenosine receptor (hA2AAR) subtype binding selectivity over the other AR subtypes. Taking this observation into account, we structurally modified an indolylpyrimidylpiperazine (IPP) scaffold, 1 (a non-selective adenosine receptors' ligand) into a modified IPP (mIPP) scaffold by switching the position of the carbonyl group, resulting in the formation of both ketone and tertiary amine groups in the new scaffold. Results showed that such modification diminished the A2A activity and instead conferred hA3AR agonistic activity. Among the new mIPP derivatives (3-6), compound 4 showed potential as a hA3AR partial agonist, with an Emax of 30% and EC50 of 2.89 ± 0.55 μM. In the cytotoxicity assays, compound 4 also exhibited higher cytotoxicity against both colorectal and liver cancer cells as compared to normal cells. Overall, this new series of compounds provide a promising starting point for further development of potent and selective hA3AR partial agonists for the treatment of gastrointestinal cancers.