Allosteric modulation of adenosine A1 receptors (A1ARs) offers a novel therapeutic approach for the treatment of numerous central and peripheral disorders; however, despite decades of research, there is a relative paucity of structural information regarding the A1AR allosteric site and mechanisms governing cooperativity with orthosteric ligands. We combined alanine-scanning mutagenesis of the A1AR second extracellular loop (ECL2) with radioligand binding and functional interaction assays to quantify effects on allosteric ligand affinity, cooperativity, and efficacy. Docking and molecular dynamics (MD) simulations were performed using an A1AR homology model based on an agonist-bound A2AAR structure. Substitution of E172ECL2 for alanine reduced the affinity of the allosteric modulators PD81723 and VCP171 for the unoccupied A1AR. Residues involved in cooperativity with the orthosteric agonist NECA were different in PD81723 and VCP171; positive cooperativity between PD81723 and NECA was reduced on alanine substitution of a number of ECL2 residues, including E170ECL2 and K173ECL2, whereas mutation of W146ECL2 and W156ECL2 decreased VCP171 cooperativity with NECA. Molecular modeling localized a likely allosteric pocket for both modulators to an extracellular vestibule that overlaps with a region used by orthosteric ligands as they transit into the canonical A1AR orthosteric site. MD simulations confirmed a key interaction between E172ECL2 and both modulators. Bound PD81723 is flanked by another residue, E170ECL2, which forms hydrogen bonds with adjacent K168ECL2 and K173ECL2. Collectively, our data suggest E172ECL2 is a key allosteric ligand-binding determinant, whereas hydrogen-bonding networks within the extracellular vestibule may facilitate the transmission of cooperativity between orthosteric and allosteric sites.

The adenosine A G protein-coupled receptor (AAR) is an important therapeutic target implicated in a wide range of cardiovascular and neuronal disorders. Although it is well established that the AAR orthosteric site is located within the receptor's transmembrane (TM) bundle, prior studies have implicated extracellular loop 2 (ECL2) as having a significant role in contributing to orthosteric ligand affinity and signaling for various G protein-coupled receptors (GPCRs). We thus performed extensive alanine scanning mutagenesis of AAR-ECL2 to explore the role of this domain on AAR orthosteric ligand pharmacology. Using quantitative analytical approaches and molecular modeling, we identified ECL2 residues that interact either directly or indirectly with orthosteric agonists and antagonists. Discrete mutations proximal to a conserved ECL2-TM3 disulfide bond selectively affected orthosteric ligand affinity, whereas a cluster of five residues near the TM4-ECL2 juncture influenced orthosteric agonist efficacy. A combination of ligand docking, molecular dynamics simulations, and mutagenesis results suggested that the orthosteric agonist 5'-N-ethylcarboxamidoadenosine binds transiently to an extracellular vestibule formed by ECL2 and the top of TM5 and TM7, prior to entry into the canonical TM bundle orthosteric site. Collectively, this study highlights a key role for ECL2 in AAR orthosteric ligand binding and receptor activation.

The binding of a small molecule ligand to its protein target is most often characterized by binding affinity and is typically viewed as an on/off switch. The more complex reality is that binding involves the ligand passing through a series of intermediate states between the solution phase and the fully bound pose. We have performed a set of 29 unbiased molecular dynamics simulations to model the binding pathways of the dopamine receptor antagonists clozapine and haloperidol binding to the D2 and D3 dopamine receptors. Through these simulations we have captured the binding pathways of clozapine and haloperidol from the extracellular vestibule to the orthosteric binding site and thereby, we also predict the bound pose of each ligand. These are the first long time scale simulations of haloperidol or clozapine binding to dopamine receptors. From these simulations, we have identified several important stages in the binding pathway, including the involvement of Tyr7.35 in a "handover" mechanism that transfers the ligand between the extracellular vestibule and Asp3.32. We have also performed interaction and cluster analyses to determine differences in binding pathways between the D2 and D3 receptors and identified metastable states that may be of use in drug design.

Biological membranes are natural barriers to the transport of molecules and drugs within human bodies. Many antibacterial agents need to cross these membranes to reach their target and elicit specific effects. Kanamycin A belongs to the family of aminoglycoside antibiotics that target cellular RNA to inhibit bacterial and viral replication. Previous studies have shown that aminoglycosides bind to mammalian but disrupt bacterial membranes. In this study, molecular dynamics (MD) simulations and infrared (IR) spectroscopy were applied to investigate the initial, first key interactions of kanamycin A, as a representative aminoglycoside, with both bacterial and mammalian lipid bilayers at the molecular level. Computational studies revealed strong hydrogen bonding interactions between the hydroxyl and amino groups of the aminoglycoside with the ester carbonyl and phosphate groups of the lipids. IR spectroscopy provided experimental verification of the important role of the lipid's ester carbonyl, phosphate and hydroxyl groups for aminoglycoside binding. The bacterial membrane became disordered upon aminoglycoside addition, whereas the mammalian membrane became stiffer and more ordered. This indicates the bacterial membrane disruption observed by previous studies.

Imatinib, a drug used for the treatment of chronic myeloid leukemia and other cancers, works by blocking the catalytic site of pathological constitutively active Abl kinase. While the binding pose is known from X-ray crystallography, the different steps leading to the formation of the complex are not well understood. The results from extensive molecular dynamics simulations show that imatinib can primarily exit the known crystallographic binding pose through the cleft of the binding site or by sliding under the αC helix. Once displaced from the crystallographic binding pose, imatinib becomes trapped in intermediate states. These intermediates are characterized by a high diversity of ligand orientations and conformations, and relaxation timescales within this region may exceed 3−4 ms. Analysis indicates that the metastable intermediate states should be spectroscopically indistinguishable from the crystallographic binding pose, in agreement with tryptophan stopped-flow fluorescence experiments.